************************************************************* 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 = 8 ***** Material: G4_Al density: 2.699 g/cm3 RadL: 8.896 cm Imean: 166.000 eV ---> Element: Al (Al) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_Si density: 2.330 g/cm3 RadL: 9.366 cm Imean: 173.000 eV ---> Element: Si (Si) Z = 14.0 N = 28.1 A = 28.09 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_Fe density: 7.874 g/cm3 RadL: 1.757 cm Imean: 286.000 eV ---> Element: Fe (Fe) Z = 26.0 N = 55.8 A = 55.85 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_Ge density: 5.323 g/cm3 RadL: 2.301 cm Imean: 350.000 eV ---> Element: Ge (Ge) Z = 32.0 N = 72.6 A = 72.61 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_W density: 19.300 g/cm3 RadL: 3.504 mm Imean: 727.000 eV ---> Element: W (W) Z = 74.0 N = 183.8 A = 183.84 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_Pb density: 11.350 g/cm3 RadL: 5.613 mm Imean: 823.000 eV ---> Element: Pb (Pb) Z = 82.0 N = 207.2 A = 207.22 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: G4_AIR density: 1.205 mg/cm3 RadL: 303.921 m Imean: 85.700 eV temperature: 273.15 K pressure: 1.00 atm ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 0.01 % ElmAbundance 0.02 % ---> Element: N (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 75.53 % ElmAbundance 78.44 % ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 23.18 % ElmAbundance 21.07 % ---> Element: Ar (Ar) Z = 18.0 N = 39.9 A = 39.95 g/mole ElmMassFraction: 1.28 % ElmAbundance 0.47 % Material: G4_WATER H_2O density: 1.000 g/cm3 RadL: 36.083 cm Imean: 75.000 eV ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.19 % ElmAbundance 66.67 % ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.81 % ElmAbundance 33.33 % Visualization Manager instantiating... Visualization Manager initialising... Registering graphics systems... You have successfully registered the following graphics systems. Current available graphics systems are: ASCIITree (ATree) DAWNFILE (DAWNFILE) G4HepRep (HepRepXML) G4HepRepFile (HepRepFile) RayTracer (RayTracer) VRML1FILE (VRML1FILE) VRML2FILE (VRML2FILE) OpenGLImmediateX (OGLIX) OpenGLStoredX (OGLSX) Registering model factories... You have successfully registered the following model factories. Registered model factories: generic drawByCharge drawByParticleID drawByOriginVolume drawByAttribute Registered filter factories: chargeFilter particleFilter originVolumeFilter attributeFilter /run/verbose 2 # /testem/det/setMat G4_WATER /testem/det/setRadius 3 cm /testem/det/setNbOfLayers 1 # /testem/phys/addPhysics standard PhysicsList::AddPhysicsList: # /run/initialize userDetector->Construct() start. --------------------------------------------------------- ---> The Absorber is a sphere of 3 cm radius of G4_WATER divided in 1 slices of 3 cm --------------------------------------------------------- Absorber is registered to the default region. physicsList->Construct() start. physicsList->Construct() start. physicsList->setCut() start. PhysicsList::SetCuts:CutLength : 1 mm # /testem/phys/setCuts 1 um # /gun/particle e- /gun/energy 4 MeV # /testem/histo/setFileType hbook /testem/histo/setFileName testem12 /testem/histo/setHisto 1 120 0. 2.4 cm ----> SetHisto 1: dE/dr (MeV/mm) along radius (cm); 120 bins from 0 cm to 2.4 cm /testem/histo/setHisto 2 180 0. 4.5 MeV ----> SetHisto 2: total Energy deposited in absorber (MeV); 180 bins from 0 MeV to 4.5 MeV /testem/histo/setHisto 3 100 0. 3. cm ----> SetHisto 3: true track length of the primary particle (cm); 100 bins from 0 cm to 3 cm /testem/histo/setHisto 4 100 0. 500. um ----> SetHisto 4: true step size of the primary particle (um); 100 bins from 0 um to 500 um /testem/histo/setHisto 5 200 -3. 3. cm ----> SetHisto 5: projected range of the primary particle (cm); 200 bins from -3 cm to 3 cm /testem/histo/setcsdaRange 2.037 cm /testem/histo/setHisto 8 120 0. 1.2 none ----> SetHisto 8: d(E/E0)/d(r/r0) along r/r0; 120 bins from 0 none to 1.2 none # /testem/event/printModulo 1000 # /run/beamOn 10000 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. msc: Model variant of multiple scattering for e- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 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= 0.2, integral: 1, fluct: 1 CSDA range table up to 1 GeV in 70 bins. 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= 0.2, integral: 1, fluct: 1 CSDA range table up to 1 GeV in 70 bins. 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. msc: Model variant of multiple scattering for proton Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 1 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= 0.2, integral: 1, fluct: 1 CSDA range table up to 1 GeV in 70 bins. msc: Model variant of multiple scattering for GenericIon LateralDisplacementFlag= 0 Skin= 0 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 1 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)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1 CSDA range table up to 1 GeV in 70 bins. 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= 0.2, integral: 1, fluct: 1 CSDA range table up to 1 GeV in 70 bins. msc: Model variant of multiple scattering for mu+ Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 1 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= 0.2, integral: 1, fluct: 1 CSDA range table up to 1 GeV in 70 bins. 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= 0.2, integral: 1, fluct: 1 CSDA range table up to 1 GeV in 70 bins. 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= 0.2, integral: 1, fluct: 1 CSDA range table up to 1 GeV in 70 bins. msc: Model variant of multiple scattering for pi- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 1 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= 0.2, integral: 1, fluct: 1 CSDA range table up to 1 GeV in 70 bins. Region -- appears in world volume Root logical volume(s) : Absorber Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] Materials : G4_WATER Production cuts : gamma 1 um e- 1 um e+ 1 um ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : G4_WATER Range cuts : gamma 1 um e- 1 um e+ 1 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV Region(s) which use this couple : DefaultRegionForTheWorld ==================================================================== Start closing geometry. G4GeometryManager::ReportVoxelStats -- Voxel Statistics Total memory consumed for geometry optimisation: 0 kByte Total CPU time elapsed for geometry optimisation: 0 seconds ### Run 0 start. --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 9876, 54321 ---------------------------------------- ---> stepMax from HistoManager = 2e+02 um Start Run processing. ---> Begin of Event: 0 ---> Begin of Event: 1000 ---> Begin of Event: 2000 ---> Begin of Event: 3000 ---> Begin of Event: 4000 ---> Begin of Event: 5000 ---> Begin of Event: 6000 ---> Begin of Event: 7000 ---> Begin of Event: 8000 ---> Begin of Event: 9000 Run terminated. Run Summary Number of events processed : 10000 User=65.4s Real=92.27s Sys=2.75s ======================== run summary ====================== The run consists of 10000 e- of 4.00 MeV through 3.00 cm of G4_WATER (density: 1.00 g/cm3 ) ============================================================ Total Energy deposited = 3.934 MeV +- 240.666 keV Track length of primary track = 2.016 cm +- 2.808 mm Range from EmCalculator = 2.021 cm (from full dE/dx) Projected range = 1.305 cm +- 4.176 mm Nb of steps of primary track = 238.40 +- 31.17 Step size= 84.562 um +- 4.576 um --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1695104246, 964521664 ---------------------------------------- Graphics systems deleted. Visualization Manager deleting... 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.