************************************************************* Geant4 version geant4-09-04-patch-01 (18-February-2011) Copyright : Geant4 Collaboration Reference : NIM A 506 (2003), 250-303 WWW : http://cern.ch/geant4 ************************************************************* /run/verbose 1 /tracking/verbose 0 # /random/setSavingFlag true /random/resetEngineFrom currentEvent.rndm RandomNumberEngineStatus restored from file: ./currentEvent.rndm --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 258382564, 915989413 ---------------------------------------- # /testhadr/TargetMat G4_WATER #/testhadr/TargetRadius 10 cm #/testhadr/TargetLength 30 cm #/testhadr/NumberDivZ 300 /testhadr/NumberDivZ 3000 /testhadr/PrintModulo 1000 # /testhadr/CutsAll 0.01 mm /testhadr/Physics PHYSLIST PhysicsList::AddPhysicsList: PhysicsList::AddPhysicsList: /testhadr/Physics emstandard_opt3 PhysicsList::AddPhysicsList: # /run/initialize ### Target consist of 3000 of G4_WATER disks with R(mm)= 100 Width(mm)= 0.1 Total Length(mm)= 300 ### ***** Table : Nb of materials = 3 ***** Material: G4_Al density: 2.699 g/cm3 RadL: 8.896 cm Nucl.Int.Length: 38.877 cm Imean: 166.000 eV ---> Element: Al (Al) Z = 13.0 N = 27.0 A = 26.98 g/mole ---> Isotope: Al27 Z = 13 N = 27 A = 26.98 g/mole abundance: 100.00 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_Galactic density: 0.000 kg/m3 RadL: 204310101.835 pc Nucl.Int.Length: 114327973.206 pc Imean: 21.800 eV temperature: 2.73 K pressure: 0.00 atm ---> 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: 100.00 % ElmAbundance 100.00 % Material: G4_WATER H_2O density: 1.000 g/cm3 RadL: 36.083 cm Nucl.Int.Length: 75.517 cm Imean: 78.000 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: 11.19 % ElmAbundance 66.67 % ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.76 % ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.04 % ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 88.81 % ElmAbundance 33.33 % ### Adding tracking cuts for neutron TimeCut(ns)= 10000 KinEnergyCut(MeV)= 0 PhysicsList::SetCuts:CutLength : 700 um # /gun/particle ion /gun/ion 6 12 6 /gun/energy 1200. MeV /run/beamOn 10000 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 ====== KleinNishina : 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 ====== BetheHeitler : Emin= 0 eV Emax= 10 TeV 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 AngularGenUrban eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 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 AngularGenUrban eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 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 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ParamICRU73 : Emin= 0 eV Emax= 10 TeV nuclearStopping: for C12[0.0] 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 msc: for kaon+ 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 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.505 keV BetheBloch : Emin= 297.505 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.505 keV BetheBloch : Emin= 297.505 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 ============================================================================================ HADRONIC PROCESSES SUMMARY (verbose level 1) Hadronic Processes for CHIPS_Inelastic Hadronic Processes for CHIPS_Inelastic Hadronic Processes for ionInelastic Models: Binary Cascade: Emin(GeV)= 0 Emax(GeV)= 20 Hadronic Processes for ElectroNuclear Models: CHIPS: Emin(GeV)= 0 Emax(GeV)= 30000 Hadronic Processes for PhotonInelastic Models: CHIPS: Emin(GeV)= 0 Emax(GeV)= 3.5 TheoFSGenerator: Emin(GeV)= 3 Emax(GeV)= 100000 Hadronic Processes for hInelastic Models: FTFP: Emin(GeV)= 4 Emax(GeV)= 100000 Bertini Cascade: Emin(GeV)= 0 Emax(GeV)= 5 Hadronic Processes for hInelastic Models: FTFP: Emin(GeV)= 4 Emax(GeV)= 100000 Bertini Cascade: Emin(GeV)= 0 Emax(GeV)= 5 Hadronic Processes for hInelastic Models: FTFP: Emin(GeV)= 4 Emax(GeV)= 100000 Bertini Cascade: Emin(GeV)= 0 Emax(GeV)= 5 Hadronic Processes for muMinusCaptureAtRest Hadronic Processes for hInelastic Models: QGSP: Emin(GeV)= 12.5 Emax(GeV)= 100000 FTFP: Emin(GeV)= 4 Emax(GeV)= 25 Bertini Cascade: Emin(GeV)= 1 Emax(GeV)= 5 Binary Cascade: Emin(GeV)= 0 Emax(GeV)= 1.5 nCapture Models: nRadCapture: Emin(GeV)= 0 Emax(GeV)= 100000 Hadronic Processes for hInelastic Models: QGSP: Emin(GeV)= 12.5 Emax(GeV)= 100000 FTFP: Emin(GeV)= 4 Emax(GeV)= 25 Bertini Cascade: Emin(GeV)= 0 Emax(GeV)= 5 Hadronic Processes for hInelastic Models: QGSP: Emin(GeV)= 12.5 Emax(GeV)= 100000 FTFP: Emin(GeV)= 4 Emax(GeV)= 25 Bertini Cascade: Emin(GeV)= 0 Emax(GeV)= 5 Hadronic Processes for hInelastic Models: QGSP: Emin(GeV)= 12.5 Emax(GeV)= 100000 FTFP: Emin(GeV)= 4 Emax(GeV)= 25 Bertini Cascade: Emin(GeV)= 1 Emax(GeV)= 5 Binary Cascade: Emin(GeV)= 0 Emax(GeV)= 1.5 ============================================================================================ ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : G4_Galactic Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 1 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes recalculation needed : No Material : G4_WATER Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um Energy thresholds : gamma 990 eV e- 14.0874 keV e+ 14.0874 keV proton 1 keV Region(s) which use this couple : DefaultRegionForTheWorld ==================================================================== ### Run 0 start ### Histo books 1 histograms in [Histo::book] File created: test44.root Start Run processing. EventAction: Event # 0 started EventAction: Event # 1000 started EventAction: Event # 2000 started G4Fragment::CalculateExcitationEnergy(): WARNING Fragment: A = 2, Z = 2, U = -2.850e+01 MeV P = (0.000e+00,-1.421e-14,1.722e+02) MeV E = 1.892e+03 MeV EventAction: Event # 3000 started EventAction: Event # 4000 started EventAction: Event # 5000 started EventAction: Event # 6000 started EventAction: Event # 7000 started EventAction: Event # 8000 started EventAction: Event # 9000 started Run terminated. Run Summary Number of events processed : 10000 User=1961.56s Real=2051.87s Sys=3.3s RunAction: End of run actions are started HistoManager: End of run actions are started ======================================================== Beam particle C12[0.0] Beam Energy(MeV) 1200 Number of events 10000 Average energy deposit (MeV) 1185 RMS(MeV) 54.94 Average number of steps of primary 4941 ======================================================== # G4 kernel has come to Quit state.