Tequila test

Reproduce STAT paper hadron RAA

General steps

• STAT group saves their hydro files in the Globus endpoint JETSCAPE-NEW.
• The data for the paper STAT paper is save in the path /Production_Runs_for_Quenching/20180405_235533_production_AuAu200_0-70/
• For the files in the corresponding folder, the generated hydro data is saved in Results/input_x_20.
• x in the folder name means the centrality of the events. For example, input_5_20 means the centrality between 5% and 6%.  My saved location on Nersc is /global/cscratch1/sd/td115/STAT_Hydro_data/20180405_235533_production_AuAu200_0-70/Results/input_x_20/event-xx. 
• The hydro folders should be renamed as event-x for JETSCAPE to read. For example, the folder event-00 should be renamed as event-0.
• In XML file, set initial_profile_path and hydro_from_file as the directory to hydro files, saying /path-to-hydro. We should put folders event-x in /path_to_hydro. The content of event-x are EventPlanesFrzout.dat, initial.hdf5, JetData.h5, results and surface.dat. Initial state information is saved in initial.hdf5 and VISHNew hydro information is saved in JetData.h5.
• JETSCAPE automatically reads hydro information from event-x one by one. We cannot specify a certain hydro file for a certain job. Set nReuseHydro to run several jet events for a certain hydro file. A possible xml setting is nEvent=100 and nReuseHydro=20, which means the first 20 jet events will be run using hydro event-0, and the second 20 jet events will be run using hydro event-1, and so on.
• If we want to run events for 0-10% centrality, we should choose same number of hydro files from 0-1% to 9-10% and put them in the same path \path-to-hydro and rename them as event-0 to event-x.
• To produce results with reasonable precision on Nersc, we can run 4000 jet events for each pTHatBin. For example, we submit one job for each pTHatBin to Nersc to run calculation in 40-50% centrality. In each job, we set nEvent=4000 and nReuseHydro=20, so we should put 200 hydro files from event-0 to event-199 in \path-to-hydro, and 20 hydro files for each 1% centrality.
• If we want to run for a better precision, we can submit 10 jobs for each pTHatBin. For each job, we should set a different value of seed!
• In XML file, we should set maxT to be 150 fm or 250 fm for JETSCAPE events with hydro input. tStart is the time to start jet calculation, but not the time to start hydro. eCMforHadronization should be set as half of center-of-mass energy.
• For AA collisions, we should use nuclear modifications to parton densities. STAT paper uses EPS09 LO results. In Pythia, we should set PDF:useHardNPDFA, PDF:useHardNPDFB as on, and also set PDF:nPDFSetA, PDF:nPDFSetB as 1.

Compare with JETSCAPE results

• JETSCAPE has published STAT paper for their AuAu 200GeV calculation using LBT and MATTER.
• Unfortunately, we cannot directly compare the public JETSCAPE results with STAT paper. STAT people actually used a LBT patch code (fixed coupling) when publishing the paper, which code is not included in public JETSCAPE repository. The default LBT code in public JETSCAPE is running coupling.
• When running LBT with vacuum MATTER using public JETSCAPE, we would better to choose $Q_0 = 2GeV$, while the coupling constant as $\alpha_s = 0.35$ to have enough energy loss.

Brick test

This brick test aims to quickly check if Tequila is running properly.

Total interaction rate