I have a question about the use of the PART bank. When I redid celeg to
produce it, I simply copied the code that produced the MCTK bank, and
changed the variables accordingly. It occurs to me that that might not
have been a good idea.
The reason is that there is no mechanism in the PART bank to tell GSIM
not to track a particular particle; with the MCTK bank, we have the
beg_vtx and end_vtx variables; if end_vtx<>0, the particle has decayed,
and we don't track it. This was done so that we could store all of the
decay chain within the bos file (except, curiously enough, for the
original resonance that was generated at the beginning of the event).
Because we don't have this flexibility within the PART bank, it seems to
me that the logical thing to do is to not store the particles that have
decayed. It should be simple enough to do, but I wanted to warn
everyone (and wait for a few return comments) before I did that out of
hand.
I will still keep the particles in the MCTK bank, because we can still
make some use of that information.
--
John Price *** *** *** *** pricej@jlab.org
Where there is no solution, there is no problem.