Hello, I simulated the radioactive decay of some radionuclides (Am241, Co60, Cs134, Cs137, Pb-210).
My tutor said that the deposited energy in the crystal by C060, Cs134 and Cs137 are good, while the Am241 and Pb210 are wrong). There is a problem in Geant4 for these radionuclides. My tutor said me that Am241 and Pb210 have radioactive sons then he think that maybe GEANT4 simulates all the decays instead of just Am241 and Pb210.
Does anyone know something about the simulation of these two radionuclides and how to fix??
You simulated these radioactive decays but didn’t mention how you specified the radionuclides. Did you use /gps/ commands in a macro?
In the past I had the opposite problem, i.e. the decay of (very) long-living states ‘polluted’ my spectrum, therefore, I think Geant4 simulates the full decay chain. But you should be able to easily check this yourself using the tracking verbosity.
I’m almost sure about Pb210 because on Table of Isotopes decay data I read only a peak at 45.5keV…
Regarding the Am241 I’ve some doubts because on Table of Isotopes decay data it gives so many peaks. Nevertheless, the Igamma% is at energy 59.6keV and I see a peak for this energy value.
I also tried to follow the radioactivedecay example, as suggested by @maire , then:
I wrote trackingaction.cc , trackingaction.hh, trackingmessenger.cc, trackingactionmessenger.cc
I am not sure if you have solved it yet. The reason why Am giving very high energy lines might be due to the spontaneous fission of 241Am. Again as mentioned by so many others in this post /grdm/nucleusLimts can solve the problem I think