Message: Same two steps repeat again and again Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Same two steps repeat again and again 

Forum: Event and Track Management
Date: 31 Aug, 2009
From: Ioannis Sechopoulos <Ioannis Sechopoulos>

A simulation I have been using for a long time now is giving me problems in that after a certain number of events (varies depending on the random number seed), it gets stuck performing the same two steps. The output of "/tracking/verbose 4" for these two steps is the following:

#Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    7     84.1     -182     44.4     0.016        0 4.64e-15       660       World Transportation

 >>AlongStepDoIt (process by process):    Process Name = Transportation

    ++G4Step Information 
      Address of G4Track    : 0x19666c0
      Step Length (mm)      : 0
      Energy Deposit (MeV)  : 0
      -----------------------------------------------------------------------
        StepPoint Information               PreStep            PostStep
      -----------------------------------------------------------------------
         Position - x (mm)   :    84.08523386931844   84.08523386931844
         Position - y (mm)   :     -181.75378512735    -181.75378512735
         Position - z (mm)   :     44.4354460094737    44.4354460094737
         Global Time (ns)    :    2.201048285186674   2.201048285186674
         Local Time (ns)     :    2.201048285186674   2.201048285186674
         Proper Time (ns)    :                    0                   0
         Momentum Direct - x :   0.4601418205088121  0.4601418205088121
         Momentum Direct - y :   0.7169432097626749  0.7169432097626749
         Momentum Direct - z :   0.5237002377257712  0.5237002377257712
         Momentum - x (MeV/c): 0.0073622691281409940.007362269128140994
         Momentum - y (MeV/c):   0.0114710913562028  0.0114710913562028
         Momentum - z (MeV/c):  0.00837920380361234 0.00837920380361234
         Total Energy (MeV)  :                0.016               0.016
         Kinetic Energy (MeV):                0.016               0.016
         Velocity (mm/ns)    :           299.792458          299.792458
         Volume Name         :                World               World
         Safety (mm)         :                5e-10               5e-10
         Polarization - x    :                    0                   0
         Polarization - y    :                    0                   0
         Polarization - Z    :                    0                   0
         Weight              :                    1                   1
         Step Status         :           Geom Limit          Geom Limit
         Process defined Step:       Transportation      Transportation
      -----------------------------------------------------------------------
          !Note! Safety of PostStep is only valid after all DoIt invocations.

    ++G4ParticleChange Information 
      -----------------------------------------------
        G4ParticleChange Information  
      -----------------------------------------------
        # of 2ndaries       :                    0
      -----------------------------------------------
        Energy Deposit (MeV):                    0
        Non-ionizing Energy Deposit (MeV):                    0
        Track Status        :                Alive
        True Path Length (mm) :                    0
        Stepping Control      :                    0
        Mass (GeV)   :                    0
        Charge (eplus)   :                    0
        MagneticMoment   :                    0
                :  =                    0*[e hbar]/[2 m]
        Position - x (mm)   :                 84.1
        Position - y (mm)   :                 -182
        Position - z (mm)   :                 44.4
        Time (ns)           :                  2.2
        Proper Time (ns)    :                    0
        Momentum Direct - x :                 0.46
        Momentum Direct - y :                0.717
        Momentum Direct - z :                0.524
        Kinetic Energy (MeV):                0.016
        Polarization - x    :                    0
        Polarization - y    :                    0
        Polarization - z    :                    0
        Touchable (pointer) :             0xc8ac20

    ++List of secondaries generated (x,y,z,kE,t,PID):  No. of secodaries = 0

 >>AlongStepDoIt (after all invocations):
    ++List of invoked processes 
      1) Transportation

    ++G4Step Information 
      Address of G4Track    : 0x19666c0
      Step Length (mm)      : 0
      Energy Deposit (MeV)  : 0
      -----------------------------------------------------------------------
        StepPoint Information               PreStep            PostStep
      -----------------------------------------------------------------------
         Position - x (mm)   :    84.08523386931844   84.08523386931844
         Position - y (mm)   :     -181.75378512735    -181.75378512735
         Position - z (mm)   :     44.4354460094737    44.4354460094737
         Global Time (ns)    :    2.201048285186674   2.201048285186674
         Local Time (ns)     :    2.201048285186674   2.201048285186674
         Proper Time (ns)    :                    0                   0
         Momentum Direct - x :   0.4601418205088121  0.4601418205088121
         Momentum Direct - y :   0.7169432097626749  0.7169432097626749
         Momentum Direct - z :   0.5237002377257712  0.5237002377257712
         Momentum - x (MeV/c): 0.0073622691281409940.007362269128140994
         Momentum - y (MeV/c):   0.0114710913562028  0.0114710913562028
         Momentum - z (MeV/c):  0.00837920380361234 0.00837920380361234
         Total Energy (MeV)  :                0.016               0.016
         Kinetic Energy (MeV):                0.016               0.016
         Velocity (mm/ns)    :           299.792458          299.792458
         Volume Name         :                World               World
         Safety (mm)         :                5e-10               5e-10
         Polarization - x    :                    0                   0
         Polarization - y    :                    0                   0
         Polarization - Z    :                    0                   0
         Weight              :                    1                   1
         Step Status         :           Geom Limit          Geom Limit
         Process defined Step:       Transportation      Transportation
      -----------------------------------------------------------------------

    ++List of secondaries generated (x,y,z,kE,t,PID):  No. of secodaries = 0

 >>PostStepDoIt (process by process):    Process Name = Transportation

    ++G4Step Information 
      Address of G4Track    : 0x19666c0
      Step Length (mm)      : 0
      Energy Deposit (MeV)  : 0
      -----------------------------------------------------------------------
        StepPoint Information               PreStep            PostStep
      -----------------------------------------------------------------------
         Position - x (mm)   :    84.08523386931844   84.08523386931844
         Position - y (mm)   :     -181.75378512735    -181.75378512735
         Position - z (mm)   :     44.4354460094737    44.4354460094737
         Global Time (ns)    :    2.201048285186674   2.201048285186674
         Local Time (ns)     :    2.201048285186674   2.201048285186674
         Proper Time (ns)    :                    0                   0
         Momentum Direct - x :   0.4601418205088121  0.4601418205088121
         Momentum Direct - y :   0.7169432097626749  0.7169432097626749
         Momentum Direct - z :   0.5237002377257712  0.5237002377257712
         Momentum - x (MeV/c): 0.0073622691281409940.007362269128140994
         Momentum - y (MeV/c):   0.0114710913562028  0.0114710913562028
         Momentum - z (MeV/c):  0.00837920380361234 0.00837920380361234
         Total Energy (MeV)  :                0.016               0.016
         Kinetic Energy (MeV):                0.016               0.016
         Velocity (mm/ns)    :           299.792458          299.792458
         Volume Name         :                World           physiSkin
         Safety (mm)         :                5e-10               5e-10
         Polarization - x    :                    0                   0
         Polarization - y    :                    0                   0
         Polarization - Z    :                    0                   0
         Weight              :                    1                   1
         Step Status         :           Geom Limit          Geom Limit
         Process defined Step:       Transportation      Transportation
      -----------------------------------------------------------------------

    ++G4ParticleChange Information 
      -----------------------------------------------
        G4ParticleChange Information  
      -----------------------------------------------
        # of 2ndaries       :                    0
      -----------------------------------------------
        Energy Deposit (MeV):                    0
        Non-ionizing Energy Deposit (MeV):                    0
        Track Status        :                Alive
        True Path Length (mm) :                    0
        Stepping Control      :                    0
        Mass (GeV)   :                    0
        Charge (eplus)   :                    0
        MagneticMoment   :                    0
                :  =                    0*[e hbar]/[2 m]
        Position - x (mm)   :                 84.1
        Position - y (mm)   :                 -182
        Position - z (mm)   :                 44.4
        Time (ns)           :                  2.2
        Proper Time (ns)    :                    0
        Momentum Direct - x :                 0.46
        Momentum Direct - y :                0.717
        Momentum Direct - z :                0.524
        Kinetic Energy (MeV):                0.016
        Polarization - x    :                    0
        Polarization - y    :                    0
        Polarization - z    :                    0
        Touchable (pointer) :             0xc8ab70

    ++List of secondaries generated (x,y,z,kE,t,PID):  No. of secodaries = 0

 **PostStepDoIt (after all invocations):
    ++List of invoked processes 
      1) Transportation
      2) LowEnPhotoElec (Forced)

    ++G4Step Information 
      Address of G4Track    : 0x19666c0
      Step Length (mm)      : 0
      Energy Deposit (MeV)  : 0
      -----------------------------------------------------------------------
        StepPoint Information               PreStep            PostStep
      -----------------------------------------------------------------------
         Position - x (mm)   :    84.08523386931844   84.08523386931844
         Position - y (mm)   :     -181.75378512735    -181.75378512735
         Position - z (mm)   :     44.4354460094737    44.4354460094737
         Global Time (ns)    :    2.201048285186674   2.201048285186674
         Local Time (ns)     :    2.201048285186674   2.201048285186674
         Proper Time (ns)    :                    0                   0
         Momentum Direct - x :   0.4601418205088121  0.4601418205088121
         Momentum Direct - y :   0.7169432097626749  0.7169432097626749
         Momentum Direct - z :   0.5237002377257712  0.5237002377257712
         Momentum - x (MeV/c): 0.0073622691281409940.007362269128140994
         Momentum - y (MeV/c):   0.0114710913562028  0.0114710913562028
         Momentum - z (MeV/c):  0.00837920380361234 0.00837920380361234
         Total Energy (MeV)  :                0.016               0.016
         Kinetic Energy (MeV):                0.016               0.016
         Velocity (mm/ns)    :           299.792458          299.792458
         Volume Name         :                World           physiSkin
         Safety (mm)         :                5e-10               5e-10
         Polarization - x    :                    0                   0
         Polarization - y    :                    0                   0
         Polarization - Z    :                    0                   0
         Weight              :                    1                   1
         Step Status         :           Geom Limit          Geom Limit
         Process defined Step:       Transportation      Transportation
      -----------------------------------------------------------------------

    ++List of secondaries generated (x,y,z,kE,t,PID):  No. of secodaries = 0
      [Note]Secondaries from AlongStepDoIt included.

    ++G4Track Information 
      -----------------------------------------------
        G4Track Information  
      -----------------------------------------------
        Step number         :                    8
        Position - x (mm)   :                 84.1
        Position - y (mm)   :                 -182
        Position - z (mm)   :                 44.4
        Global Time (ns)    :                  2.2
        Local Time (ns)     :                  2.2
        Momentum Direct - x :                 0.46
        Momentum Direct - y :                0.717
        Momentum Direct - z :                0.524
        Kinetic Energy (MeV):                0.016
        Polarization - x    :                    0
        Polarization - y    :                    0
        Polarization - z    :                    0
        Track Length        :                  660
        Track ID #          :                    1
        Parent Track ID #   :                    0
        Next Volume         :            physiSkin 
        Track Status        :                Alive
        Vertex - x (mm)     :                    0
        Vertex - y (mm)     :                 -236
        Vertex - z (mm)     :                  694
        Vertex - Px (MomDir):                0.126
        Vertex - Py (MomDir):                0.081
        Vertex - Pz (MomDir):               -0.989
        Vertex - KineE (MeV):                0.016
        Creator Process     :      Event Generator
      -----------------------------------------------

#Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    8     84.1     -182     44.4     0.016        0        0       660   physiSkin Transportation

 >>AlongStepDoIt (process by process):    Process Name = Transportation

    ++G4Step Information 
      Address of G4Track    : 0x19666c0
      Step Length (mm)      : 4.643365327944604e-15
      Energy Deposit (MeV)  : 0
      -----------------------------------------------------------------------
        StepPoint Information               PreStep            PostStep
      -----------------------------------------------------------------------
         Position - x (mm)   :    84.08523386931844   84.08523386931844
         Position - y (mm)   :     -181.75378512735    -181.75378512735
         Position - z (mm)   :     44.4354460094737    44.4354460094737
         Global Time (ns)    :    2.201048285186674   2.201048285186674
         Local Time (ns)     :    2.201048285186674   2.201048285186674
         Proper Time (ns)    :                    0                   0
         Momentum Direct - x :   0.4601418205088121  0.4601418205088121
         Momentum Direct - y :   0.7169432097626749  0.7169432097626749
         Momentum Direct - z :   0.5237002377257712  0.5237002377257712
         Momentum - x (MeV/c): 0.0073622691281409940.007362269128140994
         Momentum - y (MeV/c):   0.0114710913562028  0.0114710913562028
         Momentum - z (MeV/c):  0.00837920380361234 0.00837920380361234
         Total Energy (MeV)  :                0.016               0.016
         Kinetic Energy (MeV):                0.016               0.016
         Velocity (mm/ns)    :           299.792458          299.792458
         Volume Name         :            physiSkin           physiSkin
         Safety (mm)         :                5e-10               5e-10
         Polarization - x    :                    0                   0
         Polarization - y    :                    0                   0
         Polarization - Z    :                    0                   0
         Weight              :                    1                   1
         Step Status         :           Geom Limit          Geom Limit
         Process defined Step:       Transportation      Transportation
      -----------------------------------------------------------------------
          !Note! Safety of PostStep is only valid after all DoIt invocations.

    ++G4ParticleChange Information 
      -----------------------------------------------
        G4ParticleChange Information  
      -----------------------------------------------
        # of 2ndaries       :                    0
      -----------------------------------------------
        Energy Deposit (MeV):                    0
        Non-ionizing Energy Deposit (MeV):                    0
        Track Status        :                Alive
        True Path Length (mm) :             4.64e-15
        Stepping Control      :                    0
        Mass (GeV)   :                    0
        Charge (eplus)   :                    0
        MagneticMoment   :                    0
                :  =                    0*[e hbar]/[2 m]
        Position - x (mm)   :                 84.1
        Position - y (mm)   :                 -182
        Position - z (mm)   :                 44.4
        Time (ns)           :                  2.2
        Proper Time (ns)    :                    0
        Momentum Direct - x :                 0.46
        Momentum Direct - y :                0.717
        Momentum Direct - z :                0.524
        Kinetic Energy (MeV):                0.016
        Polarization - x    :                    0
        Polarization - y    :                    0
        Polarization - z    :                    0
        Touchable (pointer) :             0xc8ab70

    ++List of secondaries generated (x,y,z,kE,t,PID):  No. of secodaries = 0

 >>AlongStepDoIt (after all invocations):
    ++List of invoked processes 
      1) Transportation

    ++G4Step Information 
      Address of G4Track    : 0x19666c0
      Step Length (mm)      : 4.643365327944604e-15
      Energy Deposit (MeV)  : 0
      -----------------------------------------------------------------------
        StepPoint Information               PreStep            PostStep
      -----------------------------------------------------------------------
         Position - x (mm)   :    84.08523386931844   84.08523386931844
         Position - y (mm)   :     -181.75378512735    -181.75378512735
         Position - z (mm)   :     44.4354460094737    44.4354460094737
         Global Time (ns)    :    2.201048285186674   2.201048285186674
         Local Time (ns)     :    2.201048285186674   2.201048285186674
         Proper Time (ns)    :                    0                   0
         Momentum Direct - x :   0.4601418205088121  0.4601418205088121
         Momentum Direct - y :   0.7169432097626749  0.7169432097626749
         Momentum Direct - z :   0.5237002377257712  0.5237002377257712
         Momentum - x (MeV/c): 0.0073622691281409940.007362269128140994
         Momentum - y (MeV/c):   0.0114710913562028  0.0114710913562028
         Momentum - z (MeV/c):  0.00837920380361234 0.00837920380361234
         Total Energy (MeV)  :                0.016               0.016
         Kinetic Energy (MeV):                0.016               0.016
         Velocity (mm/ns)    :           299.792458          299.792458
         Volume Name         :            physiSkin           physiSkin
         Safety (mm)         :                5e-10               5e-10
         Polarization - x    :                    0                   0
         Polarization - y    :                    0                   0
         Polarization - Z    :                    0                   0
         Weight              :                    1                   1
         Step Status         :           Geom Limit          Geom Limit
         Process defined Step:       Transportation      Transportation
      -----------------------------------------------------------------------

    ++List of secondaries generated (x,y,z,kE,t,PID):  No. of secodaries = 0

 >>PostStepDoIt (process by process):    Process Name = Transportation

    ++G4Step Information 
      Address of G4Track    : 0x19666c0
      Step Length (mm)      : 4.643365327944604e-15
      Energy Deposit (MeV)  : 0
      -----------------------------------------------------------------------
        StepPoint Information               PreStep            PostStep
      -----------------------------------------------------------------------
         Position - x (mm)   :    84.08523386931844   84.08523386931844
         Position - y (mm)   :     -181.75378512735    -181.75378512735
         Position - z (mm)   :     44.4354460094737    44.4354460094737
         Global Time (ns)    :    2.201048285186674   2.201048285186674
         Local Time (ns)     :    2.201048285186674   2.201048285186674
         Proper Time (ns)    :                    0                   0
         Momentum Direct - x :   0.4601418205088121  0.4601418205088121
         Momentum Direct - y :   0.7169432097626749  0.7169432097626749
         Momentum Direct - z :   0.5237002377257712  0.5237002377257712
         Momentum - x (MeV/c): 0.0073622691281409940.007362269128140994
         Momentum - y (MeV/c):   0.0114710913562028  0.0114710913562028
         Momentum - z (MeV/c):  0.00837920380361234 0.00837920380361234
         Total Energy (MeV)  :                0.016               0.016
         Kinetic Energy (MeV):                0.016               0.016
         Velocity (mm/ns)    :           299.792458          299.792458
         Volume Name         :            physiSkin               World
         Safety (mm)         :                5e-10               5e-10
         Polarization - x    :                    0                   0
         Polarization - y    :                    0                   0
         Polarization - Z    :                    0                   0
         Weight              :                    1                   1
         Step Status         :           Geom Limit          Geom Limit
         Process defined Step:       Transportation      Transportation
      -----------------------------------------------------------------------

    ++G4ParticleChange Information 
      -----------------------------------------------
        G4ParticleChange Information  
      -----------------------------------------------
        # of 2ndaries       :                    0
      -----------------------------------------------
        Energy Deposit (MeV):                    0
        Non-ionizing Energy Deposit (MeV):                    0
        Track Status        :                Alive
        True Path Length (mm) :             4.64e-15
        Stepping Control      :                    0
        Mass (GeV)   :                    0
        Charge (eplus)   :                    0
        MagneticMoment   :                    0
                :  =                    0*[e hbar]/[2 m]
        Position - x (mm)   :                 84.1
        Position - y (mm)   :                 -182
        Position - z (mm)   :                 44.4
        Time (ns)           :                  2.2
        Proper Time (ns)    :                    0
        Momentum Direct - x :                 0.46
        Momentum Direct - y :                0.717
        Momentum Direct - z :                0.524
        Kinetic Energy (MeV):                0.016
        Polarization - x    :                    0
        Polarization - y    :                    0
        Polarization - z    :                    0
        Touchable (pointer) :             0xc8ac20

    ++List of secondaries generated (x,y,z,kE,t,PID):  No. of secodaries = 0

 **PostStepDoIt (after all invocations):
    ++List of invoked processes 
      1) Transportation
      2) LowEnPhotoElec (Forced)

    ++G4Step Information 
      Address of G4Track    : 0x19666c0
      Step Length (mm)      : 4.643365327944604e-15
      Energy Deposit (MeV)  : 0
      -----------------------------------------------------------------------
        StepPoint Information               PreStep            PostStep
      -----------------------------------------------------------------------
         Position - x (mm)   :    84.08523386931844   84.08523386931844
         Position - y (mm)   :     -181.75378512735    -181.75378512735
         Position - z (mm)   :     44.4354460094737    44.4354460094737
         Global Time (ns)    :    2.201048285186674   2.201048285186674
         Local Time (ns)     :    2.201048285186674   2.201048285186674
         Proper Time (ns)    :                    0                   0
         Momentum Direct - x :   0.4601418205088121  0.4601418205088121
         Momentum Direct - y :   0.7169432097626749  0.7169432097626749
         Momentum Direct - z :   0.5237002377257712  0.5237002377257712
         Momentum - x (MeV/c): 0.0073622691281409940.007362269128140994
         Momentum - y (MeV/c):   0.0114710913562028  0.0114710913562028
         Momentum - z (MeV/c):  0.00837920380361234 0.00837920380361234
         Total Energy (MeV)  :                0.016               0.016
         Kinetic Energy (MeV):                0.016               0.016
         Velocity (mm/ns)    :           299.792458          299.792458
         Volume Name         :            physiSkin               World
         Safety (mm)         :                5e-10               5e-10
         Polarization - x    :                    0                   0
         Polarization - y    :                    0                   0
         Polarization - Z    :                    0                   0
         Weight              :                    1                   1
         Step Status         :           Geom Limit          Geom Limit
         Process defined Step:       Transportation      Transportation
      -----------------------------------------------------------------------

    ++List of secondaries generated (x,y,z,kE,t,PID):  No. of secodaries = 0
      [Note]Secondaries from AlongStepDoIt included.

    ++G4Track Information 
      -----------------------------------------------
        G4Track Information  
      -----------------------------------------------
        Step number         :                    9
        Position - x (mm)   :                 84.1
        Position - y (mm)   :                 -182
        Position - z (mm)   :                 44.4
        Global Time (ns)    :                  2.2
        Local Time (ns)     :                  2.2
        Momentum Direct - x :                 0.46
        Momentum Direct - y :                0.717
        Momentum Direct - z :                0.524
        Kinetic Energy (MeV):                0.016
        Polarization - x    :                    0
        Polarization - y    :                    0
        Polarization - z    :                    0
        Track Length        :                  660
        Track ID #          :                    1
        Parent Track ID #   :                    0
        Next Volume         :                World 
        Track Status        :                Alive
        Vertex - x (mm)     :                    0
        Vertex - y (mm)     :                 -236
        Vertex - z (mm)     :                  694
        Vertex - Px (MomDir):                0.126
        Vertex - Py (MomDir):                0.081
        Vertex - Pz (MomDir):               -0.989
        Vertex - KineE (MeV):                0.016
        Creator Process     :      Event Generator
      -----------------------------------------------

(sorry for the length of the output, but I don't know what might be an important part of the output). The next steps are repeats of these two, with everything exactly the same but the step numbers. To the best of my knowledge, there are no geometry overlaps (I checked with DAVID) and I have been using this simulation geometry for a long time. The only thing that changed recently is that I upgraded to the new patch 2 of 9.2. Can anybody give me an idea of what I could try to avoid/fix this or let me know if you need me to try anything different and provide the output here?

Thanks! Ioannis

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 Question: Re: Same two steps repeat again and again   (Gabriele Cosmo - 01 Sep, 2009)
(_ Question: Re: Same two steps repeat again and again   (Ioannis Sechopoulos - 01 Sep, 2009)
(_ More: Re: Same two steps repeat again and again   (Gabriele Cosmo - 01 Sep, 2009)
(_ More: Re: Same two steps repeat again and again   (Ioannis Sechopoulos - 01 Sep, 2009)
(_ More: Re: Same two steps repeat again and again   (Gabriele Cosmo - 02 Sep, 2009)
(_ None: Re: Same two steps repeat again and again   (Robert Penny - 30 Sep, 2009)
(_ More: Re: Same two steps repeat again and again   (Gabriele Cosmo - 30 Sep, 2009)
(_ Question: Re: Same two steps repeat again and again   (Ron Mathis - 02 Oct, 2009)
(_ Feedback: Re: Same two steps repeat again and again   (Gabriele Cosmo - 05 Oct, 2009)
 Add Message Add Message
to: "Same two steps repeat again and again"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]