infinite loop for 50keV EM cuts
running ./bin/corsika -E 1000 -p 11 -N 1 -s 1 --emcut=0.00005 --emthin=0. --wmax=1.e4 --observation-level=8000
we get stuck in an infinite loop between proposal and particleCut:
[corsika:warning (ContinuousProcess.inl:161)] PROPOSAL::getMaxStepLength calculated a negative step length of l=-4.157991619563657 m.
Return 0_m instead.
[corsika:debug (Cascade.inl:193)] transport particle by : 0 m Medium transition after: 0.24954234943841141 m Decay after: inf m Inte
raction after: inf m Continuous limit: 0 m
[corsika:debug (ParticleCut.inl:91)] ParticleCut: checking e- (11), E_kin= 5.1227152589470494e-05 GeV, E=0.0005622271525894704 GeV,
m=0.000511 GeV
[corsika:debug (ObservationPlane.inl:32)] step limit? false
[corsika:debug (Cascade.inl:147)] total_lambda=inf g/cm2, next_interact=inf g/cm2
[corsika:debug (Decay.inl:124)] Pythia::Decay: checking if particle: e- can decay in PYTHIA? false
[corsika:debug (Cascade.inl:158)] total_lifetime=inf ns, next_decay=inf ns
[corsika:debug (TrackingLeapFrogCurved.inl:36)] TrackingLeapfrog_Curved pid: e- , E = 0.0005622271525894704 GeV
Tracking pos: (2646.53 229.165 6.40691e+06) m
Tracking p: (1.04492e-05 0.000167542 -0.000163704) [] GeV
Tracking v: (5.57174e+06 8.93372e+07 -8.72905e+07) m/s
[corsika:debug (TrackingLeapFrogCurved.inl:95)] gyroradius 12.4779 m, steplimit: 0.249542 m = 1.9959e-09 s
[corsika:debug (Intersect.inl:29)] volumeNode=0x17df02e0, numericallyInside=true
[corsika:debug (Intersect.inl:63)] search intersection times with child volume 0x17defe60:
[corsika:debug (Intersect.inl:106)] next time-intersect: 1.9959e-09 s, node 0x17df02e0
[corsika:warning (ContinuousProcess.inl:161)] PROPOSAL::getMaxStepLength calculated a negative step length of l=-4.157991619563657 m.
Return 0_m instead.
[corsika:debug (Cascade.inl:193)] transport particle by : 0 m Medium transition after: 0.24954234943841141 m Decay after: inf m Inte
raction after: inf m Continuous limit: 0 m
probably related to the fix in #563 (closed)