Dear Geant4’s users,
I’m landing 20 MeV neutron againt 500 cm of H2 (density =1.0g/cm3) divided in 50 slabs of 10cm with FTFP_BERT_HP physic list, studyng the flux of all particles resulting. I just replaced the G4IonPhysics physic by the G4IonPhysicsPHP one.
The simulation seems to led to a non conservation of mass and crash :
‘‘‘
-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : hadr01
issued by : G4ParticleHPFinalState:adjust_final_state
G4ParticleHPFinalState oneMoreSec_pd Z=2 - 0 A=4 - 4 projectile deuteronNo adjustment will be done!
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------
G4IonTable::GetIon() : illegal atomic number/mass Z =-1 A = 1 E = 0
‘‘‘
I was wondering if this could come from a non well association of HP physic list parameters, here are mine :
‘‘‘
=======================================================
====== ParticleHP Physics Parameters ========
Use only photo-evaporation 0
Skip missing isotopes 0
Neglect Doppler 0
Do not adjust final state 0
Produce fission fragments 0
Use WendtFissionModel 0
Use NRESP71Model 0
Use DBRC 0
PHP use Poisson 1
PHP check 1
CHECK HP NAMES 0
Enable DEBUG 0
Use probability tables from
=======================================================
‘‘‘
Thanks for your reading time.
Geant4 Version: 11.3.0*