############################################ !!! 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 ************************************************************* ***** Table : Nb of materials = 14 ***** Material: G4_Al density: 2.699 g/cm3 RadL: 8.896 cm Nucl.Int.Length: 38.894 cm Imean: 166.000 eV ---> Element: Al (Al) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_Si density: 2.330 g/cm3 RadL: 9.366 cm Nucl.Int.Length: 45.660 cm Imean: 173.000 eV ---> Element: Si (Si) Z = 14.0 N = 28.1 A = 28.09 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_Fe density: 7.874 g/cm3 RadL: 1.757 cm Nucl.Int.Length: 16.990 cm Imean: 286.000 eV ---> Element: Fe (Fe) Z = 26.0 N = 55.8 A = 55.85 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_Cu density: 8.960 g/cm3 RadL: 1.436 cm Nucl.Int.Length: 15.588 cm Imean: 322.000 eV ---> Element: Cu (Cu) Z = 29.0 N = 63.5 A = 63.55 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_Ge density: 5.323 g/cm3 RadL: 2.301 cm Nucl.Int.Length: 27.431 cm Imean: 350.000 eV ---> Element: Ge (Ge) Z = 32.0 N = 72.6 A = 72.61 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_Mo density: 10.220 g/cm3 RadL: 9.591 mm Nucl.Int.Length: 15.677 cm Imean: 424.000 eV ---> Element: Mo (Mo) Z = 42.0 N = 95.9 A = 95.93 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_Ta density: 16.654 g/cm3 RadL: 4.094 mm Nucl.Int.Length: 11.887 cm Imean: 718.000 eV ---> Element: Ta (Ta) Z = 73.0 N = 180.9 A = 180.95 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_W density: 19.300 g/cm3 RadL: 3.504 mm Nucl.Int.Length: 10.312 cm Imean: 727.000 eV ---> Element: W (W) Z = 74.0 N = 183.8 A = 183.84 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_Au density: 19.320 g/cm3 RadL: 3.344 mm Nucl.Int.Length: 10.540 cm Imean: 790.000 eV ---> Element: Au (Au) Z = 79.0 N = 197.0 A = 196.97 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_Pb density: 11.350 g/cm3 RadL: 5.613 mm Nucl.Int.Length: 18.248 cm Imean: 823.000 eV ---> Element: Pb (Pb) Z = 82.0 N = 207.2 A = 207.22 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_PbWO4 density: 8.280 g/cm3 RadL: 8.925 mm Nucl.Int.Length: 20.740 cm Imean: 542.741 eV ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole 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 % Material: G4_SODIUM_IODIDE density: 3.667 g/cm3 RadL: 2.588 cm Nucl.Int.Length: 42.914 cm Imean: 452.000 eV ---> Element: Na (Na) Z = 11.0 N = 23.0 A = 22.99 g/mole ElmMassFraction: 15.34 % ElmAbundance 50.00 % ---> Element: I (I) Z = 53.0 N = 126.9 A = 126.90 g/mole ElmMassFraction: 84.66 % ElmAbundance 50.00 % Material: G4_AIR density: 1.205 mg/cm3 RadL: 303.921 m Nucl.Int.Length: 710.095 m Imean: 85.700 eV temperature: 273.15 K pressure: 1.00 atm ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 0.01 % ElmAbundance 0.02 % ---> Element: N (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 75.53 % ElmAbundance 78.44 % ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 23.18 % ElmAbundance 21.07 % ---> Element: Ar (Ar) Z = 18.0 N = 39.9 A = 39.95 g/mole ElmMassFraction: 1.28 % ElmAbundance 0.47 % Material: Water density: 1.000 g/cm3 RadL: 36.083 cm Nucl.Int.Length: 75.423 cm Imean: 75.000 eV ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.19 % ElmAbundance 66.67 % ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.81 % ElmAbundance 33.33 % /run/verbose 2 # /testem/det/setMat G4_Si /testem/det/setSizeX 1 mm /testem/det/setSizeYZ 1 mm # /testem/phys/addPhysics local PhysicsList::AddPhysicsList: # /run/initialize userDetector->Construct() start. --------------------------------------------------------- ---> The Absorber is 1 mm of G4_Si divided in 1 slices of 1 mm --------------------------------------------------------- Absorber is registered to the default region. physicsList->Construct() start. physicsList->Construct() start. physicsList->CheckParticleList() start. physicsList->setCut() start. PhysicsList::SetCuts:CutLength : 1 mm # /gun/particle e- /gun/energy 500 keV # /testem/histo/setFileType root /testem/histo/setFileName testem11 /testem/histo/setHisto 1 100 0. 1. mm ----> SetHisto 1: Edep (MeV/mm) along absorber (mm); 100 bins from 0 mm to 1 mm /testem/histo/setHisto 2 102 0. 510. keV ----> SetHisto 2: total Energy deposited in absorber (keV); 102 bins from 0 keV to 510 keV /testem/histo/setHisto 3 100 0. 2. mm ----> SetHisto 3: true track length of the primary particle (mm); 100 bins from 0 mm to 2 mm /testem/histo/setHisto 4 100 0. 200. um ----> SetHisto 4: true step size of the primary particle (um); 100 bins from 0 um to 200 um /testem/histo/setHisto 5 100 0. 1. mm ----> SetHisto 5: projected range of the primary particle (mm); 100 bins from 0 mm to 1 mm # /testem/event/printModulo 10000 # /run/beamOn 20000 phot: for gamma SubType= 12 ===== EM models for the G4Region DefaultRegionForTheWorld ====== PhotoElectric : Emin= 0 eV Emax= 100 TeV compt: for gamma SubType= 13 Lambda tables from 100 eV to 100 TeV in 240 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Klein-Nishina : Emin= 0 eV Emax= 100 TeV conv: for gamma SubType= 14 Lambda tables from 1.022 MeV to 100 TeV in 240 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bethe-Heitler : Emin= 0 eV Emax= 100 TeV msc: for e- SubType= 10 Lambda tables from 100 eV to 100 TeV in 240 bins, spline: 1 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc93 : Emin= 0 eV Emax= 100 TeV eIoni: for e- SubType= 2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 100 TeV CSDA range table up to 1 GeV in 160 bins eBrem: for e- SubType= 3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 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= 100 TeV eIoni: for e+ SubType= 2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 100 TeV CSDA range table up to 1 GeV in 160 bins eBrem: for e+ SubType= 3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 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= 100 TeV annihil: for e+ SubType= 5 Lambda tables from 100 eV to 100 TeV in 240 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eplus2gg : Emin= 0 eV Emax= 100 TeV msc: for proton SubType= 10 Lambda tables from 100 eV to 100 TeV in 240 bins, spline: 1 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc90 : Emin= 0 eV Emax= 100 TeV hIoni: for proton SubType= 2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 finalRange(mm)= 0.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= 100 TeV CSDA range table up to 1 GeV in 160 bins hBrems: for proton SubType= 3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV hPairProd: for proton SubType= 4 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV msc: for GenericIon SubType= 10 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc90 : Emin= 0 eV Emax= 100 TeV ionIoni: for GenericIon SubType= 2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 Stopping Power data for 17 ion/material pairs, nuclearStopping: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ParamICRU73 : Emin= 0 eV Emax= 100 TeV CSDA range table up to 1 GeV in 160 bins nuclearStopping: for GenericIon SubType= 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 100 TeV nuclearStopping: for alpha SubType= 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 100 TeV hIoni: for anti_proton SubType= 2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 finalRange(mm)= 0.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= 100 TeV CSDA range table up to 1 GeV in 160 bins hIoni: for kaon+ SubType= 2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 finalRange(mm)= 0.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= 100 TeV CSDA range table up to 1 GeV in 160 bins hIoni: for kaon- SubType= 2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 finalRange(mm)= 0.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= 100 TeV CSDA range table up to 1 GeV in 160 bins muMsc: for mu+ SubType= 10 Lambda tables from 100 eV to 100 TeV in 240 bins, spline: 1 RangeFactor= 0.04, step limit type: 2, lateralDisplacement: 1, polarAngleLimit(deg)= 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc90 : Emin= 0 eV Emax= 100 TeV muIoni: for mu+ SubType= 2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 finalRange(mm)= 0.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= 100 TeV CSDA range table up to 1 GeV in 160 bins muBrems: for mu+ SubType= 3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 100 TeV muPairProd: for mu+ SubType= 4 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 100 TeV muIoni: for mu- SubType= 2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 finalRange(mm)= 0.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= 100 TeV CSDA range table up to 1 GeV in 160 bins muBrems: for mu- SubType= 3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 100 TeV muPairProd: for mu- SubType= 4 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 100 TeV hIoni: for pi+ SubType= 2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 finalRange(mm)= 0.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= 100 TeV CSDA range table up to 1 GeV in 160 bins hBrems: for pi+ SubType= 3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV hPairProd: for pi+ SubType= 4 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV msc: for pi- SubType= 10 Lambda tables from 100 eV to 100 TeV in 240 bins, spline: 1 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc90 : Emin= 0 eV Emax= 100 TeV hIoni: for pi- SubType= 2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 finalRange(mm)= 0.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= 100 TeV CSDA range table up to 1 GeV in 160 bins hBrems: for pi- SubType= 3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV hPairProd: for pi- SubType= 4 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV Region -- appears in world volume Root logical volume(s) : Absorber Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] Materials : G4_Si Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : G4_Si Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm Energy thresholds : gamma 6.95018 keV e- 548.291 keV e+ 526.624 keV proton 0 eV Region(s) which use this couple : DefaultRegionForTheWorld ==================================================================== Start closing geometry. G4GeometryManager::ReportVoxelStats -- Voxel Statistics Total memory consumed for geometry optimisation: 0 kByte Total CPU time elapsed for geometry optimisation: 0 seconds ### Run 0 start. --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 9876, 54321 ---------------------------------------- ---> stepMax from HistoManager = 10 um Start Run processing. ---> Begin of Event: 0 ---> Begin of Event: 10000 Run terminated. Run Summary Number of events processed : 20000 User=12.36s Real=13.63s Sys=0.52s ======================== run summary ====================== The run consists of 20000 e- of 500.00 keV through 1.00 mm of G4_Si (density: 2.33 g/cm3 ) ============================================================ Total Energy deposited = 443.666 keV +- 110.109 keV Track length of primary track = 870.965 um +- 199.370 um Range from EmCalculator = 943.428 um (from full dE/dx) Projected range = 321.497 um +- 208.092 um Nb of steps of primary track = 87.62 +- 19.93 Step size= 9.929 um +- 123.724 nm absorbed = 76.69 % transmit = 5.10 % reflected = 18.20 % --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1325532186, 417920133 ---------------------------------------- UserDetectorConstruction deleted. UserPhysicsList deleted. UserRunAction deleted. UserPrimaryGenerator deleted. G4 kernel has come to Quit state. EventManager deleted. UImanager deleted. Units table cleared. StateManager deleted. RunManagerKernel is deleted. RunManager is deleting.