Adding particle type and impact on eg: scattering angle with GBL
High momentum particles are by default used in the code (beta=1
).
For low momentum, the scattering angle formula in https://gitlab.cern.ch/corryvreckan/corryvreckan/-/blob/master/src/objects/GblTrack.cpp#L120 does not hold any more. A correction is needed as per the DPG (page 16, https://pdg.lbl.gov/2022/reviews/rpp2022-rev-passage-particles-matter.pdf) to account for beta
<1.
For now, in a separate offline branch, we use a header file with a struct containing particle type, charge number and mass is used, together with a formula to calculate the beta
and to transition if needed from kinetic energy to momentum and back.
The question is: Is the particle type needed in the code anywhere else and should this header file be sufficient? Is a more thorough implementation needed?