2542 – Particles tracking forever

Problem 2542 - Particles tracking forever

Summary: Particles tracking forever

Status:	ASSIGNED

Alias:	None

Product:	Geant4
Classification:	Unclassified
Component:	physics_lists (show other problems)
Version:	11.1
Hardware:	All All

Importance:	P4 major
Assignee:	Vladimir.Ivantchenko

URL:

Depends on:
Blocks:

Reported:	2023-05-03 16:55 CEST by Simon Spannagel
Modified:	2023-07-09 00:23 CEST (History)
CC List:	3 users (show)

See Also:	2353

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 Simon Spannagel 2023-05-03 16:55:38 CEST

This is a continuation of https://bugzilla-geant4.kek.jp/show_bug.cgi?id=2353 but now with the latest Geant4 11.1.1.

For us the changes introduced in the G4ParticleChange turned things for the worse - no exceptions are thrown when looping endlessly, hence no handle on aborting an event.

Here is the log output with tracking verbosity = 1, if you need higher levels, please let me know. The issue seems to be in the World Transportation after which the particle is outside the universe:


(E) [Geant4] **************************************************************
(E) [Geant4]  Geant4 version Name: geant4-11-01-patch-01 [MT]   (10-February-2023)
(E) [Geant4]   << in Multi-threaded mode >> 
(E) [Geant4]                        Copyright : Geant4 Collaboration
(E) [Geant4]                       References : NIM A 506 (2003), 250-303
(E) [Geant4]                                  : IEEE-TNS 53 (2006), 270-278
(E) [Geant4]                                  : NIM A 835 (2016), 186-225
(E) [Geant4]                              WWW : http://geant4.org/
(E) [Geant4] **************************************************************
(E) [Geant4] 
(E) [Geant4] Checking overlaps for volume wrapper_dut_phys:0 (G4Box) ... OK! 
(E) [Geant4] Checking overlaps for volume sensor_dut_phys:0 (G4Box) ... OK! 
(E) [Geant4] G4PhysListFactory::GetReferencePhysList <QGSP_BERT_EMZ>  EMoption= 4
(E) [Geant4] <<< Geant4 Physics List simulation engine: QGSP_BERT
(E) [Geant4] 
(E) [Geant4] G4VModularPhysicsList::ReplacePhysics: G4EmStandard with type : 2 is replaced with G4EmStandard_opt4
(E) [Geant4] 
(E) [Geant4] 
(E) [Geant4]  hInelastic QGSP_BERT  Thresholds: 
(E) [Geant4]     1) between BERT  and FTF/P over the interval 3 to 6 GeV. 
(E) [Geant4]     2) between FTF/P and QGS/P over the interval 12 to 25 GeV. 
(E) [Geant4]   -- quasiElastic: 1 for QGS  and 0 for FTF
(E) [Geant4] ### Adding tracking cuts for neutron  TimeCut(ns)= 10000  KinEnergyCut(MeV)= 0
(E) [Geant4] 
(E) (E: 1243) [Geant4] G4HadronicProcess::BuildPhysicsTable: for photonNuclear and gamma fail due to undefined pointer to the master process
(E) (E: 1243) [Geant4] 
(E) (E: 1243) [Geant4] *********************************************************************************************************
(E) (E: 1243) [Geant4] * G4Track Information:   Particle = pi+,   Track ID = 1,   Parent ID = 0
(E) (E: 1243) [Geant4] *********************************************************************************************************
(E) (E: 1243) [Geant4] 
(E) (E: 1243) [Geant4] Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
(E) (E: 1243) [Geant4]     0 -0.00217  -0.0204      -10   1.2e+05        0        0         0       World initStep
(E) (E: 1243) [Geant4]     1 -0.00217  -0.0204    -0.05   1.2e+05  9.5e-13     9.95      9.95 sensor_dut_phys Transportation
(E) (E: 1243) [Geant4]     2 -0.00217  -0.0204  -0.0495   1.2e+05        0   0.0005      9.95 sensor_dut_phys StepLimiter
(E) (E: 1243) [Geant4]     3 -0.00217  -0.0204   -0.049   1.2e+05        0   0.0005      9.95 sensor_dut_phys StepLimiter
(E) (E: 1243) [Geant4]     4 -0.00217  -0.0204  -0.0485   1.2e+05        0   0.0005      9.95 sensor_dut_phys StepLimiter

[...]

(E) (E: 1243) [Geant4]   110 -0.00217  -0.0204  0.00401   1.2e+05        0   0.0005        10 sensor_dut_phys StepLimiter
(E) (E: 1243) [Geant4]   111 -0.00217  -0.0204  0.00451   1.2e+05 0.000177   0.0005        10 sensor_dut_phys StepLimiter
(E) (E: 1243) [Geant4]   112 -0.00217  -0.0204  0.00501   1.2e+05 0.000118   0.0005        10 sensor_dut_phys StepLimiter
(E) (E: 1243) [Geant4]   113 -0.00217  -0.0204  0.00502         0        0 7.45e-06        10 sensor_dut_phys pi+Inelastic
(E) (E: 1243) [Geant4] 
(E) (E: 1243) [Geant4] *********************************************************************************************************
(E) (E: 1243) [Geant4] * G4Track Information:   Particle = gamma,   Track ID = 11,   Parent ID = 1
(E) (E: 1243) [Geant4] *********************************************************************************************************
(E) (E: 1243) [Geant4] 
(E) (E: 1243) [Geant4] Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
(E) (E: 1243) [Geant4]     0 -0.00217  -0.0204  0.00502      1.35        0        0         0 sensor_dut_phys initStep
(E) (E: 1243) [Geant4]     1  -0.0665  -0.0196    -0.05      1.35        0   0.0846    0.0846       World Transportation
(E) (E: 1243) [Geant4]     2   -0.528   -0.014   -0.445      1.35        0    0.607     0.692  OutOfWorld Transportation
(E) (E: 1243) [Geant4] 
(E) (E: 1243) [Geant4] *********************************************************************************************************
(E) (E: 1243) [Geant4] * G4Track Information:   Particle = gamma,   Track ID = 10,   Parent ID = 1
(E) (E: 1243) [Geant4] *********************************************************************************************************
(E) (E: 1243) [Geant4] 
(E) (E: 1243) [Geant4] Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
(E) (E: 1243) [Geant4]     0 -0.00217  -0.0204  0.00502      1.12        0        0         0 sensor_dut_phys initStep
(E) (E: 1243) [Geant4]     1  -0.0547  -0.0118     0.05      1.12        0   0.0697    0.0697       World Transportation
(E) (E: 1243) [Geant4]     2   -0.528   0.0657    0.455      1.12        0    0.628     0.697  OutOfWorld Transportation
(E) (E: 1243) [Geant4] 
(E) (E: 1243) [Geant4] *********************************************************************************************************
(E) (E: 1243) [Geant4] * G4Track Information:   Particle = gamma,   Track ID = 9,   Parent ID = 1
(E) (E: 1243) [Geant4] *********************************************************************************************************
(E) (E: 1243) [Geant4] 
(E) (E: 1243) [Geant4] Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
(E) (E: 1243) [Geant4]     0 -0.00217  -0.0204  0.00502      1.59        0        0         0 sensor_dut_phys initStep
(E) (E: 1243) [Geant4]     1   0.0382   0.0063     0.05      1.59        0   0.0661    0.0661       World Transportation
(E) (E: 1243) [Geant4]     2    0.528     0.33    0.596      1.59        0    0.802     0.868  OutOfWorld Transportation
(E) (E: 1243) [Geant4] 
(E) (E: 1243) [Geant4] *********************************************************************************************************
(E) (E: 1243) [Geant4] * G4Track Information:   Particle = neutron,   Track ID = 8,   Parent ID = 1
(E) (E: 1243) [Geant4] *********************************************************************************************************
(E) (E: 1243) [Geant4] 
(E) (E: 1243) [Geant4] Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
(E) (E: 1243) [Geant4]     0 -0.00217  -0.0204  0.00502     0.561        0        0         0 sensor_dut_phys initStep
(E) (E: 1243) [Geant4]     1    0.037    0.023     0.05     0.561        0   0.0738    0.0738       World Transportation
(E) (E: 1243) [Geant4]     2    0.493    0.528    0.574     0.561        0    0.858     0.932  OutOfWorld Transportation
(E) (E: 1243) [Geant4] 
(E) (E: 1243) [Geant4] *********************************************************************************************************
(E) (E: 1243) [Geant4] * G4Track Information:   Particle = alpha,   Track ID = 7,   Parent ID = 1
(E) (E: 1243) [Geant4] *********************************************************************************************************
(E) (E: 1243) [Geant4] 
(E) (E: 1243) [Geant4] Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
(E) (E: 1243) [Geant4]     0 -0.00217  -0.0204  0.00502      11.5        0        0         0 sensor_dut_phys initStep
(E) (E: 1243) [Geant4]     1 -0.00196  -0.0203  0.00458      11.5   0.0294   0.0005    0.0005 sensor_dut_phys StepLimiter

[...]

(E) (E: 1243) [Geant4]   123   0.0236 -0.00719  -0.0492      7.46   0.0329   0.0005    0.0615 sensor_dut_phys StepLimiter
(E) (E: 1243) [Geant4]   124   0.0239 -0.00709  -0.0497      7.42   0.0363   0.0005     0.062 sensor_dut_phys StepLimiter
(E) (E: 1243) [Geant4]   125    0.024   -0.007    -0.05       7.4   0.0231 0.000394    0.0624       World Transportation
(E) (E: 1243) [Geant4]   126    0.528    0.245    -1.07       7.4 5.61e-05     1.16      1.22  OutOfWorld Transportation
(E) (E: 1243) [Geant4] 
(E) (E: 1243) [Geant4] *********************************************************************************************************
(E) (E: 1243) [Geant4] * G4Track Information:   Particle = proton,   Track ID = 6,   Parent ID = 1
(E) (E: 1243) [Geant4] *********************************************************************************************************
(E) (E: 1243) [Geant4] 
(E) (E: 1243) [Geant4] Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
(E) (E: 1243) [Geant4]     0 -0.00217  -0.0204  0.00502      5.95        0        0         0 sensor_dut_phys initStep
(E) (E: 1243) [Geant4]     1 -0.00194  -0.0202   0.0046      5.94  0.00632   0.0005    0.0005 sensor_dut_phys StepLimiter

[...]

(E) (E: 1243) [Geant4]   131   0.0272 0.000847  -0.0495       5.3  0.00717   0.0005    0.0655 sensor_dut_phys StepLimiter
(E) (E: 1243) [Geant4]   132   0.0274  0.00101    -0.05      5.29  0.00501   0.0005     0.066 sensor_dut_phys StepLimiter
(E) (E: 1243) [Geant4]   133   0.0274  0.00103    -0.05      5.29 0.000375 4.56e-05     0.066       World Transportation
(E) (E: 1243) [Geant4]   134 -2.84e+21 -2.13e+21 5.47e+21      5.29        0 -6.52e+21 -6.52e+21       World hIoni
(E) (E: 1243) [Geant4]   135 -5.67e+21 -4.26e+21 1.09e+22      5.29        0 -6.52e+21  -1.3e+22       World hIoni
(E) (E: 1243) [Geant4]   136 -8.51e+21 -6.38e+21 1.64e+22      5.29        0 -6.52e+21 -1.96e+22       World hIoni
(E) (E: 1243) [Geant4]   137 -1.13e+22 -8.51e+21 2.19e+22      5.29        0 -6.52e+21 -2.61e+22       World hIoni
(E) (E: 1243) [Geant4]   138 -1.42e+22 -1.06e+22 2.73e+22      5.29        0 -6.52e+21 -3.26e+22       World hIoni
(E) (E: 1243) [Geant4]   139 -1.7e+22 -1.28e+22 3.28e+22      5.29        0 -6.52e+21 -3.91e+22       World hIoni
(E) (E: 1243) [Geant4]   140 -1.99e+22 -1.49e+22 3.83e+22      5.29        0 -6.52e+21 -4.56e+22       World hIoni
(E) (E: 1243) [Geant4]   141 -2.27e+22 -1.7e+22 4.37e+22      5.29        0 -6.52e+21 -5.21e+22       World hIoni
(E) (E: 1243) [Geant4]   142 -2.55e+22 -1.92e+22 4.92e+22      5.29        0 -6.52e+21 -5.87e+22       World hIoni
(E) (E: 1243) [Geant4]   143 -2.84e+22 -2.13e+22 5.47e+22      5.29        0 -6.52e+21 -6.52e+22       World hIoni
(E) (E: 1243) [Geant4]   144 -3.12e+22 -2.34e+22 6.01e+22      5.29        0 -6.52e+21 -7.17e+22       World hIoni
(E) (E: 1243) [Geant4]   145 -3.4e+22 -2.55e+22 6.56e+22      5.29        0 -6.52e+21 -7.82e+22       World hIoni
(E) (E: 1243) [Geant4]   146 -3.69e+22 -2.77e+22 7.11e+22      5.29        0 -6.52e+21 -8.47e+22       World hIoni

[...]

(E) (E: 1243) [Geant4]  9358 -2.62e+25 -1.96e+25 5.04e+25      5.29        0 -6.52e+21 -6.01e+25       World hIoni
(E) (E: 1243) [Geant4]  9359 -2.62e+25 -1.96e+25 5.04e+25      5.29        0 -6.52e+21 -6.01e+25       World hIoni
(E) (E: 1243) [Geant4]  9360 -2.62e+25 -1.96e+25 5.05e+25      5.29        0 -6.52e+21 -6.01e+25       World hIoni

Comment 1 Vladimir.Ivantchenko 2023-07-09 00:23:52 CEST

Hello Simon,

I am sorry for this bug. It should be fixed in two new public versions of G4: 11.1.p02 and 11.2beta. Would it be possible to check?

Vladimir