Dear Developers, I am shooting a 197Au projectile of 1 GeV/nucleon on a 0.5 mm thick 197Au target in air to study the accompanied delta-electron production along the projectile track. The program in principle works until the charge exchange reaction happens in the target with a secondary neutron. After a few hours of running the program crashes after 100 of reentries in G4HadronicProcess.cc: ........................................................................ JL reentryCount: 99 JL Process / Model: chargeExchange / Charge Exchange 0/ 1 ParticleName proton mass_pdg 938.272013 mass_dyn 938.272013 TotalEnergy 939.724490 JL G4HadronicProcess.cc !!!!!!!!!!!!!!!!!!!!!!! deltaE 183434.547811 nuclearMass 183433.330971 aPro.GetTotalEnergy() 940.941330 finalE 939.724490 nSec 1 JL G4HadronicProcess.cc !!!!!!!!!!!!!!!!!!!!!!! -------- WWWW ------- G4Exception-START -------- WWWW ------- *** G4Exception : had012 issued by : G4HadronicProcess:CheckResult() Warning: Bad energy non-conservation detected, will re-sample the interaction Process / Model: chargeExchange / Charge Exchange Primary: neutron (2112), E= 940.941, target nucleus (79, 197) E(initial - final) = 183435 MeV. *** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW ------- JL reentryCount: 100 JL Process / Model: chargeExchange / Charge Exchange 0/ 1 ParticleName proton mass_pdg 938.272013 mass_dyn 938.272013 TotalEnergy 939.726264 JL G4HadronicProcess.cc !!!!!!!!!!!!!!!!!!!!!!! deltaE 183434.546036 nuclearMass 183433.330971 aPro.GetTotalEnergy() 940.941330 finalE 939.726264 nSec 1 JL G4HadronicProcess.cc !!!!!!!!!!!!!!!!!!!!!!! -------- WWWW ------- G4Exception-START -------- WWWW ------- *** G4Exception : had012 issued by : G4HadronicProcess:CheckResult() Warning: Bad energy non-conservation detected, will re-sample the interaction Process / Model: chargeExchange / Charge Exchange Primary: neutron (2112), E= 940.941, target nucleus (79, 197) E(initial - final) = 183435 MeV. *** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW ------- -------- EEEE ------- G4Exception-START -------- EEEE ------- *** G4Exception : had006 issued by : G4HadronicProcess::PostStepDoIt Call for Charge Exchange Target element Gold197 Z= 79 A= 197 Unrecoverable error in the method ApplyYourself of chargeExchange TrackID= 122842 ParentID= 1 neutron Ekin(GeV)= 0.00137597; direction= (0.999393,-0.033417,-0.00983938) Position(mm)= (13.2471,-0.464407,599.877); material Au197 PhysicalVolume <tgt_phys> ApplyYourself does not completed after 100 attempts *** Fatal Exception *** core dump *** G4Track (0x4748fa50) - track ID = 122842, parent ID = 1 Particle type : neutron - creator process : ionInelastic, creator model : Undefined Kinetic energy : 1.37597 MeV - Momentum direction : (0.999393,-0.033417,-0.00983938) Step length : 1.33082 cm - total energy deposit : 0 eV Pre-step point : (-0.0530347,-0.0196863,600.008) - Physical volume : tgt_phys (Au197) - defined by : not available Post-step point : (13.2471,-0.464407,599.877) - Physical volume : tgt_phys (Au197) - defined by : chargeExchange - step status : 4 *** Note: Step information might not be properly updated. -------- EEEE -------- G4Exception-END --------- EEEE ------- *** G4Exception: Aborting execution *** ........................................................................ I have added some printouts starting with "JL" above in G4HadronicProcess.cc. The program terminates due to a huge energy non-conservation of 183435 MeV which roughly corresponds to the 197Au mass. It seems that in this critical event there are 2 particles in the entrance channel: the neutron and the 197Au nucleus, while in the exit channel there is only one: the proton (nSec=1). In principle there should be also a heavy companion, i.e. 197Pt in that case, so nSec should be = 2, not 1. Indeed, within 95 good events before the crash there is a similar one which ends successfully after 5 reentryCounts. I put it all below: ........................................................................ JL reentryCount: 0 JL Process / Model: chargeExchange / Charge Exchange 0/ 1 ParticleName proton mass_pdg 938.272013 mass_dyn 938.272013 TotalEnergy 957.551180 JL G4HadronicProcess.cc !!!!!!!!!!!!!!!!!!!!!!! deltaE 183434.546318 nuclearMass 183433.330971 aPro.GetTotalEnergy() 958.766528 finalE 957.551180 nSec 1 JL G4HadronicProcess.cc !!!!!!!!!!!!!!!!!!!!!!! -------- WWWW ------- G4Exception-START -------- WWWW ------- *** G4Exception : had012 issued by : G4HadronicProcess:CheckResult() Warning: Bad energy non-conservation detected, will re-sample the interaction Process / Model: chargeExchange / Charge Exchange Primary: neutron (2112), E= 958.767, target nucleus (79, 197) E(initial - final) = 183435 MeV. *** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW ------- JL reentryCount: 1 JL Process / Model: chargeExchange / Charge Exchange 0/ 1 ParticleName proton mass_pdg 938.272013 mass_dyn 938.272013 TotalEnergy 957.544826 JL G4HadronicProcess.cc !!!!!!!!!!!!!!!!!!!!!!! deltaE 183434.552673 nuclearMass 183433.330971 aPro.GetTotalEnergy() 958.766528 finalE 957.544826 nSec 1 JL G4HadronicProcess.cc !!!!!!!!!!!!!!!!!!!!!!! -------- WWWW ------- G4Exception-START -------- WWWW ------- *** G4Exception : had012 issued by : G4HadronicProcess:CheckResult() Warning: Bad energy non-conservation detected, will re-sample the interaction Process / Model: chargeExchange / Charge Exchange Primary: neutron (2112), E= 958.767, target nucleus (79, 197) E(initial - final) = 183435 MeV. *** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW ------- JL reentryCount: 2 JL Process / Model: chargeExchange / Charge Exchange 0/ 1 ParticleName proton mass_pdg 938.272013 mass_dyn 938.272013 TotalEnergy 957.551503 JL G4HadronicProcess.cc !!!!!!!!!!!!!!!!!!!!!!! deltaE 183434.545995 nuclearMass 183433.330971 aPro.GetTotalEnergy() 958.766528 finalE 957.551503 nSec 1 JL G4HadronicProcess.cc !!!!!!!!!!!!!!!!!!!!!!! -------- WWWW ------- G4Exception-START -------- WWWW ------- *** G4Exception : had012 issued by : G4HadronicProcess:CheckResult() Warning: Bad energy non-conservation detected, will re-sample the interaction Process / Model: chargeExchange / Charge Exchange Primary: neutron (2112), E= 958.767, target nucleus (79, 197) E(initial - final) = 183435 MeV. *** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW ------- JL reentryCount: 3 JL Process / Model: chargeExchange / Charge Exchange 0/ 1 ParticleName proton mass_pdg 938.272013 mass_dyn 938.272013 TotalEnergy 957.550767 JL G4HadronicProcess.cc !!!!!!!!!!!!!!!!!!!!!!! deltaE 183434.546731 nuclearMass 183433.330971 aPro.GetTotalEnergy() 958.766528 finalE 957.550767 nSec 1 JL G4HadronicProcess.cc !!!!!!!!!!!!!!!!!!!!!!! -------- WWWW ------- G4Exception-START -------- WWWW ------- *** G4Exception : had012 issued by : G4HadronicProcess:CheckResult() Warning: Bad energy non-conservation detected, will re-sample the interaction Process / Model: chargeExchange / Charge Exchange Primary: neutron (2112), E= 958.767, target nucleus (79, 197) E(initial - final) = 183435 MeV. *** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW ------- JL reentryCount: 4 JL Process / Model: chargeExchange / Charge Exchange 0/ 1 ParticleName proton mass_pdg 938.272013 mass_dyn 938.272013 TotalEnergy 957.551798 JL G4HadronicProcess.cc !!!!!!!!!!!!!!!!!!!!!!! deltaE 183434.545701 nuclearMass 183433.330971 aPro.GetTotalEnergy() 958.766528 finalE 957.551798 nSec 1 JL G4HadronicProcess.cc !!!!!!!!!!!!!!!!!!!!!!! -------- WWWW ------- G4Exception-START -------- WWWW ------- *** G4Exception : had012 issued by : G4HadronicProcess:CheckResult() Warning: Bad energy non-conservation detected, will re-sample the interaction Process / Model: chargeExchange / Charge Exchange Primary: neutron (2112), E= 958.767, target nucleus (79, 197) E(initial - final) = 183435 MeV. *** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW ------- JL reentryCount: 5 JL Process / Model: chargeExchange / Charge Exchange 0/ 2 ParticleName proton mass_pdg 938.272013 mass_dyn 938.272013 TotalEnergy 957.170432 JL Process / Model: chargeExchange / Charge Exchange 1/ 2 ParticleName Pt197 mass_pdg 183434.544684 mass_dyn 183434.544684 TotalEnergy 183434.927066 JL G4HadronicProcess.cc !!!!!!!!!!!!!!!!!!!!!!! deltaE 0.000000 nuclearMass 183433.330971 aPro.GetTotalEnergy() 958.766528 finalE 184392.097499 nSec 2 JL G4HadronicProcess.cc !!!!!!!!!!!!!!!!!!!!!!! ........................................................................ So, why the above event recovered successfully after 5 tries and the critical one did not within 100 tries? Would you have any simple workarounds on how to avoid the "Fatal Exception"? What is wrong with the chargeExchange? I have tested ModularPhysicsList and FTFP_BERT. Both crash after some hours. Best regards, Jerzy
Hello, it is not clear what Physics Lists is used and what Geant4 version, because FTFP_BERT does not have chargeExchange process. How this process is enabled? VI
Hello, because required information is not provided we will close this problem report. The problem likely exist but is not reproduced in a simple setup. VI
We close this report. If new information appears it will be opened. Vladimir
