| Summary: | infinite loops | ||
|---|---|---|---|
| Product: | Examples/Extended | Reporter: | Pronskikh <Vitali.Pronskikh> |
| Component: | electromagnetic/TestEM5 | Assignee: | Vladimir.Ivantchenko |
| Status: | RESOLVED FIXED | ||
| Severity: | major | CC: | Gabriele.Cosmo |
| Priority: | P5 | ||
| Version: | 8.2 | ||
| Hardware: | PC | ||
| OS: | Linux | ||
Well, maybe somebody can suggest at least how to switch off inelastics only for secondary ions, as this leads to the loops? I was assigned to this bug report today. Because bug was reported in 2007 I only can say that tit is not show up in recent release, so the problem is closed. VI |
PhysicsList from Hadr01 was adopted to modified TestEm7 example with rather simple geometry: two PMMA blocks and a diamond detector (tally) is irradiated by 500 AMeV C-12 beam (gauss, sigma=15 cm). The code very often gets stuck, cycling without generating new events, not occupying excess memory. Typically, the problem appears with a triton, output looks like this: ************** * G4Track Information: Particle = triton, Track ID = 752, Parent ID = 746 ********************************************************************************************************* Step# X Y Z KineE dEStep StepLeng TrakLeng Volume Process 0 9.23 cm 1.16 cm -1.79 cm 237 MeV 0 eV 0 fm 0 fm Absorber initStep 1 9.33 cm 1.15 cm -1.8 cm 236 MeV 1.04 MeV 1 mm 1 mm Absorber UserMaxStep 2 9.43 cm 1.14 cm -1.82 cm 235 MeV 925 keV 1 mm 2 mm Absorber UserMaxStep 3 9.53 cm 1.13 cm -1.83 cm 234 MeV 1.03 MeV 1 mm 3 mm Absorber UserMaxStep 4 9.63 cm 1.12 cm -1.85 cm 233 MeV 949 keV 1 mm 4 mm Absorber UserMaxStep 5 9.73 cm 1.11 cm -1.86 cm 232 MeV 965 keV 1 mm 5 mm Absorber UserMaxStep 6 9.82 cm 1.1 cm -1.88 cm 231 MeV 951 keV 1 mm 6 mm Absorber UserMaxStep 7 9.92 cm 1.09 cm -1.89 cm 230 MeV 1.02 MeV 1 mm 7 mm Absorber UserMaxStep 8 10 cm 1.07 cm -1.9 cm 229 MeV 1.03 MeV 1 mm 8 mm Absorber UserMaxStep 9 10.1 cm 1.06 cm -1.92 cm 228 MeV 1.08 MeV 1 mm 9 mm Absorber UserMaxStep 10 10.2 cm 1.05 cm -1.93 cm 227 MeV 949 keV 1 mm 1 cm Absorber UserMaxStep 11 10.3 cm 1.04 cm -1.95 cm 226 MeV 1.25 MeV 1 mm 1.1 cm Absorber UserMaxStep 12 10.4 cm 1.03 cm -1.96 cm 225 MeV 947 keV 1 mm 1.2 cm Absorber UserMaxStep 13 10.5 cm 1.02 cm -1.98 cm 224 MeV 1.02 MeV 1 mm 1.3 cm Absorber UserMaxStep 14 10.6 cm 1.01 cm -1.99 cm 223 MeV 1.06 MeV 1 mm 1.4 cm Absorber UserMaxStep 15 10.7 cm 9.96 mm -2 cm 222 MeV 1.39 MeV 1 mm 1.5 cm Absorber UserMaxStep 16 10.8 cm 9.85 mm -2.02 cm 221 MeV 984 keV 1 mm 1.6 cm Absorber UserMaxStep 17 10.9 cm 9.74 mm -2.03 cm 220 MeV 925 keV 1 mm 1.7 cm Absorber UserMaxStep 18 11 cm 9.63 mm -2.05 cm 219 MeV 1.13 MeV 1 mm 1.8 cm Absorber UserMaxStep 19 11.1 cm 9.52 mm -2.06 cm 217 MeV 1.16 MeV 1 mm 1.9 cm Absorber UserMaxStep 20 11.2 cm 9.41 mm -2.08 cm 216 MeV 957 keV 1 mm 2 cm Absorber UserMaxStep 21 11.3 cm 9.31 mm -2.09 cm 215 MeV 1.04 MeV 1 mm 2.1 cm Absorber UserMaxStep 22 11.4 cm 9.2 mm -2.1 cm 214 MeV 1.17 MeV 1 mm 2.2 cm Absorber UserMaxStep 23 11.5 cm 9.09 mm -2.12 cm 213 MeV 1.17 MeV 1 mm 2.3 cm Absorber UserMaxStep 24 11.6 cm 8.98 mm -2.13 cm 212 MeV 1 MeV 1 mm 2.4 cm Absorber UserMaxStep 25 11.7 cm 8.87 mm -2.15 cm 211 MeV 1.1 MeV 1 mm 2.5 cm Absorber UserMaxStep 26 11.8 cm 8.76 mm -2.16 cm 210 MeV 1.39 MeV 1 mm 2.6 cm Absorber UserMaxStep 27 11.9 cm 8.65 mm -2.17 cm 208 MeV 1.27 MeV 1 mm 2.7 cm Absorber UserMaxStep 28 12 cm 8.54 mm -2.19 cm 207 MeV 1.09 MeV 1 mm 2.8 cm Absorber UserMaxStep 29 12.1 cm 8.42 mm -2.2 cm 206 MeV 1.12 MeV 1 mm 2.9 cm Absorber UserMaxStep 30 12.2 cm 8.32 mm -2.22 cm 205 MeV 1.05 MeV 1 mm 3 cm Absorber UserMaxStep 31 12.3 cm 8.2 mm -2.23 cm 204 MeV 1.17 MeV 1 mm 3.1 cm Absorber UserMaxStep 32 12.4 cm 8.1 mm -2.24 cm 203 MeV 1.14 MeV 1 mm 3.2 cm Absorber UserMaxStep 33 12.5 cm 7.99 mm -2.26 cm 202 MeV 1.04 MeV 1 mm 3.3 cm Absorber UserMaxStep 34 12.6 cm 7.88 mm -2.27 cm 200 MeV 1.54 MeV 1 mm 3.4 cm Absorber UserMaxStep 35 12.7 cm 7.76 mm -2.29 cm 199 MeV 1.21 MeV 1 mm 3.5 cm Absorber UserMaxStep 36 12.8 cm 7.65 mm -2.3 cm 198 MeV 1.13 MeV 1 mm 3.6 cm Absorber UserMaxStep 37 12.9 cm 7.53 mm -2.32 cm 197 MeV 1.12 MeV 1 mm 3.7 cm Absorber UserMaxStep 38 13 cm 7.42 mm -2.33 cm 196 MeV 1.13 MeV 1 mm 3.8 cm Absorber UserMaxStep 39 13.1 cm 7.31 mm -2.34 cm 194 MeV 1.2 MeV 1 mm 3.9 cm Absorber UserMaxStep 40 13.2 cm 7.19 mm -2.36 cm 193 MeV 1.05 MeV 1 mm 4 cm Absorber UserMaxStep 41 13.3 cm 7.08 mm -2.37 cm 192 MeV 1.17 MeV 1 mm 4.1 cm Absorber UserMaxStep 42 13.4 cm 6.97 mm -2.39 cm 191 MeV 1.23 MeV 1 mm 4.2 cm Absorber UserMaxStep 43 13.5 cm 6.86 mm -2.4 cm 190 MeV 1.14 MeV 1 mm 4.3 cm Absorber UserMaxStep 44 13.6 cm 6.75 mm -2.41 cm 188 MeV 1.36 MeV 1 mm 4.4 cm Absorber UserMaxStep 45 13.7 cm 6.63 mm -2.43 cm 187 MeV 1.18 MeV 1 mm 4.5 cm Absorber UserMaxStep 46 13.8 cm 6.52 mm -2.44 cm 186 MeV 1.12 MeV 1 mm 4.6 cm Absorber UserMaxStep 47 13.9 cm 6.42 mm -2.46 cm 185 MeV 1.22 MeV 1 mm 4.7 cm Absorber UserMaxStep 48 14 cm 6.31 mm -2.47 cm 184 MeV 1.21 MeV 1 mm 4.8 cm Absorber UserMaxStep 49 14.1 cm 6.2 mm -2.48 cm 182 MeV 1.28 MeV 1 mm 4.9 cm Absorber UserMaxStep 50 14.2 cm 6.09 mm -2.5 cm 181 MeV 1.2 MeV 1 mm 5 cm Absorber UserMaxStep 51 14.3 cm 5.98 mm -2.51 cm 180 MeV 1.46 MeV 1 mm 5.1 cm Absorber UserMaxStep 52 14.3 cm 5.87 mm -2.53 cm 178 MeV 1.28 MeV 1 mm 5.2 cm Absorber UserMaxStep 53 14.4 cm 5.77 mm -2.54 cm 177 MeV 1.2 MeV 1 mm 5.3 cm Absorber UserMaxStep >>PostStepDoIt (process by process): Process Name = msc ++G4Step Information Address of G4Track : 0x9f24180 Step Length (mm) : 0.81559334423204 Energy Deposit (MeV) : 1.067205391906016 ----------------------------------------------------------------------- StepPoint Information PreStep PostStep ----------------------------------------------------------------------- Position - x (mm) : 144.4727990140769 145.2757382247816 Position - y (mm) : 5.76616381863897 5.679069959995605 Position - z (mm) : -25.39861550283754 -25.51216469773394 Global Time (ns) : 5.392123783083421 5.400139271974627 Local Time (ns) : 0.4846353561068781 0.4926508449980842 Proper Time (ns) : 0.4511834613445811 0.4587231446101036 Momentum Direct - x : 0.9844865030845913 0.9845788717671409 Momentum Direct - y : -0.1067860769443532 -0.108452735278546 Momentum Direct - z : -0.1392230548979144 -0.1372678020525772 Momentum - x (MeV/c): 997.9280990160844 994.9216426186085 Momentum - y (MeV/c): -108.2440708253236 -109.5920059061876 Momentum - z (MeV/c): -141.1239240743048 -138.7097682196617 Total Energy (MeV) : 2986.532905246621 2985.465699854715 Kinetic Energy (MeV): 177.2829052466213 176.2156998547152 Velocity (mm/ns) : 101.7519834524555 101.4721782581443 Volume Name : Absorber Absorber Safety (mm) : 29.09525750538061 28.27967634950931 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 tInelastic ----------------------------------------------------------------------- ++G4ParticleChange Information ----------------------------------------------- G4ParticleChange Information ----------------------------------------------- # of 2ndaries : 0 ----------------------------------------------- Energy Deposit (MeV): 0 Track Status : Alive True Path Length (mm) : 0.816 Stepping Control : 0 Position - x (mm) : 145 Position - y (mm) : 5.68 Position - z (mm) : -25.5 Momentum Direct - x : 0.985 Momentum Direct - y : -0.108 Momentum Direct - z : -0.137 ++List of secondaries generated (x,y,z,kE,t,PID): No. of secodaries = 0 Entering the decision point 58.9 1 1 3 1 58.9 58.9 mass table, nucleus, delta : 2.81e+03 2.81e+03 0.328 out radius - nucleus - projectile 2.71 - 0.8 mass table, nucleus, delta : 2.81e+03 2.81e+03 0.328 out radius - nucleus - projectile 2.89 - 0.8 mass table, nucleus, delta : 2.81e+03 2.81e+03 0.328 out radius - nucleus - projectile 2.81 - 0.8 mass table, nucleus, delta : 2.81e+03 2.81e+03 0.328 out radius - nucleus - projectile 2.53 - 0.8 mass table, nucleus, delta : 2.81e+03 2.81e+03 0.328 out radius - nucleus - projectile 3.09 - 0.8 mass table, nucleus, delta : 2.81e+03 2.81e+03 0.328 out radius - nucleus - projectile 2.46 - 0.8 mass table, nucleus, delta : 2.81e+03 2.81e+03 0.328 out radius - nucleus - projectile 3.08 - 0.8 mass table, nucleus, delta : 2.81e+03 2.81e+03 0.328 out radius - nucleus - projectile 2.76 - 0.8 mass table, nucleus, delta : 2.81e+03 2.81e+03 0.328 out radius - nucleus - projectile 2.23 - 0.8 mass table, nucleus, delta : 2.81e+03 2.81e+03 0.328 out radius - nucleus - projectile 2.89 - 0.8 mass table, nucleus, delta : 2.81e+03 2.81e+03 0.328 out radius - nucleus - projectile 2.19 - 0.8 mass table, nucleus, delta : 2.81e+03 2.81e+03 0.328 out radius - nucleus - projectile 2.67 - 0.8 mass table, nucleus, delta : 2.81e+03 2.81e+03 0.328 out radius - nucleus - projectile 2.56 - 0.8 ^^^^^^^^^^^^^^^^ and here it stops ^^^^^^^^^^^^^^^^^^^^ Can you please suggest a solution?