[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

StarPythia8: ” https://www.star.bnl.gov/webdata/dox/html/classStarPythia8.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 !