I am trying to model a TAC(Time to Amplitude converter) using my UserStepping.cc.
The idea is that I want to get the time difference between the detection of two specific gammas emitted from radioactive decay G4 ion which is in coincidence.
I am thinking of starting a timer as soon as fEdep == gamma1 and stopping it as soon as fEdep == gamma2, finally returning the time elapse. Also, it will be great if we scan for the second gamma (E2) only in a small time cutoff after detecting the first.
Can someone please share any suggestions on how to go about this or alternate approaches to achieving the same? Are there any relevant geant4 examples for the same?
PS: I am trying to not use Get->GlobalTIme() due to “issues” in getting the time of flight when the half-life is in order of several years.
Thanks a ton!