Dear Geant4’s users,
I’m simulating 500 MeV Triton against He4 slabs. I’m currently using the G4IonPhysicPHP in order to use the TENDL database. I already aware that TENDL doesn’t include triton data into He4 so theoritical model should be uses but that’s not the case, it leads to a crash.
Futhermore I suspect a logical issue for models HP used by ions. Models borders seems to be initialised in MeV whereas model selection apply in MeV/nuc in G4EnergyRangeManager::GetHadronicInteraction.
That can leads to an issue when an ion energy is closed to the border. The HP model can be selected where it should not be and when the process try to use it, Geant4 crashs laking of HP data.
For exemple :
The maximum border of my HP model is set to 200 MeV.
The triton arrives at 500 MeV to the the model selection, so approximtely at 166 MeV/nuc. As the kinetic energy is inside the border of HP model, the HP model is selected.
When the HP model will try to be applied, cross sections will be gathered, but the kinetic energy used for cross section gathering is in MeV. As the maximum border of the HP cross section dataset is 200 MeV Geant4 will raise an exception.
Does this logic issue (if there is one) is already known ?
Geant4 Version: 11.3.3