Add time for decay vertices #31
Conversation
src/erhic/TreeToHepMC.cxx
Outdated
|
|
||
| // recalcuate vertex time using: dt = (SV-PV).Mag()/(pMother/eMother) | ||
| // we assume primary vertex is at (0,0,0) with time = 0 | ||
| // this is done for any secondary vertices 1 um away from primary vertex |
There was a problem hiding this comment.
There is no reason to make assumptions about the primary vertex. Simply use the distance to the mother's vertex (and appropriate beta). This way, you correctly capture cascading decays.
There was a problem hiding this comment.
I agree with it. The calculation doesn't depend on (0,0,0). It is for any primary vertex.
|
Using the mother particle vertex is easy to implement. Thinking of cascading decays, we need to add the decay time for the mother to daughters as well, right? |
Yes, but that's nothing special, you just need to use t of mother instead of t=0, doesn't matter if you think of the mother as collision vertex or not, it should all recursively work out. |
…th; ii) convert to cm for BeAGLE which uses mm
stability; untabify
…s what the afterburner for beam effects expects
3 participants
Briefly, what does this PR introduce?
What kind of change does this PR introduce?
Please check if this PR fulfills the following:
Does this PR introduce breaking changes? What changes might users need to make to their code?
Does this PR change default behavior?