############################################ !!! 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/ ************************************************************** /run/verbose 1 #/tracking/verbose 1 # # select EM PhysicsList # /testem/phys/addPhysics emstandard_opt1 PhysicsList::Set emstandard_opt1 EM physics # /process/em/AddEmRegion VertexDetector G4EmStandard_opt4 /process/em/AddEmRegion MuonDetector G4EmStandard # /run/initialize Ecal is 36 cm of G4_CESIUM_IODIDE Absorber is 2 mm of G4_Al Vertex is 3 cm of 3 layers of Si of 100 um npads= 601 ### New geometry is constructed # /run/setCut 0.3 mm /run/setCutForRegion VertexDetector 0.01 mm /run/setCutForRegion MuonDetector 10 mm /testem/stepMax 100 mm # /testem/det/acceptance1 0.8271 0.02974 2.5 /testem/det/acceptance9 0.9570 0.01286 2.5 /testem/det/acceptance25 0.9793 0.00884 2.5 # /gun/particle e- /gun/energy 1 GeV /run/beamOn 100 ======================================================================= ====== Electromagnetic Physics Parameters ======== ======================================================================= LPM effect enabled 1 Enable creation and use of sampling tables 0 Apply cuts on all EM processes 1 Use combined TransportationWithMsc Enabled 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 100 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 SimplePositronium 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.8, 1 mm) Step function for muons/hadrons (0.2, 0.1 mm) Step function for light ions (0.2, 0.1 mm) Step function for general ions (0.2, 0.1 mm) Lowest e+e- kinetic energy 1 keV Lowest muon/hadron kinetic energy 1 keV Use ICRU90 data 0 Fluctuations of dE/dx are enabled 1 Type of fluctuation model for leptons and hadrons Universal Use built-in Birks saturation 0 Build CSDA range enabled 0 Use cut as a final range enabled 0 Enable angular generator interface 0 Max kinetic energy for CSDA tables 1 GeV Max kinetic energy for NIEL computation 0 eV 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+- 0 Type of msc step limit algorithm for muons/hadrons 0 Msc lateral displacement for e+- enabled 1 Msc lateral displacement for muons and hadrons 0 Urban msc model lateral displacement alg96 1 Range factor for msc step limit for e+- 0.2 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 3 regions: DefaultRegionForTheWorld 0 0 0 VertexDetector 1 0 0 MuonDetector 0 0 0 ### === Ignore cuts flag: 0 phot: for gamma applyCuts:1 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 applyCuts:1 SubType=13 BuildTable=1 Lambda table from 100 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 ====== Klein-Nishina : Emin= 0 eV Emax= 100 TeV ===== EM models for the G4Region VertexDetector ====== LowEPComptonModel : Emin= 0 eV Emax= 20 MeV Fluo KleinNishina : Emin= 20 MeV Emax= 100 TeV Fluo conv: for gamma applyCuts:1 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 ===== EM models for the G4Region VertexDetector ====== BetheHeitler5D : Emin= 0 eV Emax= 100 TeV ModifiedTsai Rayl: for gamma applyCuts:1 SubType=11 BuildTable=1 Lambda table from 100 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 ====== DummyModel : Emin= 0 eV Emax= 100 TeV StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm ===== EM models for the G4Region VertexDetector ====== LivermoreRayleigh : Emin= 0 eV Emax= 100 TeV CullenGenerator TransportationWithMsc: for e- ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 MeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm WentzelVIUni : Emin= 100 MeV Emax= 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm ===== EM models for the G4Region VertexDetector ====== GoudsmitSaunderson : Emin= 0 eV Emax= 100 MeV StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm WentzelVIUni : Emin= 100 MeV Emax= 100 TeV StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm ===== EM models for the G4Region MuonDetector ====== UrbanMsc : Emin= 0 eV Emax= 100 MeV StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm WentzelVIUni : Emin= 100 MeV Emax= 100 TeV StepLim=UseSafety 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 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.8, 1 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 100 TeV eBrem: for e- XStype:4 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 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 ModifiedTsai eBremLPM : Emin= 1 GeV Emax= 100 TeV ModifiedTsai CoulombScat: for e- XStype:1 applyCuts:1 SubType=1 BuildTable=1 Lambda table from 100 MeV to 100 TeV, 20 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 100 MeV Emax= 100 TeV ===== EM models for the G4Region VertexDetector ====== eCoulombScattering : Emin= 100 MeV Emax= 100 TeV ===== EM models for the G4Region MuonDetector ====== eCoulombScattering : Emin= 100 MeV Emax= 100 TeV TransportationWithMsc: for e+ ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 MeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm WentzelVIUni : Emin= 100 MeV Emax= 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm ===== EM models for the G4Region VertexDetector ====== GoudsmitSaunderson : Emin= 0 eV Emax= 100 MeV StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm WentzelVIUni : Emin= 100 MeV Emax= 100 TeV StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm ===== EM models for the G4Region MuonDetector ====== UrbanMsc : Emin= 0 eV Emax= 100 MeV StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm WentzelVIUni : Emin= 100 MeV Emax= 100 TeV StepLim=UseSafety 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 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.8, 1 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 100 TeV eBrem: for e+ XStype:4 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 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 ModifiedTsai eBremLPM : Emin= 1 GeV Emax= 100 TeV ModifiedTsai annihil: for e+ XStype:2 applyCuts:1 SubType=5 AtRestModel:Simple BuildTable=0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eplus2gg : Emin= 0 eV Emax= 100 TeV CoulombScat: for e+ XStype:1 applyCuts:1 SubType=1 BuildTable=1 Lambda table from 100 MeV to 100 TeV, 20 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 100 MeV Emax= 100 TeV ===== EM models for the G4Region VertexDetector ====== eCoulombScattering : Emin= 100 MeV Emax= 100 TeV ===== EM models for the G4Region MuonDetector ====== eCoulombScattering : Emin= 100 MeV Emax= 100 TeV msc: for proton SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for proton XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 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 ====== Bragg : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 100 TeV hBrems: for proton XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 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 100 eV to 100 TeV in 240 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 CoulombScat: for proton XStype:1 applyCuts:1 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 20 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV 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:0 Skin=1 Llim=1 mm ionIoni: for GenericIon XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 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.02 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 100 TeV 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:0 Skin=1 Llim=1 mm ionIoni: for alpha XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 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.02 ===== EM models for the G4Region DefaultRegionForTheWorld ====== BraggIon : Emin= 0 eV Emax=7.9452 MeV BetheBloch : Emin=7.9452 MeV Emax= 100 TeV msc: for anti_proton SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for anti_proton XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 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 ====== ICRU73QO : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 100 TeV hBrems: for anti_proton XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 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 100 eV to 100 TeV in 240 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 CoulombScat: for anti_proton XStype:1 applyCuts:1 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 20 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for kaon+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for kaon+ XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 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 ====== Bragg : Emin= 0 eV Emax=1.05231 MeV BetheBloch : Emin=1.05231 MeV Emax= 100 TeV hBrems: for kaon+ XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 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 100 eV to 100 TeV in 240 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 CoulombScat: for kaon+ XStype:1 applyCuts:1 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 20 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for kaon- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for kaon- XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 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 ====== ICRU73QO : Emin= 0 eV Emax=1.05231 MeV BetheBloch : Emin=1.05231 MeV Emax= 100 TeV hBrems: for kaon- XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 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 100 eV to 100 TeV in 240 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 CoulombScat: for kaon- XStype:1 applyCuts:1 SubType=1 BuildTable=1 Used Lambda table of kaon+ ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for mu+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm muIoni: for mu+ XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 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 ====== Bragg : Emin= 0 eV Emax= 200 keV MuBetheBloch : Emin= 200 keV Emax= 100 TeV muBrems: for mu+ XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 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 100 eV to 100 TeV in 240 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 CoulombScat: for mu+ XStype:1 applyCuts:1 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 20 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for mu- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm muIoni: for mu- XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 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 ====== ICRU73QO : Emin= 0 eV Emax= 200 keV MuBetheBloch : Emin= 200 keV Emax= 100 TeV muBrems: for mu- XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 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 100 eV to 100 TeV in 240 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 CoulombScat: for mu- XStype:1 applyCuts:1 SubType=1 BuildTable=1 Used Lambda table of mu+ ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for pi+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for pi+ XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 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 ====== Bragg : Emin= 0 eV Emax=297.505 keV BetheBloch : Emin=297.505 keV Emax= 100 TeV hBrems: for pi+ XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 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 100 eV to 100 TeV in 240 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 CoulombScat: for pi+ XStype:1 applyCuts:1 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 20 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for pi- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for pi- XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 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 ====== ICRU73QO : Emin= 0 eV Emax=297.505 keV BetheBloch : Emin=297.505 keV Emax= 100 TeV hBrems: for pi- XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 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 100 eV to 100 TeV in 240 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 CoulombScat: for pi- XStype:1 applyCuts:1 SubType=1 BuildTable=1 Used Lambda table of pi+ ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV ========= Table of registered couples ============================ Index : 0 used in the geometry : Yes Material : G4_AIR Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 30 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes Material : G4_CESIUM_IODIDE Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um Energy thresholds : gamma 18.2699 keV e- 291.571 keV e+ 281.583 keV proton 30 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes Material : G4_Al Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um Energy thresholds : gamma 3.9261 keV e- 260.03 keV e+ 253.129 keV proton 30 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 3 used in the geometry : Yes Material : G4_AIR 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 : VertexDetector Index : 4 used in the geometry : Yes Material : G4_Al Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um Energy thresholds : gamma 990 eV e- 34.1058 keV e+ 33.7225 keV proton 1 keV Region(s) which use this couple : VertexDetector Index : 5 used in the geometry : Yes Material : G4_Si 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 Region(s) which use this couple : VertexDetector Index : 6 used in the geometry : Yes Material : G4_Fe Range cuts : gamma 1 cm e- 1 cm e+ 1 cm proton 1 cm Energy thresholds : gamma 60.6431 keV e- 13.0433 MeV e+ 12.1452 MeV proton 1 MeV Region(s) which use this couple : MuonDetector Index : 7 used in the geometry : Yes Material : G4_Al Range cuts : gamma 1 cm e- 1 cm e+ 1 cm proton 1 cm Energy thresholds : gamma 20.4222 keV e- 4.56335 MeV e+ 4.27483 MeV proton 1 MeV Region(s) which use this couple : MuonDetector ================================================================== ### Run 0 starts. ### Run 0 start ### PrimaryGeneratorAction::GeneratePrimaries ########## ### e- E(MeV)= 1000 ######################################################## Run terminated. Run Summary Number of events processed : 100 User=0.310000s Real=0.311598s Sys=0.000000s RunAction: End of run actions are started HistoManager: End of run actions are started RunID# 0 ================================================================= Number of events 100 Average number of e- 403.7 Average number of gamma 635.5 Average number of e+ 47.22 Average number of steps 2946 Edep 1x1 = 0.8329 +- 0.003391 res= 4.071 % 100 Edep 3x3 = 0.959 +- 0.001318 res= 1.374 % 100 Edep 5x5 = 0.9796 +- 0.0009483 res= 0.9681 % 100 =========== Ratios without trancating =========================== E1/E9 = 0.8685 +- 0.003149 E1/E25 = 0.8503 +- 0.003311 E9/E25 = 0.979 +- 0.0008637 Beam Energy 1 GeV ================================================================== <<<<>>>> 100 events for Crystal Calorimeter Edep1x1: 0.8329 delEdep1x1= -0.1671 nrms= -1.671 Erms1x1: 0.03391 delErms1x1= -0.9661 nrms= -9.661 Edep3x3: 0.959 delEdep3x3= -0.04099 nrms= -0.4099 Erms3x3: 0.01318 delErms3x3= -0.9868 nrms= -9.868 Edep5x5: 0.9796 delEdep5x5= -0.02042 nrms= -0.2042 Erms5x5: 0.009483 delErms5x5= -0.9905 nrms= -9.905 <<<<>>>> IS ACCEPTED Z bremsstrahlung photoeffect compton conversion 13 0 0 160 0 53 261 42 0 44 55 278 48 0 48 # /testem/det/acceptance1 0.0136 0.0001 500 /testem/det/acceptance9 0.0139 0.0001 500 /testem/det/acceptance25 0.0139 0.0001 500 # /gun/particle mu- /gun/energy 20 GeV /run/beamOn 1000 ========= Table of registered couples ============================ Index : 0 used in the geometry : Yes Material : G4_AIR Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 30 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes Material : G4_CESIUM_IODIDE Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um Energy thresholds : gamma 18.27 keV e- 291.6 keV e+ 281.6 keV proton 30 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes Material : G4_Al Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um Energy thresholds : gamma 3.926 keV e- 260 keV e+ 253.1 keV proton 30 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 3 used in the geometry : Yes Material : G4_AIR 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 : VertexDetector Index : 4 used in the geometry : Yes Material : G4_Al Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um Energy thresholds : gamma 990 eV e- 34.11 keV e+ 33.72 keV proton 1 keV Region(s) which use this couple : VertexDetector Index : 5 used in the geometry : Yes Material : G4_Si Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um Energy thresholds : gamma 990 eV e- 31.7 keV e+ 31.35 keV proton 1 keV Region(s) which use this couple : VertexDetector Index : 6 used in the geometry : Yes Material : G4_Fe Range cuts : gamma 1 cm e- 1 cm e+ 1 cm proton 1 cm Energy thresholds : gamma 60.64 keV e- 13.04 MeV e+ 12.15 MeV proton 1 MeV Region(s) which use this couple : MuonDetector Index : 7 used in the geometry : Yes Material : G4_Al Range cuts : gamma 1 cm e- 1 cm e+ 1 cm proton 1 cm Energy thresholds : gamma 20.42 keV e- 4.563 MeV e+ 4.275 MeV proton 1 MeV Region(s) which use this couple : MuonDetector ================================================================== ### Run 1 starts. ### Run 1 start ### PrimaryGeneratorAction::GeneratePrimaries ########## ### mu- E(MeV)= 2e+04 ######################################################## Run terminated. Run Summary Number of events processed : 1000 User=0.470000s Real=0.472948s Sys=0.000000s RunAction: End of run actions are started HistoManager: End of run actions are started RunID# 1 ================================================================= Number of events 1000 Average number of e- 69.02 Average number of gamma 61.77 Average number of e+ 3.389 Average number of steps 421.8 Edep 1x1 = 0.0141 +- 0.0003603 Edep 3x3 = 0.01445 +- 0.0004103 Edep 5x5 = 0.01451 +- 0.0004186 =========== Ratios without trancating =========================== E1/E9 = 0.984 +- 0.0006927 E1/E25 = 0.9811 +- 0.0007849 E9/E25 = 0.997 +- 0.0001965 Beam Energy 20 GeV ================================================================== # /process/em/setSecBiasing eBrem World 0.1 2 MeV /run/physicsModified /process/em/setSecBiasing eIoni World 0.0 1 MeV /run/physicsModified /process/em/setSecBiasing compt World 0.0 1 MeV /run/physicsModified /process/em/setSecBiasing phot World 0.0 1 MeV /run/physicsModified /process/em/setSecBiasing conv World 0.0 1 MeV /run/physicsModified /run/physicsModified /gun/particle e- /gun/energy 1 GeV /run/beamOn 100 ### === Deexcitation model UAtomDeexcitation is activated for 3 regions: DefaultRegionForTheWorld 0 0 0 VertexDetector 1 0 0 MuonDetector 0 0 0 ### === Ignore cuts flag: 0 Secondary biasing is activated for e+ and eIoni inside G4Regions: DefaultRegionForTheWorld BiasingWeight= 0 Secondary biasing is activated for e+ and eBrem inside G4Regions: DefaultRegionForTheWorld BiasingWeight= 10 Secondary biasing is activated for e- and eIoni inside G4Regions: DefaultRegionForTheWorld BiasingWeight= 0 Secondary biasing is activated for e- and eBrem inside G4Regions: DefaultRegionForTheWorld BiasingWeight= 10 Secondary biasing is activated for gamma and phot inside G4Regions: DefaultRegionForTheWorld BiasingWeight= 0 Secondary biasing is activated for gamma and compt inside G4Regions: DefaultRegionForTheWorld BiasingWeight= 0 Secondary biasing is activated for gamma and conv inside G4Regions: DefaultRegionForTheWorld BiasingWeight= 0 ========= Table of registered couples ============================ Index : 0 used in the geometry : Yes Material : G4_AIR Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 30 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes Material : G4_CESIUM_IODIDE Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um Energy thresholds : gamma 18.27 keV e- 291.6 keV e+ 281.6 keV proton 30 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes Material : G4_Al Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um Energy thresholds : gamma 3.926 keV e- 260 keV e+ 253.1 keV proton 30 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 3 used in the geometry : Yes Material : G4_AIR 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 : VertexDetector Index : 4 used in the geometry : Yes Material : G4_Al Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um Energy thresholds : gamma 990 eV e- 34.11 keV e+ 33.72 keV proton 1 keV Region(s) which use this couple : VertexDetector Index : 5 used in the geometry : Yes Material : G4_Si Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um Energy thresholds : gamma 990 eV e- 31.7 keV e+ 31.35 keV proton 1 keV Region(s) which use this couple : VertexDetector Index : 6 used in the geometry : Yes Material : G4_Fe Range cuts : gamma 1 cm e- 1 cm e+ 1 cm proton 1 cm Energy thresholds : gamma 60.64 keV e- 13.04 MeV e+ 12.15 MeV proton 1 MeV Region(s) which use this couple : MuonDetector Index : 7 used in the geometry : Yes Material : G4_Al Range cuts : gamma 1 cm e- 1 cm e+ 1 cm proton 1 cm Energy thresholds : gamma 20.42 keV e- 4.563 MeV e+ 4.275 MeV proton 1 MeV Region(s) which use this couple : MuonDetector ================================================================== ### Run 2 starts. ### Run 2 start ### PrimaryGeneratorAction::GeneratePrimaries ########## ### e- E(MeV)= 1000 ######################################################## Run terminated. Run Summary Number of events processed : 100 User=0.220000s Real=0.217922s Sys=0.000000s RunAction: End of run actions are started HistoManager: End of run actions are started RunID# 2 ================================================================= Number of events 100 Average number of e- 173.2 Average number of gamma 223 Average number of e+ 48.24 Average number of steps 1838 Edep 1x1 = 0.8322 +- 0.004077 res= 4.898 % 100 Edep 3x3 = 0.9578 +- 0.003436 res= 3.587 % 100 Edep 5x5 = 0.9806 +- 0.003509 res= 3.579 % 100 =========== Ratios without trancating =========================== E1/E9 = 0.869 +- 0.003095 E1/E25 = 0.8488 +- 0.003257 E9/E25 = 0.9768 +- 0.001139 Beam Energy 1 GeV ================================================================== # /process/em/deexcitation VertexDetector true false true /run/physicsModified /process/em/fluo true /run/physicsModified /process/em/pixe true /run/physicsModified #/process/em/auger false /run/physicsModified # /run/beamOn 100 ### === Deexcitation model UAtomDeexcitation is activated for 3 regions: DefaultRegionForTheWorld 0 0 0 VertexDetector 1 0 1 MuonDetector 0 0 0 ### === G4UAtomicDeexcitation::InitialiseForNewRun() ### === Ignore cuts flag: 0 ### === PIXE model for hadrons: Empirical ### === PIXE model for e+-: Livermore Secondary biasing is activated for e+ and eIoni inside G4Regions: DefaultRegionForTheWorld BiasingWeight= 0 Secondary biasing is activated for e+ and eBrem inside G4Regions: DefaultRegionForTheWorld BiasingWeight= 10 Secondary biasing is activated for e- and eIoni inside G4Regions: DefaultRegionForTheWorld BiasingWeight= 0 Secondary biasing is activated for e- and eBrem inside G4Regions: DefaultRegionForTheWorld BiasingWeight= 10 Secondary biasing is activated for gamma and phot inside G4Regions: DefaultRegionForTheWorld BiasingWeight= 0 Secondary biasing is activated for gamma and compt inside G4Regions: DefaultRegionForTheWorld BiasingWeight= 0 Secondary biasing is activated for gamma and conv inside G4Regions: DefaultRegionForTheWorld BiasingWeight= 0 ========= Table of registered couples ============================ Index : 0 used in the geometry : Yes Material : G4_AIR Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 30 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes Material : G4_CESIUM_IODIDE Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um Energy thresholds : gamma 18.27 keV e- 291.6 keV e+ 281.6 keV proton 30 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes Material : G4_Al Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um Energy thresholds : gamma 3.926 keV e- 260 keV e+ 253.1 keV proton 30 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 3 used in the geometry : Yes Material : G4_AIR 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 : VertexDetector Index : 4 used in the geometry : Yes Material : G4_Al Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um Energy thresholds : gamma 990 eV e- 34.11 keV e+ 33.72 keV proton 1 keV Region(s) which use this couple : VertexDetector Index : 5 used in the geometry : Yes Material : G4_Si Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um Energy thresholds : gamma 990 eV e- 31.7 keV e+ 31.35 keV proton 1 keV Region(s) which use this couple : VertexDetector Index : 6 used in the geometry : Yes Material : G4_Fe Range cuts : gamma 1 cm e- 1 cm e+ 1 cm proton 1 cm Energy thresholds : gamma 60.64 keV e- 13.04 MeV e+ 12.15 MeV proton 1 MeV Region(s) which use this couple : MuonDetector Index : 7 used in the geometry : Yes Material : G4_Al Range cuts : gamma 1 cm e- 1 cm e+ 1 cm proton 1 cm Energy thresholds : gamma 20.42 keV e- 4.563 MeV e+ 4.275 MeV proton 1 MeV Region(s) which use this couple : MuonDetector ================================================================== ### Run 3 starts. ### Run 3 start ### PrimaryGeneratorAction::GeneratePrimaries ########## ### e- E(MeV)= 1000 ######################################################## Run terminated. Run Summary Number of events processed : 100 User=0.220000s Real=0.223399s Sys=0.000000s RunAction: End of run actions are started HistoManager: End of run actions are started RunID# 3 ================================================================= Number of events 100 Average number of e- 172.7 Average number of gamma 223 Average number of e+ 47.33 Average number of steps 1842 Edep 1x1 = 0.8241 +- 0.003772 res= 4.577 % 100 Edep 3x3 = 0.9537 +- 0.003438 res= 3.605 % 100 Edep 5x5 = 0.9756 +- 0.003368 res= 3.452 % 100 =========== Ratios without trancating =========================== E1/E9 = 0.8642 +- 0.00281 E1/E25 = 0.8449 +- 0.003046 E9/E25 = 0.9775 +- 0.001322 Beam Energy 1 GeV ================================================================== # ================== Deleting memory pools =================== Number of memory pools allocated: 9 of which, static: 0 Dynamic pools deleted: 9 / Total memory freed: 0.09 MB ============================================================