Problem 1728

Summary: particle erroneously propagated to the mother volume
Product: Geant4 Reporter: vitali choutko <vitali.choutko>
Component: geometry/navigationAssignee: John Apostolakis <John.Apostolakis>
Status: RESOLVED WORKSFORME    
Severity: trivial CC: asai, Gabriele.Cosmo, John.Apostolakis, Vladimir.Ivantchenko
Priority: P5    
Version: 10.1   
Hardware: All   
OS: All   

Description vitali choutko 2015-04-01 22:19:57 CEST
The particle erroneously propagated to the mother volume in a current step if

in the previous step particle exited from daughter  (G4Navigator fExiting=true)

and

in the current step G4VEnergyLossProcess::PostStepGetPhysicalInteractionLength()  returned 0   [ Most probably any process with PostStepGetPhysicalInteractionLengh() returned 0 will lead to same behaviour]



in that case

G4Transportation::AlongStepGetPhysicalInteractionLength  will not clear G4Navigator::fExiting(because of currentMinimumStep==0)  and will set  fGeometryLimitedStep=true inspite currentMinimumStep=0  /because it is on boundary(?) ,ie currentSafety==0/


G4Transportation::PostStepDoIt will call  G4Navigator::LocateGlobalPointAndSetup  which  will reduce fHistory depth /as fExiting=true/   and erroneously move particle to the mother volume.


Then one step in mother volume will be done , which  lead to miss corresponding underlying volumes.  After this step navigation resules in a normal way/




What is affected:  Any charged particles but in particularly ions with high energy (>TeV/nuc) , which naturally emit a lot of delta rays


Versions >=  4.9.6 including 4.10.01.p01 are affected.


Versions <= 4.9.4 were not affected because  G4VEnergyLossProcess::PostStepGetPhysicalInteractionLength()  never returned 0.

So trivial workarond would be just reinstate this behaviour.

The fix may be as simple as to clear fExiting in G4Transportation::AlongStepGetPhysicalInteractionLength even for  currentMinimumStep==0, however I don;t think I am qualified enough to say more.
Comment 1 Gabriele Cosmo 2015-04-02 10:01:05 CEST
Thanks for reporting this. We'll verify if this is indeed the case.
Comment 2 Gabriele Cosmo 2015-12-10 12:43:41 CET
We have been investigating this issue for quite some while but were not able to reproduce the reported problem. Can you please confirm you see the problem also in 10.1.p02 or the latest release 10.2 ?
If so, could you please provide a test case ?
Comment 3 Gabriele Cosmo 2016-12-19 11:54:08 CET
Closing issue, as not able to reproduce.