[HOW TO] run PYTHIA generator and embed at STAR
<1>Run pure pythia events
STAR provides pythia6 and pythia8 event generators :
StarPythia6: ” https://www.star.bnl.gov/webdata/dox/html/classStarPythia6.html”
Example for using StarPythia6:
- $ cvs co StRoot/StarGenerators/macros/starsim.pythia6.standalone.C
- $ ln -s StRoot/StarGenerators/macros/starsim.pythia6.standalone.C starsim.C
- $ root4star -q -b starsim.C
Modify the parameters(e.g. CMS energy, subprocess) in starsim.C depending on your analysis. Pythia has different tunes , pick the suitable one ( Please select the PDF that matches the tune you are using ) .
<P.S. Tune “perugia2012” is recommended for p-p 200 GeV collision if you use pythia6>
Relevant links or documents:
User Manual for Pythia6 : https://arxiv.org/pdf/hep-ph/0603175.pdf
Pythia6 version in StarPythia6 : $STAR/StRoot/StarGenerator/Pythia6_4_28/pythia-6.4.28.F
User Manual for Pythia8 : https://pythia.org//latest-manual/Welcome.html
Tunes for pythia8 : https://pythia.org/latest-manual/Tunes.html#anchor3
PDF index : https://lhapdf.hepforge.org/pdfsets.html
<2>Run pythia and embed
i. Copy the example to your work place
- $ cp -r /star/data01/pwg/jackw/example/example_prod_embed_pythia /YOURPATH/
ii. Set STAR software version and compile your analyzer
- $ starver SL16d_embed
- $ cons
SL16d_embed is for Run16 p-p 200GeV embedding , select the right one for your analysis.
iii. Product pythia events and do the embedding
- $ chmod +x run_embedding.sh
- $ ./run_embedding.sh
"run_pysim_pp200GeV.C " function products pythia events and return .fzd file which is needed in mixing with .daq file . Sometimes you may need specific final states . Add your own filter to “StRoot/StarGenerator/FILT” , and use AddFilter( FILTER YOUR DEFINED ) function before Init() .
Function “bfcTest.C” is for reconstructing MC event form .fzd files and real event for .daq files and do the mixing . Chain options are needed in these process ("Chain1Opt" for real events ,"Chain2Opt" for MC events and "Chain3Opt" for mixing), it is recommended to ask the people in embedding or production group to get advice depending on your own analysis.
Relevant links or documents:
Production chain option of real data : https://www.star.bnl.gov/devcgi/dbProdOptionRetrv.pl
STAR Geometry in simulation : https://drupal.star.bnl.gov/STAR/comp/prod/MCGeometry
iv. Mass production
Create your job files to run mass production . Before that , you may need to download those .daq files first .
- $ chmod +x GetDownloadList.sh
- $ chmod +x DownloadDaq.sh
- $ ./GetDownloadList.sh /YOUR/RUNLIST/ /PATH/TO/STORE/FILELIST/ /PATH/OT/STORE/YOUR/DAQFILES/ StartRun EndRun
- $ ./DownlaodDaq.sh /YOUR/RUNLIST/ /PATH/TO/STORE/FILELIST/ StartRun EndRun
To check the status of your request
- $ hpss_user.pl -w
Or in this website: http://www.star.bnl.gov/devcgi/display_accnt.cgi
Hooray !
Now enjoy your analysis !
- jackw's blog
- Login or register to post comments