Hello,
I’m trying to simulate a Cesium139 spectrum in a gamma spectroscopy simulation, but the best I can get for this isotope is the following graph:

I’m basing myself off of the rdecay02 example. Simulating Barium is no issue, but I can’t get Cesium to work. I’d appreciate any help!
Thanks
@iroisec
you are simulating Cs-139 and does not have any intense gammas leaving 627 (1% BR) and 1283 (8% BR). The spectrum you presented is from 0 to 200 keV, where only two gammas decaying via beta mode (b-) with BR less than 0.001% times.
Eg (keV) |
Ig (%) |
Decay mode |
188.88 20 |
0.0091 17 |
b- |
196.51 18 |
0.0100 17 |
b- |
230.76 9 |
0.036 3 |
b- |
Ref: - Table of Isotopes decay data
VRS
Hi, thanks for the answer! From what I understood about G4 simulations, the reactions from daughter isotopes are also tracked, so I was expecting a 166keV peak as well. I’m assuming that the issue has to do with tracking secondaries, but I’m not sure…
For tracking secondaries i suggest you to calculate the surface current of particles ( first hit of particle in tracking medium) and choose a gate on a particle to reproduce spectra. This you can achieve in steppinAction class.
VRS
@iroisec
IsFirstStepInVolume();
volume is any tracking medium.
Ref: Geant4: G4Step Class Reference
VRS