Home » Groups » STAR Protected

Software and Computing : PYTHIA: Understanding the event record

Updated on Wed, 2016-06-22 16:14. Originally created by devika on 2013-05-30 10:35.
PYTHIA: Understanding the event record 

Example event record: only first 13 lines are showing here

                            Event listing (summary)
 
    I particle/jet  KS     KF        orig    p_x      p_y      p_z                E            m
 
    1 !p+!          21    2212      0     0.000    0.000    254.998   255.000     0.938
    2 !p+!          21    2212      0     0.000    0.000   -254.998   255.000    0.938
 ==============================================================================
    3 !u!            21       2         1    -0.031    0.648    127.962   127.964    0.000
    4 !dbar!       21      -1        2     0.509   -0.635    -28.681   28.693     0.000
    5 !u!            21       2         3    5.641    3.275      62.096    62.437     0.000
    6 !dbar!       21      -1        4     0.457   -0.570    -25.754   25.764     0.000
    7 !W+!        21      24        0     6.097     2.705    36.342     88.201     80.089
    8 !e+!         21     -11       7    -11.317   24.842  -13.877   30.623    0.001
    9 !nu_e!      21      12        7    17.414   -22.137   50.219   57.578     0.000
 ==============================================================================
   10 (W+)        11      24      7     6.097     2.705       36.342     88.201    80.089
   11 e+           1     -11       8    -11.317   24.842    -13.877   30.623    0.001
   12 nu_e        1      12        9     17.414    -22.137    50.219   57.578     0.000
......

Column 1: I           ->  Particle ID
Column 2: Particle -> Particle Name
Column 3: KS        -> Particle state code
                     
                  Ground rule: codes 1-10: Currently existing partons/particles
                                      codes >10  : No longer existing particles/partons or other kind of event information

                                      eg:   1 - An undecayed particle or unfragmened parton
                                            11 - A decayed particle or fragmened parton
                                             
 Column 4: KF       -> Particle code
                                     eg:  2212 -> P+
                                            -1     ->dbar                                        

Column 5: orig     -> Particle ID of the Particle which originated from

Column 6,7,8 : Px,Py,Pz  -> Paricle momentum in X,Y,and Z directions in GeV/c.
Column 9:  E        -> Energy of the particle in GeV.
Column 10: m     -> Mass of the particle in GeV/c2

 


More details and description : Pythia manual 6.4: Chapter 5

 

How the event record works,

A few examples may help clarify the picture. For a single diffractive event ppbar -> pdiffrpbar,

the event record will start with
 

I     K(I,1) (KF)         K(I,2)(KS)           K(I,3) (orig)           comment

1    21                        2212                 0                  incoming p

2    21                       -2212                 0                 incoming p

========================= not part of record; appears in listings

3    21                        9902210            1              outgoing pdiffr

4    21                            -2212              2                outgoing p

========================= again not part of record

The typical QCD 2 ! 2 process would be
 

I      K(I,1)  (KF)          K(I,2)(KS)        K(I,3)(orig)        comment

1      21                       2212                  0               incoming p

2      21                         -2212              0            incoming p

=========================

3      21                                 2                    1                 u picked from incoming p

4      21                                -1                    2               d picked from incoming p

5      21                                21                   3              u evolved to g at hard scattering

6      21                                -1                    4              still d at hard scattering

7      21                                21                   0              outgoing g from hard scattering

 



Groups: