############################################ !!! 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 2 # /testem/det/setAbsMat Silicon /testem/det/setAbsThick 20 um /testem/det/setAbsYZ 50 um # /testem/phys/addPhysics emstandard_opt3 PhysicsList::AddPhysicsList: # /testem/phys/setCuts 7 um # /run/initialize userDetector->Construct() start. Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 113804112.800 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Silicon density: 2.330 g/cm3 RadL: 9.368 cm Nucl.Int.Length: 45.663 cm Imean: 173.000 eV ---> Element: Silicon ( ) Z = 14.0 N = 28.1 A = 28.09 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % The WORLD is made of 30 um of Galactic. The transverse size (YZ) of the world is 60 um The ABSORBER is made of 20 um of Silicon. The transverse size (YZ) is 50 um X position of the middle of the absorber 0 fm World is registered to the default region. physicsList->Construct() start. physicsList->Construct() start. physicsList->CheckParticleList() start. physicsList->setCut() start. PhysicsList::SetCuts:CutLength : 1 mm # /testem/gun/setDefault /gun/particle pi+ /gun/energy 5 GeV # /testem/histo/setFileName testem5 /testem/histo/setFileType root /testem/histo/setHisto 13 100 -0.25 0.25 mrad ----> SetHisto 13: (transmit, charged) : projected angle at exit (mrad); 100 bins from -0.25 mrad to 0.25 mrad # /run/beamOn 100000 phot: for gamma SubType= 12 ===== EM models for the G4Region DefaultRegionForTheWorld ====== PhotoElectric : Emin= 0 eV Emax= 10 TeV compt: for gamma SubType= 13 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Klein-Nishina : Emin= 0 eV Emax= 10 TeV conv: for gamma SubType= 14 Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bethe-Heitler : Emin= 0 eV Emax= 10 TeV Rayl: for gamma SubType= 11 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV msc: for e- SubType= 10 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5 ===== 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 220 bins Lambda tables from threshold to 10 TeV in 220 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= 10 TeV eBrem: for e- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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= 10 TeV eBrem: for e+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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+ SubType= 5 Lambda tables from 100 eV to 10 TeV in 220 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 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.01, 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 ====== ParamICRU73 : Emin= 0 eV Emax= 10 TeV nuclearStopping: for GenericIon SubType= 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV nuclearStopping: for alpha SubType= 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV hIoni: for anti_proton SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins, spline: 1 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 10 TeV muIoni: for mu+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 10 TeV CoulombScat: for mu+ SubType= 1 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 10 TeV muIoni: for mu- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 10 TeV CoulombScat: for mu- SubType= 1 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 10 TeV hIoni: for pi+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 10 TeV Region -- appears in world volume Root logical volume(s) : World Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] Materials : Galactic Silicon Production cuts : gamma 7 um e- 7 um e+ 7 um proton 0 fm ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : Galactic Range cuts : gamma 7 um e- 7 um e+ 7 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 : Silicon Range cuts : gamma 7 um e- 7 um e+ 7 um proton 0 fm Energy thresholds : gamma 990 eV e- 21.4886 keV e+ 20.8954 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 ---------------------------------------- Start Run processing. ---> Begin of Event: 0 ---> Begin of Event: 10000 ---> Begin of Event: 20000 ---> Begin of Event: 30000 ---> Begin of Event: 40000 ---> Begin of Event: 50000 ---> Begin of Event: 60000 ---> Begin of Event: 70000 ---> Begin of Event: 80000 ---> Begin of Event: 90000 Run terminated. Run Summary Number of events processed : 100000 User=4.22s Real=11.65s Sys=2.4s ======================== run summary ====================== The run was 100000 pi+ of 5 GeV through 20 um of Silicon (density: 2.33 g/cm3 ) Total energy deposit in absorber per event = 6.329 keV +- 15.51 eV -----> Mean dE/dx = 3.165 MeV/cm (1.358 MeV*cm2/g) From formulas : restricted dEdx = 2.963 MeV/cm (1.272 MeV*cm2/g) full dEdx = 4.726 MeV/cm (2.028 MeV*cm2/g) Leakage : primary = 5 GeV +- 294.7 eV secondaries = 1.778 keV +- 293.3 eV Energy balance : edep + eleak = 5 GeV Total track length (charged) in absorber per event = 20.26 um +- 7.916 nm Total track length (neutral) in absorber per event = 2.089 Ang +- 1.085 Ang Number of steps (charged) in absorber per event = 1.156 +- 8.897e-06 Number of steps (neutral) in absorber per event = 6e-05 +- 8.897e-06 Number of secondaries per event : Gammas = 6e-05; electrons = 0.01646; positrons = 0 Number of events with the primary particle transmitted = 100 % Number of events with at least 1 particle transmitted (same charge as primary) = 100 % Number of events with the primary particle reflected = 0 % Number of events with at least 1 particle reflected (same charge as primary) = 0 % MultipleScattering: rms proj angle of transmit primary particle = 0.02452 mrad (central part only) computed theta0 (Highland formula) = 0.02627 mrad central part defined as +- 0.0788 mrad; Tail ratio = 1.573 % --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 206812152, 1002169913 ---------------------------------------- # /testem/det/setAbsMat Gold /testem/det/setAbsThick 9.658 um /testem/det/setAbsYZ 1 mm /testem/det/update Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 113804112.800 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Gold density: 19.320 g/cm3 RadL: 3.344 mm Nucl.Int.Length: 10.540 cm Imean: 790.000 eV ---> Element: Gold ( ) Z = 79.0 N = 197.0 A = 196.97 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % The WORLD is made of 14.487 um of Galactic. The transverse size (YZ) of the world is 1.2 mm The ABSORBER is made of 9.658 um of Gold. The transverse size (YZ) is 1 mm X position of the middle of the absorber 0 fm Obsolete world logical volume is removed from the default region. World is registered to the default region. /testem/gun/setDefault /gun/particle e- /gun/energy 15.7 MeV /run/beamOn 50000 phot: for gamma SubType= 12 ===== EM models for the G4Region DefaultRegionForTheWorld ====== PhotoElectric : Emin= 0 eV Emax= 10 TeV compt: for gamma SubType= 13 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Klein-Nishina : Emin= 0 eV Emax= 10 TeV conv: for gamma SubType= 14 Lambda tables from 1.0219978 MeV to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bethe-Heitler : Emin= 0 eV Emax= 10 TeV Rayl: for gamma SubType= 11 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV msc: for e- SubType= 10 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5 ===== 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 220 bins Lambda tables from threshold to 10 TeV in 220 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= 10 TeV eBrem: for e- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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= 10 TeV eBrem: for e+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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+ SubType= 5 Lambda tables from 100 eV to 10 TeV in 220 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 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.01, 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 ====== ParamICRU73 : Emin= 0 eV Emax= 10 TeV nuclearStopping: for GenericIon SubType= 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV nuclearStopping: for alpha SubType= 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV hIoni: for anti_proton SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 1.052311 MeV BetheBloch : Emin= 1.052311 MeV Emax= 10 TeV hBrems: for kaon+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 1.052311 MeV BetheBloch : Emin= 1.052311 MeV Emax= 10 TeV hBrems: for kaon- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins, spline: 1 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.459156 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 10 TeV muIoni: for mu+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 10 TeV CoulombScat: for mu+ SubType= 1 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 11.459156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 10 TeV muIoni: for mu- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 10 TeV CoulombScat: for mu- SubType= 1 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 11.459156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 10 TeV hIoni: for pi+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 297.50434 keV BetheBloch : Emin= 297.50434 keV Emax= 10 TeV hBrems: for pi+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 297.50434 keV BetheBloch : Emin= 297.50434 keV Emax= 10 TeV hBrems: for pi- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 10 TeV Region -- appears in world volume Root logical volume(s) : World Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] Materials : Galactic Gold Production cuts : gamma 7 um e- 7 um e+ 7 um proton 0 fm ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : Galactic Range cuts : gamma 7 um e- 7 um e+ 7 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 : No recalculation needed : No Material : Silicon Range cuts : gamma 7 um e- 7 um e+ 7 um proton 0 fm Energy thresholds : gamma 990 eV e- 21.488579 keV e+ 20.895427 keV proton 0 eV Index : 2 used in the geometry : Yes recalculation needed : No Material : Gold Range cuts : gamma 7 um e- 7 um e+ 7 um proton 0 fm Energy thresholds : gamma 6.5422867 keV e- 65.140571 keV e+ 63.840474 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 1 start. --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 206812152, 1002169913 ---------------------------------------- Start Run processing. ---> Begin of Event: 0 ---> Begin of Event: 10000 ---> Begin of Event: 20000 ---> Begin of Event: 30000 ---> Begin of Event: 40000 Run terminated. Run Summary Number of events processed : 50000 User=3.27s Real=4.5s Sys=1.22s ======================== run summary ====================== The run was 50000 e- of 15.7 MeV through 9.66 um of Gold (density: 19.3 g/cm3 ) Total energy deposit in absorber per event = 18.6 keV +- 64.68 eV -----> Mean dE/dx = 19.26 MeV/cm (0.9971 MeV*cm2/g) From formulas : restricted dEdx = 17.9 MeV/cm (0.9266 MeV*cm2/g) full dEdx = 59.7 MeV/cm (3.09 MeV*cm2/g) Leakage : primary = 15.64 MeV +- 2.154 keV secondaries = 40.1 keV +- 2.145 keV Energy balance : edep + eleak = 15.7 MeV Total track length (charged) in absorber per event = 9.95 um +- 12.57 nm Total track length (neutral) in absorber per event = 133.2 nm +- 6.153 nm Number of steps (charged) in absorber per event = 5.071 +- 1.585e-05 Number of steps (neutral) in absorber per event = 0.02956 +- 1.585e-05 Number of secondaries per event : Gammas = 0.02944; electrons = 0.02198; positrons = 0 Number of events with the primary particle transmitted = 99.97 % Number of events with at least 1 particle transmitted (same charge as primary) = 99.97 % Number of events with the primary particle reflected = 0.028 % Number of events with at least 1 particle reflected (same charge as primary) = 0.276 % MultipleScattering: rms proj angle of transmit primary particle = 32.92 mrad (central part only) computed theta0 (Highland formula) = 35.1 mrad central part defined as +- 105.3 mrad; Tail ratio = 2.702 % --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1822484465, 1275230702 ---------------------------------------- # /testem/phys/setCuts 10 um # /testem/det/setAbsMat G4_POLYETHYLENE /testem/det/setAbsThick 4.74 mm /testem/det/setAbsYZ 50 mm /testem/det/update Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 113804112.800 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_POLYETHYLENE (C_2H_4)_N-Polyethylene density: 940.000 mg/cm3 RadL: 47.632 cm Nucl.Int.Length: 72.090 cm Imean: 57.400 eV ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 14.37 % ElmAbundance 66.67 % ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 85.63 % ElmAbundance 33.33 % The WORLD is made of 7.11 mm of Galactic. The transverse size (YZ) of the world is 6 cm The ABSORBER is made of 4.74 mm of G4_POLYETHYLENE. The transverse size (YZ) is 5 cm X position of the middle of the absorber 0 fm Obsolete world logical volume is removed from the default region. World is registered to the default region. /testem/gun/setDefault /gun/particle mu+ /gun/energy 96.2 MeV /run/beamOn 50000 phot: for gamma SubType= 12 ===== EM models for the G4Region DefaultRegionForTheWorld ====== PhotoElectric : Emin= 0 eV Emax= 10 TeV compt: for gamma SubType= 13 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Klein-Nishina : Emin= 0 eV Emax= 10 TeV conv: for gamma SubType= 14 Lambda tables from 1.0219978 MeV to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bethe-Heitler : Emin= 0 eV Emax= 10 TeV Rayl: for gamma SubType= 11 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV msc: for e- SubType= 10 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5 ===== 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 220 bins Lambda tables from threshold to 10 TeV in 220 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= 10 TeV eBrem: for e- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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= 10 TeV eBrem: for e+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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+ SubType= 5 Lambda tables from 100 eV to 10 TeV in 220 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 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.01, 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 ====== ParamICRU73 : Emin= 0 eV Emax= 10 TeV nuclearStopping: for GenericIon SubType= 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV nuclearStopping: for alpha SubType= 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV hIoni: for anti_proton SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 1.052311 MeV BetheBloch : Emin= 1.052311 MeV Emax= 10 TeV hBrems: for kaon+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 1.052311 MeV BetheBloch : Emin= 1.052311 MeV Emax= 10 TeV hBrems: for kaon- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins, spline: 1 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.459156 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 10 TeV muIoni: for mu+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 10 TeV CoulombScat: for mu+ SubType= 1 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 11.459156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 10 TeV muIoni: for mu- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 10 TeV CoulombScat: for mu- SubType= 1 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 11.459156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 10 TeV hIoni: for pi+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 297.50434 keV BetheBloch : Emin= 297.50434 keV Emax= 10 TeV hBrems: for pi+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 297.50434 keV BetheBloch : Emin= 297.50434 keV Emax= 10 TeV hBrems: for pi- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 10 TeV Region -- appears in world volume Root logical volume(s) : World Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] Materials : Galactic G4_POLYETHYLENE Production cuts : gamma 10 um e- 10 um e+ 10 um proton 0 fm ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : Galactic Range cuts : gamma 10 um e- 10 um e+ 10 um 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 : No recalculation needed : No Material : Silicon Range cuts : gamma 10 um e- 10 um e+ 10 um proton 0 fm Energy thresholds : gamma 990 eV e- 31.951563 keV e+ 31.737567 keV proton 0 eV Index : 2 used in the geometry : No recalculation needed : No Material : Gold Range cuts : gamma 10 um e- 10 um e+ 10 um proton 0 fm Energy thresholds : gamma 8.2770687 keV e- 80.768699 keV e+ 79.156691 keV proton 0 eV Index : 3 used in the geometry : Yes recalculation needed : No Material : G4_POLYETHYLENE Range cuts : gamma 10 um e- 10 um e+ 10 um proton 0 fm Energy thresholds : gamma 990 eV e- 14.174755 keV e+ 13.754766 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 2 start. --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1822484465, 1275230702 ---------------------------------------- Start Run processing. ---> Begin of Event: 0 ---> Begin of Event: 10000 ---> Begin of Event: 20000 ---> Begin of Event: 30000 ---> Begin of Event: 40000 Run terminated. Run Summary Number of events processed : 50000 User=6.73s Real=9.07s Sys=1.29s ======================== run summary ====================== The run was 50000 mu+ of 96.2 MeV through 4.74 mm of G4_POLYETHYLENE (density: 940 mg/cm3) Total energy deposit in absorber per event = 1.033 MeV +- 1.01 keV -----> Mean dE/dx = 2.18 MeV/cm (2.319 MeV*cm2/g) From formulas : restricted dEdx = 1.751 MeV/cm (1.863 MeV*cm2/g) full dEdx = 2.263 MeV/cm (2.407 MeV*cm2/g) Leakage : primary = 95.14 MeV +- 2.335 keV secondaries = 30.95 keV +- 3.642 keV Energy balance : edep + eleak = 96.2 MeV Total track length (charged) in absorber per event = 5.106 mm +- 3.914 um Total track length (neutral) in absorber per event = 31.28 um +- 2.947 um Number of steps (charged) in absorber per event = 11.22 +- 0.0002798 Number of steps (neutral) in absorber per event = 0.01258 +- 0.0002798 Number of secondaries per event : Gammas = 0.01166; electrons = 3.988; positrons = 2e-05 Number of events with the primary particle transmitted = 100 % Number of events with at least 1 particle transmitted (same charge as primary) = 100 % Number of events with the primary particle reflected = 0 % Number of events with at least 1 particle reflected (same charge as primary) = 0 % MultipleScattering: rms proj angle of transmit primary particle = 7.29 mrad (central part only) computed theta0 (Highland formula) = 7.636 mrad central part defined as +- 22.91 mrad; Tail ratio = 1.284 % --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1953603820, 877017596 ---------------------------------------- # /testem/phys/setCuts 1 mm # /testem/stack/killSecondaries # /testem/det/setAbsMat Aluminium /testem/det/setAbsThick 8.004 mm /testem/det/update Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 113804112.800 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Nucl.Int.Length: 38.879 cm Imean: 166.000 eV ---> Element: Aluminium ( ) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % The WORLD is made of 1.2006 cm of Galactic. The transverse size (YZ) of the world is 6 cm The ABSORBER is made of 8.004 mm of Aluminium. The transverse size (YZ) is 5 cm X position of the middle of the absorber 0 fm Obsolete world logical volume is removed from the default region. World is registered to the default region. /testem/gun/setDefault /gun/particle proton /gun/energy 174.664 GeV /run/beamOn 50000 phot: for gamma SubType= 12 ===== EM models for the G4Region DefaultRegionForTheWorld ====== PhotoElectric : Emin= 0 eV Emax= 10 TeV compt: for gamma SubType= 13 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Klein-Nishina : Emin= 0 eV Emax= 10 TeV conv: for gamma SubType= 14 Lambda tables from 1.0219978 MeV to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bethe-Heitler : Emin= 0 eV Emax= 10 TeV Rayl: for gamma SubType= 11 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV msc: for e- SubType= 10 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5 ===== 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 220 bins Lambda tables from threshold to 10 TeV in 220 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= 10 TeV eBrem: for e- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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= 10 TeV eBrem: for e+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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+ SubType= 5 Lambda tables from 100 eV to 10 TeV in 220 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 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.01, 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 ====== ParamICRU73 : Emin= 0 eV Emax= 10 TeV nuclearStopping: for GenericIon SubType= 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV nuclearStopping: for alpha SubType= 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV hIoni: for anti_proton SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 1.052311 MeV BetheBloch : Emin= 1.052311 MeV Emax= 10 TeV hBrems: for kaon+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 1.052311 MeV BetheBloch : Emin= 1.052311 MeV Emax= 10 TeV hBrems: for kaon- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins, spline: 1 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.459156 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 10 TeV muIoni: for mu+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 10 TeV CoulombScat: for mu+ SubType= 1 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 11.459156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 10 TeV muIoni: for mu- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 10 TeV CoulombScat: for mu- SubType= 1 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 11.459156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 10 TeV hIoni: for pi+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 297.50434 keV BetheBloch : Emin= 297.50434 keV Emax= 10 TeV hBrems: for pi+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 297.50434 keV BetheBloch : Emin= 297.50434 keV Emax= 10 TeV hBrems: for pi- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 10 TeV Region -- appears in world volume Root logical volume(s) : World Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] Materials : Galactic Aluminium 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 : Galactic Range cuts : gamma 1 mm e- 1 mm e+ 1 mm 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 : No recalculation needed : No Material : Silicon Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm Energy thresholds : gamma 6.9501782 keV e- 548.29137 keV e+ 526.62382 keV proton 0 eV Index : 2 used in the geometry : No recalculation needed : No Material : Gold Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm Energy thresholds : gamma 112.64401 keV e- 2.3097622 MeV e+ 2.1308135 MeV proton 0 eV Index : 3 used in the geometry : No recalculation needed : No Material : G4_POLYETHYLENE Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm Energy thresholds : gamma 2.2626058 keV e- 351.8774 keV e+ 342.54485 keV proton 0 eV Index : 4 used in the geometry : Yes recalculation needed : No Material : Aluminium Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm Energy thresholds : gamma 6.9036292 keV e- 598.34507 keV e+ 570.85042 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 3 start. --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1953603820, 877017596 ---------------------------------------- Start Run processing. ---> Begin of Event: 0 ---> Begin of Event: 10000 ---> Begin of Event: 20000 ---> Begin of Event: 30000 ---> Begin of Event: 40000 Run terminated. Run Summary Number of events processed : 50000 User=2.22s Real=3.49s Sys=1.25s ======================== run summary ====================== The run was 50000 proton of 175 GeV through 8 mm of Aluminium (density: 2.7 g/cm3 ) Total energy deposit in absorber per event = 4.999 MeV +- 290.4 keV -----> Mean dE/dx = 6.245 MeV/cm (2.313 MeV*cm2/g) From formulas : restricted dEdx = 4.024 MeV/cm (1.491 MeV*cm2/g) full dEdx = 6.049 MeV/cm (2.24 MeV*cm2/g) Leakage : primary = 174.7 GeV +- 290.4 keV secondaries = 0 eV +- 0 eV Energy balance : edep + eleak = 174.7 GeV Total track length (charged) in absorber per event = 8.004 mm +- 0 fm Total track length (neutral) in absorber per event = 0 fm +- 0 fm Number of steps (charged) in absorber per event = 1.268 +- 0 Number of steps (neutral) in absorber per event = 0 +- 0 Number of secondaries per event : Gammas = 0; electrons = 0.268; positrons = 0.0009 Number of events with the primary particle transmitted = 100 % Number of events with at least 1 particle transmitted (same charge as primary) = 100 % Number of events with the primary particle reflected = 0 % Number of events with at least 1 particle reflected (same charge as primary) = 0 % MultipleScattering: rms proj angle of transmit primary particle = 0.02108 mrad (central part only) computed theta0 (Highland formula) = 0.02111 mrad central part defined as +- 0.06333 mrad; Tail ratio = 1.836 % --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1162542043, 318101973 ---------------------------------------- # /testem/det/setAbsMat Iron /testem/det/setAbsThick 1 m /testem/det/update Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 113804112.800 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Iron density: 7.870 g/cm3 RadL: 1.759 cm Nucl.Int.Length: 16.999 cm Imean: 286.000 eV ---> Element: Iron ( ) Z = 26.0 N = 55.8 A = 55.85 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % The WORLD is made of 1.5 m of Galactic. The transverse size (YZ) of the world is 6 cm The ABSORBER is made of 1 m of Iron. The transverse size (YZ) is 5 cm X position of the middle of the absorber 0 fm Obsolete world logical volume is removed from the default region. World is registered to the default region. /testem/gun/setDefault /gun/particle mu+ /gun/energy 100 GeV /run/beamOn 50000 phot: for gamma SubType= 12 ===== EM models for the G4Region DefaultRegionForTheWorld ====== PhotoElectric : Emin= 0 eV Emax= 10 TeV compt: for gamma SubType= 13 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Klein-Nishina : Emin= 0 eV Emax= 10 TeV conv: for gamma SubType= 14 Lambda tables from 1.0219978 MeV to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bethe-Heitler : Emin= 0 eV Emax= 10 TeV Rayl: for gamma SubType= 11 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV msc: for e- SubType= 10 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5 ===== 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 220 bins Lambda tables from threshold to 10 TeV in 220 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= 10 TeV eBrem: for e- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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= 10 TeV eBrem: for e+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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+ SubType= 5 Lambda tables from 100 eV to 10 TeV in 220 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 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.01, 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 ====== ParamICRU73 : Emin= 0 eV Emax= 10 TeV nuclearStopping: for GenericIon SubType= 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV nuclearStopping: for alpha SubType= 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV hIoni: for anti_proton SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 1.052311 MeV BetheBloch : Emin= 1.052311 MeV Emax= 10 TeV hBrems: for kaon+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 1.052311 MeV BetheBloch : Emin= 1.052311 MeV Emax= 10 TeV hBrems: for kaon- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins, spline: 1 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.459156 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 10 TeV muIoni: for mu+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 10 TeV CoulombScat: for mu+ SubType= 1 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 11.459156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 10 TeV muIoni: for mu- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 10 TeV CoulombScat: for mu- SubType= 1 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1 11.459156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 10 TeV hIoni: for pi+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 297.50434 keV BetheBloch : Emin= 297.50434 keV Emax= 10 TeV hBrems: for pi+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 297.50434 keV BetheBloch : Emin= 297.50434 keV Emax= 10 TeV hBrems: for pi- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 220 bins Lambda tables from threshold to 10 TeV in 220 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 220 bins Lambda tables from threshold to 10 TeV in 220 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 10 TeV Region -- appears in world volume Root logical volume(s) : World Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] Materials : Galactic Iron 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 : Galactic Range cuts : gamma 1 mm e- 1 mm e+ 1 mm 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 : No recalculation needed : No Material : Silicon Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm Energy thresholds : gamma 6.9501782 keV e- 548.29137 keV e+ 526.62382 keV proton 0 eV Index : 2 used in the geometry : No recalculation needed : No Material : Gold Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm Energy thresholds : gamma 112.64401 keV e- 2.3097622 MeV e+ 2.1308135 MeV proton 0 eV Index : 3 used in the geometry : No recalculation needed : No Material : G4_POLYETHYLENE Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm Energy thresholds : gamma 2.2626058 keV e- 351.8774 keV e+ 342.54485 keV proton 0 eV Index : 4 used in the geometry : No recalculation needed : No Material : Aluminium Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm Energy thresholds : gamma 6.9036292 keV e- 598.34507 keV e+ 570.85042 keV proton 0 eV Index : 5 used in the geometry : Yes recalculation needed : No Material : Iron Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm Energy thresholds : gamma 20.643844 keV e- 1.295918 MeV e+ 1.2116932 MeV 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 4 start. --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1162542043, 318101973 ---------------------------------------- Start Run processing. ---> Begin of Event: 0 ---> Begin of Event: 10000 ---> Begin of Event: 20000 ---> Begin of Event: 30000 ---> Begin of Event: 40000 Run terminated. Run Summary Number of events processed : 50000 User=26.38s Real=44.83s Sys=1.28s ======================== run summary ====================== The run was 50000 mu+ of 100 GeV through 1 m of Iron (density: 7.87 g/cm3 ) Total energy deposit in absorber per event = 2.126 GeV +- 14.13 MeV -----> Mean dE/dx = 21.26 MeV/cm (2.702 MeV*cm2/g) From formulas : restricted dEdx = 11.22 MeV/cm (1.425 MeV*cm2/g) full dEdx = 21.6 MeV/cm (2.745 MeV*cm2/g) Leakage : primary = 97.87 GeV +- 14.15 MeV secondaries = 0 eV +- 0 eV Energy balance : edep + eleak = 100 GeV Total track length (charged) in absorber per event = 99 cm +- 52.13 um Total track length (neutral) in absorber per event = 0 fm +- 0 fm Number of steps (charged) in absorber per event = 47.44 +- 0.001021 Number of steps (neutral) in absorber per event = 0 +- 0 Number of secondaries per event : Gammas = 0.04906; electrons = 45.84; positrons = 1.251 Number of events with the primary particle transmitted = 100 % Number of events with at least 1 particle transmitted (same charge as primary) = 100 % Number of events with the primary particle reflected = 0 % Number of events with at least 1 particle reflected (same charge as primary) = 0 % MultipleScattering: rms proj angle of transmit primary particle = 1.121 mrad (central part only) computed theta0 (Highland formula) = 1.182 mrad central part defined as +- 3.545 mrad; Tail ratio = 0.289 % --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1305857416, 170091693 ---------------------------------------- # 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.