############################################ !!! WARNING - FPE detection is activated !!! ############################################ ************************************************************* Geant4 version Name: geant4-10-03-ref-03 [MT] (31-March-2017) Copyright : Geant4 Collaboration Reference : NIM A 506 (2003), 250-303 WWW : http://cern.ch/geant4 ************************************************************* G4PhysListFactory::GetReferencePhysList EMoption= 0 <<< Geant4 Physics List simulation engine: FTFP_BERT 2.0 <<< Reference Physics List FTFP_BERT is built /run/verbose 1 /tracking/verbose 0 # #/tile/setX1 -9.128 mm #/tile/setX2 9.128 mm /tile/setX1 -15 mm /tile/setX2 15 mm #/run/printProgress 100 # /run/setCut 1 mm # /run/initialize ***** Table : Nb of materials = 3 ***** Material: G4_Fe density: 7.874 g/cm3 RadL: 1.757 cm Nucl.Int.Length: 16.990 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: G4_POLYSTYRENE density: 1.060 g/cm3 RadL: 41.313 cm Nucl.Int.Length: 68.750 cm Imean: 68.700 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: 92.26 % ElmAbundance 50.00 % ---> 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: 7.74 % ElmAbundance 50.00 % Material: G4_AIR density: 1.205 mg/cm3 RadL: 303.921 m Nucl.Int.Length: 710.095 m Imean: 85.700 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: 0.01 % ElmAbundance 0.02 % ---> Element: N (N) Z = 7.0 N = 14 A = 14.007 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: 75.53 % ElmAbundance 78.44 % ---> Element: O (O) Z = 8.0 N = 16 A = 15.999 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: 23.18 % ElmAbundance 21.07 % ---> 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: 1.28 % ElmAbundance 0.47 % FTFP_BERT : new threshold between BERT and FTFP is over the interval for pions : 3 to 12 GeV for kaons : 3 to 12 GeV for proton : 3 to 12 GeV for neutron : 3 to 12 GeV ### Adding tracking cuts for neutron TimeCut(ns)= 10000 KinEnergyCut(MeV)= 0 /testhadr/histo/fileName e-_20GeV_FTFP_BERT_th-21 # /gun/particle e- /gun/direction 0.358367949545 0 0.933580426497 /gun/energy 20 GeV /run/beamOn 10000 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 AngularGenSauterGavrila FluoActive compt: for gamma SubType= 13 BuildTable= 1 Lambda table from 100 eV to 1 MeV, 7 bins per 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 per decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== BetheHeitler : Emin= 0 eV Emax= 80 GeV BetheHeitlerLPM : Emin= 80 GeV Emax= 100 TeV msc: for e- SubType= 10 RangeFactor= 0.04, stepLimitType: 1, latDisplacement: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 MeV Table with 42 bins Emin= 100 eV Emax= 100 MeV WentzelVIUni : Emin= 100 MeV Emax= 100 TeV Table with 42 bins Emin= 100 MeV Emax= 100 TeV eIoni: for e- 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 per decade, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 100 TeV eBrem: for e- 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 per 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 DipBustGen eBremLPM : Emin= 1 GeV Emax= 100 TeV DipBustGen CoulombScat: for e-, integral: 1 SubType= 1 BuildTable= 1 Lambda table from 100 MeV to 100 TeV, 7 bins per decade, spline: 1 180 < 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 RangeFactor= 0.04, stepLimitType: 1, latDisplacement: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 MeV Table with 42 bins Emin= 100 eV Emax= 100 MeV WentzelVIUni : Emin= 100 MeV Emax= 100 TeV Table with 42 bins Emin= 100 MeV Emax= 100 TeV eIoni: for e+ 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 per decade, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 100 TeV eBrem: for e+ 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 per 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 DipBustGen eBremLPM : Emin= 1 GeV Emax= 100 TeV DipBustGen annihil: for e+, integral: 1 SubType= 5 BuildTable= 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eplus2gg : Emin= 0 eV Emax= 100 TeV CoulombScat: for e+, integral: 1 SubType= 1 BuildTable= 1 Lambda table from 100 MeV to 100 TeV, 7 bins per decade, spline: 1 180 < 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 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 0, polarAngleLimit(deg)= 180 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Table with 84 bins Emin= 100 eV Emax= 100 TeV hIoni: for proton 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 per decade, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 100 TeV hBrems: for proton 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 per decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV hPairProd: for proton 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 per decade, spline: 1 Sampling table 17x1001; from 7.50618 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV CoulombScat: for proton, integral: 1 SubType= 1 BuildTable= 1 Used Lambda table of anti_proton 180 < 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 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV ionIoni: for GenericIon 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 per decade, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.02 Stopping Power data for 17 ion/material pairs ===== EM models for the G4Region DefaultRegionForTheWorld ====== BraggIon : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 100 TeV msc: for alpha SubType= 10 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Table with 84 bins Emin= 100 eV Emax= 100 TeV ionIoni: for alpha 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 per decade, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.02 ===== EM models for the G4Region DefaultRegionForTheWorld ====== BraggIon : Emin= 0 eV Emax= 7.9452 MeV BetheBloch : Emin= 7.9452 MeV Emax= 100 TeV msc: for anti_proton SubType= 10 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 0, polarAngleLimit(deg)= 180 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Table with 84 bins Emin= 100 eV Emax= 100 TeV hIoni: for anti_proton 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 per decade, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 100 TeV hBrems: for anti_proton 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 per decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV hPairProd: for anti_proton 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 per decade, spline: 1 Sampling table 17x1001; from 7.50618 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV CoulombScat: for anti_proton, integral: 1 SubType= 1 BuildTable= 1 Lambda table from threshold to 100 TeV, 7 bins per decade, spline: 1 180 < 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 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 0, polarAngleLimit(deg)= 180 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Table with 84 bins Emin= 100 eV Emax= 100 TeV hIoni: for kaon+ 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 per decade, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 1.05231 MeV BetheBloch : Emin= 1.05231 MeV Emax= 100 TeV hBrems: for kaon+ 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 per decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV hPairProd: for kaon+ 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 per decade, spline: 1 Sampling table 18x1001; from 3.94942 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV CoulombScat: for kaon+, integral: 1 SubType= 1 BuildTable= 1 Lambda table from threshold to 100 TeV, 7 bins per decade, spline: 1 180 < 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 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 0, polarAngleLimit(deg)= 180 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Table with 84 bins Emin= 100 eV Emax= 100 TeV hIoni: for kaon- 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 per decade, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV BetheBloch : Emin= 1.05231 MeV Emax= 100 TeV hBrems: for kaon- 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 per decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV hPairProd: for kaon- 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 per decade, spline: 1 Sampling table 18x1001; from 3.94942 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV CoulombScat: for kaon-, integral: 1 SubType= 1 BuildTable= 1 Used Lambda table of kaon+ 180 < 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 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 0, polarAngleLimit(deg)= 180 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Table with 84 bins Emin= 100 eV Emax= 100 TeV muIoni: for mu+ 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 per decade, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 200 keV BetheBloch : Emin= 200 keV Emax= 1 GeV MuBetheBloch : Emin= 1 GeV Emax= 100 TeV muBrems: for mu+ 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 per decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 100 TeV muPairProd: for mu+ 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 per decade, spline: 1 Sampling table 21x1001; from 1 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 100 TeV CoulombScat: for mu+, integral: 1 SubType= 1 BuildTable= 1 Lambda table from threshold to 100 TeV, 7 bins per decade, spline: 1 180 < 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 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 0, polarAngleLimit(deg)= 180 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Table with 84 bins Emin= 100 eV Emax= 100 TeV muIoni: for mu- 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 per decade, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 200 keV BetheBloch : Emin= 200 keV Emax= 1 GeV MuBetheBloch : Emin= 1 GeV Emax= 100 TeV muBrems: for mu- 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 per decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 100 TeV muPairProd: for mu- 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 per decade, spline: 1 Sampling table 21x1001; from 1 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 100 TeV CoulombScat: for mu-, integral: 1 SubType= 1 BuildTable= 1 Used Lambda table of mu+ 180 < 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 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 0, polarAngleLimit(deg)= 180 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Table with 84 bins Emin= 100 eV Emax= 100 TeV hIoni: for pi+ 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 per decade, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 297.505 keV BetheBloch : Emin= 297.505 keV Emax= 100 TeV hBrems: for pi+ 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 per decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV hPairProd: for pi+ 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 per decade, spline: 1 Sampling table 20x1001; from 1.11656 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV CoulombScat: for pi+, integral: 1 SubType= 1 BuildTable= 1 Lambda table from threshold to 100 TeV, 7 bins per decade, spline: 1 180 < 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 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 0, polarAngleLimit(deg)= 180 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Table with 84 bins Emin= 100 eV Emax= 100 TeV hIoni: for pi- 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 per decade, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 297.505 keV BetheBloch : Emin= 297.505 keV Emax= 100 TeV hBrems: for pi- 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 per decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV hPairProd: for pi- 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 per decade, spline: 1 Sampling table 20x1001; from 1.11656 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV CoulombScat: for pi-, integral: 1 SubType= 1 BuildTable= 1 Used Lambda table of pi+ 180 < Theta(degree) < 180; pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV ==================================================================== HADRONIC PROCESSES SUMMARY (verbose level 1) --------------------------------------------------- Hadronic Processes for neutron Process: hadElastic Model: hElasticCHIPS: 0 eV ---> 100 TeV Cr_sctns: G4NeutronElasticXS: 0 eV ---> 100 TeV Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: neutronInelastic Model: FTFP: 3 GeV ---> 100 TeV Model: BertiniCascade: 0 eV ---> 12 GeV Cr_sctns: G4NeutronInelasticXS: 0 eV ---> 100 TeV Cr_sctns: Barashenkov-Glauber: 0 eV ---> 100 TeV Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV Process: nCapture Model: nRadCapture: 0 eV ---> 100 TeV Cr_sctns: G4NeutronCaptureXS: 0 eV ---> 100 TeV Cr_sctns: GheishaCaptureXS: 0 eV ---> 100 TeV Process: nKiller --------------------------------------------------- Hadronic Processes for GenericIon Process: ionInelastic Model: Binary Light Ion Cascade: 0 eV /n ---> 4 GeV/n Model: FTFP: 2 GeV/n ---> 100 TeV/n Cr_sctns: Glauber-Gribov nucleus nucleus: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV --------------------------------------------------- Hadronic Processes for He3 Process: hadElastic Model: hElasticLHEP: 0 eV /n ---> 100 TeV/n Cr_sctns: Glauber-Gribov nucleus nucleus: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: He3Inelastic Model: Binary Light Ion Cascade: 0 eV /n ---> 4 GeV/n Model: FTFP: 2 GeV/n ---> 100 TeV/n Cr_sctns: Glauber-Gribov nucleus nucleus: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV --------------------------------------------------- Hadronic Processes for alpha Process: hadElastic Model: hElasticLHEP: 0 eV /n ---> 100 TeV/n Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: alphaInelastic Model: Binary Light Ion Cascade: 0 eV /n ---> 4 GeV/n Model: FTFP: 2 GeV/n ---> 100 TeV/n Cr_sctns: Glauber-Gribov nucleus nucleus: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV --------------------------------------------------- Hadronic Processes for anti_He3 Process: hadElastic Model: hElasticLHEP: 0 eV /n ---> 100.1 MeV/n Model: AntiAElastic: 100 MeV/n ---> 100 TeV/n Cr_sctns: AntiAGlauber: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: anti_He3Inelastic Model: FTFP: 0 eV /n ---> 100 TeV/n Cr_sctns: AntiAGlauber: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV Process: hFritiofCaptureAtRest --------------------------------------------------- Hadronic Processes for anti_alpha Process: hadElastic Model: hElasticLHEP: 0 eV /n ---> 100.1 MeV/n Model: AntiAElastic: 100 MeV/n ---> 100 TeV/n Cr_sctns: AntiAGlauber: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: anti_alphaInelastic Model: FTFP: 0 eV /n ---> 100 TeV/n Cr_sctns: AntiAGlauber: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV Process: hFritiofCaptureAtRest --------------------------------------------------- Hadronic Processes for anti_deuteron Process: hadElastic Model: hElasticLHEP: 0 eV /n ---> 100.1 MeV/n Model: AntiAElastic: 100 MeV/n ---> 100 TeV/n Cr_sctns: AntiAGlauber: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: anti_deuteronInelastic Model: FTFP: 0 eV /n ---> 100 TeV/n Cr_sctns: AntiAGlauber: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV Process: hFritiofCaptureAtRest --------------------------------------------------- Hadronic Processes for anti_neutron Process: hadElastic Model: hElasticLHEP: 0 eV ---> 100 TeV Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: anti_neutronInelastic Model: FTFP: 0 eV ---> 100 TeV Cr_sctns: AntiAGlauber: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV --------------------------------------------------- Hadronic Processes for anti_proton Process: hadElastic Model: hElasticLHEP: 0 eV ---> 100.1 MeV Model: AntiAElastic: 100 MeV ---> 100 TeV Cr_sctns: AntiAGlauber: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: anti_protonInelastic Model: FTFP: 0 eV ---> 100 TeV Cr_sctns: AntiAGlauber: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV Process: hFritiofCaptureAtRest --------------------------------------------------- Hadronic Processes for anti_triton Process: hadElastic Model: hElasticLHEP: 0 eV /n ---> 100.1 MeV/n Model: AntiAElastic: 100 MeV/n ---> 100 TeV/n Cr_sctns: AntiAGlauber: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: anti_tritonInelastic Model: FTFP: 0 eV /n ---> 100 TeV/n Cr_sctns: AntiAGlauber: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV Process: hFritiofCaptureAtRest --------------------------------------------------- Hadronic Processes for deuteron Process: hadElastic Model: hElasticLHEP: 0 eV /n ---> 100 TeV/n Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: dInelastic Model: Binary Light Ion Cascade: 0 eV /n ---> 4 GeV/n Model: FTFP: 2 GeV/n ---> 100 TeV/n Cr_sctns: Glauber-Gribov nucleus nucleus: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV --------------------------------------------------- Hadronic Processes for e+ Process: positronNuclear Model: G4ElectroVDNuclearModel: 0 eV ---> 1 PeV Cr_sctns: ElectroNuclearXS: 0 eV ---> 100 TeV Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV --------------------------------------------------- Hadronic Processes for e- Process: electronNuclear Model: G4ElectroVDNuclearModel: 0 eV ---> 1 PeV Cr_sctns: ElectroNuclearXS: 0 eV ---> 100 TeV Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV --------------------------------------------------- Hadronic Processes for gamma Process: photonNuclear Model: BertiniCascade: 0 eV ---> 3.5 GeV Model: TheoFSGenerator: 3 GeV ---> 100 TeV Cr_sctns: PhotoNuclearXS: 0 eV ---> 100 TeV Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV --------------------------------------------------- Hadronic Processes for kaon+ Process: hadElastic Model: hElasticLHEP: 0 eV ---> 100 TeV Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: kaon+Inelastic Model: FTFP: 3 GeV ---> 100 TeV Model: BertiniCascade: 0 eV ---> 12 GeV Cr_sctns: Glauber-Gribov: 0 eV ---> 2.88022e+295 J Cr_sctns: ChipsKaonPlusInelasticXS: 0 eV ---> 100 TeV Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV --------------------------------------------------- Hadronic Processes for kaon- Process: hadElastic Model: hElasticLHEP: 0 eV ---> 100 TeV Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: kaon-Inelastic Model: FTFP: 3 GeV ---> 100 TeV Model: BertiniCascade: 0 eV ---> 12 GeV Cr_sctns: Glauber-Gribov: 0 eV ---> 2.88022e+295 J Cr_sctns: ChipsKaonMinusInelasticXS: 0 eV ---> 100 TeV Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV Process: hBertiniCaptureAtRest --------------------------------------------------- Hadronic Processes for lambda Process: hadElastic Model: hElasticLHEP: 0 eV ---> 100 TeV Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: lambdaInelastic Model: BertiniCascade: 0 eV ---> 6 GeV Model: FTFP: 2 GeV ---> 100 TeV Cr_sctns: ChipsHyperonInelasticXS: 0 eV ---> 100 TeV Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV --------------------------------------------------- Hadronic Processes for mu+ Process: muonNuclear Model: G4MuonVDNuclearModel: 0 eV ---> 1 PeV Cr_sctns: KokoulinMuonNuclearXS: 0 eV ---> 100 TeV --------------------------------------------------- Hadronic Processes for mu- Process: muonNuclear Model: G4MuonVDNuclearModel: 0 eV ---> 1 PeV Cr_sctns: KokoulinMuonNuclearXS: 0 eV ---> 100 TeV Process: muMinusCaptureAtRest --------------------------------------------------- Hadronic Processes for pi+ Process: hadElastic Model: hElasticLHEP: 0 eV ---> 1.0001 GeV Model: hElasticGlauber: 1 GeV ---> 100 TeV Cr_sctns: Barashenkov-Glauber: 0 eV ---> 100 TeV Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: pi+Inelastic Model: FTFP: 3 GeV ---> 100 TeV Model: BertiniCascade: 0 eV ---> 12 GeV Cr_sctns: G4CrossSectionPairGG: 0 eV ---> 100 TeV G4CrossSectionPairGG: G4PiNuclearCrossSection cross sections below 91 GeV, Glauber-Gribov above Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV --------------------------------------------------- Hadronic Processes for pi- Process: hadElastic Model: hElasticLHEP: 0 eV ---> 1.0001 GeV Model: hElasticGlauber: 1 GeV ---> 100 TeV Cr_sctns: Barashenkov-Glauber: 0 eV ---> 100 TeV Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: pi-Inelastic Model: FTFP: 3 GeV ---> 100 TeV Model: BertiniCascade: 0 eV ---> 12 GeV Cr_sctns: G4CrossSectionPairGG: 0 eV ---> 100 TeV G4CrossSectionPairGG: G4PiNuclearCrossSection cross sections below 91 GeV, Glauber-Gribov above Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV Process: hBertiniCaptureAtRest --------------------------------------------------- Hadronic Processes for proton Process: hadElastic Model: hElasticCHIPS: 0 eV ---> 100 TeV Cr_sctns: ChipsProtonElasticXS: 0 eV ---> 100 TeV Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: protonInelastic Model: FTFP: 3 GeV ---> 100 TeV Model: BertiniCascade: 0 eV ---> 12 GeV Cr_sctns: Barashenkov-Glauber: 0 eV ---> 100 TeV Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV --------------------------------------------------- Hadronic Processes for triton Process: hadElastic Model: hElasticLHEP: 0 eV /n ---> 100 TeV/n Cr_sctns: GheishaElastic: 0 eV ---> 100 TeV Process: tInelastic Model: Binary Light Ion Cascade: 0 eV /n ---> 4 GeV/n Model: FTFP: 2 GeV/n ---> 100 TeV/n Cr_sctns: Glauber-Gribov nucleus nucleus: 0 eV ---> 2.88022e+295 J Cr_sctns: GheishaInelastic: 0 eV ---> 100 TeV ================================================================ ======================================================================= ====== Pre-compound/De-excitation Physics Parameters ======== ======================================================================= Type of pre-compound inverse x-section 3 Type of de-excitation inverse x-section 3 Min excitation energy (keV) 0.01 Level density (1/MeV) 0.1 Time limit for long lived isomeres (ns) 1e+08 Use new data files 1 Use complete data files 0 Correlated gamma emission flag 0 ======================================================================= Number of de-excitation channels 8 ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes Material : G4_AIR 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 : G4_Fe Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 20.6438 keV e- 1.29592 MeV e+ 1.21169 MeV proton 100 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes Material : G4_POLYSTYRENE Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 2.4362 keV e- 361.464 keV e+ 351.877 keV proton 100 keV Region(s) which use this couple : DefaultRegionForTheWorld ==================================================================== ### Run 0 starts. ### Run 0 start ### Histo::Save: Opended file for 26 histograms EventAction: Event # 0 started EventAction: Event # 100 started EventAction: Event # 200 started EventAction: Event # 300 started EventAction: Event # 400 started EventAction: Event # 500 started EventAction: Event # 600 started EventAction: Event # 700 started EventAction: Event # 800 started EventAction: Event # 900 started EventAction: Event # 1000 started EventAction: Event # 1100 started EventAction: Event # 1200 started EventAction: Event # 1300 started EventAction: Event # 1400 started EventAction: Event # 1500 started EventAction: Event # 1600 started EventAction: Event # 1700 started EventAction: Event # 1800 started EventAction: Event # 1900 started EventAction: Event # 2000 started EventAction: Event # 2100 started EventAction: Event # 2200 started EventAction: Event # 2300 started EventAction: Event # 2400 started EventAction: Event # 2500 started EventAction: Event # 2600 started EventAction: Event # 2700 started EventAction: Event # 2800 started EventAction: Event # 2900 started EventAction: Event # 3000 started EventAction: Event # 3100 started EventAction: Event # 3200 started EventAction: Event # 3300 started EventAction: Event # 3400 started EventAction: Event # 3500 started EventAction: Event # 3600 started EventAction: Event # 3700 started EventAction: Event # 3800 started EventAction: Event # 3900 started EventAction: Event # 4000 started EventAction: Event # 4100 started EventAction: Event # 4200 started EventAction: Event # 4300 started EventAction: Event # 4400 started EventAction: Event # 4500 started EventAction: Event # 4600 started EventAction: Event # 4700 started EventAction: Event # 4800 started EventAction: Event # 4900 started EventAction: Event # 5000 started EventAction: Event # 5100 started EventAction: Event # 5200 started EventAction: Event # 5300 started EventAction: Event # 5400 started EventAction: Event # 5500 started EventAction: Event # 5600 started EventAction: Event # 5700 started EventAction: Event # 5800 started EventAction: Event # 5900 started EventAction: Event # 6000 started EventAction: Event # 6100 started EventAction: Event # 6200 started EventAction: Event # 6300 started EventAction: Event # 6400 started EventAction: Event # 6500 started EventAction: Event # 6600 started EventAction: Event # 6700 started EventAction: Event # 6800 started EventAction: Event # 6900 started EventAction: Event # 7000 started EventAction: Event # 7100 started EventAction: Event # 7200 started EventAction: Event # 7300 started EventAction: Event # 7400 started EventAction: Event # 7500 started EventAction: Event # 7600 started EventAction: Event # 7700 started EventAction: Event # 7800 started EventAction: Event # 7900 started EventAction: Event # 8000 started EventAction: Event # 8100 started EventAction: Event # 8200 started EventAction: Event # 8300 started EventAction: Event # 8400 started EventAction: Event # 8500 started EventAction: Event # 8600 started EventAction: Event # 8700 started EventAction: Event # 8800 started EventAction: Event # 8900 started EventAction: Event # 9000 started EventAction: Event # 9100 started EventAction: Event # 9200 started EventAction: Event # 9300 started EventAction: Event # 9400 started EventAction: Event # 9500 started EventAction: Event # 9600 started EventAction: Event # 9700 started EventAction: Event # 9800 started EventAction: Event # 9900 started Run terminated. Run Summary Number of events processed : 10000 User=5784.7s Real=5791.4s Sys=1.11s RunAction: End of run actions are started HistoManager: End of run actions are started ======================================================== Beam particle e- Beam Energy(MeV) 20000 Number of events 10000 Average energy deposit (MeV) 641.8 RMS(MeV) 49.04 Average number of steps 3.203e+04 Average number of gamma 7348 Average number of e- 2.512e+04 Average number of e+ 621.3 Average number of neutrons 7.88 Average number of protons 6.691 Average number of antiprotons 0.0002 Average number of pi+ & pi- 0.0519 Average number of pi0 0.0322 Average number of kaons 0.0009 Average number of muons 0.0168 Average number of deuterons+tritons 0.0788 Average number of He3+alpha 0.1477 Average number of ions 11.69 Average number of forward neutrons 0 Average number of reflected neutrons 0 Average number of leaked neutrons 0 Average number of proton leak 0 Average number of pion leak 0 ======================================================== # G4 kernel has come to Quit state. ================== Deleting memory pools =================== Number of memory pools allocated: 12; of which, static: 0 Dynamic pools deleted: 12 / Total memory freed: 1.2 MB ============================================================ RunManagerKernel is deleted. Good bye :)