Hello everyone,
I’m trying to use Geant4 to get the beta spectrum of Tb161 beta decay. I modified the example/extended/radioactivedecay\rdecay01\src\TrackingAction.cc to fit my need. Here is my code:
// energy spectrum
//
G4int ih = 0;
if (particle == G4Electron::Electron() || particle == G4Positron::Positron())
{
int creator_id = track->GetCreatorModelID();
if(25520==creator_id || 25530==creator_id) ih = 1;
}
The e-/e+ decay creator ID are 25520, 25530 respectively.
And this is my macro for Tb161:
# macro for rdecay01
#
/control/cout/ignoreThreadsExcept 0
/control/verbose 2
/run/verbose 1
#
# Set a very high time threshold to allow all decays to happen
/process/had/rdm/thresholdForVeryLongDecayTime 1.0e+60 year
#
/gun/particle ion
/gun/ion 65 161
#
/tracking/verbose 2
/run/beamOn 1
/tracking/verbose 0
#
/analysis/setFileName Tb161
/analysis/h1/set 1 600 0. 600 keV #e+ e-
/analysis/h1/set 2 600 0. 600 keV #neutrino
/analysis/h1/set 3 600 0. 600 keV #gamma
/analysis/h1/set 6 100 0. 2500 keV #EkinTot (Q)
/analysis/h1/set 7 150 0. 15e3 keV #P balance
/analysis/h1/set 8 100 0. 100. year #time of life
/analysis/h1/set 9 100 1. 3. MeV #EvisTot
#
/run/printProgress 1000
/run/beamOn 10000
The result of e- spectrum:
The mean energy is 140 keV. According to the NNDC database, the mean energy should be 154keV(NNDC Tb161).
I also tried many other nuclide like I131 which showed a consistent result with NNDC. There is only one discrepancy so far.
Dose anyone know what’s wrong here?
TrackingAction.cc (7.0 KB)
Geant4 Version: 11.3.2
Operating System: Ubuntu 22.04
Compiler/Version: gcc version 11.4.0
CMake Version: 3.30.1
