************************************************************* Geant4 version Name: global-V09-00-03 (9-May-2008) Copyright : Geant4 Collaboration Reference : NIM A 506 (2003), 250-303 WWW : http://cern.ch/geant4 ************************************************************* ***** Table : Nb of materials = 15 ***** Material: Air density: 1.290 mg/cm3 RadL: 285.161 m Imean: 85.684 eV temperature: 273.15 K pressure: 1.00 atm ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 70.00 % ElmAbundance 72.71 % ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 30.00 % ElmAbundance 27.29 % Material: H2liquid density: 70.800 mg/cm3 RadL: 8.923 m Imean: 21.800 eV ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Water H_2O density: 1.000 g/cm3 RadL: 36.092 cm Imean: 75.000 eV ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.21 % ElmAbundance 66.67 % ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.79 % ElmAbundance 33.33 % Material: liquidArgon density: 1.390 g/cm3 RadL: 14.065 cm Imean: 187.200 eV ---> Element: liquidArgon ( ) Z = 18.0 N = 40.0 A = 39.95 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Carbon density: 2.267 g/cm3 RadL: 18.833 cm Imean: 78.000 eV ---> Element: Carbon ( ) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Imean: 166.400 eV ---> Element: Aluminium ( ) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Silicon density: 2.330 g/cm3 RadL: 9.368 cm Imean: 173.600 eV ---> Element: Silicon ( ) Z = 14.0 N = 28.1 A = 28.09 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Germanium density: 5.323 g/cm3 RadL: 2.301 cm Imean: 332.800 eV ---> Element: Germanium ( ) Z = 32.0 N = 72.6 A = 72.61 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: NaI density: 3.670 g/cm3 RadL: 2.586 cm Imean: 452.000 eV ---> Element: Sodium (Na) Z = 11.0 N = 23.0 A = 22.99 g/mole ElmMassFraction: 15.34 % ElmAbundance 50.00 % ---> Element: Iodine (I) Z = 53.0 N = 126.9 A = 126.90 g/mole ElmMassFraction: 84.66 % ElmAbundance 50.00 % Material: Iodine density: 4.930 g/cm3 RadL: 1.720 cm Imean: 492.900 eV ---> Element: Iodine (I) Z = 53.0 N = 126.9 A = 126.90 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: BGO density: 7.100 g/cm3 RadL: 1.123 cm Imean: 469.142 eV ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 15.41 % ElmAbundance 63.16 % ---> Element: Germanium (Ge) Z = 32.0 N = 72.6 A = 72.59 g/mole ElmMassFraction: 17.48 % ElmAbundance 15.79 % ---> Element: Bismuth (Bi) Z = 83.0 N = 209.0 A = 208.98 g/mole ElmMassFraction: 67.10 % ElmAbundance 21.05 % Material: Iron density: 7.870 g/cm3 RadL: 1.759 cm Imean: 286.000 eV ---> Element: Iron ( ) Z = 26.0 N = 55.9 A = 55.85 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Tungsten density: 19.300 g/cm3 RadL: 3.504 mm Imean: 725.200 eV ---> Element: Tungsten ( ) Z = 74.0 N = 183.8 A = 183.85 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Lead density: 11.350 g/cm3 RadL: 5.612 mm Imean: 820.000 eV ---> Element: Lead ( ) Z = 82.0 N = 207.2 A = 207.19 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Uranium density: 18.950 g/cm3 RadL: 3.166 mm Imean: 892.400 eV ---> Element: Uranium ( ) Z = 92.0 N = 238.0 A = 238.03 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % /run/verbose 2 # /testem/det/setMat Water #### /testem/det/setSize 100 m # #### /testem/phys/addPhysics livermore #### /testem/phys/addPhysics penelope # /run/initialize userDetector->Construct() start. The Box is 100 m of Water Water is registered to the default region. physicsList->Construct() start. physicsList->Construct() start. physicsList->setCut() start. # /gun/particle gamma /gun/energy 100 keV # /run/beamOn 100000 phot: Total cross sections from Sandia parametrisation. Sampling according PhotoElectric model compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV Sampling according Klein-Nishina model tables are built for gamma Lambda tables from 100 eV to 100 GeV in 90 bins. conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; sampling secondary e+e- according Bethe-Heitler model tables are built for gamma Lambda tables from 1.022 MeV to 100 GeV in 100 bins. eIoni: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0 eBrem: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1 eIoni: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0 eBrem: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1 annihil: Sampling according eplus2gg model tables are built for e+ Lambda tables from 100 eV to 100 TeV in 120 bins. hIoni: tables are built for proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. NuclearStopping= 1 Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0 ionIoni: tables are built for GenericIon dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from BraggIon for protons below. NuclearStopping= 1 Stopping Power data for 8 ion/material pairs are used. Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0 hIoni: tables are built for anti_proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. NuclearStopping= 1 Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0 muIoni: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0 muBrems: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model muPairProd: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model muIoni: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0 muBrems: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model muPairProd: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model hIoni: tables are built for pi+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV Parametrisation from Bragg for protons below. NuclearStopping= 1 Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0 hIoni: tables are built for pi- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV Parametrisation from Bragg for protons below. NuclearStopping= 1 Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0 Region -- appears in world volume Root logical volume(s) : Water Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] Materials : Water Production cuts : gamma 1 mm e- 1 mm e+ 1 mm ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : Water Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 2.90186 keV e- 347.138 keV e+ 338.695 keV 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=1.4s Real=1.4s Sys=0s The run consists of 100000 gamma of 100 keV through 100 m of Water (density: 1 g/cm3 ) Process calls frequency ---> compt = 98335 phot = 1665 MeanFreePath: 6.0496 cm +- 6.0503 cm massic: 6.0496 g/cm2 CrossSection: 0.1653 cm^-1 massic: 0.1653 cm2/g mean energy of charged secondaries: 15.209 keV mass_energy_transfer coef: 0.025141 cm2/g Verification : crossSections from G4EmCalculator compt= 0.1623 cm2/g phot= 0.0026853 cm2/g total= 0.16498 cm2/g --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1982964064, 678286457 ---------------------------------------- # /gun/particle e- /gun/energy 100 MeV # /run/beamOn 100000 Region -- appears in world volume Root logical volume(s) : Water Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] Materials : Water Production cuts : gamma 1 mm e- 1 mm e+ 1 mm ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : Water Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 2.9 keV e- 3.5e+02 keV e+ 3.4e+02 keV Region(s) which use this couple : DefaultRegionForTheWorld ==================================================================== ### Run 1 start. --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1982964064, 678286457 ---------------------------------------- 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=1.4s Real=1.4s Sys=0s The run consists of 100000 e- of 100 MeV through 100 m of Water (density: 1 g/cm3 ) Process calls frequency ---> eBrem = 61697 eIoni = 38303 MeanFreePath: 1.5679 cm +- 1.5695 cm massic: 1.5679 g/cm2 CrossSection: 0.63781 cm^-1 massic: 0.63781 cm2/g mean energy of charged secondaries: 705.83 keV mass_energy_transfer coef: 0.0045019 cm2/g Verification : crossSections from G4EmCalculator eBrem= 0.39397 cm2/g eIoni= 0.24497 cm2/g total= 0.63894 cm2/g --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 679237908, 105871640 ---------------------------------------- UserDetectorConstruction deleted. UserPhysicsList deleted. UserRunAction deleted. UserPrimaryGenerator deleted. G4 kernel has come to Quit state. EventManager deleted. Default detector region deleted. UImanager deleted. Units table cleared. StateManager deleted. RunManagerKernel is deleted. RunManager is deleting.