I ran these two commands and looked at the output:
gevgen -s -n 100 -p 14 -t 1000060120 -e 1
gevdump -f gntp.0.ghep.root
Many events do not conserve baryon number, and some do not conserve charge.
Attached is a screenshot of two puzzling events. The initial state is always a 1
GeV nu_mu plus C12.
The first has a final state consisting of:
C11, mu-, 2 neutrons, 3 protons
That conserves neither baryons nor charge.
The second has a final state that includes a "NucBindE" with PDGid 2000000101 --
why does that have Ist=1 ?? -- I thought I could select particles with Ist==1 to
send to my other simulation program, but that obviously is not enough -- how
should I determine the final particles (plus recoil nucleus)?
Also: what are the notations off to the right? In particular, what is the
3-vector P with each final-state lepton? and how does it relate to (Px,Py,Pz) ?
Obviously I am missing something. Can anyone please explain....
Tom Roberts
|