Pythia - Event generation
More details on the first lessons
The code we used to generate elementary proton+proton -> JPsi in the first lessons is based on the generator Pythia 6.2. To understand what we did, have a look in the $TUTORIAL/event_gen directory to :
- The gen_pythia script which launch an event production.
- Two arguments : $EVENTS $NAME (by default : 1000 jpsi)
- One input file : $NAME_pyt.par (for instance jpsi_pyt.par)
If the file exists in your local directory the script will use it, otherwise it will take it from $TUTORIAL/event_gen/parpythia
- Its outputs are :
- $NAME_pyt.log contains the cout and cerr output of the generation.
- $NAME_pyt.dat contains the quadri-momenta of the stable particles going out of the process (in a computer format)
- $NAME_pyt.root contains an ntuple with the accepted events (pid and momentum of each particle and its parent and grand parent).
- The jpsi_pyt.par file which contains the cards used for our JPsi simulation. The comments in the file inform you on the meaning of the variables.
Running your own stuff
Thus, if you want to generate 999 events of your stuff, then you just have to build your own mystuff_pyt.par in your directory and to launch :
> $TUTORIAL/event_gen/gen_pythia 999 mystuff
Then you can execute the other steps the same way we did in the first lessons :
> $TUTORIAL/simu/simulate 999 mystuff
> $TUTORIAL/resp/respond 999 mystuff
> $TUTORIAL/preco/reconstruct 999 mystuff
Even more details !
The code we use in the tutorial to generate JPsi events was developped by Frédéric.
Please refer to his pythia web page for more details, including :
Page maintained by Raphaël /
Last modified: Thu Feb 21 11:45:05 CET 2002
