Created attachment 212 [details] My modified version of example B1 which produce in Event 0 an inelastic event with energy not conserved Dear Geant4 developers, I am doing Geant4 simulations, where a 0.15 mm thick steel foil is bombarded with 7 TeV protons. I observed that the tail of the energy deposit distribution of protons goes to negative, Ordo(-10 GeV). I tried to find the root of the problem, and I found that during a protonInelastic process the energy seems not to be conserved correctly. Below I show the relevant output of my simplified code, where "Event 0" already contains such protonInelastic calculation. Using 938.2, 139.5, 134.9 MeV /c^2 as the proton, pi+/-, pi0 rest energy, the secondaries carry more energy than the incoming proton. I have attached my example code, which is basically example B1, but: - the 2 shapes are from G4_Al. The World and Envelope is G4_Galactic. - particle is 7 TeV proton - I calculate the energy deposit in the Step function (for 1 event only). The 0. event already contains such process, can be seen with /tracking/verbose 2 ********************************************************************************************************* * G4Track Information: Particle = proton, Track ID = 1, Parent ID = 0 ********************************************************************************************************* Step# X(mm) Y(mm) Z(mm) KinE(MeV) dE(MeV) StepLeng TrackLeng NextVolume ProcName 0 47.4 -53 -150 7e+06 0 0 0 Envelope initStep 1 47.4 -53 40 7e+06 9.1e-24 190 190 Shape2 Transportation 7000.938272013001 2 47.4 -53 47.1 7e+06 1.97 7.09 197 Shape2 hIoni :----- List of 2ndaries - #SpawnInStep= 1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 1 --------------- : 47.4 -53 47.1 0.603 e- :----------------------------------------------------------------- EndOf2ndaries Info --------------- 7000.938272013001 3 47.4 -53 52.5 0 1.53 5.45 203 Shape2 protonInelastic :----- List of 2ndaries - #SpawnInStep= 59(Rest= 0,Along= 0,Post=59), #SpawnTotal= 60 --------------- : 47.4 -53 52.5 4.47e+05 pi- : 47.4 -53 52.5 9.7e+05 proton : 47.4 -53 52.5 7.95e+05 pi+ : 47.4 -53 52.5 1.45e+06 eta_prime : 47.4 -53 52.5 6.25e+05 eta_prime : 47.4 -53 52.5 2.72e+04 eta_prime : 47.4 -53 52.5 3.59e+04 pi+ : 47.4 -53 52.5 1.15e+04 pi0 : 47.4 -53 52.5 1.24e+04 pi- : 47.4 -53 52.5 34.5 neutron : 47.4 -53 52.5 29 neutron : 47.4 -53 52.5 3.28e+04 pi- : 47.4 -53 52.5 6.75e+04 pi0 : 47.4 -53 52.5 1.61e+04 pi+ : 47.4 -53 52.5 22.6 proton : 47.4 -53 52.5 1e+03 pi0 : 47.4 -53 52.5 3.33e+03 pi0 : 47.4 -53 52.5 2.67e+03 pi0 : 47.4 -53 52.5 1.26e+03 pi- : 47.4 -53 52.5 6.24e+03 proton : 47.4 -53 52.5 1.93e+04 pi+ : 47.4 -53 52.5 1.7e+05 pi- : 47.4 -53 52.5 6.39e+05 eta_prime : 47.4 -53 52.5 1.12e+05 pi+ : 47.4 -53 52.5 7.99e+04 pi- : 47.4 -53 52.5 200 pi0 : 47.4 -53 52.5 73.8 neutron : 47.4 -53 52.5 408 proton : 47.4 -53 52.5 4.27e+05 pi- : 47.4 -53 52.5 5.63e+05 pi+ : 47.4 -53 52.5 2.82e+03 pi- : 47.4 -53 52.5 1.88e+04 proton : 47.4 -53 52.5 596 pi+ : 47.4 -53 52.5 3.52e+03 anti_neutron : 47.4 -53 52.5 2.16e+04 pi0 : 47.4 -53 52.5 7.47e+03 gamma : 47.4 -53 52.5 118 pi- : 47.4 -53 52.5 1.28e+03 pi0 : 47.4 -53 52.5 490 pi+ : 47.4 -53 52.5 1.22e+03 pi- : 47.4 -53 52.5 1.29e+03 neutron : 47.4 -53 52.5 1.19e+03 pi0 : 47.4 -53 52.5 1.92e+04 lambda : 47.4 -53 52.5 1.43e+03 pi0 : 47.4 -53 52.5 1.41e+03 kaon+ : 47.4 -53 52.5 1.3e+04 pi- : 47.4 -53 52.5 2.63e+03 pi0 : 47.4 -53 52.5 1.24e+03 pi+ : 47.4 -53 52.5 1.58e+04 pi0 : 47.4 -53 52.5 1.1e+04 pi+ : 47.4 -53 52.5 7.73e+04 pi- : 47.4 -53 52.5 1.78e+05 pi+ : 47.4 -53 52.5 3.9e+04 pi- : 47.4 -53 52.5 1.87e+04 pi0 : 47.4 -53 52.5 2.99e+04 pi+ : 47.4 -53 52.5 397 proton : 47.4 -53 52.5 138 proton : 47.4 -53 52.5 98.2 triton : 47.4 -53 52.5 109 triton :----------------------------------------------------------------- EndOf2ndaries Info ---------------
Sorry, it seems that I made at least one serious mistake.
Hi, I'll close this for now. If the problem persists, please open a new report. Dennis
