############################################ !!! 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 = 7 ***** 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 % 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: CaCO3 density: 2.800 g/cm3 RadL: 8.581 cm Nucl.Int.Length: 34.749 cm Imean: 123.236 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Calcium (Ca) Z = 20.0 N = 40 A = 40.080 g/mole ---> Isotope: Ca40 Z = 20 N = 40 A = 39.96 g/mole abundance: 96.941 % ---> Isotope: Ca42 Z = 20 N = 42 A = 41.96 g/mole abundance: 0.647 % ---> Isotope: Ca43 Z = 20 N = 43 A = 42.96 g/mole abundance: 0.135 % ---> Isotope: Ca44 Z = 20 N = 44 A = 43.96 g/mole abundance: 2.086 % ---> Isotope: Ca46 Z = 20 N = 46 A = 45.95 g/mole abundance: 0.004 % ---> Isotope: Ca48 Z = 20 N = 48 A = 47.95 g/mole abundance: 0.187 % ElmMassFraction: 40.04 % ElmAbundance 20.00 % ---> Element: Carbon (C) Z = 6.0 N = 12 A = 12.010 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: 12.00 % ElmAbundance 20.00 % ---> 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: 47.96 % ElmAbundance 60.00 % Material: Carbon density: 2.265 g/cm3 RadL: 18.851 cm Nucl.Int.Length: 35.388 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: 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: Tin density: 7.310 g/cm3 RadL: 1.206 cm Nucl.Int.Length: 23.531 cm Imean: 488.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Sn (Sn) Z = 50.0 N = 119 A = 118.710 g/mole ---> Isotope: Sn112 Z = 50 N = 112 A = 111.91 g/mole abundance: 0.970 % ---> Isotope: Sn114 Z = 50 N = 114 A = 113.90 g/mole abundance: 0.660 % ---> Isotope: Sn115 Z = 50 N = 115 A = 114.90 g/mole abundance: 0.340 % ---> Isotope: Sn116 Z = 50 N = 116 A = 115.90 g/mole abundance: 14.540 % ---> Isotope: Sn117 Z = 50 N = 117 A = 116.90 g/mole abundance: 7.680 % ---> Isotope: Sn118 Z = 50 N = 118 A = 117.90 g/mole abundance: 24.220 % ---> Isotope: Sn119 Z = 50 N = 119 A = 118.90 g/mole abundance: 8.590 % ---> Isotope: Sn120 Z = 50 N = 120 A = 119.90 g/mole abundance: 32.580 % ---> Isotope: Sn122 Z = 50 N = 122 A = 121.90 g/mole abundance: 4.630 % ---> Isotope: Sn124 Z = 50 N = 124 A = 123.91 g/mole abundance: 5.790 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: HI density: 2.800 g/cm3 RadL: 3.049 cm Nucl.Int.Length: 60.880 cm Imean: 462.393 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: 0.79 % 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: 99.21 % ElmAbundance 50.00 % /run/verbose 2 # /testem/det/setMat Iron /testem/det/setSize 1 m # /testem/phys/addPhysics emstandard_opt0 PhysicsList::AddPhysicsList: /run/setCut 1 mm # /run/initialize userDetector->Construct() start. The Box is 1 m of Iron Iron is registered to the default region. physicsList->Construct() start. physicsList->CheckParticleList() start. physicsList->setCut() start. # /gun/particle mu+ /gun/energy 10 TeV # /run/printProgress 1000 # /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 1 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 1000 PeV 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+- (0.2, 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 1 Type of fluctuation model for leptons and hadrons Urban 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 1000 PeV in 89 bins ===== EM models for the G4Region DefaultRegionForTheWorld ====== LivermorePhElectric : Emin= 0 eV Emax= 1000 PeV SauterGavrila Fluo 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 1000 PeV in 84 bins ===== EM models for the G4Region DefaultRegionForTheWorld ====== Klein-Nishina : Emin= 0 eV Emax= 1000 PeV conv: for gamma SubType=14 BuildTable=1 Lambda table from 1.022 MeV to 1000 PeV, 14 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== BetheHeitlerLPM : Emin= 0 eV Emax= 1000 PeV ModifiedTsai Rayl: for gamma SubType=11 BuildTable=1 Lambda table from 100 eV to 150 keV, 7 bins/decade, spline: 0 LambdaPrime table from 150 keV to 1000 PeV in 90 bins ===== EM models for the G4Region DefaultRegionForTheWorld ====== LivermoreRayleigh : Emin= 0 eV Emax= 100 TeV CullenGenerator msc: for e- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 MeV Nbins=42 100 eV - 100 MeV StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm WentzelVIUni : Emin= 100 MeV Emax= 100 TeV Nbins=42 100 MeV - 100 TeV StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm eIoni: for e- XStype:1 SubType=2 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 1000 PeV eBrem: for e- XStype:4 SubType=3 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 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= 1000 PeV ModifiedTsai CoulombScat: for e- XStype:1 SubType=1 BuildTable=1 Lambda table from 100 MeV to 1000 PeV, 7 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 100 MeV Emax= 100 TeV msc: for e+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 MeV Nbins=42 100 eV - 100 MeV StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm WentzelVIUni : Emin= 100 MeV Emax= 100 TeV Nbins=42 100 MeV - 100 TeV StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm eIoni: for e+ XStype:1 SubType=2 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 1000 PeV eBrem: for e+ XStype:4 SubType=3 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 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= 1000 PeV ModifiedTsai annihil: for e+ XStype:2 SubType=5 AtRestModel:Simple BuildTable=0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eplus2gg : Emin= 0 eV Emax= 100 TeV CoulombScat: for e+ XStype:1 SubType=1 BuildTable=1 Lambda table from 100 MeV to 1000 PeV, 7 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 100 MeV Emax= 100 TeV msc: for proton SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : 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:1 SubType=2 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 1000 PeV hBrems: for proton XStype:1 SubType=3 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 1000 PeV ModifiedMephi hPairProd: for proton XStype:1 SubType=4 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 Sampling table 33x1001, from 7.50618 GeV to 1e+06 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 1000 PeV ModifiedMephi CoulombScat: for proton XStype:1 SubType=1 BuildTable=1 Lambda table from threshold to 1000 PeV, 7 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV 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:1 SubType=2 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.02 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 1000 PeV 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 ionIoni: for alpha XStype:1 SubType=2 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.02 ===== EM models for the G4Region DefaultRegionForTheWorld ====== BraggIon : Emin= 0 eV Emax=7.9452 MeV BetheBloch : Emin=7.9452 MeV Emax= 1000 PeV msc: for anti_proton SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : 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:1 SubType=2 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 1000 PeV hBrems: for anti_proton XStype:1 SubType=3 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 1000 PeV ModifiedMephi hPairProd: for anti_proton XStype:1 SubType=4 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 Sampling table 33x1001, from 7.50618 GeV to 1e+06 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 1000 PeV ModifiedMephi CoulombScat: for anti_proton XStype:1 SubType=1 BuildTable=1 Lambda table from threshold to 1000 PeV, 7 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for kaon+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : 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:1 SubType=2 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax=1.05231 MeV BetheBloch : Emin=1.05231 MeV Emax= 1000 PeV hBrems: for kaon+ XStype:1 SubType=3 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 1000 PeV ModifiedMephi hPairProd: for kaon+ XStype:1 SubType=4 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 Sampling table 34x1001, from 3.94942 GeV to 1e+06 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 1000 PeV ModifiedMephi CoulombScat: for kaon+ XStype:1 SubType=1 BuildTable=1 Lambda table from threshold to 1000 PeV, 7 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for kaon- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : 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:1 SubType=2 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax=1.05231 MeV BetheBloch : Emin=1.05231 MeV Emax= 1000 PeV hBrems: for kaon- XStype:1 SubType=3 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 1000 PeV ModifiedMephi hPairProd: for kaon- XStype:1 SubType=4 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 Sampling table 34x1001, from 3.94942 GeV to 1e+06 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 1000 PeV ModifiedMephi CoulombScat: for kaon- XStype:1 SubType=1 BuildTable=1 Used Lambda table of kaon+ ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for mu+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : 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:1 SubType=2 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 200 keV MuBetheBloch : Emin= 200 keV Emax= 1000 PeV muBrems: for mu+ XStype:1 SubType=3 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 1000 PeV ModifiedMephi muPairProd: for mu+ XStype:1 SubType=4 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 Sampling table 37x1001, from 0.85 GeV to 1e+06 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 1000 PeV ModifiedMephi CoulombScat: for mu+ XStype:1 SubType=1 BuildTable=1 Lambda table from threshold to 1000 PeV, 7 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for mu- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : 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:1 SubType=2 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 200 keV MuBetheBloch : Emin= 200 keV Emax= 1000 PeV muBrems: for mu- XStype:1 SubType=3 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 1000 PeV ModifiedMephi muPairProd: for mu- XStype:1 SubType=4 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 Sampling table 37x1001, from 0.85 GeV to 1e+06 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 1000 PeV ModifiedMephi CoulombScat: for mu- XStype:1 SubType=1 BuildTable=1 Used Lambda table of mu+ ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for pi+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : 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:1 SubType=2 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax=297.505 keV BetheBloch : Emin=297.505 keV Emax= 1000 PeV hBrems: for pi+ XStype:1 SubType=3 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 1000 PeV ModifiedMephi hPairProd: for pi+ XStype:1 SubType=4 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 Sampling table 36x1001, from 1.11656 GeV to 1e+06 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 1000 PeV ModifiedMephi CoulombScat: for pi+ XStype:1 SubType=1 BuildTable=1 Lambda table from threshold to 1000 PeV, 7 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for pi- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : 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:1 SubType=2 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax=297.505 keV BetheBloch : Emin=297.505 keV Emax= 1000 PeV hBrems: for pi- XStype:1 SubType=3 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 1000 PeV ModifiedMephi hPairProd: for pi- XStype:1 SubType=4 dE/dx and range tables from 100 eV to 1000 PeV in 112 bins Lambda tables from threshold to 1000 PeV, 7 bins/decade, spline: 1 Sampling table 36x1001, from 1.11656 GeV to 1e+06 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 1000 PeV ModifiedMephi CoulombScat: for pi- XStype:1 SubType=1 BuildTable=1 Used Lambda table of pi+ ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV Region -- -- appears in world volume This region is in the mass world. Root logical volume(s) : Iron Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] Materials : Iron 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 : Iron Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 20.7519 keV e- 1.28828 MeV e+ 1.21724 MeV 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. ### Run 0 start. ------- 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 --------------------------------------- --> Event 0 starts. --> Event 1000 starts. --> Event 2000 starts. --> Event 3000 starts. --> Event 4000 starts. --> Event 5000 starts. --> Event 6000 starts. --> Event 7000 starts. --> Event 8000 starts. --> Event 9000 starts. Run terminated. Run Summary Number of events processed : 10000 User=0.680000s Real=0.682784s Sys=0.000000s The run consists of 10000 mu+ of 10 TeV through 1 m of Iron (density: 7.9 g/cm3 ) Number of process calls ---> CoulombScat : 5817 muIoni : 465476 muPairProd : 63643 muBrems : 651 Simulation: total CrossSection = 0.53559 /cm MeanFreePath = 1.8671 cm massicCrossSection = 0.068054 cm2/g Theory: total CrossSection = 0.52947 /cm MeanFreePath = 1.8887 cm massicCrossSection = 0.067277 cm2/g # /gun/particle pi+ /run/beamOn 10000 Region -- -- appears in world volume This region is in the mass world. Root logical volume(s) : Iron Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] Materials : Iron 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 : Iron Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 20.7519 keV e- 1.28828 MeV e+ 1.21724 MeV proton 100 keV Region(s) which use this couple : DefaultRegionForTheWorld ================================================================== ### Run 1 starts. ### Run 1 start. ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={882454084924701664, 2270259229220857311, 818577390581940897, 1426604042526396198, 480266226005173583, 188530774845791518, 1049894952900439173, 2108254491617653549, 67514774349258322, 175092158049147372, 393893130671613424, 1529173506857517301, 1747618713700092971, 2151849549381702022, 721140998242719756, 579078769140434288, 1816672379963618080} counter= 17sumtot= 2265974108483199772 --------------------------------------- --> Event 0 starts. --> Event 1000 starts. --> Event 2000 starts. --> Event 3000 starts. --> Event 4000 starts. --> Event 5000 starts. --> Event 6000 starts. --> Event 7000 starts. --> Event 8000 starts. --> Event 9000 starts. Run terminated. Run Summary Number of events processed : 10000 User=0.620000s Real=0.630539s Sys=0.000000s The run consists of 10000 pi+ of 10 TeV through 1 m of Iron (density: 7.9 g/cm3 ) Number of process calls ---> hIoni : 436874 hPairProd : 59172 CoulombScat : 5569 hBrems : 393 Simulation: total CrossSection = 0.50201 /cm MeanFreePath = 1.992 cm massicCrossSection = 0.063788 cm2/g # /gun/particle proton /run/beamOn 10000 Region -- -- appears in world volume This region is in the mass world. Root logical volume(s) : Iron Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] Materials : Iron 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 : Iron Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 20.7519 keV e- 1.28828 MeV e+ 1.21724 MeV proton 100 keV Region(s) which use this couple : DefaultRegionForTheWorld ================================================================== ### Run 2 starts. ### Run 2 start. ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={2103257823577124208, 469380737387581810, 386042433805554138, 1802430644607842505, 804448371540535249, 1779466859180359812, 1937382105945851007, 365603651364239426, 1236294566231758049, 530694216919414896, 1036938576786736727, 1786903812069656131, 966308926261532388, 907406917814257267, 1632242837280314955, 454243073969002515, 212028608884676175} counter= 2sumtot= 2270173099130579601 --------------------------------------- --> Event 0 starts. --> Event 1000 starts. --> Event 2000 starts. --> Event 3000 starts. --> Event 4000 starts. --> Event 5000 starts. --> Event 6000 starts. --> Event 7000 starts. --> Event 8000 starts. --> Event 9000 starts. Run terminated. Run Summary Number of events processed : 10000 User=0.590000s Real=0.589387s Sys=0.000000s The run consists of 10000 proton of 10 TeV through 1 m of Iron (density: 7.9 g/cm3 ) Number of process calls ---> hIoni : 436614 hPairProd : 33957 CoulombScat : 5719 hBrems : 11 Simulation: total CrossSection = 0.4763 /cm MeanFreePath = 2.0995 cm massicCrossSection = 0.060521 cm2/g # G4 kernel has come to Quit state. Deleting G4Run (id:2) UserDetectorConstruction deleted 0x15c8ca0 UserPhysicsList deleted 0x1609098 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 '17G4DynamicParticle', size : 0.0154 MB Pool ID '7G4Event', size : 0.000961 MB Pool ID '15G4PrimaryVertex', size : 0.000961 MB Pool ID '17G4PrimaryParticle', size : 0.000961 MB Pool ID '7G4Track', size : 0.0308 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.059 MB ============================================================ G4Allocator objects are deleted. UImanager deleted. StateManager deleted. RunManagerKernel is deleted. Good bye :) RunManager is deleted.