############################################ !!! WARNING - FPE detection is activated !!! ############################################ ************************************************************* Geant4 version Name: geant4-09-03-ref-06 (25-June-2010) Copyright : Geant4 Collaboration Reference : NIM A 506 (2003), 250-303 WWW : http://cern.ch/geant4 ************************************************************* /run/verbose 1 #/tracking/verbose 1 # # select EM PhysicsList # /testem/phys/addPhysics emstandard_opt1 PhysicsList::Set emstandard_opt1 EM physics # /testem/phys/setCuts 0.3 mm /testem/phys/VertexCuts 0.01 mm /testem/phys/MuonCuts 10.0 mm /testem/stepMax 100 mm # /testem/det/acceptance1 0.8271 0.02974 2.5 /testem/det/acceptance9 0.9570 0.01286 2.5 /testem/det/acceptance25 0.9793 0.00884 2.5 # #/process/eLoss/preciseRange true #/process/eLoss/integral true /run/initialize Ecal is 36 cm of G4_CESIUM_IODIDE Absorber is 2 mm of G4_Al Vertex is 3 cm of 3 layers of Si of 100 um npads= 601 ### New geometry is constructed # #/testem/histo/fileName e_pr.root #/testem/histo/fileType root /gun/particle e- /gun/energy 1 GeV /run/beamOn 100 phot: for gamma, applyCuts: 1 SubType= 12 ===== EM models for the G4Region DefaultRegionForTheWorld ====== PhotoElectric : Emin= 0 eV Emax= 10 TeV compt: for gamma, applyCuts: 1 SubType= 13 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Klein-Nishina : Emin= 0 eV Emax= 10 TeV conv: for gamma, applyCuts: 1 SubType= 14 Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bethe-Heitler : Emin= 0 eV Emax= 10 TeV msc: for e- SubType= 10 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc93 : Emin= 0 eV Emax= 10 TeV eIoni: for e- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.8, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 10 TeV eBrem: for e- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 LPM flag: 1 for E > 1 GeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== eBrem : Emin= 0 eV Emax= 1 GeV eBremRel : Emin= 1 GeV Emax= 10 TeV eIoni: for e+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.8, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 10 TeV eBrem: for e+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 LPM flag: 1 for E > 1 GeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== eBrem : Emin= 0 eV Emax= 1 GeV eBremRel : Emin= 1 GeV Emax= 10 TeV annihil: for e+, applyCuts: 1 SubType= 5 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eplus2gg : Emin= 0 eV Emax= 10 TeV msc: for proton SubType= 10 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc90 : Emin= 0 eV Emax= 10 TeV hIoni: for proton SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 10 TeV hBrems: for proton SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 10 TeV hPairProd: for proton SubType= 4 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 10 TeV msc: for GenericIon SubType= 10 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc90 : Emin= 0 eV Emax= 10 TeV ionIoni: for GenericIon SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 Stopping Power data for 17 ion/material pairs, nuclearStopping: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== BraggIon : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 10 TeV hIoni: for anti_proton SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 10 TeV hIoni: for kaon+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 1.05231 MeV BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV hBrems: for kaon+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 10 TeV hPairProd: for kaon+ SubType= 4 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 10 TeV hIoni: for kaon- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV hBrems: for kaon- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 10 TeV hPairProd: for kaon- SubType= 4 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 10 TeV muMsc: for mu+ SubType= 10 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc90 : Emin= 0 eV Emax= 10 TeV muIoni: for mu+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 200 keV BetheBloch : Emin= 200 keV Emax= 1 GeV MuBetheBloch : Emin= 1 GeV Emax= 10 TeV muBrems: for mu+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 10 TeV muPairProd: for mu+ SubType= 4 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 10 TeV muIoni: for mu- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 200 keV BetheBloch : Emin= 200 keV Emax= 1 GeV MuBetheBloch : Emin= 1 GeV Emax= 10 TeV muBrems: for mu- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 10 TeV muPairProd: for mu- SubType= 4 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 10 TeV hIoni: for pi+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 297.504 keV BetheBloch : Emin= 297.504 keV Emax= 10 TeV hBrems: for pi+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 10 TeV hPairProd: for pi+ SubType= 4 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 10 TeV msc: for pi- SubType= 10 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc90 : Emin= 0 eV Emax= 10 TeV hIoni: for pi- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 297.504 keV BetheBloch : Emin= 297.504 keV Emax= 10 TeV hBrems: for pi- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 10 TeV hPairProd: for pi- SubType= 4 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 10 TeV ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : G4_AIR Range cuts : gamma 300 um e- 300 um e+ 300 um proton 0 fm Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 0 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes recalculation needed : No Material : G4_CESIUM_IODIDE Range cuts : gamma 300 um e- 300 um e+ 300 um proton 0 fm Energy thresholds : gamma 18.2919 keV e- 291.524 keV e+ 281.891 keV proton 0 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : G4_Al Range cuts : gamma 300 um e- 300 um e+ 300 um proton 0 fm Energy thresholds : gamma 3.92577 keV e- 258.31 keV e+ 253.154 keV proton 0 eV 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 proton 0 fm Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 0 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 proton 0 fm Energy thresholds : gamma 990 eV e- 34.1725 keV e+ 33.9436 keV proton 0 eV 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 proton 0 fm Energy thresholds : gamma 990 eV e- 31.9516 keV e+ 31.7376 keV proton 0 eV 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 proton 0 fm Energy thresholds : gamma 60.499 keV e- 13.0779 MeV e+ 12.2279 MeV proton 0 eV 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 proton 0 fm Energy thresholds : gamma 20.3682 keV e- 4.52289 MeV e+ 4.28616 MeV proton 0 eV Region(s) which use this couple : MuonDetector ==================================================================== ### Run 0 start HistoManager: Histograms are booked and run has been started Start Run processing. Run terminated. Run Summary Number of events processed : 100 User=1.85s Real=2.86s Sys=0s RunAction: End of run actions are started HistoManager: End of run actions are started ================================================================= Number of events 100 Average number of e- 404.1 Average number of gamma 680.3 Average number of e+ 48.44 Average number of steps 3032 Edep 1x1 = 0.8268 +- 0.003033 res= 3.576 % Edep 3x3 = 0.9564 +- 0.00132 res= 1.36 % Edep 5x5 = 0.98 +- 0.0007972 res= 0.8052 % =========== Mean values without trancating ===================== Edep 1x1 = 0.8215 +- 0.003784 res= 4.607 % Edep 3x3 = 0.9552 +- 0.001452 res= 1.52 % Edep 5x5 = 0.9793 +- 0.0009604 res= 0.9807 % =========== Ratios without trancating =========================== E1/E9 = 0.8599 +- 0.003448 E1/E25 = 0.8388 +- 0.003691 E9/E25 = 0.9754 +- 0.001015 Beam Energy 1 GeV ================================================================== <<<<>>>> 100 events for Crystal Calorimeter Edep1x1: 0.8268 delEdep1x1= -0.0003204 nrms= -0.1077 Erms1x1: 0.02957 delErms1x1= -0.0001739 nrms= -0.05847 Edep3x3: 0.9564 delEdep3x3= -0.0006442 nrms= -0.5009 Erms3x3: 0.01301 delErms3x3= 0.0001451 nrms= 0.1128 Edep5x5: 0.98 delEdep5x5= 0.0007303 nrms= 0.8261 Erms5x5: 0.007892 delErms5x5= -0.0009484 nrms= -1.073 <<<<>>>> IS ACCEPTED Z bremsstrahlung photoeffect compton conversion 13 0 0 158 0 53 284 41 0 45 55 298 48 0 49 # G4 kernel has come to Quit state.