I have implemented radioactive decays of certain isotopes. I am able to export relevant particle-related information (see attached). I am wondering how to determine the generation process of emitted particle(s)? For example, in the attached case, how to determine that 129.359 keV electron was produced from internal conversion. Similarly, discriminating beta-decay electrons from internal conversion for example?
Can we set up something in SteppingAction to export this information at each step to a ROOT branch/console? Any assistance is appreciated.
Hello, I dont know if it can help a lot but for my application i wanted to know generation process for gamma to discriminate neutron induced gamma from other generation process and i did this in StyeppingAction.cc :
@Geoffrey_Varignier : Thanks, I just tried. The GetProcessName() gives me âRadioactiveDecayBaseâ Is there a way to look deeper which decay channel lead to this? For example internal conversion?
Be careful with your language. The decay process (in Geant4 talk), is RadioactiveDecayBase. Full stop. What youâre asking for is which âdecay channelâ was selected by the RadioactiveDecayBase process for the particular nuclear decay instance.
I am not sure that RDecay has a public interface to register or report which channel was chosen. It might be that you have to âinferâ it in your SteppingActon by looking at the secondaries, but thatâs a terrible approach.
@mkelsey : Thanks for noting that. I have updated the words; I am still new to G4 and appreciate you clarifying this. Can you please shed some more light on how that information can be âinferredâ using secondaries? I guess any way via which I can discriminate beta-decay and internal conversion electrons is good for me?
To be a bit more concrete. What youâd do in your stepping action is look to see if the particle has decayed (status fStopAndKill, secondaries in the current step). Loop over the secondaries (use GetSecondaryInCurrentStep() to avoid picking up other secondaries from past steps), and look for the kind of decay you expect (if thereâs an alpha, if thereâs a neutrino, and so on).