I simulated U, Th decay series and other natural radionuclides with Geant4 example ; examples/extended/radioactivedecay/rdecay01.
In the case of the U decay chain, the gamma-ray yield values match well with the BIPM table of nuclides, but notably, the yield values for the thorium decay chain appear to be consistently lower. For instance, in the case of Pb-212, one of daughter isotopes of Th-232, which should have a 238.632 keV gamma-ray yield of around 43.6%, it is instead approximately 34.48%. It appears that in other Monte Carlo codes, the simulations work as expected, but in Geant4, the results are consistently underestimated for Th decay chain. Please explain why this is the case…
_Operating System:Window 10 OS