I am simulating a gamma camera/scintillator, and need to separate the total number of produced photons by gamma-scintillator interaction location and process(photoelectric, Compton…).
So for example, if a photoelectric interaction occurs, the gamma interacts and produces a secondary electron, which in turn produces photons in later steps. So every secondary particle must be tracked for several steps until all possible photon production has been recorded.
How would I count these photons? Currently, the interaction locations are recorded from SteppingAction, but given the example above, a simple count of secondary photons is insufficient (as in the photoelectric example above, most photons are produced by the secondary electron).
I appreciate any suggestions you may have.