************************************************************* Geant4 version Name: geant4-09-02-ref-05 (5-June-2009) Copyright : Geant4 Collaboration Reference : NIM A 506 (2003), 250-303 WWW : http://cern.ch/geant4 ************************************************************* /run/verbose 1 /tracking/verbose 0 # /testem/det/setNbOfLayers 66 /testem/det/setNbOfAbsor 2 /testem/det/setAbsor 1 Lead 0.2 cm /testem/det/setAbsor 2 Scintillator 0.4 cm /testem/det/setSizeYZ 20 cm # /testem/phys/addPhysics emstandard_opt2 # #/testem/run/setHisto 0 100 1150. 1850. MeV #/testem/run/setHisto 1 100 27450. 28150. MeV # /gun/particle e- /gun/energy 10 GeV /gun/position -198 0 0 mm # #/random/setSavingFlag true #/random/resetEngineFrom first.rndm /testem/event/printModulo 1000 # /testem/phys/setCuts 0.01 mm /run/initialize ------------------------------------------------------------- ---> The calorimeter is 66 layers of: Lead: 2 mm Scintillator: 4 mm ------------------------------------------------------------- Material: Galactic density: 0.000 kg/m3 RadL: 204322111.300 pc Nucl.Int.Length: 113728944.943 pc Imean: 21.800 eV temperature: 2.73 K pressure: 0.00 atm ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Lead density: 11.350 g/cm3 RadL: 5.613 mm Nucl.Int.Length: 18.247 cm Imean: 820.000 eV ---> Element: Pb (Pb) Z = 82.0 N = 207.2 A = 207.22 g/mole ---> Isotope: Pb204 Z = 82 N = 204 A = 203.97 g/mole abundance: 1.40 % ---> Isotope: Pb206 Z = 82 N = 206 A = 205.97 g/mole abundance: 24.10 % ---> Isotope: Pb207 Z = 82 N = 207 A = 206.98 g/mole abundance: 22.10 % ---> Isotope: Pb208 Z = 82 N = 208 A = 207.98 g/mole abundance: 52.40 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Scintillator density: 1.032 g/cm3 RadL: 42.544 cm Nucl.Int.Length: 70.066 cm Imean: 63.913 eV ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 91.47 % ElmAbundance 47.37 % ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 8.53 % ElmAbundance 52.63 % ------------------------------------------------------------- PhysicsList::SetCuts:CutLength : 1 mm /run/beamOn 10000 phot: for gamma, applyCuts: 1 SubType= 12 ===== EM models for the G4Region DefaultRegionForTheWorld ====== PhotoElectric : Emin= 0 eV Emax= 10 TeV compt: for gamma, applyCuts: 1 SubType= 13 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Klein-Nishina : Emin= 0 eV Emax= 10 TeV conv: for gamma, applyCuts: 1 SubType= 14 Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bethe-Heitler : Emin= 0 eV Emax= 10 TeV muMsc: for e- SubType= 10 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 RangeFactor= 0.2, step limit type: 1, lateralDisplacement: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 10 TeV eIoni: for e- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, 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= 10 TeV eBrem: for e- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 LPM flag: 1 for E > 1 GeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== eBrem : Emin= 0 eV Emax= 1 GeV eBremRel : Emin= 1 GeV Emax= 10 TeV CoulombScat: for e-, applyCuts: 1 SubType= 1 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 11.4592 < Theta(degree) < 180, Eth(MeV)= 10000; q2Max(GeV^2)= 1e+06 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 10 GeV CoulombScattering : Emin= 10 GeV Emax= 10 TeV eIoni: for e+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, 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= 10 TeV eBrem: for e+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 LPM flag: 1 for E > 1 GeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== eBrem : Emin= 0 eV Emax= 1 GeV eBremRel : Emin= 1 GeV Emax= 10 TeV annihil: for e+, applyCuts: 1 SubType= 5 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eplus2gg : Emin= 0 eV Emax= 10 TeV CoulombScat: for e+, applyCuts: 1 SubType= 1 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 11.4592 < Theta(degree) < 180, Eth(MeV)= 10000; q2Max(GeV^2)= 1e+06 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 10 GeV CoulombScattering : Emin= 10 GeV Emax= 10 TeV muMsc: for proton SubType= 10 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 RangeFactor= 0.2, step limit type: 1, lateralDisplacement: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 10 TeV hIoni: for proton SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 NuclearStopping= 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 10 TeV hBrems: for proton SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, 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 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV msc: for GenericIon SubType= 10 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 0, skin= 3 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMscUni90 : Emin= 0 eV Emax= 100 TeV ionIoni: for GenericIon SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.15 Stopping Power data for 17 ion/material pairs, nuclearStopping: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== BraggIon : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 10 TeV hIoni: for anti_proton SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 NuclearStopping= 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 10 TeV hIoni: for kaon+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 NuclearStopping= 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 1.05231 MeV BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV hBrems: for kaon+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, 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 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV hIoni: for kaon- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 NuclearStopping= 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 1.05231 MeV BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV hBrems: for kaon- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, 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 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV muMsc: for mu+ SubType= 10 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 RangeFactor= 0.2, step limit type: 1, lateralDisplacement: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 10 TeV muIoni: for mu+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 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= 10 TeV muBrems: for mu+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 10 TeV muPairProd: for mu+ SubType= 4 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 10 TeV CoulombScat: for mu+, applyCuts: 1 SubType= 1 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; q2Max(GeV^2)= 1e+06 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 10 MeV CoulombScattering : Emin= 10 MeV Emax= 10 TeV muIoni: for mu- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 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= 10 TeV muBrems: for mu- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 10 TeV muPairProd: for mu- SubType= 4 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 10 TeV CoulombScat: for mu-, applyCuts: 1 SubType= 1 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; q2Max(GeV^2)= 1e+06 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 10 MeV CoulombScattering : Emin= 10 MeV Emax= 10 TeV hIoni: for pi+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 NuclearStopping= 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 297.504 keV BetheBloch : Emin= 297.504 keV Emax= 10 TeV hBrems: for pi+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, 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 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV muMsc: for pi- SubType= 10 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 RangeFactor= 0.2, step limit type: 1, lateralDisplacement: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 10 TeV hIoni: for pi- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 NuclearStopping= 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 297.504 keV BetheBloch : Emin= 297.504 keV Emax= 10 TeV hBrems: for pi- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, 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 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : Galactic Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes recalculation needed : No Material : Lead Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 5.97389 keV e- 57.5352 keV e+ 56.8313 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : Scintillator Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 14.4914 keV e+ 14.112 keV Region(s) which use this couple : DefaultRegionForTheWorld ==================================================================== ### Run 0 start. --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 9876, 54321 ---------------------------------------- Start Run processing. ---> Begin Of Event: 0 ---> Begin Of Event: 1000 ---> Begin Of Event: 2000 ---> Begin Of Event: 3000 ---> Begin Of Event: 4000 ---> Begin Of Event: 5000 ---> Begin Of Event: 6000 ---> Begin Of Event: 7000 ---> Begin Of Event: 8000 ---> Begin Of Event: 9000 Run terminated. Run Summary Number of events processed : 10000 User=21546.12s Real=21914.98s Sys=16.32s ------------------------------------------------------------ material Edep RMS sqrt(E0(GeV))*rmsE/Emean total tracklen Lead: 8.2958 GeV : 65.36 MeV 2.492 +- 0.02492 % 6.22 m +- 5.13 cm Scintillator: 1.4456 GeV : 50.42 MeV 11.03 +- 0.1103 % 7.57 m +- 27.3 cm ------------------------------------------------------------ Beam particle e- E = 10 GeV Energy deposition from Energy flow balance : material Total Edep Lead: 8.29579 GeV Scintillator: 1.44559 GeV ------------------------------------------------------------ --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1347107054, 153013207 ---------------------------------------- # G4 kernel has come to Quit state.