I am running a simulation using a carbon beam, which irradiates a simple cubic phantom of different materials (water, Bismuth, Gadolinium, etc), and retrieving the secundary gammas produced.
When running it I came across this error, which appeared very often for the simulated events:
"G4BinaryLightIonReaction.cc: mom check:
G4WT20 > initial - final (-249.546,-171.35,2308.67;5621.6) 3.mag 2328.43
G4WT20 > .. pInitialState/pFinalState/spectators
G4WT20 > (0,0,4330.94;26879.8)
G4WT20 > (249.546,171.35,2022.27;21258.2)
G4WT20 > (-249.546,-171.35,2308.67;6128.79)
G4WT20 > .. A,Z 6 4
G4WT20 > G4BinaryLightIonReaction invalid final state for:
G4WT20 > Primary 0x62d835e0, (A,Z)=(12,6) , kinetic energy 809.902
G4WT20 > Target nucleus (A,Z)=(16,8)
G4WT20 > if frequent, please submit above information as bug report"
Another error also pops up, way less frequent but still appears:
-------- WWWW ------- G4Exception-START -------- WWWW -------
G4Exception : HAD_FANCY3DNUCLEUS_001
issued by : G4Fancy3DNucleus::ChooseFermiMomenta(): difficulty finding proton momentum, set it to (0,0,0)
Nucleus Z A 6 12
proton with eMax=938.186
This is just a warning message.
-------- WWWW -------- G4Exception-END --------- WWWW -------
I am defining the carbon beam through a macro in the following way:
/gps/particle ion
/gps/ion 6 12
/gps/pos/type Beam
/gps/pos/centre 0.0 0.0 400. mm
/gps/pos/sigma_x 4.756 mm
/gps/pos/sigma_y 4.756 mm
/gps/ang/type beam1d
/gps/ang/sigma_x 0.0032 rad
/gps/ang/sigma_y 0.0032 rad
/gps/ene/type Gauss
/gps/ene/mono 3421.2 MeV
/gps/ene/sigma 50 keV
And I have added to my simulation the following modular physics:
G4EmStandardPhysics_option4(); G4EmExtraPhysics(); G4DecayPhysics(); G4HadronElasticPhysicsHP(); G4HadronPhysicsQGSP_BIC_HP(); G4StoppingPhysics(); G4IonPhysics(); G4RadioactiveDecayPhysics();
Additionally I have turned on/off the following electromagnetic flags in the macro file:
/process/em/fluo true
/process/em/pixe true
/process/em/pixeXSmodel ECPSSR_ANSTO
/process/em/auger true
/process/em/augerCascade true
/process/em/deexcitationIgnoreCut false
Could someone shed some ligth into the meaning of these errors and how could I fix it?
Geant4 Version: 11.3.2
_Operating System:_Linux
_Compiler/Version:_GCC14.1.0
_CMake Version:_3.30.0