Environment variable "G4FORCE_RUN_MANAGER_TYPE" enabled with value == Serial. Forcing G4RunManager type... ############################################ !!! WARNING - FPE detection is activated !!! ############################################ ################################ !!! G4Backtrace is activated !!! ################################ ************************************************************** Geant4 version Name: geant4-11-04-patch-01 (13-March-2026) Copyright : Geant4 Collaboration References : NIM A 506 (2003), 250-303 : IEEE-TNS 53 (2006), 270-278 : NIM A 835 (2016), 186-225 WWW : http://geant4.org/ ************************************************************** /control/cout/ignoreThreadsExcept 0 /run/numberOfThreads 2 *** /run/numberOfThreads command is issued in sequential mode. Command is ignored. /run/verbose 1 # /testem/det/setAbsMat Silicon /testem/det/setAbsThick 20 um /testem/det/setAbsYZ 50 um # /testem/phys/addPhysics emstandard_opt3 PhysicsList::AddPhysicsList: /process/em/verbose 1 /process/em/workerVerbose 0 # /run/setCut 7 um # /run/initialize Material: G4_Galactic density: 0.000 mg/cm3 RadL: 204310098.490 pc Nucl.Int.Length: 3240.779 pc Imean: 21.800 eV temperature: 2.73 K pressure: 0.00 atm ---> Element: H (H) Z = 1.0 N = 1 A = 1.008 g/mole ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 % ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Silicon density: 2.330 g/cm3 RadL: 9.366 cm Nucl.Int.Length: 45.660 cm Imean: 173.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Si (Si) Z = 14.0 N = 28 A = 28.085 g/mole ---> Isotope: Si28 Z = 14 N = 28 A = 27.98 g/mole abundance: 92.230 % ---> Isotope: Si29 Z = 14 N = 29 A = 28.98 g/mole abundance: 4.683 % ---> Isotope: Si30 Z = 14 N = 30 A = 29.97 g/mole abundance: 3.087 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % The WORLD is made of 24 um of G4_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 # /testem/gun/setDefault /gun/particle pi+ /gun/energy 5 GeV # /run/printProgress 10000 /run/beamOn 100000 ======================================================================= ====== Electromagnetic Physics Parameters ======== ======================================================================= LPM effect enabled 1 Enable creation and use of sampling tables 0 Apply cuts on all EM processes 0 Use combined TransportationWithMsc Disabled Use general process 1 Enable linear polarisation for gamma 0 Enable photoeffect sampling below K-shell 1 Enable sampling of quantum entanglement 0 X-section factor for integral approach 0.8 Min kinetic energy for tables 10 eV Max kinetic energy for tables 100 TeV Number of bins per decade of a table 20 Verbose level 1 Verbose level for worker thread 0 Bremsstrahlung energy threshold above which primary e+- is added to the list of secondary 100 TeV Bremsstrahlung energy threshold above which primary muon/hadron is added to the list of secondary 100 TeV Positron annihilation at rest model AllisonPositronium Enable 3 gamma annihilation on fly 0 Lowest triplet kinetic energy 1 MeV Enable sampling of gamma linear polarisation 0 5D gamma conversion model type 0 5D gamma conversion model on isolated ion 0 Use RiGe 5D e+e- pair production model by muons 0 Livermore data directory epics_2017 ======================================================================= ====== Ionisation Parameters ======== ======================================================================= Step function for e+- (0.2, 0.1 mm) Step function for muons/hadrons (0.2, 0.05 mm) Step function for light ions (0.1, 0.02 mm) Step function for general ions (0.1, 0.001 mm) Lowest e+e- kinetic energy 100 eV Lowest muon/hadron kinetic energy 1 keV Use ICRU90 data 1 Fluctuations of dE/dx are enabled 1 Type of fluctuation model for leptons and hadrons Urban Use built-in Birks saturation 0 Build CSDA range enabled 0 Use cut as a final range enabled 0 Enable angular generator interface 1 Max kinetic energy for CSDA tables 1 GeV Max kinetic energy for NIEL computation 1 MeV Linear loss limit 0.01 Read data from file for e+e- pair production by mu 0 ======================================================================= ====== Multiple Scattering Parameters ======== ======================================================================= Type of msc step limit algorithm for e+- 3 Type of msc step limit algorithm for muons/hadrons 0 Msc lateral displacement for e+- enabled 1 Msc lateral displacement for muons and hadrons 1 Urban msc model lateral displacement alg96 1 Range factor for msc step limit for e+- 0.04 Range factor for msc step limit for muons/hadrons 0.2 Geometry factor for msc step limitation of e+- 2.5 Safety factor for msc step limit for e+- 0.6 Skin parameter for msc step limitation of e+- 1 Lambda limit for msc step limit for e+- 1 mm Use Mott correction for e- scattering 1 Factor used for dynamic computation of angular limit between single and multiple scattering 1 Fixed angular limit between single and multiple scattering 3.1416 rad Upper energy limit for e+- multiple scattering 100 MeV Type of electron single scattering model 0 Type of nuclear form-factor 1 Screening factor 1 ======================================================================= ====== Atomic Deexcitation Parameters ======== ======================================================================= Fluorescence enabled 1 Directory in G4LEDATA for fluorescence data files fluor Auger electron cascade enabled 0 PIXE atomic de-excitation enabled 0 De-excitation module ignores cuts 0 Type of PIXE cross section for hadrons Empirical Type of PIXE cross section for e+- Livermore ======================================================================= ### === Deexcitation model UAtomDeexcitation is activated for 1 region: DefaultRegionForTheWorld 1 0 0 ### === Ignore cuts flag: 0 phot: for gamma SubType=12 BuildTable=0 LambdaPrime table from 200 keV to 100 TeV in 174 bins ===== EM models for the G4Region DefaultRegionForTheWorld ====== LivermorePhElectric : Emin= 0 eV Emax= 100 TeV SauterGavrila Fluo compt: for gamma SubType=13 BuildTable=1 Lambda table from 10 eV to 1 MeV, 20 bins/decade, spline: 1 LambdaPrime table from 1 MeV to 100 TeV in 160 bins ===== EM models for the G4Region DefaultRegionForTheWorld ====== KleinNishina : Emin= 0 eV Emax= 100 TeV Fluo conv: for gamma SubType=14 BuildTable=1 Lambda table from 1.022 MeV to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== BetheHeitlerLPM : Emin= 0 eV Emax= 100 TeV ModifiedTsai Rayl: for gamma SubType=11 BuildTable=1 Lambda table from 10 eV to 150 keV, 20 bins/decade, spline: 0 LambdaPrime table from 150 keV to 100 TeV in 176 bins ===== EM models for the G4Region DefaultRegionForTheWorld ====== LivermoreRayleigh : Emin= 0 eV Emax= 100 TeV CullenGenerator msc: for e- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=DistanceToBoundary Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm eIoni: for e- XStype:3 SubType=2 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 100 TeV deltaVI eBrem: for e- XStype:4 SubType=3 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eBremSB : Emin= 0 eV Emax= 1 GeV AngularGen2BS eBremLPM : Emin= 1 GeV Emax= 100 TeV AngularGen2BS ePairProd: for e- XStype:1 SubType=4 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 0 Sampling table 25x1001 from 0.1 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== ePairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi msc: for e+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=DistanceToBoundary Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm eIoni: for e+ XStype:3 SubType=2 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 100 TeV deltaVI eBrem: for e+ XStype:4 SubType=3 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eBremSB : Emin= 0 eV Emax= 1 GeV AngularGen2BS eBremLPM : Emin= 1 GeV Emax= 100 TeV AngularGen2BS ePairProd: for e+ XStype:1 SubType=4 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 0 Sampling table 25x1001 from 0.1 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== ePairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi annihil: for e+ XStype:2 SubType=5 AtRestModel:Allison BuildTable=0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eplus2gg : Emin= 0 eV Emax= 100 TeV msc: for proton SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm hIoni: for proton XStype:3 SubType=2 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 2 MeV deltaVI BetheBloch : Emin= 2 MeV Emax= 100 TeV deltaVI hBrems: for proton XStype:1 SubType=3 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi hPairProd: for proton XStype:1 SubType=4 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 Sampling table 17x1001, from 7.50618 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi nuclearStopping: for proton SubType=8 BuildTable=0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV msc: for GenericIon SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm ionIoni: for GenericIon XStype:3 SubType=2 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.1, 0.001 mm), integ: 3, fluct: 1, linLossLim= 0.02 ===== EM models for the G4Region DefaultRegionForTheWorld ====== LindhardSorensen : Emin= 0 eV Emax= 100 TeV deltaVI nuclearStopping: for GenericIon SubType=8 BuildTable=0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV msc: for alpha SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm ionIoni: for alpha XStype:3 SubType=2 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.1, 0.02 mm), integ: 3, fluct: 1, linLossLim= 0.02 ===== EM models for the G4Region DefaultRegionForTheWorld ====== BraggIon : Emin= 0 eV Emax=7.9452 MeV deltaVI BetheBloch : Emin=7.9452 MeV Emax= 100 TeV deltaVI nuclearStopping: for alpha SubType=8 BuildTable=0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV msc: for anti_proton SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm hIoni: for anti_proton XStype:3 SubType=2 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 2 MeV deltaVI BetheBloch : Emin= 2 MeV Emax= 100 TeV deltaVI hBrems: for anti_proton XStype:1 SubType=3 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi hPairProd: for anti_proton XStype:1 SubType=4 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 Sampling table 17x1001, from 7.50618 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi msc: for kaon+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm hIoni: for kaon+ XStype:3 SubType=2 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax=1.05231 MeV deltaVI BetheBloch : Emin=1.05231 MeV Emax= 100 TeV deltaVI hBrems: for kaon+ XStype:1 SubType=3 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi hPairProd: for kaon+ XStype:1 SubType=4 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 Sampling table 18x1001, from 3.94942 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi msc: for kaon- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm hIoni: for kaon- XStype:3 SubType=2 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax=1.05231 MeV deltaVI BetheBloch : Emin=1.05231 MeV Emax= 100 TeV deltaVI hBrems: for kaon- XStype:1 SubType=3 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi hPairProd: for kaon- XStype:1 SubType=4 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 Sampling table 18x1001, from 3.94942 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi msc: for mu+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm muIoni: for mu+ XStype:3 SubType=2 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 200 keV deltaVI MuBetheBloch : Emin= 200 keV Emax= 100 TeV deltaVI muBrems: for mu+ XStype:1 SubType=3 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi muPairProd: for mu+ XStype:1 SubType=4 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 Sampling table 21x1001, from 0.85 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi msc: for mu- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm muIoni: for mu- XStype:3 SubType=2 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 200 keV deltaVI MuBetheBloch : Emin= 200 keV Emax= 100 TeV deltaVI muBrems: for mu- XStype:1 SubType=3 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi muPairProd: for mu- XStype:1 SubType=4 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 Sampling table 21x1001, from 0.85 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi msc: for pi+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm hIoni: for pi+ XStype:3 SubType=2 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax=297.505 keV deltaVI BetheBloch : Emin=297.505 keV Emax= 100 TeV deltaVI hBrems: for pi+ XStype:1 SubType=3 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi hPairProd: for pi+ XStype:1 SubType=4 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 Sampling table 20x1001, from 1.11656 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi msc: for pi- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm hIoni: for pi- XStype:3 SubType=2 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax=297.505 keV deltaVI BetheBloch : Emin=297.505 keV Emax= 100 TeV deltaVI hBrems: for pi- XStype:1 SubType=3 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi hPairProd: for pi- XStype:1 SubType=4 dE/dx and range tables from 10 eV to 100 TeV in 260 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 Sampling table 20x1001, from 1.11656 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi ========= Table of registered couples ============================ Index : 0 used in the geometry : Yes Material : G4_Galactic Range cuts : gamma 7 um e- 7 um e+ 7 um proton 7 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 700 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes Material : Silicon Range cuts : gamma 7 um e- 7 um e+ 7 um proton 7 um Energy thresholds : gamma 990 eV e- 21.0638 keV e+ 20.5957 keV proton 700 eV Region(s) which use this couple : DefaultRegionForTheWorld ================================================================== ### Run 0 starts. ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={906770732717044781, 629165745432651234, 1235682547346241386, 68420008233404568, 2262190991329242458, 2266470399991071809, 1976726662926872232, 245458862506414172, 1955974201201518530, 2155248512522080758, 604170912935414061, 1116171330120743511, 1861018313684488333, 1296715403254578286, 1549011045957234151, 370819759640195970, 2230139271784837643} counter= 17sumtot= 1977567618660788324 --------------------------------------- --> Event 0 starts. --> Event 10000 starts. --> Event 20000 starts. --> Event 30000 starts. --> Event 40000 starts. --> Event 50000 starts. --> Event 60000 starts. --> Event 70000 starts. --> Event 80000 starts. --> Event 90000 starts. Run terminated. Run Summary Number of events processed : 100000 User=0.430000s Real=0.434872s Sys=0.000000s ======================== 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.353 keV +- 16.34 eV -----> Mean dE/dx = 3.177 MeV/cm (1.363 MeV*cm2/g) From formulas : restricted dEdx = 2.96 MeV/cm (1.27 MeV*cm2/g) full dEdx = 4.726 MeV/cm (2.028 MeV*cm2/g) Leakage : primary = 5 GeV +- 725.3 eV secondaries = 2.672 keV +- 724.6 eV Energy balance : edep + eleak = 5 GeV Total track length (charged) in absorber per event = 20.26 um +- 7.771 nm Total track length (neutral) in absorber per event = 7.453e+04 fm +- 7.453e+04 fm Number of steps (charged) in absorber per event = 1.179 +- 0.005528 Number of steps (neutral) in absorber per event = 1e-05 +- 1e-05 Number of secondaries per event : Gammas = 1e-05 electrons = 0.017 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.02504 mrad (central part only) computed theta0 (Highland formula) = 0.02627 mrad central part defined as +- 0.07881 mrad Tail ratio = 1.951 % Gamma process counts: Photoeffect 0 Compton 0 Conversion 0 Rayleigh 0 ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={1378456096184970916, 645125664873386088, 289313246081823319, 692181631145386720, 1554795592127445348, 429523111405020466, 1767267687826005316, 887948550155982521, 745800083929026403, 309730682002155907, 2175523263042675480, 1439287206468397773, 1375105445691351298, 782374144723864041, 1715829858399311388, 655882193793759711, 1581270085403564578} counter= 8sumtot= 2284513478758269616 --------------------------------------- # /testem/det/setAbsMat Gold /testem/det/setAbsThick 9.658 um /testem/det/setAbsYZ 1 mm /testem/gun/setDefault # Field Should be issued after /run/initialize #/testem/det/setField 1 tesla /gun/particle e- /gun/energy 15.7 MeV /run/beamOn 50000 ### === Deexcitation model UAtomDeexcitation is activated for 1 region: DefaultRegionForTheWorld 1 0 0 ### === Ignore cuts flag: 0 ========= Table of registered couples ============================ Index : 0 used in the geometry : Yes Material : G4_Galactic Range cuts : gamma 7 um e- 7 um e+ 7 um proton 7 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 700 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : No Material : Silicon Range cuts : gamma 7 um e- 7 um e+ 7 um proton 7 um Energy thresholds : gamma 990 eV e- 21.0638 keV e+ 20.5957 keV proton 700 eV Index : 2 used in the geometry : Yes Material : Gold Range cuts : gamma 7 um e- 7 um e+ 7 um proton 7 um Energy thresholds : gamma 6.52544 keV e- 65.1638 keV e+ 63.9139 keV proton 700 eV Region(s) which use this couple : DefaultRegionForTheWorld ================================================================== ### Run 1 starts. ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={1378456096184970916, 645125664873386088, 289313246081823319, 692181631145386720, 1554795592127445348, 429523111405020466, 1767267687826005316, 887948550155982521, 745800083929026403, 309730682002155907, 2175523263042675480, 1439287206468397773, 1375105445691351298, 782374144723864041, 1715829858399311388, 655882193793759711, 1581270085403564578} counter= 8sumtot= 2284513478758269616 --------------------------------------- --> Event 0 starts. --> Event 10000 starts. --> Event 20000 starts. --> Event 30000 starts. --> Event 40000 starts. Run terminated. Run Summary Number of events processed : 50000 User=0.300000s Real=0.298352s Sys=0.000000s ======================== 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.38 keV +- 64.44 eV -----> Mean dE/dx = 19.03 MeV/cm (0.9848 MeV*cm2/g) From formulas : restricted dEdx = 17.91 MeV/cm (0.927 MeV*cm2/g) full dEdx = 61.41 MeV/cm (3.179 MeV*cm2/g) Leakage : primary = 15.64 MeV +- 2.11 keV secondaries = 37.48 keV +- 2.102 keV Energy balance : edep + eleak = 15.7 MeV Total track length (charged) in absorber per event = 9.911 um +- 11.24 nm Total track length (neutral) in absorber per event = 141.4 nm +- 13.24 nm Number of steps (charged) in absorber per event = 2.602 +- 0.009024 Number of steps (neutral) in absorber per event = 0.02698 +- 0.0007528 Number of secondaries per event : Gammas = 0.0269 electrons = 0.02112 positrons = 0 Number of events with the primary particle transmitted = 99.98 % Number of events with at least 1 particle transmitted (same charge as primary) = 99.98 % Number of events with the primary particle reflected = 0.02 % Number of events with at least 1 particle reflected (same charge as primary) = 0.31 % MultipleScattering: rms proj angle of transmit primary particle = 34.71 mrad (central part only) computed theta0 (Highland formula) = 35.1 mrad central part defined as +- 105.3 mrad Tail ratio = 2.948 % Gamma process counts: Photoeffect 198 Compton 0 Conversion 0 Rayleigh 4 ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={187465526861526162, 1919470187343084000, 1200600146555496649, 2254698630955541145, 1010695742965213971, 386425232931658995, 1773838893809632618, 1545813557997812675, 1661416122218187049, 1838452107617613130, 852398828744869876, 2030305144316100, 2300250181027420193, 126900021668951897, 823556986218677264, 1235812847093715199, 1130805818374825746} counter= 11sumtot= 1803887063818991061 --------------------------------------- # /run/setCut 10 um # /testem/det/setAbsMat G4_POLYETHYLENE /testem/det/setAbsThick 4.74 mm /testem/det/setAbsYZ 50 mm /testem/gun/setDefault /gun/particle mu+ /gun/energy 96.2 MeV /run/beamOn 50000 ### === Deexcitation model UAtomDeexcitation is activated for 1 region: DefaultRegionForTheWorld 1 0 0 ### === Ignore cuts flag: 0 ========= Table of registered couples ============================ Index : 0 used in the geometry : Yes 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 : No Material : Silicon Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um Energy thresholds : gamma 990 eV e- 31.6994 keV e+ 31.3501 keV proton 1 keV Index : 2 used in the geometry : No Material : Gold Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um Energy thresholds : gamma 8.32988 keV e- 80.7452 keV e+ 79.0938 keV proton 1 keV Index : 3 used in the geometry : Yes Material : G4_POLYETHYLENE Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um Energy thresholds : gamma 990 eV e- 13.7224 keV e+ 13.4607 keV proton 1 keV Region(s) which use this couple : DefaultRegionForTheWorld ================================================================== ### Run 2 starts. ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={187465526861526162, 1919470187343084000, 1200600146555496649, 2254698630955541145, 1010695742965213971, 386425232931658995, 1773838893809632618, 1545813557997812675, 1661416122218187049, 1838452107617613130, 852398828744869876, 2030305144316100, 2300250181027420193, 126900021668951897, 823556986218677264, 1235812847093715199, 1130805818374825746} counter= 11sumtot= 1803887063818991061 --------------------------------------- --> Event 0 starts. --> Event 10000 starts. --> Event 20000 starts. --> Event 30000 starts. --> Event 40000 starts. Run terminated. Run Summary Number of events processed : 50000 User=0.910000s Real=0.936989s Sys=0.000000s ======================== 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.048 MeV +- 993.7 eV -----> Mean dE/dx = 2.21 MeV/cm (2.351 MeV*cm2/g) From formulas : restricted dEdx = 1.749 MeV/cm (1.861 MeV*cm2/g) full dEdx = 2.265 MeV/cm (2.409 MeV*cm2/g) Leakage : primary = 95.13 MeV +- 1.332 keV secondaries = 26.01 keV +- 701.3 eV Energy balance : edep + eleak = 96.2 MeV Total track length (charged) in absorber per event = 5.17 mm +- 3.775 um Total track length (neutral) in absorber per event = 22.35 um +- 2.273 um Number of steps (charged) in absorber per event = 12.17 +- 0.03721 Number of steps (neutral) in absorber per event = 0.00988 +- 0.0004846 Number of secondaries per event : Gammas = 0.00914 electrons = 4.181 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 = 7.373 mrad (central part only) computed theta0 (Highland formula) = 7.636 mrad central part defined as +- 22.91 mrad Tail ratio = 1.214 % Gamma process counts: Photoeffect 305 Compton 23 Conversion 0 Rayleigh 14 ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={2083361880407413973, 558014575596605443, 155743980189950201, 2182957000861726697, 1838655335096136537, 1946843552728816280, 1681406852924851701, 474130547566224162, 1018240738768540513, 1872432096900372193, 472993835442224762, 10222201860268061, 1779133025925577121, 2263978871037324573, 1473228915816656057, 35068858945786700, 1174357358231345229} counter= 2sumtot= 268182545376574644 --------------------------------------- # /run/setCut 1 mm # /testem/stack/killSecondaries # /testem/det/setAbsMat Aluminium /testem/det/setAbsThick 8.004 mm /testem/gun/setDefault /gun/particle proton /gun/energy 174.664 GeV /run/beamOn 50000 ### === Deexcitation model UAtomDeexcitation is activated for 1 region: DefaultRegionForTheWorld 1 0 0 ### === Ignore cuts flag: 0 ========= Table of registered couples ============================ Index : 0 used in the geometry : Yes Material : G4_Galactic Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : No Material : Silicon Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 6.93672 keV e- 546.311 keV e+ 525.593 keV proton 100 keV Index : 2 used in the geometry : No Material : Gold Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 112.405 keV e- 2.28837 MeV e+ 2.12466 MeV proton 100 keV Index : 3 used in the geometry : No Material : G4_POLYETHYLENE Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 2.25978 keV e- 351.127 keV e+ 341.776 keV proton 100 keV Index : 4 used in the geometry : Yes Material : Aluminium Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 6.91691 keV e- 597.603 keV e+ 574.29 keV proton 100 keV Region(s) which use this couple : DefaultRegionForTheWorld ================================================================== ### Run 3 starts. ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={2083361880407413973, 558014575596605443, 155743980189950201, 2182957000861726697, 1838655335096136537, 1946843552728816280, 1681406852924851701, 474130547566224162, 1018240738768540513, 1872432096900372193, 472993835442224762, 10222201860268061, 1779133025925577121, 2263978871037324573, 1473228915816656057, 35068858945786700, 1174357358231345229} counter= 2sumtot= 268182545376574644 --------------------------------------- --> Event 0 starts. --> Event 10000 starts. --> Event 20000 starts. --> Event 30000 starts. --> Event 40000 starts. Run terminated. Run Summary Number of events processed : 50000 User=0.220000s Real=0.217912s Sys=0.000000s ======================== 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 = 5.362 MeV +- 381.9 keV -----> Mean dE/dx = 6.699 MeV/cm (2.481 MeV*cm2/g) From formulas : restricted dEdx = 4.024 MeV/cm (1.49 MeV*cm2/g) full dEdx = 6.049 MeV/cm (2.24 MeV*cm2/g) Leakage : primary = 174.7 GeV +- 381.9 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.002316 Number of steps (neutral) in absorber per event = 0 +- 0 Number of secondaries per event : Gammas = 0 electrons = 0.268 positrons = 0.00088 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.02062 mrad (central part only) computed theta0 (Highland formula) = 0.02111 mrad central part defined as +- 0.06333 mrad Tail ratio = 1.327 % Gamma process counts: Photoeffect 0 Compton 0 Conversion 0 Rayleigh 0 ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={864120598200464168, 189420282748189000, 2125101903705486721, 1863466081725669111, 1383960382039097673, 550892631529147885, 2244069534527202395, 933979936320874605, 791647544461653037, 286412522794523854, 1616800613386878212, 478416970401061220, 1548587089976077386, 917812281992815632, 2284924600166281660, 94168853074452663, 1674934093723585259} counter= 15sumtot= 1401971847063908873 --------------------------------------- # /testem/det/setAbsMat Iron /testem/det/setAbsThick 1 m /testem/gun/setDefault /gun/particle mu+ /gun/energy 100 GeV /run/beamOn 50000 ### === Deexcitation model UAtomDeexcitation is activated for 1 region: DefaultRegionForTheWorld 1 0 0 ### === Ignore cuts flag: 0 ========= Table of registered couples ============================ Index : 0 used in the geometry : Yes Material : G4_Galactic Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : No Material : Silicon Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 6.93672 keV e- 546.311 keV e+ 525.593 keV proton 100 keV Index : 2 used in the geometry : No Material : Gold Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 112.405 keV e- 2.28837 MeV e+ 2.12466 MeV proton 100 keV Index : 3 used in the geometry : No Material : G4_POLYETHYLENE Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 2.25978 keV e- 351.127 keV e+ 341.776 keV proton 100 keV Index : 4 used in the geometry : No Material : Aluminium Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 6.91691 keV e- 597.603 keV e+ 574.29 keV proton 100 keV Index : 5 used in the geometry : Yes Material : Iron Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 20.7519 keV e- 1.28828 MeV e+ 1.21724 MeV proton 100 keV Region(s) which use this couple : DefaultRegionForTheWorld ================================================================== ### Run 4 starts. ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={864120598200464168, 189420282748189000, 2125101903705486721, 1863466081725669111, 1383960382039097673, 550892631529147885, 2244069534527202395, 933979936320874605, 791647544461653037, 286412522794523854, 1616800613386878212, 478416970401061220, 1548587089976077386, 917812281992815632, 2284924600166281660, 94168853074452663, 1674934093723585259} counter= 15sumtot= 1401971847063908873 --------------------------------------- --> Event 0 starts. --> Event 10000 starts. --> Event 20000 starts. --> Event 30000 starts. --> Event 40000 starts. Run terminated. Run Summary Number of events processed : 50000 User=3.840000s Real=3.908897s Sys=0.000000s ======================== 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.163 GeV +- 15.06 MeV -----> Mean dE/dx = 21.63 MeV/cm (2.748 MeV*cm2/g) From formulas : restricted dEdx = 11.21 MeV/cm (1.425 MeV*cm2/g) full dEdx = 21.61 MeV/cm (2.746 MeV*cm2/g) Leakage : primary = 97.84 GeV +- 15.06 MeV secondaries = 0 eV +- 0 eV Energy balance : edep + eleak = 100 GeV Total track length (charged) in absorber per event = 1 m +- 7.06 nm Total track length (neutral) in absorber per event = 0 fm +- 0 fm Number of steps (charged) in absorber per event = 47.58 +- 0.03044 Number of steps (neutral) in absorber per event = 0 +- 0 Number of secondaries per event : Gammas = 0.05044 electrons = 46.52 positrons = 1.254 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.178 mrad (central part only) computed theta0 (Highland formula) = 1.182 mrad central part defined as +- 3.546 mrad Tail ratio = 0.579 % Gamma process counts: Photoeffect 0 Compton 0 Conversion 0 Rayleigh 0 ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={1046585659322464193, 1700724626221482459, 1710146587566228662, 2274062821915192324, 1093598942862882934, 859975426035031537, 1845450543346184012, 1347610440888355163, 1747669745252223770, 1589127766525861706, 1326010817860615296, 367647672722605028, 58019342554185087, 2022010207749596709, 2304966030133393681, 135872583825542285, 522422170598193944} counter= 6sumtot= 1199314302456793231 --------------------------------------- # ================== Deleting memory pools =================== Number of memory pools allocated: 9 of which, static: 0 Dynamic pools deleted: 9 / Total memory freed: 0.056 MB ============================================================