In addition to the bug reported by me some days ago, the same problem appeares without any changes even in TestEm7 (code gets stuck with infinite loops). In TestEm7 input file set: /addPhysics Livermore, binary_ion, elastic, /gun/ion 6 12 6 /gun/energy 6.0 GeV /testem/gun/rndm 1.5 cm anything else remained as in the initial example. After 3933 events the program gets stuck with the final output: >>PostStepDoIt (process by process): Process Name = msc ++G4Step Information Address of G4Track : 0xa3b9310 Step Length (mm) : 0.5538911774433177 Energy Deposit (MeV) : 2.519489280684367 ----------------------------------------------------------------------- StepPoint Information PreStep PostStep ----------------------------------------------------------------------- Position - x (mm) : 165.3819604203123 165.841243621768 Position - y (mm) : 1.239854125826569 1.443312155913882 Position - z (mm) : 2.815707542980468 3.049070631041495 Global Time (ns) : 2.108259809932339 2.113985357902642 Local Time (ns) : 0.01020952153379051 0.0159350695040934 Proper Time (ns) : 0.009649242706264067 0.01506850000927662 Momentum Direct - x : 0.8291951130232393 0.8291943207775567 Momentum Direct - y : 0.3673254404232669 0.3661178565826357 Momentum Direct - z : 0.4213163720486426 0.4223677230584539 Momentum - x (MeV/c): 794.158857916648 787.6601128810146 Momentum - y (MeV/c): 351.8047172114673 347.779073032261 Momentum - z (MeV/c): 403.5143521623218 401.2113928970887 Total Energy (MeV) : 2968.004696081307 2965.485206800623 Kinetic Energy (MeV): 158.7746960813073 156.255206800623 Velocity (mm/ns) : 96.74015293471828 96.0301208593886 Volume Name : Absorber Absorber Safety (mm) : 34.61803957968769 34.06414919340551 Polarization - x : 0 0 Polarization - y : 0 0 Polarization - Z : 0 0 Weight : 1 1 Step Status : PostStep Proc PostStep Proc Process defined Step: UserMaxStep He3Inelastic ----------------------------------------------------------------------- ++G4ParticleChange Information ----------------------------------------------- G4ParticleChange Information ----------------------------------------------- # of 2ndaries : 0 ----------------------------------------------- Energy Deposit (MeV): 0 Track Status : Alive True Path Length (mm) : 0.554 Stepping Control : 0 Position - x (mm) : 166 Position - y (mm) : 1.44 Position - z (mm) : 3.05 Momentum Direct - x : 0.829 Momentum Direct - y : 0.366 Momentum Direct - z : 0.422 ++List of secondaries generated (x,y,z,kE,t,PID): No. of secodaries = 0 ^^^^^^^^^^^^^^^ here it stops ^^^^^^^^^^^^^^^^^^^^^^^ If this is a model general problem, is it possible to set some cuts (which) to prevent these light ion problems, before the next release ?
One more detail: to reproduce the bug, the length of the absorber should be chosen: testem/det/SetSizeX 40 cm
