Warning, /geant4/examples/extended/electromagnetic/TestEm5/TestEm5-pixe_ANSTO.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 /run/verbose 1
0023 #
0024 /testem/det/setAbsMat Gold
0025 /testem/det/setAbsThick 50 um
0026 /testem/det/setAbsYZ 50 um
0027 #
0028 /testem/phys/addPhysics emstandard_opt4
0029 PhysicsList::AddPhysicsList: <emstandard_opt4>
0030 ###/testem/phys/addPhysics emlivermore
0031 #
0032 /run/setCut 2 um
0033 /run/setCutForAGivenParticle gamma 0.5 um
0034 #
0035 # to (de)activate pixe or atomic deexcitation
0036 /process/em/pixe true
0037 /run/physicsModified
0038 /process/em/augerCascade true
0039 /run/physicsModified
0040 /process/em/deexcitationIgnoreCut true
0041 /run/physicsModified
0042 #
0043 # to activate the ANSTO cross sections
0044 /process/em/pixeXSmodel ECPSSR_ANSTO
0045 /run/physicsModified
0046 #
0047 # to activate the ANSTO fluorescence
0048 # radiation yields
0049 /process/em/fluoANSTO true
0050 /run/physicsModified
0051 #
0052 /process/em/verbose 0
0053 /process/em/printParameters
0054 =======================================================================
0055 ====== Electromagnetic Physics Parameters ========
0056 =======================================================================
0057 LPM effect enabled 1
0058 Enable creation and use of sampling tables 0
0059 Apply cuts on all EM processes 0
0060 Use combined TransportationWithMsc Disabled
0061 Use general process 1
0062 Enable linear polarisation for gamma 0
0063 Enable photoeffect sampling below K-shell 1
0064 Enable sampling of quantum entanglement 0
0065 X-section factor for integral approach 0.8
0066 Min kinetic energy for tables 100 eV
0067 Max kinetic energy for tables 100 TeV
0068 Number of bins per decade of a table 20
0069 Verbose level 0
0070 Verbose level for worker thread 0
0071 Bremsstrahlung energy threshold above which
0072 primary e+- is added to the list of secondary 100 TeV
0073 Bremsstrahlung energy threshold above which primary
0074 muon/hadron is added to the list of secondary 100 TeV
0075 Lowest triplet kinetic energy 1 MeV
0076 Enable sampling of gamma linear polarisation 0
0077 5D gamma conversion model type 0
0078 5D gamma conversion model on isolated ion 0
0079 Livermore data directory epics_2017
0080 =======================================================================
0081 ====== Ionisation Parameters ========
0082 =======================================================================
0083 Step function for e+- (0.2, 0.01 mm)
0084 Step function for muons/hadrons (0.1, 0.05 mm)
0085 Step function for light ions (0.1, 0.02 mm)
0086 Step function for general ions (0.1, 0.001 mm)
0087 Lowest e+e- kinetic energy 100 eV
0088 Lowest muon/hadron kinetic energy 1 keV
0089 Use ICRU90 data 1
0090 Fluctuations of dE/dx are enabled 1
0091 Type of fluctuation model for leptons and hadrons Urban
0092 Use built-in Birks satuaration 0
0093 Build CSDA range enabled 0
0094 Use cut as a final range enabled 0
0095 Enable angular generator interface 1
0096 Max kinetic energy for CSDA tables 1 GeV
0097 Max kinetic energy for NIEL computation 1 MeV
0098 Linear loss limit 0.01
0099 Read data from file for e+e- pair production by mu 0
0100 =======================================================================
0101 ====== Multiple Scattering Parameters ========
0102 =======================================================================
0103 Type of msc step limit algorithm for e+- 2
0104 Type of msc step limit algorithm for muons/hadrons 0
0105 Msc lateral displacement for e+- enabled 1
0106 Msc lateral displacement for muons and hadrons 1
0107 Urban msc model lateral displacement alg96 1
0108 Range factor for msc step limit for e+- 0.08
0109 Range factor for msc step limit for muons/hadrons 0.2
0110 Geometry factor for msc step limitation of e+- 2.5
0111 Safety factor for msc step limit for e+- 0.6
0112 Skin parameter for msc step limitation of e+- 3
0113 Lambda limit for msc step limit for e+- 1 mm
0114 Use Mott correction for e- scattering 1
0115 Factor used for dynamic computation of angular
0116 limit between single and multiple scattering 1
0117 Fixed angular limit between single
0118 and multiple scattering 3.1416 rad
0119 Upper energy limit for e+- multiple scattering 100 MeV
0120 Type of electron single scattering model 0
0121 Type of nuclear form-factor 1
0122 Screening factor 1
0123 =======================================================================
0124 ====== Atomic Deexcitation Parameters ========
0125 =======================================================================
0126 Fluorescence enabled 1
0127 Directory in G4LEDATA for fluorescence data files fluor_ANSTO
0128 Auger electron cascade enabled 1
0129 PIXE atomic de-excitation enabled 1
0130 De-excitation module ignores cuts 1
0131 Type of PIXE cross section for hadrons ECPSSR_ANSTO
0132 Type of PIXE cross section for e+- Livermore
0133 =======================================================================
0134 #
0135 /run/initialize
0136
0137 Material: G4_Galactic density: 0.000 mg/cm3 RadL: 204310098.490 pc Nucl.Int.Length: 113427284.261 pc
0138 Imean: 21.800 eV temperature: 2.73 K pressure: 0.00 atm
0139
0140 ---> Element: H (H) Z = 1.0 N = 1 A = 1.008 g/mole
0141 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 %
0142 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 %
0143 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0144
0145
0146 Material: Gold density: 19.320 g/cm3 RadL: 3.344 mm Nucl.Int.Length: 10.540 cm
0147 Imean: 790.000 eV temperature: 293.15 K pressure: 1.00 atm
0148
0149 ---> Element: Au (Au) Z = 79.0 N = 197 A = 196.967 g/mole
0150 ---> Isotope: Au197 Z = 79 N = 197 A = 196.97 g/mole abundance: 100.000 %
0151 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0152
0153
0154 The WORLD is made of 60 um of G4_Galactic. The transverse size (YZ) of the world is 60 um
0155 The ABSORBER is made of 50 um of Gold. The transverse size (YZ) is 50 um
0156 X position of the middle of the absorber 0 fm
0157 #
0158 /testem/gun/setDefault
0159 /gun/particle proton
0160 /gun/energy 3 MeV
0161 #
0162 # to be run interactively
0163 ###/control/execute vis.mac
0164 ###/tracking/verbose 2
0165 ###/run/beamOn 1
0166 #
0167 /analysis/setFileName pixe
0168 /analysis/h1/set 3 1000 0 100 keV
0169 /analysis/h1/set 51 1000 0 100 keV
0170 /analysis/h1/set 55 1000 0 100 keV
0171 #
0172 /analysis/h1/set 2 500 0 50 keV
0173 /analysis/h1/set 50 200 0 20 keV
0174 /analysis/h1/set 54 200 0 20 keV
0175 #
0176 /analysis/h1/set 4 500 0.01 50 keV log10
0177 /analysis/h1/set 52 200 0.01 20 keV log10
0178 /analysis/h1/set 56 200 0.01 20 keV log10
0179 #
0180 /tracking/verbose 0
0181 /run/printProgress 10000
0182 /run/beamOn 100000
0183
0184 ### === G4UAtomicDeexcitation::InitialiseForNewRun()
0185 Using ANSTO K Cross Sections!
0186 Using ANSTO L Cross Sections!
0187 Using ANSTO M Cross Sections!
0188 ======================================================================
0189 The G4PenelopeIonisationModel is being used with the PIXE flag ON.
0190 Atomic de-excitation will be produced statistically by the PIXE
0191 interface by using the shell cross section --> Livermore
0192 The built-in model procedure for atomic de-excitation is disabled.
0193 *Please be sure this is intended*, or disable PIXE by
0194 /process/em/pixe false
0195 ======================================================================
0196
0197 ========= Table of registered couples ============================
0198
0199 Index : 0 used in the geometry : Yes
0200 Material : G4_Galactic
0201 Range cuts : gamma 500 nm e- 2 um e+ 2 um proton 2 um
0202 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 200 eV
0203 Region(s) which use this couple :
0204 DefaultRegionForTheWorld
0205
0206 Index : 1 used in the geometry : Yes
0207 Material : Gold
0208 Range cuts : gamma 500 nm e- 2 um e+ 2 um proton 2 um
0209 Energy thresholds : gamma 1.07128 keV e- 30.6738 keV e+ 30.1802 keV proton 200 eV
0210 Region(s) which use this couple :
0211 DefaultRegionForTheWorld
0212
0213 ==================================================================
0214
0215 ### Run 0 starts.
0216
0217 ------- MixMaxRng engine status -------
0218 Current state vector is:
0219 mixmax state, file version 1.0
0220 N=17 V[N]={906770732717044781, 629165745432651234, 1235682547346241386, 68420008233404568, 2262190991329242458, 2266470399991071809, 1976726662926872232, 245458862506414172, 1955974201201518530, 2155248512522080758, 604170912935414061, 1116171330120743511, 1861018313684488333, 1296715403254578286, 1549011045957234151, 370819759640195970, 2230139271784837643} counter= 17sumtot= 1977567618660788324
0221 ---------------------------------------
0222 ... set ntuple merging row mode : row-wise - done
0223 ... create file : pixe.root - done
0224 ... open analysis file : pixe.root - done
0225 ... open analysis file : pixe.root - done
0226 --> Event 0 starts.
0227 --> Event 10000 starts.
0228 --> Event 20000 starts.
0229 --> Event 30000 starts.
0230 --> Event 40000 starts.
0231 --> Event 50000 starts.
0232 --> Event 60000 starts.
0233 --> Event 70000 starts.
0234 --> Event 80000 starts.
0235 --> Event 90000 starts.
0236 Run terminated.
0237 Run Summary
0238 Number of events processed : 100000
0239 User=73.480000s Real=93.646424s Sys=0.160000s
0240
0241 ======================== run summary ======================
0242
0243 The run was 100000 proton of 3 MeV through 50 um of Gold (density: 19.3 g/cm3 )
0244
0245 Total energy deposit in absorber per event = 3 MeV +- 36.45 eV
0246
0247 -----> Mean dE/dx = 600 MeV/cm (31.05 MeV*cm2/g)
0248
0249 From formulas :
0250 restricted dEdx = 676.4 MeV/cm (35.01 MeV*cm2/g)
0251 full dEdx = 676.5 MeV/cm (35.01 MeV*cm2/g)
0252
0253 Leakage : primary = 65.74 eV +- 35.92 eV secondaries = 14.07 eV +- 1.179 eV
0254 Energy balance : edep + eleak = 3 MeV
0255
0256 Total track length (charged) in absorber per event = 31.06 um +- 1.083 nm
0257 Total track length (neutral) in absorber per event = 362.7 nm +- 3.184 nm
0258
0259 Number of steps (charged) in absorber per event = 141.8 +- 0.1263
0260 Number of steps (neutral) in absorber per event = 0.5861 +- 0.002446
0261
0262 Number of secondaries per event : Gammas = 0.5839 electrons = 109.6 positrons = 0
0263
0264 Number of events with the primary particle transmitted = 0 %
0265 Number of events with at least 1 particle transmitted (same charge as primary) = 0 %
0266
0267 Number of events with the primary particle reflected = 0.005 %
0268 Number of events with at least 1 particle reflected (same charge as primary) = 0.005 %
0269
0270 MultipleScattering:
0271 rms proj angle of transmit primary particle = 0 mrad (central part only)
0272 computed theta0 (Highland formula) = 233.3 mrad
0273 central part defined as +- 699.8 mrad Tail ratio = 0 %
0274
0275 Gamma process counts:
0276 Photoeffect 58109
0277 Compton 3
0278 Conversion 0
0279 Rayleigh 218
0280 ... write file : pixe.root - done
0281 ... close file : pixe.root - done
0282
0283 ------- MixMaxRng engine status -------
0284 Current state vector is:
0285 mixmax state, file version 1.0
0286 N=17 V[N]={944671353633884230, 300194246912256983, 1338926093308238266, 1441599044978111444, 559828476273747537, 976665045002376518, 1655586542061687849, 1507013437582523992, 2052944524787068341, 2173254900581545545, 733849232472135636, 379249519059016924, 877093410379007237, 1488918075031648810, 1263787258307413205, 589725557209658928, 749180548014423557} counter= 14sumtot= 585743191885193394
0287 ---------------------------------------
0288 ================== Deleting memory pools ===================
0289 Number of memory pools allocated: 9 of which, static: 0
0290 Dynamic pools deleted: 9 / Total memory freed: 0.16 MB
0291 ============================================================
