939 – infinite loops

Problem 939 - infinite loops

Summary: infinite loops

Status:	RESOLVED FIXED

Alias:	None

Product:	Examples/Extended
Classification:	Unclassified
Component:	electromagnetic/TestEM5 (show other problems)
Version:	8.2
Hardware:	PC Linux

Importance:	P5 major
Assignee:	Vladimir.Ivantchenko

URL:

Depends on:
Blocks:

Reported:	2007-05-02 13:58 CEST by Pronskikh
Modified:	2009-01-20 19:47 CET (History)
CC List:	1 user (show)

See Also:

Attachments
Add an attachment (proposed patch, testcase, etc.)

Note You need to log in before you can comment on or make changes to this problem.

Description Pronskikh 2007-05-02 13:58:27 CEST

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?

Comment 1 Pronskikh 2007-05-08 14:08:14 CEST

Well, maybe somebody can suggest at least how to
switch off inelastics only for secondary ions, as this
leads to the loops?

Comment 2 Vladimir.Ivantchenko 2009-01-20 19:47:07 CET

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