I’m working on a simulation using Geant4 where I want to pick the decays that are happening in the simulation based on lifetime of the particle. If the lifetime is longer than a given value I want to stop the decay chain.
For Radioactive particles this worked fine by using for the G4ParticleDefinition::
SetPDGLifetime(-1). This stopped all decays of the radioactive particles.
This however did not work for the exited states in the simulation.
I simulated Cs137 decay which first goes to Ba137[661.659] before releasing its characteristic 662KeV photon and going to Ba137. The only way I could find to stop the exited state from decaying is to stop the track from propagating by
I would however like to propagate these states without them “relaxing” to their non-exited state.
Does anyone have an idea how I can do this?
(note: I already tried setting
SetPDGStable(true) or deleting its decay table but neither worked)