************************************************************* Geant4 version Name: geant4-09-01-patch-02-ref (9-May-2008) 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 mg/cm3 RadL: 204322175.595 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 Imean: 820.000 eV ---> Element: Pb (Pb) Z = 82.0 N = 207.2 A = 207.22 g/mole ---> Isotope: Pb Z = 82 N = 204 A = 203.97 g/mole abundance: 1.40 % ---> Isotope: Pb Z = 82 N = 206 A = 205.97 g/mole abundance: 24.10 % ---> Isotope: Pb Z = 82 N = 207 A = 206.98 g/mole abundance: 22.10 % ---> Isotope: Pb 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 Imean: 63.913 eV ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole ---> Isotope: C Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % ---> Isotope: C 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: H Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % ---> Isotope: H 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: Total cross sections from Sandia parametrisation. Sampling according PhotoElectric model compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV Sampling according Klein-Nishina model tables are built for gamma Lambda tables from 100 eV to 10 TeV in 220 bins. conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; sampling secondary e+e- according Bethe-Heitler model tables are built for gamma Lambda tables from 1.022 MeV to 10 TeV in 220 bins. msc: Model variant of multiple scattering for e- Lambda tables from 100 eV to 10 TeV in 220 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 1 ### The process eBrem is added to the list of collaborative processes of eIoni ### The process eBrem is added to the list of collaborative processes of eIoni ### The process muBrems is added to the list of collaborative processes of muIoni ### The process muPairProd is added to the list of collaborative processes of muIoni ### The process muBrems is added to the list of collaborative processes of muIoni ### The process muPairProd is added to the list of collaborative processes of muIoni eIoni: tables are built for e- dE/dx and range tables from 100 eV to 10 TeV in 220 bins. Lambda tables from threshold to 10 TeV in 220 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 Subcutoff sampling in 1 regions eBrem: tables are built for e- dE/dx and range tables from 100 eV to 10 TeV in 220 bins. Lambda tables from threshold to 10 TeV in 220 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1 Subcutoff sampling in 1 regions eIoni: tables are built for e+ dE/dx and range tables from 100 eV to 10 TeV in 220 bins. Lambda tables from threshold to 10 TeV in 220 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 Subcutoff sampling in 1 regions eBrem: tables are built for e+ dE/dx and range tables from 100 eV to 10 TeV in 220 bins. Lambda tables from threshold to 10 TeV in 220 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1 Subcutoff sampling in 1 regions annihil: Sampling according eplus2gg model tables are built for e+ Lambda tables from 100 eV to 100 TeV in 120 bins. msc: Model variant of multiple scattering for proton Lambda tables from 100 eV to 10 TeV in 220 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 0 hIoni: tables are built for proton dE/dx and range tables from 100 eV to 10 TeV in 220 bins. Lambda tables from threshold to 10 TeV in 220 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. NuclearStopping= 1 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for GenericIon LateralDisplacementFlag= 0 Skin= 0 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 1 ionIoni: tables are built for GenericIon dE/dx and range tables from 100 eV to 10 TeV in 220 bins. Lambda tables from threshold to 10 TeV in 220 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from BraggIon for protons below. NuclearStopping= 1 Stopping Power data for 8 ion/material pairs are used. Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1 Subcutoff sampling in 1 regions hIoni: tables are built for anti_proton dE/dx and range tables from 100 eV to 10 TeV in 220 bins. Lambda tables from threshold to 10 TeV in 220 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. NuclearStopping= 1 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for mu+ Lambda tables from 100 eV to 10 TeV in 220 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 0 muIoni: tables are built for mu+ dE/dx and range tables from 100 eV to 10 TeV in 220 bins. Lambda tables from threshold to 10 TeV in 220 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 Subcutoff sampling in 1 regions muBrems: tables are built for mu+ dE/dx and range tables from 100 eV to 10 TeV in 220 bins. Lambda tables from threshold to 10 TeV in 220 bins. Parametrised model Subcutoff sampling in 1 regions muPairProd: tables are built for mu+ dE/dx and range tables from 100 eV to 10 TeV in 220 bins. Lambda tables from threshold to 10 TeV in 220 bins. Parametrised model Subcutoff sampling in 1 regions muIoni: tables are built for mu- dE/dx and range tables from 100 eV to 10 TeV in 220 bins. Lambda tables from threshold to 10 TeV in 220 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 Subcutoff sampling in 1 regions muBrems: tables are built for mu- dE/dx and range tables from 100 eV to 10 TeV in 220 bins. Lambda tables from threshold to 10 TeV in 220 bins. Parametrised model Subcutoff sampling in 1 regions muPairProd: tables are built for mu- dE/dx and range tables from 100 eV to 10 TeV in 220 bins. Lambda tables from threshold to 10 TeV in 220 bins. Parametrised model Subcutoff sampling in 1 regions hIoni: tables are built for pi+ dE/dx and range tables from 100 eV to 10 TeV in 220 bins. Lambda tables from threshold to 10 TeV in 220 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV Parametrisation from Bragg for protons below. NuclearStopping= 1 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for pi- Lambda tables from 100 eV to 10 TeV in 220 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 0 hIoni: tables are built for pi- dE/dx and range tables from 100 eV to 10 TeV in 220 bins. Lambda tables from threshold to 10 TeV in 220 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV Parametrisation from Bragg for protons below. NuclearStopping= 1 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 Subcutoff sampling in 1 regions ========= 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=7754.58s Real=8087.67s Sys=4.4s ------------------------------------------------------------ material Total Edep sqrt(E0(GeV))*rmsE/Emean total tracklen Lead: 8.1015 GeV +- 56.04 MeV 2.187 +- 0.02187 % 6.08 m +- 4.37 cm Scintillator: 1.6426 GeV +- 40.09 MeV 7.718 +- 0.07718 % 8.81 m +- 21.8 cm ------------------------------------------------------------ Energy deposition from Energy flow balance : material Total Edep Lead: 17.7645 GeV Scintillator: -8.02465 GeV ------------------------------------------------------------ --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1955402620, 1872830150 ---------------------------------------- # G4 kernel has come to Quit state.