Message: Re: endless call of UserSteppingAction Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Re: endless call of UserSteppingAction 

Keywords: UserSteppingAction
Forum: Event and Track Management
Re: None endless call of UserSteppingAction (Konstantin)
Re: None Re: endless call of UserSteppingAction (Makoto Asai)
Date: 09 May, 2006
From: Konstantin <Konstantin>

    I've found the source of the problem reported. In my simulations I used physics
list based on the one from examples/advanced/radioprotection. When process of decay
is being registered there in RemSimDecay.cc, all short-lived particles are excluded:

       if (theDecayProcess -> IsApplicable(*particle) && !particle->IsShortLived())
         {...
         }

So when a rho0 (a short-lived meson) appears during simulation it has no processes
registered ("Transportation" is also not registered for short-lived in this example).
So it is endlessly transported with step ~1E+308.
   Does anybody know what was the reason of such exclusion?


   When I changed the code line above to:

      if (theDecayProcess -> IsApplicable(*particle))

my program started to work OK. But still several times my program has similar
problems: seldom it starts running idle (CPU is busy, no real simulation).
SteppinAction and UserAction are not called during this idle working. Geometry is
very simple (a several tens GeV proton impinges scintillator, divided into strips).
The last lines of tracking info are given below.
   What is the problem?


*********************************************************************************************************
* G4Track Information:   Particle = pi+,   Track ID = 19,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0    -9.84     31.1    -1.65       640        0        0         0 scint_division initStep
    1    -9.92     31.2    -1.09       640    0.428    0.578     0.578 scint_division msc
    2      -10     31.3   -0.519       640    0.455    0.578      1.16 scint_division msc
    3      -10     31.3   -0.511       639   0.0159  0.00885      1.16 scint_division Transportation
    4    -10.1     31.4     0.37       639    0.733      0.9      2.06 scint_division msc
    5    -10.2     31.5     1.25       638    0.785      0.9      2.96 scint_division msc
    6    -10.3     31.6     1.85       637    0.457     0.61      3.57 scint_division hLowEIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  1 ---------------
    :     -10.3      31.6      1.85     0.557                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------


         Thanks in advance,
         Konstantin

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

1 Note: Re: endless call of UserSteppingAction   (Konstantin - 09 May, 2006)
1 None: Re: endless call of UserSteppingAction   (Makoto Asai - 18 May, 2006)
3 None: Re: endless call of UserSteppingAction   (Dennis H. Wright - 19 May, 2006)
5 Question: Re: endless call of UserSteppingAction   (Alexey Lubashevskiy - 08 Nov, 2006)
 Add Message Add Message
to: "Re: endless call of UserSteppingAction"

 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 ]