Warning, /geant4/examples/extended/electromagnetic/TestEm16/TestEm16.out is written in an unsupported language. File is not indexed.
0001 Environment variable "G4FORCE_RUN_MANAGER_TYPE" enabled with value == Serial. Forcing G4RunManager type...
0002
0003 ############################################
0004 !!! WARNING - FPE detection is activated !!!
0005 ############################################
0006
0007
0008 ################################
0009 !!! G4Backtrace is activated !!!
0010 ################################
0011
0012
0013 **************************************************************
0014 Geant4 version Name: geant4-11-03-ref-06 (30-June-2025)
0015 Copyright : Geant4 Collaboration
0016 References : NIM A 506 (2003), 250-303
0017 : IEEE-TNS 53 (2006), 270-278
0018 : NIM A 835 (2016), 186-225
0019 WWW : http://geant4.org/
0020 **************************************************************
0021
0022
0023 ***** Table : Nb of materials = 4 *****
0024
0025 Material: Beryllium density: 1.848 g/cm3 RadL: 35.276 cm Nucl.Int.Length: 39.413 cm
0026 Imean: 63.700 eV temperature: 293.15 K pressure: 1.00 atm
0027
0028 ---> Element: Be (Be) Z = 4.0 N = 9 A = 9.012 g/mole
0029 ---> Isotope: Be9 Z = 4 N = 9 A = 9.01 g/mole abundance: 100.000 %
0030 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0031
0032
0033 Material: Carbon density: 2.265 g/cm3 RadL: 18.851 cm Nucl.Int.Length: 35.388 cm
0034 Imean: 81.000 eV temperature: 293.15 K pressure: 1.00 atm
0035
0036 ---> Element: C (C) Z = 6.0 N = 12 A = 12.011 g/mole
0037 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.930 %
0038 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.070 %
0039 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0040
0041
0042 Material: Iron density: 7.870 g/cm3 RadL: 1.758 cm Nucl.Int.Length: 16.999 cm
0043 Imean: 286.000 eV temperature: 293.15 K pressure: 1.00 atm
0044
0045 ---> Element: Fe (Fe) Z = 26.0 N = 56 A = 55.845 g/mole
0046 ---> Isotope: Fe54 Z = 26 N = 54 A = 53.94 g/mole abundance: 5.845 %
0047 ---> Isotope: Fe56 Z = 26 N = 56 A = 55.93 g/mole abundance: 91.754 %
0048 ---> Isotope: Fe57 Z = 26 N = 57 A = 56.94 g/mole abundance: 2.119 %
0049 ---> Isotope: Fe58 Z = 26 N = 58 A = 57.93 g/mole abundance: 0.282 %
0050 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0051
0052
0053 Material: Vacuum density: 0.000 mg/cm3 RadL: 25045177323.588 km Nucl.Int.Length: 55623331308.548 km
0054 Imean: 82.000 eV temperature: 296.15 K pressure: 0.00 atm
0055
0056 ---> Element: N (N) Z = 7.0 N = 14 A = 14.007 g/mole
0057 ---> Isotope: N14 Z = 7 N = 14 A = 14.00 g/mole abundance: 99.632 %
0058 ---> Isotope: N15 Z = 7 N = 15 A = 15.00 g/mole abundance: 0.368 %
0059 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0060
0061
0062
0063 /run/numberOfThreads 2
0064 *** /run/numberOfThreads command is issued in sequential mode.
0065 Command is ignored.
0066 /control/cout/ignoreThreadsExcept 0
0067 /run/verbose 1
0068 #
0069 /testem/det/setMat Vacuum
0070 #
0071 /run/initialize
0072
0073 The Box is 500 m of Vacuum
0074 #
0075 # inactivate other processes to make SynchrotronRadiation more visible
0076 /process/inactivate msc
0077 /run/physicsModified
0078 /process/inactivate eIoni
0079 /run/physicsModified
0080 /process/inactivate eBrem
0081 /run/physicsModified
0082 /process/inactivate annihil
0083 /run/physicsModified
0084 /process/inactivate phot
0085 /run/physicsModified
0086 /process/inactivate compt
0087 /run/physicsModified
0088 /process/inactivate conv
0089 /run/physicsModified
0090 /process/inactivate Rayl
0091 /run/physicsModified
0092 #
0093 /globalField/setValue 0 0 1 tesla
0094 G4ChordFinder: stepperDriverId: 2
0095 #
0096 /gun/particle e-
0097 /gun/energy 9.9994890009 GeV
0098 #
0099 # expected results is:
0100 # Emean = 20.48 keV
0101 # E_rms = 37.18 keV
0102 #
0103 /run/beamOn 100
0104 =======================================================================
0105 ====== Electromagnetic Physics Parameters ========
0106 =======================================================================
0107 LPM effect enabled 1
0108 Enable creation and use of sampling tables 0
0109 Apply cuts on all EM processes 0
0110 Use combined TransportationWithMsc Disabled
0111 Use general process 0
0112 Enable linear polarisation for gamma 0
0113 Enable photoeffect sampling below K-shell 1
0114 Enable sampling of quantum entanglement 0
0115 X-section factor for integral approach 0.8
0116 Min kinetic energy for tables 100 eV
0117 Max kinetic energy for tables 100 TeV
0118 Number of bins per decade of a table 7
0119 Verbose level 1
0120 Verbose level for worker thread 0
0121 Bremsstrahlung energy threshold above which
0122 primary e+- is added to the list of secondary 100 TeV
0123 Bremsstrahlung energy threshold above which primary
0124 muon/hadron is added to the list of secondary 100 TeV
0125 Positron annihilation at rest model SimplePositronium
0126 Enable 3 gamma annihilation on fly 0
0127 Lowest triplet kinetic energy 1 MeV
0128 Enable sampling of gamma linear polarisation 0
0129 5D gamma conversion model type 0
0130 5D gamma conversion model on isolated ion 0
0131 Use RiGe 5D e+e- pair production model by muons 0
0132 Livermore data directory epics_2017
0133 =======================================================================
0134 ====== Ionisation Parameters ========
0135 =======================================================================
0136 Step function for e+- (0.2, 1 mm)
0137 Step function for muons/hadrons (0.2, 0.1 mm)
0138 Step function for light ions (0.2, 0.1 mm)
0139 Step function for general ions (0.2, 0.1 mm)
0140 Lowest e+e- kinetic energy 1 keV
0141 Lowest muon/hadron kinetic energy 1 keV
0142 Use ICRU90 data 0
0143 Fluctuations of dE/dx are enabled 1
0144 Type of fluctuation model for leptons and hadrons Universal
0145 Use built-in Birks satuaration 0
0146 Build CSDA range enabled 0
0147 Use cut as a final range enabled 0
0148 Enable angular generator interface 0
0149 Max kinetic energy for CSDA tables 1 GeV
0150 Max kinetic energy for NIEL computation 0 eV
0151 Linear loss limit 0.01
0152 Read data from file for e+e- pair production by mu 0
0153 =======================================================================
0154 ====== Multiple Scattering Parameters ========
0155 =======================================================================
0156 Type of msc step limit algorithm for e+- 1
0157 Type of msc step limit algorithm for muons/hadrons 0
0158 Msc lateral displacement for e+- enabled 1
0159 Msc lateral displacement for muons and hadrons 0
0160 Urban msc model lateral displacement alg96 1
0161 Range factor for msc step limit for e+- 0.04
0162 Range factor for msc step limit for muons/hadrons 0.2
0163 Geometry factor for msc step limitation of e+- 2.5
0164 Safety factor for msc step limit for e+- 0.6
0165 Skin parameter for msc step limitation of e+- 1
0166 Lambda limit for msc step limit for e+- 1 mm
0167 Use Mott correction for e- scattering 0
0168 Factor used for dynamic computation of angular
0169 limit between single and multiple scattering 1
0170 Fixed angular limit between single
0171 and multiple scattering 3.1416 rad
0172 Upper energy limit for e+- multiple scattering 100 MeV
0173 Type of electron single scattering model 0
0174 Type of nuclear form-factor 1
0175 Screening factor 1
0176 =======================================================================
0177
0178 phot: for gamma SubType=12 BuildTable=0
0179 LambdaPrime table from 200 keV to 100 TeV in 61 bins
0180 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0181 PhotoElectric : Emin= 0 eV Emax= 100 TeV SauterGavrila
0182
0183 compt: for gamma SubType=13 BuildTable=1
0184 Lambda table from 100 eV to 1 MeV, 7 bins/decade, spline: 1
0185 LambdaPrime table from 1 MeV to 100 TeV in 56 bins
0186 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0187 Klein-Nishina : Emin= 0 eV Emax= 100 TeV
0188
0189 conv: for gamma SubType=14 BuildTable=1
0190 Lambda table from 1.022 MeV to 100 TeV, 18 bins/decade, spline: 1
0191 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0192 BetheHeitlerLPM : Emin= 0 eV Emax= 100 TeV ModifiedTsai
0193
0194 Rayl: for gamma SubType=11 BuildTable=1
0195 Lambda table from 100 eV to 150 keV, 7 bins/decade, spline: 0
0196 LambdaPrime table from 150 keV to 100 TeV in 62 bins
0197 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0198 LivermoreRayleigh : Emin= 0 eV Emax= 100 TeV CullenGenerator
0199
0200 XrayReflection: Gamma specular reflection for energies > 30 eV.
0201
0202 msc: for e- SubType= 10
0203 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0204 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV
0205 StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0206
0207 eIoni: for e- XStype:3 SubType=2
0208 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0209 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0210 StepFunction=(0.2, 1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0211 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0212 MollerBhabha : Emin= 0 eV Emax= 100 TeV
0213
0214 eBrem: for e- XStype:4 SubType=3
0215 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0216 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0217 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000
0218 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0219 eBremSB : Emin= 0 eV Emax= 1 GeV ModifiedTsai
0220 eBremLPM : Emin= 1 GeV Emax= 100 TeV ModifiedTsai
0221 SynRad: Incoherent Synchrotron Radiation
0222 Good description for long magnets at all energies.
0223
0224 msc: for e+ SubType= 10
0225 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0226 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV
0227 StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0228
0229 eIoni: for e+ XStype:3 SubType=2
0230 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0231 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0232 StepFunction=(0.2, 1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0233 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0234 MollerBhabha : Emin= 0 eV Emax= 100 TeV
0235
0236 eBrem: for e+ XStype:4 SubType=3
0237 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0238 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0239 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000
0240 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0241 eBremSB : Emin= 0 eV Emax= 1 GeV ModifiedTsai
0242 eBremLPM : Emin= 1 GeV Emax= 100 TeV ModifiedTsai
0243
0244 annihil: for e+ XStype:2 SubType=5 AtRestModel:Simple BuildTable=0
0245 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0246 eplus2gg : Emin= 0 eV Emax= 100 TeV
0247
0248 msc: for proton SubType= 10
0249 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0250 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV
0251 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0252
0253 hIoni: for proton XStype:3 SubType=2
0254 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0255 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0256 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0257 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0258 Bragg : Emin= 0 eV Emax= 2 MeV
0259 BetheBloch : Emin= 2 MeV Emax= 100 TeV
0260
0261 hBrems: for proton XStype:1 SubType=3
0262 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0263 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0264 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0265 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0266
0267 hPairProd: for proton XStype:1 SubType=4
0268 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0269 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0270 Sampling table 17x1001 from 7.50618 GeV to 100 TeV
0271 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0272 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0273
0274 msc: for mu+ SubType= 10
0275 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0276 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV
0277 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0278
0279 muIoni: for mu+ XStype:3 SubType=2
0280 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0281 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0282 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0283 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0284 Bragg : Emin= 0 eV Emax= 200 keV
0285 MuBetheBloch : Emin= 200 keV Emax= 100 TeV
0286
0287 muBrems: for mu+ XStype:1 SubType=3
0288 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0289 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0290 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0291 MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0292
0293 muPairProd: for mu+ XStype:1 SubType=4
0294 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0295 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0296 Sampling table 21x1001 from 0.85 GeV to 100 TeV
0297 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0298 muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0299
0300 msc: for mu- SubType= 10
0301 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0302 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV
0303 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0304
0305 muIoni: for mu- XStype:3 SubType=2
0306 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0307 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0308 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0309 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0310 ICRU73QO : Emin= 0 eV Emax= 200 keV
0311 MuBetheBloch : Emin= 200 keV Emax= 100 TeV
0312
0313 ========= Table of registered couples ============================
0314
0315 Index : 0 used in the geometry : Yes
0316 Material : Vacuum
0317 Range cuts : gamma 1 km e- 1 km e+ 1 km proton 1 km
0318 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 GeV
0319 Region(s) which use this couple :
0320 DefaultRegionForTheWorld
0321
0322 ==================================================================
0323
0324 ### Run 0 starts.
0325
0326 ------- MixMaxRng engine status -------
0327 Current state vector is:
0328 mixmax state, file version 1.0
0329 N=17 V[N]={906770732717044781, 629165745432651234, 1235682547346241386, 68420008233404568, 2262190991329242458, 2266470399991071809, 1976726662926872232, 245458862506414172, 1955974201201518530, 2155248512522080758, 604170912935414061, 1116171330120743511, 1861018313684488333, 1296715403254578286, 1549011045957234151, 370819759640195970, 2230139271784837643} counter= 17sumtot= 1977567618660788324
0330 ---------------------------------------
0331 G4SynchrotronRadiation::GetMeanFreePath for particle e-:
0332 MeanFreePath = 16.1828 cm
0333 G4SynchrotronRadiation::GetRandomEnergySR :
0334 Ecr = 66.5 keV
0335 Emean = 20.48 keV
0336 E_rms = 37.18 keV
0337 Run terminated.
0338 Run Summary
0339 Number of events processed : 100
0340 User=0.120000s Real=0.124169s Sys=0.000000s
0341 Summary for synchrotron radiation :
0342 Number of photons = 64566
0343 Emean = 20.42 +/- 0.1465 keV
0344 E_rms = 37.22 keV
0345 Energy Max / Mean = 35.69
0346 MeanFreePath = 16.2 cm
0347
0348 ------- MixMaxRng engine status -------
0349 Current state vector is:
0350 mixmax state, file version 1.0
0351 N=17 V[N]={2257445806913462277, 632527905652818191, 631474175170306566, 1914435244802913599, 1306329175570610787, 762881766574265544, 719637845945592025, 1589777450882882221, 1702736512761014061, 890200915960239409, 384368303897083373, 453040775896457888, 307076404477827112, 1945929717505046236, 544360434503550427, 215256680747272403, 182941379877301700} counter= 3sumtot= 299519432642786162
0352 ---------------------------------------
0353 ================== Deleting memory pools ===================
0354 Number of memory pools allocated: 9 of which, static: 0
0355 Dynamic pools deleted: 9 / Total memory freed: 0.35 MB
0356 ============================================================
