source: trunk/examples/extended/electromagnetic/TestEm5/TestEm5.out @ 1333

*************************************************************
 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
*************************************************************

/run/verbose 2
#
/testem/det/setAbsMat Silicon
/testem/det/setAbsThick 20 um
/testem/det/setAbsYZ 50 um
#
/testem/phys/addPhysics standard
PhysicsList::AddPhysicsList: <standard>
#
/testem/phys/setGCut 2.13 mm
/testem/phys/setECut 7.0 um
/testem/phys/setPCut 7.0 um
#
/run/initialize
userDetector->Construct() start.

 Material: Galactic     density:  0.000 mg/cm3  RadL: 204727576.737 pc   Imean:  21.800 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   Imean: 173.600 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->setCut() start.
PhysicsList::SetCuts:CutLength : 1 mm 
#
/testem/gun/setDefault
/gun/particle pi+
/gun/energy 5 GeV
#
/testem/histo/setFileName testem5
/testem/histo/setFileType hbook
/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:  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 TeV in 1200 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 TeV in 1200 bins.

msc:  Model variant of multiple scattering for e-
      Lambda tables from 100 eV  to 100 TeV in 1200 bins.
      LateralDisplacementFlag=  1   Skin= 4
      Boundary/stepping algorithm is active with RangeFactor= 0.02  Step limit type 2

eIoni:   tables are built for  e-
      dE/dx and range tables from 100 eV  to 100 TeV in 1200 bins.
      Lambda tables from threshold to 100 TeV in 1200 bins.
      Delta cross sections and sampling from MollerBhabha model
      Good description from 1 KeV to 100 GeV.
      Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1

eBrem:   tables are built for  e-
      dE/dx and range tables from 100 eV  to 100 TeV in 1200 bins.
      Lambda tables from threshold to 100 TeV in 1200 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 1200 bins.
      Lambda tables from threshold to 100 TeV in 1200 bins.
      Delta cross sections and sampling from MollerBhabha model
      Good description from 1 KeV to 100 GeV.
      Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1

eBrem:   tables are built for  e+
      dE/dx and range tables from 100 eV  to 100 TeV in 1200 bins.
      Lambda tables from threshold to 100 TeV in 1200 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 1200 bins.
      LateralDisplacementFlag=  1   Skin= 4
      Boundary/stepping algorithm is active with RangeFactor= 0.02  Step limit type 2

hIoni:   tables are built for  proton
      dE/dx and range tables from 100 eV  to 100 TeV in 1200 bins.
      Lambda tables from threshold to 100 TeV in 1200 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)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1

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 1200 bins.
      Lambda tables from threshold to 100 TeV in 1200 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.2, integral: 1, fluct: 1

hIoni:   tables are built for  anti_proton
      dE/dx and range tables from 100 eV  to 100 TeV in 1200 bins.
      Lambda tables from threshold to 100 TeV in 1200 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)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1

msc:  Model variant of multiple scattering for mu+
      Lambda tables from 100 eV  to 100 TeV in 1200 bins.
      LateralDisplacementFlag=  1   Skin= 4
      Boundary/stepping algorithm is active with RangeFactor= 0.02  Step limit type 2

muIoni:   tables are built for  mu+
      dE/dx and range tables from 100 eV  to 100 TeV in 1200 bins.
      Lambda tables from threshold to 100 TeV in 1200 bins.
      Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, 
      radiative corrections for E > 1 GeV
      Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1

muBrems:   tables are built for  mu+
      dE/dx and range tables from 100 eV  to 100 TeV in 1200 bins.
      Lambda tables from threshold to 100 TeV in 1200 bins.
      Parametrised model 

muPairProd:   tables are built for  mu+
      dE/dx and range tables from 100 eV  to 100 TeV in 1200 bins.
      Lambda tables from threshold to 100 TeV in 1200 bins.
      Parametrised model 

muIoni:   tables are built for  mu-
      dE/dx and range tables from 100 eV  to 100 TeV in 1200 bins.
      Lambda tables from threshold to 100 TeV in 1200 bins.
      Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, 
      radiative corrections for E > 1 GeV
      Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1

muBrems:   tables are built for  mu-
      dE/dx and range tables from 100 eV  to 100 TeV in 1200 bins.
      Lambda tables from threshold to 100 TeV in 1200 bins.
      Parametrised model 

muPairProd:   tables are built for  mu-
      dE/dx and range tables from 100 eV  to 100 TeV in 1200 bins.
      Lambda tables from threshold to 100 TeV in 1200 bins.
      Parametrised model 

hIoni:   tables are built for  pi+
      dE/dx and range tables from 100 eV  to 100 TeV in 1200 bins.
      Lambda tables from threshold to 100 TeV in 1200 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)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1

msc:  Model variant of multiple scattering for pi-
      Lambda tables from 100 eV  to 100 TeV in 1200 bins.
      LateralDisplacementFlag=  1   Skin= 4
      Boundary/stepping algorithm is active with RangeFactor= 0.02  Step limit type 2

hIoni:   tables are built for  pi-
      dE/dx and range tables from 100 eV  to 100 TeV in 1200 bins.
      Lambda tables from threshold to 100 TeV in 1200 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)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1

Region <DefaultRegionForTheWorld> -- appears in <World> world volume
 Root logical volume(s) : World 
 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
 Materials : Galactic Silicon 
 Production cuts :  gamma 2.13 mm     e- 7 um     e+ 7 um 

========= Table of registered couples ==============================

Index : 0     used in the geometry : Yes     recalculation needed : No 
 Material : Galactic
 Range cuts        :  gamma 2.13 mm     e- 7 um     e+ 7 um 
 Energy thresholds :  gamma 990 eV     e- 990 eV     e+ 990 eV 
 Region(s) which use this couple : 
    DefaultRegionForTheWorld

Index : 1     used in the geometry : Yes     recalculation needed : No 
 Material : Silicon
 Range cuts        :  gamma 2.13 mm     e- 7 um     e+ 7 um 
 Energy thresholds :  gamma 9.93829 keV    e- 20.7856 keV    e+ 20.268 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=18.6s Real=260.66s Sys=23.45s

 ======================== 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.331 keV +- 16.59 eV 

 -----> Mean dE/dx = 3.165 MeV/cm       (1.359 MeV*cm2/g)

 From formulas :
   restricted dEdx = 2.936 MeV/cm       (1.26 MeV*cm2/g)
   full dEdx       = 4.697 MeV/cm       (2.016 MeV*cm2/g)

 Leakage :  primary = 5 GeV +- 414.7 eV    secondaries = 1.972 keV +- 413.7 eV 
 Energy balance :  edep + eleak = 5 GeV

 Total track length (charged) in absorber per event = 20.28 um  +- 8.08 nm 
 Total track length (neutral) in absorber per event = 6.125 Ang +- 4.382 Ang

 Number of steps (charged) in absorber per event = 10.15 +- 8.989e-06
 Number of steps (neutral) in absorber per event = 2e-05 +- 8.989e-06

 Number of secondaries per event : Gammas = 2e-05;   electrons = 0.01764;   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.02411 mrad (central part only)
  computed theta0 (Highland formula)          = 0.02627 mrad
  central part defined as +- 0.0788 mrad;   Tail ratio = 1.011 %

--------- Ranecu engine status ---------
 Initial seed (index) = 0
 Current couple of seeds = 921430354, 552442089
----------------------------------------
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.