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-02-ref-06 (28-June-2024)
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 Lowest triplet kinetic energy 1 MeV
0126 Enable sampling of gamma linear polarisation 0
0127 5D gamma conversion model type 0
0128 5D gamma conversion model on isolated ion 0
0129 Livermore data directory epics_2017
0130 =======================================================================
0131 ====== Ionisation Parameters ========
0132 =======================================================================
0133 Step function for e+- (0.2, 1 mm)
0134 Step function for muons/hadrons (0.2, 0.1 mm)
0135 Step function for light ions (0.2, 0.1 mm)
0136 Step function for general ions (0.2, 0.1 mm)
0137 Lowest e+e- kinetic energy 1 keV
0138 Lowest muon/hadron kinetic energy 1 keV
0139 Use ICRU90 data 0
0140 Fluctuations of dE/dx are enabled 1
0141 Type of fluctuation model for leptons and hadrons Universal
0142 Use built-in Birks satuaration 0
0143 Build CSDA range enabled 0
0144 Use cut as a final range enabled 0
0145 Enable angular generator interface 0
0146 Max kinetic energy for CSDA tables 1 GeV
0147 Max kinetic energy for NIEL computation 0 eV
0148 Linear loss limit 0.01
0149 Read data from file for e+e- pair production by mu 0
0150 =======================================================================
0151 ====== Multiple Scattering Parameters ========
0152 =======================================================================
0153 Type of msc step limit algorithm for e+- 1
0154 Type of msc step limit algorithm for muons/hadrons 0
0155 Msc lateral displacement for e+- enabled 1
0156 Msc lateral displacement for muons and hadrons 0
0157 Urban msc model lateral displacement alg96 1
0158 Range factor for msc step limit for e+- 0.04
0159 Range factor for msc step limit for muons/hadrons 0.2
0160 Geometry factor for msc step limitation of e+- 2.5
0161 Safety factor for msc step limit for e+- 0.6
0162 Skin parameter for msc step limitation of e+- 1
0163 Lambda limit for msc step limit for e+- 1 mm
0164 Use Mott correction for e- scattering 0
0165 Factor used for dynamic computation of angular
0166 limit between single and multiple scattering 1
0167 Fixed angular limit between single
0168 and multiple scattering 3.1416 rad
0169 Upper energy limit for e+- multiple scattering 100 MeV
0170 Type of electron single scattering model 0
0171 Type of nuclear form-factor 1
0172 Screening factor 1
0173 =======================================================================
0174
0175 phot: for gamma SubType=12 BuildTable=0
0176 LambdaPrime table from 200 keV to 100 TeV in 61 bins
0177 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0178 PhotoElectric : Emin= 0 eV Emax= 100 TeV SauterGavrila
0179
0180 compt: for gamma SubType=13 BuildTable=1
0181 Lambda table from 100 eV to 1 MeV, 7 bins/decade, spline: 1
0182 LambdaPrime table from 1 MeV to 100 TeV in 56 bins
0183 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0184 Klein-Nishina : Emin= 0 eV Emax= 100 TeV
0185
0186 conv: for gamma SubType=14 BuildTable=1
0187 Lambda table from 1.022 MeV to 100 TeV, 18 bins/decade, spline: 1
0188 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0189 BetheHeitlerLPM : Emin= 0 eV Emax= 100 TeV ModifiedTsai
0190
0191 Rayl: for gamma SubType=11 BuildTable=1
0192 Lambda table from 100 eV to 150 keV, 7 bins/decade, spline: 0
0193 LambdaPrime table from 150 keV to 100 TeV in 62 bins
0194 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0195 LivermoreRayleigh : Emin= 0 eV Emax= 100 TeV CullenGenerator
0196
0197 XrayReflection: Gamma specular reflection for energies > 30 eV.
0198
0199 msc: for e- SubType= 10
0200 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0201 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV
0202 StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0203
0204 eIoni: for e- XStype:3 SubType=2
0205 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0206 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0207 StepFunction=(0.2, 1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0208 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0209 MollerBhabha : Emin= 0 eV Emax= 100 TeV
0210
0211 eBrem: for e- XStype:4 SubType=3
0212 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0213 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0214 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000
0215 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0216 eBremSB : Emin= 0 eV Emax= 1 GeV ModifiedTsai
0217 eBremLPM : Emin= 1 GeV Emax= 100 TeV ModifiedTsai
0218 SynRad: Incoherent Synchrotron Radiation
0219 Good description for long magnets at all energies.
0220
0221 msc: for e+ SubType= 10
0222 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0223 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV
0224 StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0225
0226 eIoni: for e+ XStype:3 SubType=2
0227 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0228 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0229 StepFunction=(0.2, 1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0230 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0231 MollerBhabha : Emin= 0 eV Emax= 100 TeV
0232
0233 eBrem: for e+ XStype:4 SubType=3
0234 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0235 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0236 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000
0237 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0238 eBremSB : Emin= 0 eV Emax= 1 GeV ModifiedTsai
0239 eBremLPM : Emin= 1 GeV Emax= 100 TeV ModifiedTsai
0240
0241 annihil: for e+ XStype:2 SubType=5 AtRestModel:Simple BuildTable=0
0242 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0243 eplus2gg : Emin= 0 eV Emax= 100 TeV
0244
0245 msc: for proton SubType= 10
0246 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0247 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV
0248 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0249
0250 hIoni: for proton XStype:3 SubType=2
0251 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0252 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0253 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0254 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0255 Bragg : Emin= 0 eV Emax= 2 MeV
0256 BetheBloch : Emin= 2 MeV Emax= 100 TeV
0257
0258 hBrems: for proton XStype:1 SubType=3
0259 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0260 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0261 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0262 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0263
0264 hPairProd: for proton XStype:1 SubType=4
0265 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0266 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0267 Sampling table 17x1001 from 7.50618 GeV to 100 TeV
0268 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0269 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0270
0271 msc: for mu+ SubType= 10
0272 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0273 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV
0274 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0275
0276 muIoni: for mu+ XStype:3 SubType=2
0277 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0278 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0279 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0280 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0281 Bragg : Emin= 0 eV Emax= 200 keV
0282 MuBetheBloch : Emin= 200 keV Emax= 100 TeV
0283
0284 muBrems: for mu+ XStype:1 SubType=3
0285 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0286 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0287 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0288 MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0289
0290 muPairProd: for mu+ XStype:1 SubType=4
0291 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0292 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0293 Sampling table 21x1001 from 0.85 GeV to 100 TeV
0294 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0295 muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0296
0297 msc: for mu- SubType= 10
0298 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0299 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV
0300 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0301
0302 muIoni: for mu- XStype:3 SubType=2
0303 dE/dx and range tables from 100 eV to 100 TeV in 84 bins
0304 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
0305 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0306 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0307 ICRU73QO : Emin= 0 eV Emax= 200 keV
0308 MuBetheBloch : Emin= 200 keV Emax= 100 TeV
0309
0310 ========= Table of registered couples ============================
0311
0312 Index : 0 used in the geometry : Yes
0313 Material : Vacuum
0314 Range cuts : gamma 1 km e- 1 km e+ 1 km proton 1 km
0315 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 GeV
0316 Region(s) which use this couple :
0317 DefaultRegionForTheWorld
0318
0319 ==================================================================
0320
0321 ### Run 0 starts.
0322
0323 ------- MixMaxRng engine status -------
0324 Current state vector is:
0325 mixmax state, file version 1.0
0326 N=17 V[N]={906770732717044781, 629165745432651234, 1235682547346241386, 68420008233404568, 2262190991329242458, 2266470399991071809, 1976726662926872232, 245458862506414172, 1955974201201518530, 2155248512522080758, 604170912935414061, 1116171330120743511, 1861018313684488333, 1296715403254578286, 1549011045957234151, 370819759640195970, 2230139271784837643} counter= 17sumtot= 1977567618660788324
0327 ---------------------------------------
0328 G4SynchrotronRadiation::GetMeanFreePath for particle e-:
0329 MeanFreePath = 16.1828 cm
0330 G4SynchrotronRadiation::GetRandomEnergySR :
0331 Ecr = 66.5 keV
0332 Emean = 20.48 keV
0333 E_rms = 37.18 keV
0334 Run terminated.
0335 Run Summary
0336 Number of events processed : 100
0337 User=0.200000s Real=0.233569s Sys=0.000000s
0338 Summary for synchrotron radiation :
0339 Number of photons = 64566
0340 Emean = 20.42 +/- 0.1465 keV
0341 E_rms = 37.22 keV
0342 Energy Max / Mean = 35.69
0343 MeanFreePath = 16.2 cm
0344
0345 ------- MixMaxRng engine status -------
0346 Current state vector is:
0347 mixmax state, file version 1.0
0348 N=17 V[N]={2257445806913462277, 632527905652818191, 631474175170306566, 1914435244802913599, 1306329175570610787, 762881766574265544, 719637845945592025, 1589777450882882221, 1702736512761014061, 890200915960239409, 384368303897083373, 453040775896457888, 307076404477827112, 1945929717505046236, 544360434503550427, 215256680747272403, 182941379877301700} counter= 3sumtot= 299519432642786162
0349 ---------------------------------------
0350 ================== Deleting memory pools ===================
0351 Number of memory pools allocated: 9 of which, static: 0
0352 Dynamic pools deleted: 9 / Total memory freed: 0.35 MB
0353 ============================================================
