************************************************************* Geant4 version Name: (10-February-2006) Copyright : Geant4 Collaboration Reference : NIM A 506 (2003), 250-303 WWW : http://cern.ch/geant4 ************************************************************* Visualization Manager instantiating... Visualization Manager initialising... Registering graphics systems... You have successfully registered the following graphics systems. Current available graphics systems are: ASCIITree (ATree) DAWNFILE (DAWNFILE) GAGTree (GAGTree) G4HepRep (HepRepXML) G4HepRepFile (HepRepFile) RayTracer (RayTracer) VRML1FILE (VRML1FILE) VRML2FILE (VRML2FILE) FukuiRenderer (DAWN) OpenGLImmediateX (OGLIX) OpenGLStoredX (OGLSX) VRML1 (VRML1) VRML2 (VRML2) Registering model factories... You have successfully registered the following model factories. Registered model factories: drawByCharge drawByParticleID drawByOriginVolume Registered models: None Registered filter factories: particleFilter Registered Filters: None /run/verbose 1 /tracking/verbose 0 # /testem/det/setNbOfLayers 30 /testem/det/setNbOfAbsor 3 /testem/det/setAbsor 1 LeadSb 0.5 cm /testem/det/setAbsor 2 Scintillator 0.5 cm /testem/det/setAbsor 3 G4_AIR 0.3 cm /testem/det/setSizeYZ 60 cm # #/testem/phys/addPhysics g4v52 /testem/phys/addPhysics standard # #/testem/run/setHisto 0 100 1150. 1850. MeV #/testem/run/setHisto 1 100 27450. 28150. MeV # /gun/particle e- /gun/energy 10 GeV # #/process/eLoss/StepFunction 0.2 1 mm #/process/eLoss/preciseRange true # #/testem/tracking/stepMax 1 mum # #/random/setSavingFlag true #/random/resetEngineFrom first.rndm /testem/event/printModulo 1000 # /testem/phys/setCuts 0.01 mm /run/initialize ------------------------------------------------------------- ---> The calorimeter is 30 layers of: LeadSb: 5 mm Scintillator: 5 mm G4_AIR: 3 mm ------------------------------------------------------------- Material: Galactic density: 0.000 kg/m3 temperature: 2.73 K pressure: 0.00 atm RadLength: 204322175.595 pc ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: LeadSb density: 11.350 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 5.674 mm ---> Element: Sb (Sb) Z = 51.0 N = 121.9 A = 121.76 g/mole ---> Isotope: Sb Z = 51 N = 121 A = 120.90 g/mole abundance: 57.21 % ---> Isotope: Sb Z = 51 N = 123 A = 122.90 g/mole abundance: 42.79 % ElmMassFraction: 4.00 % ElmAbundance 6.62 % ---> 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: 96.00 % ElmAbundance 93.38 % Material: Scintillator density: 1.032 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 42.544 cm ---> 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 % Material: G4_AIR density: 1.205 mg/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 303.921 m ---> 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: 0.01 % ElmAbundance 0.02 % ---> Element: N (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ---> Isotope: N Z = 7 N = 14 A = 14.00 g/mole abundance: 99.63 % ---> Isotope: N Z = 7 N = 15 A = 15.00 g/mole abundance: 0.37 % ElmMassFraction: 75.53 % ElmAbundance 78.44 % ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ---> Isotope: O Z = 8 N = 16 A = 15.99 g/mole abundance: 99.76 % ---> Isotope: O Z = 8 N = 17 A = 17.00 g/mole abundance: 0.04 % ---> Isotope: O Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 23.18 % ElmAbundance 21.07 % ---> Element: Ar (Ar) Z = 18.0 N = 40.0 A = 39.95 g/mole ---> Isotope: Ar Z = 18 N = 36 A = 35.97 g/mole abundance: 0.34 % ---> Isotope: Ar Z = 18 N = 38 A = 37.96 g/mole abundance: 0.06 % ---> Isotope: Ar Z = 18 N = 40 A = 39.96 g/mole abundance: 99.60 % ElmMassFraction: 1.28 % ElmAbundance 0.47 % ------------------------------------------------------------- PhysicsList::SetCuts:CutLength : 1 mm /run/beamOn 10000 phot: Total cross sections from Sandia parametrisation. 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 100 GeV in 90 bins. conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; sampling secondary e+e- according to the Bethe-Heitler model tables are built for gamma Lambda tables from 1.022 MeV to 100 GeV in 100 bins. msc: Model variant of multiple scattering for e- Lambda tables from 100 eV to 100 TeV in 120 bins. Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 eIoni: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections from Moller+Bhabha, good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 eBrem: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections from a parametrisation based on the EEDL data library. Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. eIoni: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections from Moller+Bhabha, good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 eBrem: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections from a parametrisation based on the EEDL data library. Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. annihil: Heilter model of formula of annihilation into 2 photons 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 100 TeV in 120 bins. Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 hIoni: tables are built for proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used to proton dE/dx and range Bether-Bloch model for Escaled > 2 MeV, ICRU49 parametrisation for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 msc: Model variant of multiple scattering for GenericIon Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 ionIoni: tables are built for GenericIon dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used to proton dE/dx and range Bether-Bloch model for Escaled > 2 MeV, ICRU49 parametrisation for alpha particles below. Stopping Power data for 8 ion/material pairs are used. Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1 hIoni: tables are built for anti_proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used to proton dE/dx and range Bether-Bloch model for Escaled > 2 MeV, ICRU49 parametrisation for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 msc: Model variant of multiple scattering for mu+ Lambda tables from 100 eV to 100 TeV in 120 bins. Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 muIoni: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 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 muBrems: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model muPairProd: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model muIoni: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 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 muBrems: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model muPairProd: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model hIoni: tables are built for pi+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used to proton dE/dx and range Bether-Bloch model for Escaled > 0.297504 MeV, ICRU49 parametrisation for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 msc: Model variant of multiple scattering for pi- Lambda tables from 100 eV to 100 TeV in 120 bins. Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 hIoni: tables are built for pi- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used to proton dE/dx and range Bether-Bloch model for Escaled > 0.297504 MeV, ICRU49 parametrisation for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 ========= 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 : LeadSb Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 5.92678 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- 6.67753 keV e+ 6.59583 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 3 used in the geometry : Yes recalculation needed : No Material : G4_AIR 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 ==================================================================== ### 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=15453.47s Real=15682.25s Sys=4.34s ------------------------------------------------------------ material Total Edep sqrt(E0(GeV))*rmsE/Emean total tracklen LeadSb: 9.0875 GeV +- 43.06 MeV 1.498 +- 0.01498 % 6.8 m +- 3.54 cm Scintillator: 859.54 MeV +- 36.1 MeV 13.28 +- 0.1328 % 4.61 m +- 19.8 cm G4_AIR: 776.28 keV +- 84.44 keV 34.4 +- 0.344 % 3.56 m +- 33 cm ------------------------------------------------------------ Energy deposition from Energy flow balance : material Total Edep LeadSb: 9.08764 GeV Scintillator: 821.65 MeV G4_AIR: 38.3455 MeV ------------------------------------------------------------ --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1886945797, 1516082406 ---------------------------------------- # /testem/phys/setCuts 0.001 mm /run/initialize /run/beamOn 10000 phot: Total cross sections from Sandia parametrisation. 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 100 GeV in 90 bins. conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; sampling secondary e+e- according to the Bethe-Heitler model tables are built for gamma Lambda tables from 1.0219981 MeV to 100 GeV in 100 bins. msc: Model variant of multiple scattering for e- Lambda tables from 100 eV to 100 TeV in 120 bins. Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 eIoni: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections from Moller+Bhabha, good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 eBrem: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections from a parametrisation based on the EEDL data library. Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. eIoni: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections from Moller+Bhabha, good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 eBrem: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections from a parametrisation based on the EEDL data library. Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. annihil: Heilter model of formula of annihilation into 2 photons 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 100 TeV in 120 bins. Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 hIoni: tables are built for proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used to proton dE/dx and range Bether-Bloch model for Escaled > 2 MeV, ICRU49 parametrisation for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 msc: Model variant of multiple scattering for GenericIon Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 ionIoni: tables are built for GenericIon dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used to proton dE/dx and range Bether-Bloch model for Escaled > 2 MeV, ICRU49 parametrisation for alpha particles below. Stopping Power data for 8 ion/material pairs are used. Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1 hIoni: tables are built for anti_proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used to proton dE/dx and range Bether-Bloch model for Escaled > 2 MeV, ICRU49 parametrisation for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 msc: Model variant of multiple scattering for mu+ Lambda tables from 100 eV to 100 TeV in 120 bins. Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 muIoni: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 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 muBrems: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model muPairProd: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model muIoni: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 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 muBrems: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model muPairProd: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model hIoni: tables are built for pi+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used to proton dE/dx and range Bether-Bloch model for Escaled > 0.29750425 MeV, ICRU49 parametrisation for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 msc: Model variant of multiple scattering for pi- Lambda tables from 100 eV to 100 TeV in 120 bins. Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 hIoni: tables are built for pi- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used to proton dE/dx and range Bether-Bloch model for Escaled > 0.29750425 MeV, ICRU49 parametrisation for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : Galactic Range cuts : gamma 1 um e- 1 um e+ 1 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 : LeadSb Range cuts : gamma 1 um e- 1 um e+ 1 um Energy thresholds : gamma 1.1972034 keV e- 4.4697838 keV e+ 4.3610775 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : Scintillator Range cuts : gamma 1 um e- 1 um e+ 1 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 3 used in the geometry : Yes recalculation needed : No Material : G4_AIR Range cuts : gamma 1 um e- 1 um e+ 1 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV Region(s) which use this couple : DefaultRegionForTheWorld ==================================================================== ### Run 1 start. --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1886945797, 1516082406 ---------------------------------------- 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=68767.56s Real=69393.39s Sys=12.09s ------------------------------------------------------------ material Total Edep sqrt(E0(GeV))*rmsE/Emean total tracklen LeadSb: 9.0795 GeV +- 41.43 MeV 1.443 +- 0.01443 % 6.94 m +- 3.47 cm Scintillator: 868.48 MeV +- 35.78 MeV 13.03 +- 0.1303 % 4.65 m +- 19.6 cm G4_AIR: 806.71 keV +- 87.68 keV 34.37 +- 0.3437 % 3.57 m +- 32 cm ------------------------------------------------------------ Energy deposition from Energy flow balance : material Total Edep LeadSb: 9.07943 GeV Scintillator: 859.777 MeV G4_AIR: 9.47616 MeV ------------------------------------------------------------ --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 461617910, 273966158 ---------------------------------------- # Graphics systems deleted. Visualization Manager deleting... G4 kernel has come to Quit state.