############################################ !!! WARNING - FPE detection is activated !!! ############################################ ################################ !!! G4Backtrace is activated !!! ################################ ************************************************************** Geant4 version Name: geant4-11-04-patch-01 (13-March-2026) Copyright : Geant4 Collaboration References : NIM A 506 (2003), 250-303 : IEEE-TNS 53 (2006), 270-278 : NIM A 835 (2016), 186-225 WWW : http://geant4.org/ ************************************************************** ***** Table : Nb of materials = 16 ***** Material: Air density: 1.290 mg/cm3 RadL: 285.161 m Nucl.Int.Length: 662.680 m Imean: 85.703 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Nitrogen (N) Z = 7.0 N = 14 A = 14.010 g/mole ---> Isotope: N14 Z = 7 N = 14 A = 14.00 g/mole abundance: 99.632 % ---> Isotope: N15 Z = 7 N = 15 A = 15.00 g/mole abundance: 0.368 % ElmMassFraction: 70.00 % ElmAbundance 72.71 % ---> Element: Oxygen (O) Z = 8.0 N = 16 A = 16.000 g/mole ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 % ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 % ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 % ElmMassFraction: 30.00 % ElmAbundance 27.29 % Material: H2liquid density: 70.800 mg/cm3 RadL: 8.923 m Nucl.Int.Length: 4.944 m Imean: 19.200 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Hydrogen (H) Z = 1.0 N = 1 A = 1.010 g/mole ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 % ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Water H_2O density: 1.000 g/cm3 RadL: 36.092 cm Nucl.Int.Length: 75.356 cm Imean: 75.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Hydrogen (H) Z = 1.0 N = 1 A = 1.010 g/mole ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 % ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 % ElmMassFraction: 11.21 % ElmAbundance 66.67 % ---> Element: Oxygen (O) Z = 8.0 N = 16 A = 16.000 g/mole ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 % ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 % ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 % ElmMassFraction: 88.79 % ElmAbundance 33.33 % Material: liquidArgon density: 1.390 g/cm3 RadL: 14.064 cm Nucl.Int.Length: 86.076 cm Imean: 188.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Ar (Ar) Z = 18.0 N = 40 A = 39.948 g/mole ---> Isotope: Ar36 Z = 18 N = 36 A = 35.97 g/mole abundance: 0.337 % ---> Isotope: Ar38 Z = 18 N = 38 A = 37.96 g/mole abundance: 0.063 % ---> Isotope: Ar40 Z = 18 N = 40 A = 39.96 g/mole abundance: 99.600 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Carbon density: 2.267 g/cm3 RadL: 18.834 cm Nucl.Int.Length: 35.357 cm Imean: 81.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: C (C) Z = 6.0 N = 12 A = 12.011 g/mole ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.930 % ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.070 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Nucl.Int.Length: 38.880 cm Imean: 166.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Al (Al) Z = 13.0 N = 27 A = 26.982 g/mole ---> Isotope: Al27 Z = 13 N = 27 A = 26.98 g/mole abundance: 100.000 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Silicon density: 2.330 g/cm3 RadL: 9.366 cm Nucl.Int.Length: 45.660 cm Imean: 173.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Si (Si) Z = 14.0 N = 28 A = 28.085 g/mole ---> Isotope: Si28 Z = 14 N = 28 A = 27.98 g/mole abundance: 92.230 % ---> Isotope: Si29 Z = 14 N = 29 A = 28.98 g/mole abundance: 4.683 % ---> Isotope: Si30 Z = 14 N = 30 A = 29.97 g/mole abundance: 3.087 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Germanium density: 5.323 g/cm3 RadL: 2.301 cm Nucl.Int.Length: 27.431 cm Imean: 350.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Ge (Ge) Z = 32.0 N = 73 A = 72.613 g/mole ---> Isotope: Ge70 Z = 32 N = 70 A = 69.92 g/mole abundance: 20.840 % ---> Isotope: Ge72 Z = 32 N = 72 A = 71.92 g/mole abundance: 27.540 % ---> Isotope: Ge73 Z = 32 N = 73 A = 72.92 g/mole abundance: 7.730 % ---> Isotope: Ge74 Z = 32 N = 74 A = 73.92 g/mole abundance: 36.280 % ---> Isotope: Ge76 Z = 32 N = 76 A = 75.92 g/mole abundance: 7.610 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: NaI density: 3.670 g/cm3 RadL: 2.586 cm Nucl.Int.Length: 42.878 cm Imean: 452.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Sodium (Na) Z = 11.0 N = 23 A = 22.990 g/mole ---> Isotope: Na23 Z = 11 N = 23 A = 22.99 g/mole abundance: 100.000 % ElmMassFraction: 15.34 % ElmAbundance 50.00 % ---> Element: Iodine (I) Z = 53.0 N = 127 A = 126.900 g/mole ---> Isotope: I127 Z = 53 N = 127 A = 126.90 g/mole abundance: 100.000 % ElmMassFraction: 84.66 % ElmAbundance 50.00 % Material: Iodine density: 4.930 g/cm3 RadL: 1.720 cm Nucl.Int.Length: 35.676 cm Imean: 491.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Iodine (I) Z = 53.0 N = 127 A = 126.900 g/mole ---> Isotope: I127 Z = 53 N = 127 A = 126.90 g/mole abundance: 100.000 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: BGO density: 7.100 g/cm3 RadL: 1.123 cm Nucl.Int.Length: 22.806 cm Imean: 473.785 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Oxygen (O) Z = 8.0 N = 16 A = 16.000 g/mole ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 % ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 % ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 % ElmMassFraction: 15.41 % ElmAbundance 63.16 % ---> Element: Germanium (Ge) Z = 32.0 N = 73 A = 72.590 g/mole ---> Isotope: Ge70 Z = 32 N = 70 A = 69.92 g/mole abundance: 20.840 % ---> Isotope: Ge72 Z = 32 N = 72 A = 71.92 g/mole abundance: 27.540 % ---> Isotope: Ge73 Z = 32 N = 73 A = 72.92 g/mole abundance: 7.730 % ---> Isotope: Ge74 Z = 32 N = 74 A = 73.92 g/mole abundance: 36.280 % ---> Isotope: Ge76 Z = 32 N = 76 A = 75.92 g/mole abundance: 7.610 % ElmMassFraction: 17.48 % ElmAbundance 15.79 % ---> Element: Bismuth (Bi) Z = 83.0 N = 209 A = 208.980 g/mole ---> Isotope: Bi209 Z = 83 N = 209 A = 208.98 g/mole abundance: 100.000 % ElmMassFraction: 67.10 % ElmAbundance 21.05 % Material: Iron density: 7.870 g/cm3 RadL: 1.758 cm Nucl.Int.Length: 16.999 cm Imean: 286.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Fe (Fe) Z = 26.0 N = 56 A = 55.845 g/mole ---> Isotope: Fe54 Z = 26 N = 54 A = 53.94 g/mole abundance: 5.845 % ---> Isotope: Fe56 Z = 26 N = 56 A = 55.93 g/mole abundance: 91.754 % ---> Isotope: Fe57 Z = 26 N = 57 A = 56.94 g/mole abundance: 2.119 % ---> Isotope: Fe58 Z = 26 N = 58 A = 57.93 g/mole abundance: 0.282 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Tungsten density: 19.300 g/cm3 RadL: 3.504 mm Nucl.Int.Length: 10.312 cm Imean: 727.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: W (W) Z = 74.0 N = 184 A = 183.842 g/mole ---> Isotope: W180 Z = 74 N = 180 A = 179.95 g/mole abundance: 0.120 % ---> Isotope: W182 Z = 74 N = 182 A = 181.95 g/mole abundance: 26.500 % ---> Isotope: W183 Z = 74 N = 183 A = 182.95 g/mole abundance: 14.310 % ---> Isotope: W184 Z = 74 N = 184 A = 183.95 g/mole abundance: 30.640 % ---> Isotope: W186 Z = 74 N = 186 A = 185.95 g/mole abundance: 28.430 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Lead density: 11.350 g/cm3 RadL: 5.613 mm Nucl.Int.Length: 18.248 cm Imean: 823.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Pb (Pb) Z = 82.0 N = 207 A = 207.217 g/mole ---> Isotope: Pb204 Z = 82 N = 204 A = 203.97 g/mole abundance: 1.400 % ---> Isotope: Pb206 Z = 82 N = 206 A = 205.97 g/mole abundance: 24.100 % ---> Isotope: Pb207 Z = 82 N = 207 A = 206.98 g/mole abundance: 22.100 % ---> Isotope: Pb208 Z = 82 N = 208 A = 207.98 g/mole abundance: 52.400 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Uranium density: 18.950 g/cm3 RadL: 3.166 mm Nucl.Int.Length: 11.446 cm Imean: 890.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: U (U) Z = 92.0 N = 238 A = 238.029 g/mole ---> Isotope: U234 Z = 92 N = 234 A = 234.04 g/mole abundance: 0.005 % ---> Isotope: U235 Z = 92 N = 235 A = 235.04 g/mole abundance: 0.720 % ---> Isotope: U238 Z = 92 N = 238 A = 238.05 g/mole abundance: 99.275 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Galactic density: 0.000 mg/cm3 RadL: 204310098.490 pc Nucl.Int.Length: 3240.779 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm ---> Element: H (H) Z = 1.0 N = 1 A = 1.008 g/mole ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 % ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % /run/verbose 2 # /testem/det/setMat Water # /run/initialize userDetector->Construct() start. The Box is 100 m of Water World is registered to the default region. physicsList->Construct() start. physicsList->CheckParticleList() start. physicsList->setCut() start. # /process/inactivate eIoni /run/physicsModified /process/inactivate eBrem /run/physicsModified # /gun/particle e- # ##/tracking/verbose 1 ##/run/beamOn 1 # /gun/energy 5 MeV /run/beamOn 10000 ======================================================================= ====== Electromagnetic Physics Parameters ======== ======================================================================= LPM effect enabled 1 Enable creation and use of sampling tables 0 Apply cuts on all EM processes 0 Use combined TransportationWithMsc Disabled Use general process 0 Enable linear polarisation for gamma 0 Enable photoeffect sampling below K-shell 1 Enable sampling of quantum entanglement 0 X-section factor for integral approach 0.8 Min kinetic energy for tables 100 eV Max kinetic energy for tables 100 TeV Number of bins per decade of a table 7 Verbose level 1 Verbose level for worker thread 0 Bremsstrahlung energy threshold above which primary e+- is added to the list of secondary 100 TeV Bremsstrahlung energy threshold above which primary muon/hadron is added to the list of secondary 100 TeV Positron annihilation at rest model SimplePositronium Enable 3 gamma annihilation on fly 0 Lowest triplet kinetic energy 1 MeV Enable sampling of gamma linear polarisation 0 5D gamma conversion model type 0 5D gamma conversion model on isolated ion 0 Use RiGe 5D e+e- pair production model by muons 0 Livermore data directory epics_2017 ======================================================================= ====== Ionisation Parameters ======== ======================================================================= Step function for e+- (1, 1 mm) Step function for muons/hadrons (0.2, 0.1 mm) Step function for light ions (0.2, 0.1 mm) Step function for general ions (0.2, 0.1 mm) Lowest e+e- kinetic energy 1 keV Lowest muon/hadron kinetic energy 1 keV Use ICRU90 data 0 Fluctuations of dE/dx are enabled 0 Type of fluctuation model for leptons and hadrons Universal Use built-in Birks saturation 0 Build CSDA range enabled 0 Use cut as a final range enabled 0 Enable angular generator interface 0 Max kinetic energy for CSDA tables 1 GeV Max kinetic energy for NIEL computation 0 eV Linear loss limit 0.01 Read data from file for e+e- pair production by mu 0 ======================================================================= ====== Multiple Scattering Parameters ======== ======================================================================= Type of msc step limit algorithm for e+- 1 Type of msc step limit algorithm for muons/hadrons 0 Msc lateral displacement for e+- enabled 1 Msc lateral displacement for muons and hadrons 0 Urban msc model lateral displacement alg96 1 Range factor for msc step limit for e+- 0.04 Range factor for msc step limit for muons/hadrons 0.2 Geometry factor for msc step limitation of e+- 2.5 Safety factor for msc step limit for e+- 0.6 Skin parameter for msc step limitation of e+- 1 Lambda limit for msc step limit for e+- 1 mm Use Mott correction for e- scattering 0 Factor used for dynamic computation of angular limit between single and multiple scattering 1 Fixed angular limit between single and multiple scattering 3.1416 rad Upper energy limit for e+- multiple scattering 100 MeV Type of electron single scattering model 0 Type of nuclear form-factor 1 Screening factor 1 ======================================================================= phot: for gamma SubType=12 BuildTable=0 LambdaPrime table from 200 keV to 100 TeV in 61 bins ===== EM models for the G4Region DefaultRegionForTheWorld ====== PhotoElectric : Emin= 0 eV Emax= 100 TeV SauterGavrila compt: for gamma SubType=13 BuildTable=1 Lambda table from 100 eV to 1 MeV, 7 bins/decade, spline: 1 LambdaPrime table from 1 MeV to 100 TeV in 56 bins ===== EM models for the G4Region DefaultRegionForTheWorld ====== Klein-Nishina : Emin= 0 eV Emax= 100 TeV conv: for gamma SubType=14 BuildTable=1 Lambda table from 1.022 MeV to 100 TeV, 18 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== BetheHeitlerLPM : Emin= 0 eV Emax= 100 TeV ModifiedTsai msc: for e- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm eIoni: for e- XStype:0 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(1, 1 mm), integ: 0, fluct: 0, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 100 TeV eBrem: for e- XStype:0 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eBremSB : Emin= 0 eV Emax= 1 GeV ModifiedTsai eBremLPM : Emin= 1 GeV Emax= 100 TeV ModifiedTsai msc: for e+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm eIoni: for e+ XStype:0 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(1, 1 mm), integ: 0, fluct: 0, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 100 TeV eBrem: for e+ XStype:0 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eBremSB : Emin= 0 eV Emax= 1 GeV ModifiedTsai eBremLPM : Emin= 1 GeV Emax= 100 TeV ModifiedTsai annihil: for e+ SubType=5 AtRestModel:Simple BuildTable=0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eplus2gg : Emin= 0 eV Emax= 100 TeV msc: for proton SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for proton XStype:0 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 0, fluct: 0, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 100 TeV msc: for GenericIon SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm ionIoni: for GenericIon XStype:0 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 0, fluct: 0, linLossLim= 0.02 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 100 TeV msc: for alpha SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for alpha XStype:0 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 0, fluct: 0, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax=7.9452 MeV BetheBloch : Emin=7.9452 MeV Emax= 100 TeV msc: for anti_proton SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for anti_proton XStype:0 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 0, fluct: 0, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 100 TeV msc: for kaon+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for kaon+ XStype:0 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 0, fluct: 0, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax=1.05231 MeV BetheBloch : Emin=1.05231 MeV Emax= 100 TeV msc: for kaon- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for kaon- XStype:0 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 0, fluct: 0, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax=1.05231 MeV BetheBloch : Emin=1.05231 MeV Emax= 100 TeV msc: for mu+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm muIoni: for mu+ XStype:0 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 0, fluct: 0, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 200 keV MuBetheBloch : Emin= 200 keV Emax= 100 TeV muBrems: for mu+ XStype:0 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi muPairProd: for mu+ XStype:0 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 Sampling table 21x1001, from 0.85 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi msc: for mu- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm muIoni: for mu- XStype:0 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 0, fluct: 0, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 200 keV MuBetheBloch : Emin= 200 keV Emax= 100 TeV muBrems: for mu- XStype:0 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi muPairProd: for mu- XStype:0 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 Sampling table 21x1001, from 0.85 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi msc: for pi+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for pi+ XStype:0 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 0, fluct: 0, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax=297.505 keV BetheBloch : Emin=297.505 keV Emax= 100 TeV msc: for pi- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for pi- XStype:0 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 0, fluct: 0, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax=297.505 keV BetheBloch : Emin=297.505 keV Emax= 100 TeV Region -- -- appears in world volume This region is in the mass world. Root logical volume(s) : World Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] Materials : Galactic Water Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Region -- -- is not associated to any world. Root logical volume(s) : Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] Materials : Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm ========= Table of registered couples ============================ Index : 0 used in the geometry : Yes Material : Galactic Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes Material : Water Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 2.93515 keV e- 351.291 keV e+ 341.598 keV proton 100 keV Region(s) which use this couple : DefaultRegionForTheWorld ================================================================== Start closing geometry. -------------------------------------------------------------------------------- G4VoxelisationHelper::ReportVoxelStats -- Voxel Statistics Total memory consumed for geometry optimisation: 0 kByte Total CPU time elapsed for geometry optimisation: 0 seconds -------------------------------------------------------------------------------- ### Run 0 starts. ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={906770732717044781, 629165745432651234, 1235682547346241386, 68420008233404568, 2262190991329242458, 2266470399991071809, 1976726662926872232, 245458862506414172, 1955974201201518530, 2155248512522080758, 604170912935414061, 1116171330120743511, 1861018313684488333, 1296715403254578286, 1549011045957234151, 370819759640195970, 2230139271784837643} counter= 17sumtot= 1977567618660788324 --------------------------------------- Run terminated. Run Summary Number of events processed : 10000 User=0.030000s Real=0.023743s Sys=0.000000s The run consists of 10000 e- of 5 MeV through 100 m of Water (density: 1 g/cm3 ) Process calls frequency ---> msc = 10000 truePathLength : 2.7531 cm +- 6.8356 nm geomPathLength : 2.0688 cm +- 15.962 nm lateralDisplac : 0 fm +- 0 fm Psi : 0 mrad +- 0 mrad (0 deg +- 0 deg) Theta_plane : 0 mrad (0 deg) phi correlation: 0 +- 0 (std::cos(phi_pos - phi_dir)) Verification from G4EmCalculator. transport mean free path : 8.3239 cm range from restrict dE/dx: 2.7531 cm ---> effective facRange : 1 compute theta0 from Highland : 620.25 mrad (35.538 deg) ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={906770732717044781, 629165745432651234, 1235682547346241386, 68420008233404568, 2262190991329242458, 2266470399991071809, 1976726662926872232, 245458862506414172, 1955974201201518530, 2155248512522080758, 604170912935414061, 1116171330120743511, 1861018313684488333, 1296715403254578286, 1549011045957234151, 370819759640195970, 2230139271784837643} counter= 17sumtot= 1977567618660788324 --------------------------------------- # /gun/energy 100 keV /run/beamOn 10000 Region -- -- appears in world volume This region is in the mass world. Root logical volume(s) : World Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] Materials : Galactic Water Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Region -- -- is not associated to any world. Root logical volume(s) : Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] Materials : Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm ========= Table of registered couples ============================ Index : 0 used in the geometry : Yes Material : Galactic Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes Material : Water Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 2.93515 keV e- 351.291 keV e+ 341.598 keV proton 100 keV Region(s) which use this couple : DefaultRegionForTheWorld ================================================================== ### Run 1 starts. ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={906770732717044781, 629165745432651234, 1235682547346241386, 68420008233404568, 2262190991329242458, 2266470399991071809, 1976726662926872232, 245458862506414172, 1955974201201518530, 2155248512522080758, 604170912935414061, 1116171330120743511, 1861018313684488333, 1296715403254578286, 1549011045957234151, 370819759640195970, 2230139271784837643} counter= 17sumtot= 1977567618660788324 --------------------------------------- Run terminated. Run Summary Number of events processed : 10000 User=0.030000s Real=0.029904s Sys=0.000000s The run consists of 10000 e- of 100 keV through 100 m of Water (density: 1 g/cm3 ) Process calls frequency ---> msc = 10000 truePathLength : 6.1496 um +- 591.23 nm geomPathLength : 6.0268 um +- 567.72 nm lateralDisplac : 891.71 nm +- 128.01 nm Psi : 146.21 mrad +- 7.2125 mrad (8.377 deg +- 0.41325 deg) Theta_plane : 236.21 mrad (13.534 deg) phi correlation: 0.17482 +- 0.1395 (std::cos(phi_pos - phi_dir)) Verification from G4EmCalculator. transport mean free path : 153.44 um range from restrict dE/dx: 143.25 um ---> effective facRange : 0.040077 compute theta0 from Highland : 178.15 mrad (10.207 deg) ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={310976880687426508, 1713548588833330960, 1561521303267605489, 588569209103310290, 1954856900276614597, 724140034554777249, 713551730147133754, 402013335515644345, 971644034657465189, 193536399066193402, 893980623575290694, 766203625014548430, 498977489459909191, 1292557791726455116, 505565288099017699, 1782657595941510752, 964103503856943415} counter= 14sumtot= 2003346278501013374 --------------------------------------- G4 kernel has come to Quit state. Deleting G4Run (id:1) UserDetectorConstruction deleted 0x101a410 UserPhysicsList deleted 0x105dee0 UserActionInitialization deleted 0 UserWorkerInitialization deleted 0 UserWorkerThreadInitialization deleted 0 UserRunAction deleted. UserPrimaryGenerator deleted. RunManager is deleting RunManagerKernel. G4SDManager deleted. EventManager deleted. Units table cleared. TransportationManager deleted. Total navigation history collections cleaned: 6 G4RNGHelper object is deleted. ================== Deleting memory pools =================== Pool ID '20G4NavigationLevelRep', size : 0.00673 MB Pool ID '24G4ReferenceCountedHandleIvE', size : 0.000961 MB Pool ID '7G4Event', size : 0.000961 MB Pool ID '15G4PrimaryVertex', size : 0.000961 MB Pool ID '17G4PrimaryParticle', size : 0.000961 MB Pool ID '17G4DynamicParticle', size : 0.000961 MB Pool ID '7G4Track', size : 0.000961 MB Pool ID '18G4TouchableHistory', size : 0.000961 MB Pool ID '15G4CountedObjectIvE', size : 0.000961 MB Number of memory pools allocated: 9 of which, static: 0 Dynamic pools deleted: 9 / Total memory freed: 0.014 MB ============================================================ G4Allocator objects are deleted. UImanager deleted. StateManager deleted. RunManagerKernel is deleted. Good bye :) RunManager is deleted.