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/setMat G4_PbWO4 /testem/det/setLbin 40 0.5 /testem/det/setRbin 50 0.1 /testem/run/acceptance 0.949 0.0082 3.0 /testem/phys/addPhysics emstandard_opt0 /run/setCut 1.0 mm /process/em/verbose 1 /process/em/workerVerbose 0 /process/em/printParameters ======================================================================= ====== 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 100 eV Max kinetic energy for tables 100 TeV Number of bins per decade of a table 7 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.2, 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 Urban 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+- 1 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.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 0 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 ======================================================================= /run/initialize Absorber is 17.8491 cm of G4_PbWO4 R= 4.46227 cm Material: G4_PbWO4 density: 8.280 g/cm3 RadL: 8.925 mm Nucl.Int.Length: 20.740 cm Imean: 542.741 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: O (O) Z = 8.0 N = 16 A = 15.999 g/mole ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 % ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 % ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 % ElmMassFraction: 14.06 % ElmAbundance 66.67 % ---> Element: Pb (Pb) Z = 82.0 N = 207 A = 207.217 g/mole ---> Isotope: Pb204 Z = 82 N = 204 A = 203.97 g/mole abundance: 1.400 % ---> Isotope: Pb206 Z = 82 N = 206 A = 205.97 g/mole abundance: 24.100 % ---> Isotope: Pb207 Z = 82 N = 207 A = 206.98 g/mole abundance: 22.100 % ---> Isotope: Pb208 Z = 82 N = 208 A = 207.98 g/mole abundance: 52.400 % ElmMassFraction: 45.54 % ElmAbundance 16.67 % ---> Element: W (W) Z = 74.0 N = 184 A = 183.842 g/mole ---> Isotope: W180 Z = 74 N = 180 A = 179.95 g/mole abundance: 0.120 % ---> Isotope: W182 Z = 74 N = 182 A = 181.95 g/mole abundance: 26.500 % ---> Isotope: W183 Z = 74 N = 183 A = 182.95 g/mole abundance: 14.310 % ---> Isotope: W184 Z = 74 N = 184 A = 183.95 g/mole abundance: 30.640 % ---> Isotope: W186 Z = 74 N = 186 A = 185.95 g/mole abundance: 28.430 % ElmMassFraction: 40.40 % ElmAbundance 16.67 % /gun/particle e- /gun/energy 5 GeV /run/printProgress 20 /testem/run/verbose 1 /run/beamOn 100 ======================================================================= ====== 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 100 eV Max kinetic energy for tables 100 TeV Number of bins per decade of a table 7 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.2, 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 Urban 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+- 1 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.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 0 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 ======================================================================= phot: for gamma SubType=12 BuildTable=0 LambdaPrime table from 200 keV to 100 TeV in 61 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 100 eV to 1 MeV, 7 bins/decade, spline: 1 LambdaPrime table from 1 MeV to 100 TeV in 56 bins ===== EM models for the G4Region DefaultRegionForTheWorld ====== Klein-Nishina : Emin= 0 eV Emax= 100 TeV conv: for gamma SubType=14 BuildTable=1 Lambda table from 1.022 MeV to 100 TeV, 18 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 100 eV to 150 keV, 7 bins/decade, spline: 0 LambdaPrime table from 150 keV to 100 TeV in 62 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 MeV Nbins=42 100 eV - 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 Nbins=42 100 MeV - 100 TeV StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm eIoni: for e- XStype:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 100 TeV eBrem: for e- XStype:3 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 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 SubType=1 BuildTable=1 Lambda table from 100 MeV to 100 TeV, 7 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 msc: for e+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 MeV Nbins=42 100 eV - 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 Nbins=42 100 MeV - 100 TeV StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm eIoni: for e+ XStype:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 100 TeV eBrem: for e+ XStype:3 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 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 SubType=5 AtRestModel:Simple BuildTable=0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eplus2gg : Emin= 0 eV Emax= 100 TeV CoulombScat: for e+ XStype:1 SubType=1 BuildTable=1 Lambda table from 100 MeV to 100 TeV, 7 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 msc: for proton SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=84 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:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, 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 84 bins Lambda tables from threshold to 100 TeV, 7 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 84 bins Lambda tables from threshold to 100 TeV, 7 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 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 7 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:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, 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:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, 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=84 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:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, 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 84 bins Lambda tables from threshold to 100 TeV, 7 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 84 bins Lambda tables from threshold to 100 TeV, 7 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 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 7 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=84 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:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, 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 84 bins Lambda tables from threshold to 100 TeV, 7 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 84 bins Lambda tables from threshold to 100 TeV, 7 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 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 7 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=84 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:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, 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 84 bins Lambda tables from threshold to 100 TeV, 7 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 84 bins Lambda tables from threshold to 100 TeV, 7 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 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=84 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:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, 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 84 bins Lambda tables from threshold to 100 TeV, 7 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 84 bins Lambda tables from threshold to 100 TeV, 7 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 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 7 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=84 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:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, 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 84 bins Lambda tables from threshold to 100 TeV, 7 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 84 bins Lambda tables from threshold to 100 TeV, 7 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 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=84 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:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, 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 84 bins Lambda tables from threshold to 100 TeV, 7 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 84 bins Lambda tables from threshold to 100 TeV, 7 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 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 7 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=84 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:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, 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 84 bins Lambda tables from threshold to 100 TeV, 7 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 84 bins Lambda tables from threshold to 100 TeV, 7 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 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_PbWO4 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 89.5666 keV e- 1.13293 MeV e+ 1.06705 MeV proton 100 keV 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 --------------------------------------- ... set ntuple merging row mode : row-wise - done ... create file : testem2.root - done ... open analysis file : testem2.root - done ... open analysis file : testem2.root - done ----> Histogram file is opened in testem2.root. --> Event 0 starts. --> Event 20 starts. --> Event 40 starts. --> Event 60 starts. --> Event 80 starts. Run terminated. Run Summary Number of events processed : 100 User=1.050000s Real=1.054195s Sys=0.000000s LOGITUDINAL PROFILE CUMULATIVE LOGITUDINAL PROFILE bin Mean rms bin Mean rms 0.00-> 0.50 radl: 0.19% 0.10% 0-> 0.50 radl: 0.19% 0.10% 0.50-> 1.00 radl: 0.44% 0.28% 0-> 1.00 radl: 0.64% 0.35% 1.00-> 1.50 radl: 0.90% 0.55% 0-> 1.50 radl: 1.53% 0.80% 1.50-> 2.00 radl: 1.50% 0.87% 0-> 2.00 radl: 3.03% 1.57% 2.00-> 2.50 radl: 2.17% 1.09% 0-> 2.50 radl: 5.20% 2.54% 2.50-> 3.00 radl: 3.01% 1.34% 0-> 3.00 radl: 8.21% 3.71% 3.00-> 3.50 radl: 3.70% 1.48% 0-> 3.50 radl: 11.91% 4.92% 3.50-> 4.00 radl: 4.31% 1.52% 0-> 4.00 radl: 16.23% 6.16% 4.00-> 4.50 radl: 4.92% 1.59% 0-> 4.50 radl: 21.15% 7.28% 4.50-> 5.00 radl: 5.16% 1.66% 0-> 5.00 radl: 26.31% 8.29% 5.00-> 5.50 radl: 5.44% 1.58% 0-> 5.50 radl: 31.75% 9.38% 5.50-> 6.00 radl: 5.41% 1.56% 0-> 6.00 radl: 37.16% 10.30% 6.00-> 6.50 radl: 5.23% 1.30% 0-> 6.50 radl: 42.39% 11.03% 6.50-> 7.00 radl: 5.18% 1.26% 0-> 7.00 radl: 47.57% 11.50% 7.00-> 7.50 radl: 5.06% 1.33% 0-> 7.50 radl: 52.63% 11.61% 7.50-> 8.00 radl: 4.72% 1.17% 0-> 8.00 radl: 57.35% 11.52% 8.00-> 8.50 radl: 4.42% 1.15% 0-> 8.50 radl: 61.78% 11.32% 8.50-> 9.00 radl: 4.20% 1.17% 0-> 9.00 radl: 65.97% 10.92% 9.00-> 9.50 radl: 3.73% 1.15% 0-> 9.50 radl: 69.70% 10.42% 9.50->10.00 radl: 3.30% 0.96% 0->10.00 radl: 73.01% 9.99% 10.00->10.50 radl: 2.99% 1.05% 0->10.50 radl: 75.99% 9.40% 10.50->11.00 radl: 2.62% 1.01% 0->11.00 radl: 78.61% 8.84% 11.00->11.50 radl: 2.40% 0.98% 0->11.50 radl: 81.02% 8.27% 11.50->12.00 radl: 2.08% 0.92% 0->12.00 radl: 83.10% 7.67% 12.00->12.50 radl: 1.80% 0.91% 0->12.50 radl: 84.90% 7.05% 12.50->13.00 radl: 1.56% 0.78% 0->13.00 radl: 86.46% 6.44% 13.00->13.50 radl: 1.32% 0.73% 0->13.50 radl: 87.77% 5.89% 13.50->14.00 radl: 1.14% 0.61% 0->14.00 radl: 88.92% 5.44% 14.00->14.50 radl: 1.00% 0.65% 0->14.50 radl: 89.92% 4.98% 14.50->15.00 radl: 0.85% 0.54% 0->15.00 radl: 90.77% 4.66% 15.00->15.50 radl: 0.76% 0.52% 0->15.50 radl: 91.54% 4.37% 15.50->16.00 radl: 0.65% 0.44% 0->16.00 radl: 92.18% 4.08% 16.00->16.50 radl: 0.50% 0.40% 0->16.50 radl: 92.68% 3.77% 16.50->17.00 radl: 0.46% 0.45% 0->17.00 radl: 93.14% 3.41% 17.00->17.50 radl: 0.39% 0.41% 0->17.50 radl: 93.53% 3.09% 17.50->18.00 radl: 0.33% 0.38% 0->18.00 radl: 93.86% 2.78% 18.00->18.50 radl: 0.32% 0.33% 0->18.50 radl: 94.18% 2.50% 18.50->19.00 radl: 0.25% 0.30% 0->19.00 radl: 94.43% 2.28% 19.00->19.50 radl: 0.20% 0.29% 0->19.50 radl: 94.63% 2.02% 19.50->20.00 radl: 0.18% 0.28% 0->20.00 radl: 94.81% 1.79% RADIAL PROFILE CUMULATIVE RADIAL PROFILE bin Mean rms bin Mean rms 0.00-> 0.10 radl: 17.53% 2.26% 0-> 0.10 radl: 17.53% 2.26% 0.10-> 0.20 radl: 13.19% 1.36% 0-> 0.20 radl: 30.72% 2.67% 0.20-> 0.30 radl: 9.91% 1.01% 0-> 0.30 radl: 40.62% 2.81% 0.30-> 0.40 radl: 7.63% 0.82% 0-> 0.40 radl: 48.25% 2.72% 0.40-> 0.50 radl: 6.06% 0.74% 0-> 0.50 radl: 54.31% 2.67% 0.50-> 0.60 radl: 4.85% 0.64% 0-> 0.60 radl: 59.16% 2.48% 0.60-> 0.70 radl: 3.98% 0.67% 0-> 0.70 radl: 63.14% 2.26% 0.70-> 0.80 radl: 3.26% 0.55% 0-> 0.80 radl: 66.41% 2.14% 0.80-> 0.90 radl: 2.79% 0.52% 0-> 0.90 radl: 69.20% 2.02% 0.90-> 1.00 radl: 2.30% 0.39% 0-> 1.00 radl: 71.49% 2.03% 1.00-> 1.10 radl: 2.05% 0.34% 0-> 1.10 radl: 73.55% 1.98% 1.10-> 1.20 radl: 1.80% 0.39% 0-> 1.20 radl: 75.35% 1.98% 1.20-> 1.30 radl: 1.59% 0.31% 0-> 1.30 radl: 76.94% 1.95% 1.30-> 1.40 radl: 1.41% 0.28% 0-> 1.40 radl: 78.35% 1.91% 1.40-> 1.50 radl: 1.26% 0.29% 0-> 1.50 radl: 79.61% 1.90% 1.50-> 1.60 radl: 1.19% 0.31% 0-> 1.60 radl: 80.80% 1.87% 1.60-> 1.70 radl: 1.03% 0.25% 0-> 1.70 radl: 81.83% 1.85% 1.70-> 1.80 radl: 0.94% 0.24% 0-> 1.80 radl: 82.77% 1.85% 1.80-> 1.90 radl: 0.85% 0.22% 0-> 1.90 radl: 83.61% 1.83% 1.90-> 2.00 radl: 0.80% 0.22% 0-> 2.00 radl: 84.41% 1.80% 2.00-> 2.10 radl: 0.76% 0.22% 0-> 2.10 radl: 85.17% 1.80% 2.10-> 2.20 radl: 0.72% 0.21% 0-> 2.20 radl: 85.89% 1.84% 2.20-> 2.30 radl: 0.62% 0.18% 0-> 2.30 radl: 86.52% 1.84% 2.30-> 2.40 radl: 0.62% 0.18% 0-> 2.40 radl: 87.13% 1.85% 2.40-> 2.50 radl: 0.60% 0.16% 0-> 2.50 radl: 87.73% 1.86% 2.50-> 2.60 radl: 0.53% 0.17% 0-> 2.60 radl: 88.27% 1.82% 2.60-> 2.70 radl: 0.51% 0.17% 0-> 2.70 radl: 88.78% 1.82% 2.70-> 2.80 radl: 0.46% 0.16% 0-> 2.80 radl: 89.24% 1.81% 2.80-> 2.90 radl: 0.45% 0.16% 0-> 2.90 radl: 89.69% 1.79% 2.90-> 3.00 radl: 0.41% 0.13% 0-> 3.00 radl: 90.11% 1.79% 3.00-> 3.10 radl: 0.38% 0.14% 0-> 3.10 radl: 90.49% 1.78% 3.10-> 3.20 radl: 0.37% 0.14% 0-> 3.20 radl: 90.86% 1.75% 3.20-> 3.30 radl: 0.34% 0.14% 0-> 3.30 radl: 91.20% 1.76% 3.30-> 3.40 radl: 0.32% 0.12% 0-> 3.40 radl: 91.52% 1.76% 3.40-> 3.50 radl: 0.32% 0.13% 0-> 3.50 radl: 91.84% 1.76% 3.50-> 3.60 radl: 0.29% 0.12% 0-> 3.60 radl: 92.13% 1.77% 3.60-> 3.70 radl: 0.27% 0.11% 0-> 3.70 radl: 92.40% 1.77% 3.70-> 3.80 radl: 0.26% 0.12% 0-> 3.80 radl: 92.66% 1.77% 3.80-> 3.90 radl: 0.26% 0.11% 0-> 3.90 radl: 92.92% 1.76% 3.90-> 4.00 radl: 0.22% 0.09% 0-> 4.00 radl: 93.14% 1.77% 4.00-> 4.10 radl: 0.23% 0.10% 0-> 4.10 radl: 93.37% 1.76% 4.10-> 4.20 radl: 0.20% 0.09% 0-> 4.20 radl: 93.57% 1.76% 4.20-> 4.30 radl: 0.21% 0.10% 0-> 4.30 radl: 93.77% 1.77% 4.30-> 4.40 radl: 0.19% 0.09% 0-> 4.40 radl: 93.96% 1.78% 4.40-> 4.50 radl: 0.17% 0.08% 0-> 4.50 radl: 94.13% 1.79% 4.50-> 4.60 radl: 0.15% 0.07% 0-> 4.60 radl: 94.28% 1.79% 4.60-> 4.70 radl: 0.16% 0.09% 0-> 4.70 radl: 94.44% 1.78% 4.70-> 4.80 radl: 0.14% 0.07% 0-> 4.80 radl: 94.58% 1.78% 4.80-> 4.90 radl: 0.12% 0.07% 0-> 4.90 radl: 94.71% 1.78% 4.90-> 5.00 radl: 0.10% 0.05% 0-> 5.00 radl: 94.81% 1.79% ===== SUMMARY ===== Total number of events: 100 Mean number of charged steps: 7242.32 Mean number of neutral steps: 4101.23 energy deposit : 94.81 % E0 +- 1.79 % E0 charged traklen: 483.38 radl +- 9.33 radl neutral traklen: 4206.39 radl +- 143.52 radl 90.00 % confinement: radius = 2.97 radl (2.65 cm ) <<<>>> 100 events for Total Energy in Absorber Edep: 0.948093 delEdep= -0.000906621 nrms= -1.10564 Erms: 0.0178703 delErms= 0.00967029 nrms= 11.793 <<<>>> IS NOT ACCEPTED ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={2217242971309508852, 1772600625655739516, 1083784801228937735, 1376652001525196202, 923881048880974269, 391157941837653129, 120038680055995158, 1455996709604147955, 444702727052845783, 2056176360790095207, 12388529709462456, 506878535865677272, 1486203893460794122, 1678230504276076655, 1948481978275914262, 1381823691342543834, 677819068304920907} counter= 14sumtot= 1087315995466931706 --------------------------------------- ... write file : testem2.root - done ... close file : testem2.root - done ... clear all data - done ================== Deleting memory pools =================== Number of memory pools allocated: 9 of which, static: 0 Dynamic pools deleted: 9 / Total memory freed: 0.059 MB ============================================================