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-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 /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 Positron annihilation at rest model AllisonPositronium
0076 Enable 3 gamma annihilation on fly 1
0077 Lowest triplet kinetic energy 1 MeV
0078 Enable sampling of gamma linear polarisation 0
0079 5D gamma conversion model type 0
0080 5D gamma conversion model on isolated ion 0
0081 Use RiGe 5D e+e- pair production model by muons 0
0082 Livermore data directory epics_2017
0083 =======================================================================
0084 ====== Ionisation Parameters ========
0085 =======================================================================
0086 Step function for e+- (0.2, 0.01 mm)
0087 Step function for muons/hadrons (0.1, 0.05 mm)
0088 Step function for light ions (0.1, 0.02 mm)
0089 Step function for general ions (0.1, 0.001 mm)
0090 Lowest e+e- kinetic energy 100 eV
0091 Lowest muon/hadron kinetic energy 1 keV
0092 Use ICRU90 data 1
0093 Fluctuations of dE/dx are enabled 1
0094 Type of fluctuation model for leptons and hadrons Urban
0095 Use built-in Birks satuaration 0
0096 Build CSDA range enabled 0
0097 Use cut as a final range enabled 0
0098 Enable angular generator interface 1
0099 Max kinetic energy for CSDA tables 1 GeV
0100 Max kinetic energy for NIEL computation 1 MeV
0101 Linear loss limit 0.01
0102 Read data from file for e+e- pair production by mu 0
0103 =======================================================================
0104 ====== Multiple Scattering Parameters ========
0105 =======================================================================
0106 Type of msc step limit algorithm for e+- 2
0107 Type of msc step limit algorithm for muons/hadrons 0
0108 Msc lateral displacement for e+- enabled 1
0109 Msc lateral displacement for muons and hadrons 1
0110 Urban msc model lateral displacement alg96 1
0111 Range factor for msc step limit for e+- 0.08
0112 Range factor for msc step limit for muons/hadrons 0.2
0113 Geometry factor for msc step limitation of e+- 2.5
0114 Safety factor for msc step limit for e+- 0.6
0115 Skin parameter for msc step limitation of e+- 3
0116 Lambda limit for msc step limit for e+- 1 mm
0117 Use Mott correction for e- scattering 1
0118 Factor used for dynamic computation of angular
0119 limit between single and multiple scattering 1
0120 Fixed angular limit between single
0121 and multiple scattering 3.1416 rad
0122 Upper energy limit for e+- multiple scattering 100 MeV
0123 Type of electron single scattering model 0
0124 Type of nuclear form-factor 1
0125 Screening factor 1
0126 =======================================================================
0127 ====== Atomic Deexcitation Parameters ========
0128 =======================================================================
0129 Fluorescence enabled 1
0130 Directory in G4LEDATA for fluorescence data files fluor_ANSTO
0131 Auger electron cascade enabled 1
0132 PIXE atomic de-excitation enabled 1
0133 De-excitation module ignores cuts 1
0134 Type of PIXE cross section for hadrons ECPSSR_ANSTO
0135 Type of PIXE cross section for e+- Livermore
0136 =======================================================================
0137 #
0138 /run/initialize
0139
0140 Material: G4_Galactic density: 0.000 mg/cm3 RadL: 204310098.490 pc Nucl.Int.Length: 3240.779 pc
0141 Imean: 21.800 eV temperature: 2.73 K pressure: 0.00 atm
0142
0143 ---> Element: H (H) Z = 1.0 N = 1 A = 1.008 g/mole
0144 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 %
0145 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 %
0146 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0147
0148
0149 Material: Gold density: 19.320 g/cm3 RadL: 3.344 mm Nucl.Int.Length: 10.540 cm
0150 Imean: 790.000 eV temperature: 293.15 K pressure: 1.00 atm
0151
0152 ---> Element: Au (Au) Z = 79.0 N = 197 A = 196.967 g/mole
0153 ---> Isotope: Au197 Z = 79 N = 197 A = 196.97 g/mole abundance: 100.000 %
0154 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0155
0156
0157 The WORLD is made of 60 um of G4_Galactic. The transverse size (YZ) of the world is 60 um
0158 The ABSORBER is made of 50 um of Gold. The transverse size (YZ) is 50 um
0159 X position of the middle of the absorber 0 fm
0160 #
0161 /testem/gun/setDefault
0162 /gun/particle proton
0163 /gun/energy 3 MeV
0164 #
0165 # to be run interactively
0166 ###/control/execute vis.mac
0167 ###/tracking/verbose 2
0168 ###/run/beamOn 1
0169 #
0170 /analysis/setFileName pixe
0171 /analysis/h1/set 3 1000 0 100 keV
0172 /analysis/h1/set 51 1000 0 100 keV
0173 /analysis/h1/set 55 1000 0 100 keV
0174 #
0175 /analysis/h1/set 2 500 0 50 keV
0176 /analysis/h1/set 50 200 0 20 keV
0177 /analysis/h1/set 54 200 0 20 keV
0178 #
0179 /analysis/h1/set 4 500 0.01 50 keV log10
0180 /analysis/h1/set 52 200 0.01 20 keV log10
0181 /analysis/h1/set 56 200 0.01 20 keV log10
0182 #
0183 /tracking/verbose 0
0184 /run/printProgress 10000
0185 /run/beamOn 100000
0186
0187 ### === G4UAtomicDeexcitation::InitialiseForNewRun()
0188 Using ANSTO K Cross Sections!
0189 Using ANSTO L Cross Sections!
0190 Using ANSTO M Cross Sections!
0191 ======================================================================
0192 The G4PenelopeIonisationModel is being used with the PIXE flag ON.
0193 Atomic de-excitation will be produced statistically by the PIXE
0194 interface by using the shell cross section --> Livermore
0195 The built-in model procedure for atomic de-excitation is disabled.
0196 *Please be sure this is intended*, or disable PIXE by
0197 /process/em/pixe false
0198 ======================================================================
0199
0200 ========= Table of registered couples ============================
0201
0202 Index : 0 used in the geometry : Yes
0203 Material : G4_Galactic
0204 Range cuts : gamma 500 nm e- 2 um e+ 2 um proton 2 um
0205 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 200 eV
0206 Region(s) which use this couple :
0207 DefaultRegionForTheWorld
0208
0209 Index : 1 used in the geometry : Yes
0210 Material : Gold
0211 Range cuts : gamma 500 nm e- 2 um e+ 2 um proton 2 um
0212 Energy thresholds : gamma 1.07128 keV e- 30.6738 keV e+ 30.1802 keV proton 200 eV
0213 Region(s) which use this couple :
0214 DefaultRegionForTheWorld
0215
0216 ==================================================================
0217
0218 ### Run 0 starts.
0219
0220 ------- MixMaxRng engine status -------
0221 Current state vector is:
0222 mixmax state, file version 1.0
0223 N=17 V[N]={906770732717044781, 629165745432651234, 1235682547346241386, 68420008233404568, 2262190991329242458, 2266470399991071809, 1976726662926872232, 245458862506414172, 1955974201201518530, 2155248512522080758, 604170912935414061, 1116171330120743511, 1861018313684488333, 1296715403254578286, 1549011045957234151, 370819759640195970, 2230139271784837643} counter= 17sumtot= 1977567618660788324
0224 ---------------------------------------
0225 ... set ntuple merging row mode : row-wise - done
0226 ... create file : pixe.root - done
0227 ... open analysis file : pixe.root - done
0228 ... open analysis file : pixe.root - done
0229 --> Event 0 starts.
0230 --> Event 10000 starts.
0231 --> Event 20000 starts.
0232 --> Event 30000 starts.
0233 --> Event 40000 starts.
0234 --> Event 50000 starts.
0235 --> Event 60000 starts.
0236 --> Event 70000 starts.
0237 --> Event 80000 starts.
0238 --> Event 90000 starts.
0239 Run terminated.
0240 Run Summary
0241 Number of events processed : 100000
0242 User=35.380000s Real=35.465336s Sys=0.000000s
0243
0244 ======================== run summary ======================
0245
0246 The run was 100000 proton of 3 MeV through 50 um of Gold (density: 19.3 g/cm3 )
0247
0248 Total energy deposit in absorber per event = 3 MeV +- 76.71 eV
0249
0250 -----> Mean dE/dx = 599.9 MeV/cm (31.05 MeV*cm2/g)
0251
0252 From formulas :
0253 restricted dEdx = 676.4 MeV/cm (35.01 MeV*cm2/g)
0254 full dEdx = 676.5 MeV/cm (35.01 MeV*cm2/g)
0255
0256 Leakage : primary = 278.1 eV +- 76.3 eV secondaries = 14.72 eV +- 955.4 meV
0257 Energy balance : edep + eleak = 3 MeV
0258
0259 Total track length (charged) in absorber per event = 31.06 um +- 1.249 nm
0260 Total track length (neutral) in absorber per event = 355.2 nm +- 3.118 nm
0261
0262 Number of steps (charged) in absorber per event = 141.8 +- 0.1266
0263 Number of steps (neutral) in absorber per event = 0.5817 +- 0.002434
0264
0265 Number of secondaries per event : Gammas = 0.5798 electrons = 109.6 positrons = 0
0266
0267 Number of events with the primary particle transmitted = 0 %
0268 Number of events with at least 1 particle transmitted (same charge as primary) = 0 %
0269
0270 Number of events with the primary particle reflected = 0.016 %
0271 Number of events with at least 1 particle reflected (same charge as primary) = 0.016 %
0272
0273 MultipleScattering:
0274 rms proj angle of transmit primary particle = 0 mrad (central part only)
0275 computed theta0 (Highland formula) = 233.3 mrad
0276 central part defined as +- 699.8 mrad Tail ratio = 0 %
0277
0278 Gamma process counts:
0279 Photoeffect 57671
0280 Compton 2
0281 Conversion 0
0282 Rayleigh 195
0283 ... write file : pixe.root - done
0284 ... close file : pixe.root - done
0285
0286 ------- MixMaxRng engine status -------
0287 Current state vector is:
0288 mixmax state, file version 1.0
0289 N=17 V[N]={1865847532392181060, 99825406610688157, 359444420570531604, 1595240065015446652, 893119202495921306, 1324979107203650855, 852295147626347754, 826571437727771892, 2049915322115205500, 757520833902555854, 348265068551329111, 748267151128822368, 1950628009764165982, 1535417471323711813, 1798702300717679741, 1338321521261254999, 2232042557627417187} counter= 15sumtot= 2129658482325130227
0290 ---------------------------------------
0291 ================== Deleting memory pools ===================
0292 Number of memory pools allocated: 9 of which, static: 0
0293 Dynamic pools deleted: 9 / Total memory freed: 0.15 MB
0294 ============================================================
