************************************************************* Geant4 version Name: (10-February-2006) Copyright : Geant4 Collaboration Reference : NIM A 506 (2003), 250-303 WWW : http://cern.ch/geant4 ************************************************************* ### NIST DataBase for Materials is used G4Material: Mean excitation energy is changed for G4_AIR Iold= 85.627eV; Inew= 85.7 eV; G4Material: Mean excitation energy is changed for G4_Al Iold= 166.4eV; Inew= 166 eV; G4Material: Mean excitation energy is changed for G4_Si Iold= 173.6eV; Inew= 173 eV; G4Material: Mean excitation energy is changed for G4_Fe Iold= 286eV; Inew= 286 eV; G4Material: Mean excitation energy is changed for G4_CESIUM_IODIDE Iold= 491.166eV; Inew= 553.1 eV; 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/CalMat G4_PbWO4 /testem/det/EcalLength 230 mm /testem/det/EcalWidth 22 mm # #/testem/phys/addPhysics g4v52 /testem/phys/addPhysics standard PhysicsList::AddPhysicsList: /testem/phys/mscStepLimit false # /testem/phys/VertexCuts 0.01 mm /testem/phys/MuonCuts 10.0 mm #/testem/stepMax 10 mm # # /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/phys/setCuts 100 mm /run/initialize ***** Table : Nb of materials = 6 ***** 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 % Material: G4_Al density: 2.699 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 8.896 cm ---> Element: Al (Al) Z = 13.0 N = 27.0 A = 26.98 g/mole ---> Isotope: Al Z = 13 N = 27 A = 26.98 g/mole abundance: 100.00 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_Si density: 2.330 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 9.366 cm ---> Element: Si (Si) Z = 14.0 N = 28.1 A = 28.09 g/mole ---> Isotope: Si Z = 14 N = 28 A = 27.98 g/mole abundance: 92.23 % ---> Isotope: Si Z = 14 N = 29 A = 28.98 g/mole abundance: 4.68 % ---> Isotope: Si Z = 14 N = 30 A = 29.97 g/mole abundance: 3.09 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_Fe density: 7.874 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 1.757 cm ---> Element: Fe (Fe) Z = 26.0 N = 55.9 A = 55.85 g/mole ---> Isotope: Fe Z = 26 N = 54 A = 53.94 g/mole abundance: 5.84 % ---> Isotope: Fe Z = 26 N = 56 A = 55.93 g/mole abundance: 91.75 % ---> Isotope: Fe Z = 26 N = 57 A = 56.94 g/mole abundance: 2.12 % ---> Isotope: Fe Z = 26 N = 58 A = 57.93 g/mole abundance: 0.28 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_CESIUM_IODIDE density: 4.510 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 1.860 cm ---> Element: I (I) Z = 53.0 N = 127.0 A = 126.90 g/mole ---> Isotope: I Z = 53 N = 127 A = 126.90 g/mole abundance: 100.00 % ElmMassFraction: 48.85 % ElmAbundance 50.00 % ---> Element: Cs (Cs) Z = 55.0 N = 133.0 A = 132.91 g/mole ---> Isotope: Cs Z = 55 N = 133 A = 132.91 g/mole abundance: 100.00 % ElmMassFraction: 51.15 % ElmAbundance 50.00 % Material: G4_PbWO4 density: 8.280 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 8.925 mm ---> 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: 14.06 % ElmAbundance 66.67 % ---> Element: Pb (Pb) Z = 82.0 N = 207.2 A = 207.22 g/mole ElmMassFraction: 45.54 % ElmAbundance 16.67 % ---> Element: W (W) Z = 74.0 N = 183.8 A = 183.84 g/mole ElmMassFraction: 40.40 % ElmAbundance 16.67 % Ecal is 23 cm of G4_PbWO4 Absorber is 2 mm of G4_Al Vertex is 3 cm of 3 layers of Si of 100 um npads= 219 ### New geometry is constructed PhysicsList: world cuts are set cutG= 100 mm cutE= 100 mm Vertex cuts are set Muon cuts are set /testem/event/printModulo 1000 /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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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 : G4_AIR Range cuts : gamma 10 cm e- 10 cm e+ 10 cm Energy thresholds : gamma 1.20612 keV e- 86.3701 keV e+ 85.3134 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes recalculation needed : No Material : G4_PbWO4 Range cuts : gamma 10 cm e- 10 cm e+ 10 cm Energy thresholds : gamma 1.15638 MeV e- 230.876 MeV e+ 219.783 MeV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 10 cm e- 10 cm e+ 10 cm Energy thresholds : gamma 79.141 keV e- 56.7407 MeV e+ 54.0144 MeV 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 : VertexDetector Index : 4 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 33.8876 keV e+ 33.4729 keV Region(s) which use this couple : VertexDetector Index : 5 used in the geometry : Yes recalculation needed : No Material : G4_Si Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 31.4747 keV e+ 31.0896 keV Region(s) which use this couple : VertexDetector Index : 6 used in the geometry : Yes recalculation needed : No Material : G4_Fe Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 60.4019 keV e- 12.9518 MeV e+ 12.0297 MeV Region(s) which use this couple : MuonDetector Index : 7 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 20.4949 keV e- 4.54874 MeV e+ 4.2772 MeV Region(s) which use this couple : MuonDetector ==================================================================== ### Run 0 start HistoManager: Histograms are booked and run has been started /vis/scene/notifyHandlers WARNING: Scene "refresh" not found. /vis/scene/list to see scenes. Start Run processing. Run terminated. Run Summary Number of events processed : 10000 User=2724.35s Real=3220.49s Sys=0.47s RunAction: End of run actions are started HistoManager: End of run actions are started ======================================================== Number of events 10000 Average number of e- 5120 Average number of gamma 2656 Average number of e+ 583.8 Average number of steps 1.982e+04 Edep 1x1 = 0.7683 +- 0.0001331 res= 5.479 % Edep 3x3 = 0.9312 +- 6.525e-05 res= 2.216 % Edep 5x5 = 0.9708 +- 4.441e-05 res= 1.447 % ======================================================== /vis/viewer/update ERROR: G4VisCommandsViewerUpdate::SetNewValue: no current viewer. # /testem/phys/setCuts 10 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.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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2975 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.2 Step limitation 0 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.2975 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 : G4_AIR Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 990 eV e- 7.432 keV e+ 7.341 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes recalculation needed : No Material : G4_PbWO4 Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 219.3 keV e- 11.45 MeV e+ 10.64 MeV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 20.49 keV e- 4.549 MeV e+ 4.277 MeV 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 : VertexDetector Index : 4 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 33.89 keV e+ 33.47 keV Region(s) which use this couple : VertexDetector Index : 5 used in the geometry : Yes recalculation needed : No Material : G4_Si Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 31.47 keV e+ 31.09 keV Region(s) which use this couple : VertexDetector Index : 6 used in the geometry : Yes recalculation needed : No Material : G4_Fe Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 60.4 keV e- 12.95 MeV e+ 12.03 MeV Region(s) which use this couple : MuonDetector Index : 7 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 20.49 keV e- 4.549 MeV e+ 4.277 MeV Region(s) which use this couple : MuonDetector ==================================================================== ### Run 1 start HistoManager: Histograms are booked and run has been started /vis/scene/notifyHandlers WARNING: Scene "refresh" not found. /vis/scene/list to see scenes. Start Run processing. Run terminated. Run Summary Number of events processed : 10000 User=3759.69s Real=3916.46s Sys=1.12s RunAction: End of run actions are started HistoManager: End of run actions are started ======================================================== Number of events 10000 Average number of e- 7643 Average number of gamma 4208 Average number of e+ 551.2 Average number of steps 2.672e+04 Edep 1x1 = 0.7652 +- 0.0001325 res= 5.475 % Edep 3x3 = 0.935 +- 6.395e-05 res= 2.163 % Edep 5x5 = 0.9734 +- 4.62e-05 res= 1.501 % ======================================================== /vis/viewer/update ERROR: G4VisCommandsViewerUpdate::SetNewValue: no current viewer. # /testem/phys/setCuts 3 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.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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2975 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.2 Step limitation 0 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.2975 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 : G4_AIR Range cuts : gamma 3 mm e- 3 mm e+ 3 mm Energy thresholds : gamma 990 eV e- 1.396 keV e+ 1.379 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes recalculation needed : No Material : G4_PbWO4 Range cuts : gamma 3 mm e- 3 mm e+ 3 mm Energy thresholds : gamma 114.1 keV e- 3.144 MeV e+ 2.92 MeV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 3 mm e- 3 mm e+ 3 mm Energy thresholds : gamma 11.49 keV e- 1.448 MeV e+ 1.378 MeV 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 : VertexDetector Index : 4 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 33.89 keV e+ 33.47 keV Region(s) which use this couple : VertexDetector Index : 5 used in the geometry : Yes recalculation needed : No Material : G4_Si Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 31.47 keV e+ 31.09 keV Region(s) which use this couple : VertexDetector Index : 6 used in the geometry : Yes recalculation needed : No Material : G4_Fe Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 60.4 keV e- 12.95 MeV e+ 12.03 MeV Region(s) which use this couple : MuonDetector Index : 7 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 20.49 keV e- 4.549 MeV e+ 4.277 MeV Region(s) which use this couple : MuonDetector ==================================================================== ### Run 2 start HistoManager: Histograms are booked and run has been started /vis/scene/notifyHandlers WARNING: Scene "refresh" not found. /vis/scene/list to see scenes. Start Run processing. Run terminated. Run Summary Number of events processed : 10000 User=4189.29s Real=4767.96s Sys=0.93s RunAction: End of run actions are started HistoManager: End of run actions are started ======================================================== Number of events 10000 Average number of e- 8505 Average number of gamma 5052 Average number of e+ 546.5 Average number of steps 2.909e+04 Edep 1x1 = 0.7667 +- 0.0001316 res= 5.427 % Edep 3x3 = 0.9357 +- 6.246e-05 res= 2.111 % Edep 5x5 = 0.9738 +- 4.157e-05 res= 1.35 % ======================================================== /vis/viewer/update ERROR: G4VisCommandsViewerUpdate::SetNewValue: no current viewer. # /testem/phys/setCuts 1 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.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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2975 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.2 Step limitation 0 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.2975 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 : G4_AIR Range cuts : gamma 1 mm e- 1 mm e+ 1 mm 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 : G4_PbWO4 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 84.78 keV e- 1.132 MeV e+ 1.064 MeV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 6.887 keV e- 596.7 keV e+ 568 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 : VertexDetector Index : 4 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 33.89 keV e+ 33.47 keV Region(s) which use this couple : VertexDetector Index : 5 used in the geometry : Yes recalculation needed : No Material : G4_Si Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 31.47 keV e+ 31.09 keV Region(s) which use this couple : VertexDetector Index : 6 used in the geometry : Yes recalculation needed : No Material : G4_Fe Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 60.4 keV e- 12.95 MeV e+ 12.03 MeV Region(s) which use this couple : MuonDetector Index : 7 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 20.49 keV e- 4.549 MeV e+ 4.277 MeV Region(s) which use this couple : MuonDetector ==================================================================== ### Run 3 start HistoManager: Histograms are booked and run has been started /vis/scene/notifyHandlers WARNING: Scene "refresh" not found. /vis/scene/list to see scenes. Start Run processing. Run terminated. Run Summary Number of events processed : 10000 User=4477.38s Real=5540.04s Sys=0.63s RunAction: End of run actions are started HistoManager: End of run actions are started ======================================================== Number of events 10000 Average number of e- 9079 Average number of gamma 5511 Average number of e+ 543.7 Average number of steps 3.069e+04 Edep 1x1 = 0.767 +- 0.0001341 res= 5.529 % Edep 3x3 = 0.936 +- 6.308e-05 res= 2.131 % Edep 5x5 = 0.9739 +- 4.225e-05 res= 1.372 % ======================================================== /vis/viewer/update ERROR: G4VisCommandsViewerUpdate::SetNewValue: no current viewer. # /testem/phys/setCuts 0.7 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.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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2975 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.2 Step limitation 0 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.2975 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 : G4_AIR Range cuts : gamma 700 um e- 700 um e+ 700 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 : G4_PbWO4 Range cuts : gamma 700 um e- 700 um e+ 700 um Energy thresholds : gamma 74.69 keV e- 842.2 keV e+ 792 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 700 um e- 700 um e+ 700 um Energy thresholds : gamma 5.845 keV e- 460.8 keV e+ 444 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 : VertexDetector Index : 4 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 33.89 keV e+ 33.47 keV Region(s) which use this couple : VertexDetector Index : 5 used in the geometry : Yes recalculation needed : No Material : G4_Si Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 31.47 keV e+ 31.09 keV Region(s) which use this couple : VertexDetector Index : 6 used in the geometry : Yes recalculation needed : No Material : G4_Fe Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 60.4 keV e- 12.95 MeV e+ 12.03 MeV Region(s) which use this couple : MuonDetector Index : 7 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 20.49 keV e- 4.549 MeV e+ 4.277 MeV Region(s) which use this couple : MuonDetector ==================================================================== ### Run 4 start HistoManager: Histograms are booked and run has been started /vis/scene/notifyHandlers WARNING: Scene "refresh" not found. /vis/scene/list to see scenes. Start Run processing. Run terminated. Run Summary Number of events processed : 10000 User=4604.13s Real=5552.97s Sys=0.59s RunAction: End of run actions are started HistoManager: End of run actions are started ======================================================== Number of events 10000 Average number of e- 9343 Average number of gamma 5704 Average number of e+ 542.7 Average number of steps 3.136e+04 Edep 1x1 = 0.7673 +- 0.0001366 res= 5.629 % Edep 3x3 = 0.9361 +- 6.543e-05 res= 2.21 % Edep 5x5 = 0.974 +- 4.172e-05 res= 1.355 % ======================================================== /vis/viewer/update ERROR: G4VisCommandsViewerUpdate::SetNewValue: no current viewer. # /testem/phys/setCuts 0.3 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.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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2975 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.2 Step limitation 0 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.2975 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 : G4_AIR Range cuts : gamma 300 mum e- 300 mum e+ 300 mum 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 : G4_PbWO4 Range cuts : gamma 300 mum e- 300 mum e+ 300 mum Energy thresholds : gamma 42.23 keV e- 438.6 keV e+ 422.7 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 300 mum e- 300 mum e+ 300 mum Energy thresholds : gamma 3.897 keV e- 258.3 keV e+ 252.1 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 : VertexDetector Index : 4 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 33.89 keV e+ 33.47 keV Region(s) which use this couple : VertexDetector Index : 5 used in the geometry : Yes recalculation needed : No Material : G4_Si Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 31.47 keV e+ 31.09 keV Region(s) which use this couple : VertexDetector Index : 6 used in the geometry : Yes recalculation needed : No Material : G4_Fe Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 60.4 keV e- 12.95 MeV e+ 12.03 MeV Region(s) which use this couple : MuonDetector Index : 7 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 20.49 keV e- 4.549 MeV e+ 4.277 MeV Region(s) which use this couple : MuonDetector ==================================================================== ### Run 5 start HistoManager: Histograms are booked and run has been started /vis/scene/notifyHandlers WARNING: Scene "refresh" not found. /vis/scene/list to see scenes. Start Run processing. Run terminated. Run Summary Number of events processed : 10000 User=5142s Real=7039.73s Sys=0.98s RunAction: End of run actions are started HistoManager: End of run actions are started ======================================================== Number of events 10000 Average number of e- 1.051e+04 Average number of gamma 6561 Average number of e+ 539.4 Average number of steps 3.431e+04 Edep 1x1 = 0.7686 +- 0.0001322 res= 5.441 % Edep 3x3 = 0.9367 +- 6.365e-05 res= 2.149 % Edep 5x5 = 0.9743 +- 4.22e-05 res= 1.37 % ======================================================== /vis/viewer/update ERROR: G4VisCommandsViewerUpdate::SetNewValue: no current viewer. # /testem/phys/setCuts 0.1 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.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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2975 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.2 Step limitation 0 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.2975 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 : G4_AIR Range cuts : gamma 100 um e- 100 um e+ 100 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 : G4_PbWO4 Range cuts : gamma 100 um e- 100 um e+ 100 um Energy thresholds : gamma 20.06 keV e- 209.6 keV e+ 204.5 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 100 um e- 100 um e+ 100 um Energy thresholds : gamma 2.314 keV e- 129.7 keV e+ 128.1 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 : VertexDetector Index : 4 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 33.89 keV e+ 33.47 keV Region(s) which use this couple : VertexDetector Index : 5 used in the geometry : Yes recalculation needed : No Material : G4_Si Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 31.47 keV e+ 31.09 keV Region(s) which use this couple : VertexDetector Index : 6 used in the geometry : Yes recalculation needed : No Material : G4_Fe Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 60.4 keV e- 12.95 MeV e+ 12.03 MeV Region(s) which use this couple : MuonDetector Index : 7 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 20.49 keV e- 4.549 MeV e+ 4.277 MeV Region(s) which use this couple : MuonDetector ==================================================================== ### Run 6 start HistoManager: Histograms are booked and run has been started /vis/scene/notifyHandlers WARNING: Scene "refresh" not found. /vis/scene/list to see scenes. Start Run processing. Run terminated. Run Summary Number of events processed : 10000 User=6028.05s Real=8454.82s Sys=0.84s RunAction: End of run actions are started HistoManager: End of run actions are started ======================================================== Number of events 10000 Average number of e- 1.256e+04 Average number of gamma 7736 Average number of e+ 536.2 Average number of steps 3.929e+04 Edep 1x1 = 0.7688 +- 0.0001454 res= 5.98 % Edep 3x3 = 0.9371 +- 6.734e-05 res= 2.272 % Edep 5x5 = 0.9745 +- 4.255e-05 res= 1.381 % ======================================================== /vis/viewer/update ERROR: G4VisCommandsViewerUpdate::SetNewValue: no current viewer. # /testem/phys/setCuts 0.07 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.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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2975 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.2 Step limitation 0 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.2975 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 : G4_AIR Range cuts : gamma 70 um e- 70 um e+ 70 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 : G4_PbWO4 Range cuts : gamma 70 um e- 70 um e+ 70 um Energy thresholds : gamma 15.57 keV e- 167.9 keV e+ 161.8 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 70 um e- 70 um e+ 70 um Energy thresholds : gamma 1.964 keV e- 105.2 keV e+ 103.9 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 : VertexDetector Index : 4 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 33.89 keV e+ 33.47 keV Region(s) which use this couple : VertexDetector Index : 5 used in the geometry : Yes recalculation needed : No Material : G4_Si Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 31.47 keV e+ 31.09 keV Region(s) which use this couple : VertexDetector Index : 6 used in the geometry : Yes recalculation needed : No Material : G4_Fe Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 60.4 keV e- 12.95 MeV e+ 12.03 MeV Region(s) which use this couple : MuonDetector Index : 7 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 20.49 keV e- 4.549 MeV e+ 4.277 MeV Region(s) which use this couple : MuonDetector ==================================================================== ### Run 7 start HistoManager: Histograms are booked and run has been started /vis/scene/notifyHandlers WARNING: Scene "refresh" not found. /vis/scene/list to see scenes. Start Run processing. Run terminated. Run Summary Number of events processed : 10000 User=6386.85s Real=6458.66s Sys=1.49s RunAction: End of run actions are started HistoManager: End of run actions are started ======================================================== Number of events 10000 Average number of e- 1.343e+04 Average number of gamma 8138 Average number of e+ 535.2 Average number of steps 4.135e+04 Edep 1x1 = 0.7694 +- 0.0001356 res= 5.573 % Edep 3x3 = 0.9374 +- 6.207e-05 res= 2.094 % Edep 5x5 = 0.9746 +- 3.962e-05 res= 1.285 % ======================================================== /vis/viewer/update ERROR: G4VisCommandsViewerUpdate::SetNewValue: no current viewer. # /testem/phys/setCuts 0.03 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.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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2 Step limitation 0 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.2975 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.2 Step limitation 0 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.2975 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 : G4_AIR Range cuts : gamma 30 mum e- 30 mum e+ 30 mum 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 : G4_PbWO4 Range cuts : gamma 30 mum e- 30 mum e+ 30 mum Energy thresholds : gamma 8.806 keV e- 98.89 keV e+ 97.68 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 30 mum e- 30 mum e+ 30 mum Energy thresholds : gamma 1.322 keV e- 64.28 keV e+ 63.49 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 : VertexDetector Index : 4 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 33.89 keV e+ 33.47 keV Region(s) which use this couple : VertexDetector Index : 5 used in the geometry : Yes recalculation needed : No Material : G4_Si Range cuts : gamma 10 um e- 10 um e+ 10 um Energy thresholds : gamma 990 eV e- 31.47 keV e+ 31.09 keV Region(s) which use this couple : VertexDetector Index : 6 used in the geometry : Yes recalculation needed : No Material : G4_Fe Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 60.4 keV e- 12.95 MeV e+ 12.03 MeV Region(s) which use this couple : MuonDetector Index : 7 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 20.49 keV e- 4.549 MeV e+ 4.277 MeV Region(s) which use this couple : MuonDetector ==================================================================== ### Run 8 start HistoManager: Histograms are booked and run has been started /vis/scene/notifyHandlers WARNING: Scene "refresh" not found. /vis/scene/list to see scenes. Start Run processing. Run terminated. Run Summary Number of events processed : 10000 User=7471.44s Real=7534.93s Sys=1.9s RunAction: End of run actions are started HistoManager: End of run actions are started ======================================================== Number of events 10000 Average number of e- 1.619e+04 Average number of gamma 9053 Average number of e+ 533.4 Average number of steps 4.76e+04 Edep 1x1 = 0.7692 +- 0.000139 res= 5.712 % Edep 3x3 = 0.9374 +- 6.205e-05 res= 2.093 % Edep 5x5 = 0.9746 +- 4.005e-05 res= 1.3 % ======================================================== /vis/viewer/update ERROR: G4VisCommandsViewerUpdate::SetNewValue: no current viewer. # Graphics systems deleted. Visualization Manager deleting... G4 kernel has come to Quit state.