hepmc - Blame information for rev 432
Subversion Repositories:
| Rev | Author | Line No. | Line |
|---|---|---|---|
| 76 | garren | 1 | ////////////////////////////////////////////////////////////////////////// |
| 2 | // garren@fnal.gov, April 2007 | ||
| 3 | // | ||
| 4 | // Initialize pythia | ||
| 5 | // Using a separate cc file allows us to test issues with PythiaWrapper.h | ||
| 6 | // | ||
| 7 | ////////////////////////////////////////////////////////////////////////// | ||
| 8 | |||
| 9 | #include "HepMC/PythiaWrapper.h" | ||
| 432 | garren | 10 | #include "PythiaHelper.h" |
| 76 | garren | 11 | |
| 12 | void initPythia() | ||
| 13 | { | ||
| 14 | // (Some platforms may require the initialization of pythia PYDATA block | ||
| 15 | // data as external - if you get pythia initialization errors try | ||
| 16 | // commenting in/out the below call to initpydata() ) | ||
| 17 | // initpydata(); | ||
| 18 | // | ||
| 19 | // Select W+gamma process (process number 20) | ||
| 20 | // (here we have to be careful of C/F77 differences: arrays in C | ||
| 21 | // start at 0, F77 at 1, so we need to subtract 1 from the process #) | ||
| 22 | pysubs.msel=0; | ||
| 23 | pysubs.msub[20-1] = 1; | ||
| 24 | // set random number seed (mandatory!) | ||
| 25 | pydatr.mrpy[0]=55122 ; | ||
| 26 | // Tell Pythia not to write multiple copies of particles in event record. | ||
| 27 | pypars.mstp[128-1] = 2; | ||
| 28 | // Example of setting a Pythia parameter: set the top mass | ||
| 29 | pydat2.pmas[1-1][6-1]= 175; | ||
| 30 | // | ||
| 31 | // Call pythia initialization | ||
| 32 | call_pyinit( "CMS", "p", "p", 14000. ); | ||
| 33 | } |
