Hi,
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 SetTrackStatus(fStopAndKill)
.
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)
Cheers,
Koen