Can't get proper Cesium spectrum

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:
Screen Shot 2024-07-02 at 12.43.28
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