Warning, /geant4/examples/extended/medical/fanoCavity/fanoCavity.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 = 8 *****
0024
0025 Material: Water density: 1.000 g/cm3 RadL: 36.092 cm Nucl.Int.Length: 75.356 cm
0026 Imean: 78.000 eV temperature: 293.15 K pressure: 1.00 atm
0027
0028 ---> Element: Hydrogen (H) Z = 1.0 N = 1 A = 1.010 g/mole
0029 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 %
0030 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 %
0031 ElmMassFraction: 11.21 % ElmAbundance 66.67 %
0032
0033 ---> Element: Oxygen (O) Z = 8.0 N = 16 A = 16.000 g/mole
0034 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 %
0035 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 %
0036 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 %
0037 ElmMassFraction: 88.79 % ElmAbundance 33.33 %
0038
0039
0040 Material: Water_gas density: 1.000 mg/cm3 RadL: 360.924 m Nucl.Int.Length: 753.560 m
0041 Imean: 78.000 eV temperature: 293.15 K pressure: 1.00 atm
0042
0043 ---> Element: Hydrogen (H) Z = 1.0 N = 1 A = 1.010 g/mole
0044 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 %
0045 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 %
0046 ElmMassFraction: 11.21 % ElmAbundance 66.67 %
0047
0048 ---> Element: Oxygen (O) Z = 8.0 N = 16 A = 16.000 g/mole
0049 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 %
0050 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 %
0051 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 %
0052 ElmMassFraction: 88.79 % ElmAbundance 33.33 %
0053
0054
0055 Material: Graphite density: 2.265 g/cm3 RadL: 18.851 cm Nucl.Int.Length: 35.388 cm
0056 Imean: 81.000 eV temperature: 293.15 K pressure: 1.00 atm
0057
0058 ---> Element: C (C) Z = 6.0 N = 12 A = 12.011 g/mole
0059 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.930 %
0060 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.070 %
0061 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0062
0063
0064 Material: Graphite_gas density: 2.265 mg/cm3 RadL: 188.508 m Nucl.Int.Length: 353.880 m
0065 Imean: 81.000 eV temperature: 293.15 K pressure: 1.00 atm
0066
0067 ---> Element: C (C) Z = 6.0 N = 12 A = 12.011 g/mole
0068 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.930 %
0069 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.070 %
0070 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0071
0072
0073 Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Nucl.Int.Length: 38.880 cm
0074 Imean: 166.000 eV temperature: 293.15 K pressure: 1.00 atm
0075
0076 ---> Element: Al (Al) Z = 13.0 N = 27 A = 26.982 g/mole
0077 ---> Isotope: Al27 Z = 13 N = 27 A = 26.98 g/mole abundance: 100.000 %
0078 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0079
0080
0081 Material: Aluminium_gas density: 2.700 mg/cm3 RadL: 88.930 m Nucl.Int.Length: 388.800 m
0082 Imean: 166.000 eV temperature: 293.15 K pressure: 1.00 atm
0083
0084 ---> Element: Al (Al) Z = 13.0 N = 27 A = 26.982 g/mole
0085 ---> Isotope: Al27 Z = 13 N = 27 A = 26.98 g/mole abundance: 100.000 %
0086 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0087
0088
0089 Material: G4_WATER H_2O density: 1.000 g/cm3 RadL: 36.083 cm Nucl.Int.Length: 75.375 cm
0090 Imean: 78.000 eV temperature: 293.15 K pressure: 1.00 atm
0091
0092 ---> Element: H (H) Z = 1.0 N = 1 A = 1.008 g/mole
0093 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 %
0094 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 %
0095 ElmMassFraction: 11.19 % ElmAbundance 66.67 %
0096
0097 ---> Element: O (O) Z = 8.0 N = 16 A = 15.999 g/mole
0098 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 %
0099 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 %
0100 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 %
0101 ElmMassFraction: 88.81 % ElmAbundance 33.33 %
0102
0103
0104 Material: g4Water_gas H_2O density: 1.000 mg/cm3 RadL: 360.830 m Nucl.Int.Length: 753.748 m
0105 Imean: 78.000 eV temperature: 293.15 K pressure: 1.00 atm
0106
0107 ---> Element: H (H) Z = 1.0 N = 1 A = 1.008 g/mole
0108 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 %
0109 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 %
0110 ElmMassFraction: 11.19 % ElmAbundance 66.67 %
0111
0112 ---> Element: O (O) Z = 8.0 N = 16 A = 15.999 g/mole
0113 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 %
0114 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 %
0115 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 %
0116 ElmMassFraction: 88.81 % ElmAbundance 33.33 %
0117
0118
0119
0120 /control/cout/ignoreThreadsExcept 0
0121 /run/verbose 1
0122 /testem/phys/addPhysics standard_opt3
0123 PhysicsList::AddPhysicsList: <standard_opt3>
0124 /process/eLoss/StepFunction 0.2 10 um
0125 /run/physicsModified
0126 /run/initialize
0127
0128 ---------------------------------------------------------
0129 ---> The Wall is 5 mm of G4_WATER ( 1 g/cm3 )
0130 The Cavity is 2 mm of g4Water_gas ( 1 mg/cm3 )
0131 ---------------------------------------------------------
0132
0133 PhysicsList::SetCuts:CutLength : 10 km
0134 /run/printProgress 1000
0135 /run/beamOn 4000
0136 =======================================================================
0137 ====== Electromagnetic Physics Parameters ========
0138 =======================================================================
0139 LPM effect enabled 1
0140 Enable creation and use of sampling tables 0
0141 Apply cuts on all EM processes 0
0142 Use combined TransportationWithMsc Disabled
0143 Use general process 0
0144 Enable linear polarisation for gamma 0
0145 Enable photoeffect sampling below K-shell 1
0146 Enable sampling of quantum entanglement 0
0147 X-section factor for integral approach 0.8
0148 Min kinetic energy for tables 10 eV
0149 Max kinetic energy for tables 10 GeV
0150 Number of bins per decade of a table 20
0151 Verbose level 1
0152 Verbose level for worker thread 0
0153 Bremsstrahlung energy threshold above which
0154 primary e+- is added to the list of secondary 100 TeV
0155 Bremsstrahlung energy threshold above which primary
0156 muon/hadron is added to the list of secondary 100 TeV
0157 Lowest triplet kinetic energy 1 MeV
0158 Enable sampling of gamma linear polarisation 0
0159 5D gamma conversion model type 0
0160 5D gamma conversion model on isolated ion 0
0161 Livermore data directory epics_2017
0162 =======================================================================
0163 ====== Ionisation Parameters ========
0164 =======================================================================
0165 Step function for e+- (0.2, 0.01 mm)
0166 Step function for muons/hadrons (0.2, 0.1 mm)
0167 Step function for light ions (0.2, 0.1 mm)
0168 Step function for general ions (0.2, 0.1 mm)
0169 Lowest e+e- kinetic energy 100 eV
0170 Lowest muon/hadron kinetic energy 1 keV
0171 Use ICRU90 data 0
0172 Fluctuations of dE/dx are enabled 1
0173 Type of fluctuation model for leptons and hadrons Urban
0174 Use built-in Birks satuaration 0
0175 Build CSDA range enabled 1
0176 Use cut as a final range enabled 0
0177 Enable angular generator interface 0
0178 Max kinetic energy for CSDA tables 10 GeV
0179 Max kinetic energy for NIEL computation 0 eV
0180 Linear loss limit 0.01
0181 Read data from file for e+e- pair production by mu 0
0182 =======================================================================
0183 ====== Multiple Scattering Parameters ========
0184 =======================================================================
0185 Type of msc step limit algorithm for e+- 3
0186 Type of msc step limit algorithm for muons/hadrons 0
0187 Msc lateral displacement for e+- enabled 1
0188 Msc lateral displacement for muons and hadrons 0
0189 Urban msc model lateral displacement alg96 1
0190 Range factor for msc step limit for e+- 0.04
0191 Range factor for msc step limit for muons/hadrons 0.2
0192 Geometry factor for msc step limitation of e+- 2.5
0193 Safety factor for msc step limit for e+- 0.6
0194 Skin parameter for msc step limitation of e+- 1
0195 Lambda limit for msc step limit for e+- 1 mm
0196 Use Mott correction for e- scattering 0
0197 Factor used for dynamic computation of angular
0198 limit between single and multiple scattering 1
0199 Fixed angular limit between single
0200 and multiple scattering 3.1416 rad
0201 Upper energy limit for e+- multiple scattering 100 MeV
0202 Type of electron single scattering model 0
0203 Type of nuclear form-factor 1
0204 Screening factor 1
0205 =======================================================================
0206
0207 phot: for gamma SubType=12 BuildTable=0
0208 LambdaPrime table from 200 keV to 10 GeV in 94 bins
0209 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0210 PhotoElectric : Emin= 0 eV Emax= 10 GeV SauterGavrila
0211
0212 compt: for gamma SubType=13 BuildTable=1
0213 Lambda table from 10 eV to 1 MeV, 20 bins/decade, spline: 1
0214 LambdaPrime table from 1 MeV to 10 GeV in 80 bins
0215 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0216 myKlein-Nishina : Emin= 0 eV Emax= 10 GeV
0217
0218 conv: for gamma SubType=14 BuildTable=1
0219 Lambda table from 1.022 MeV to 10 GeV, 25 bins/decade, spline: 1
0220 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0221 BetheHeitlerLPM : Emin= 0 eV Emax= 10 GeV ModifiedTsai
0222
0223 msc: for e- SubType= 10
0224 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0225 UrbanMsc : Emin= 0 eV Emax= 10 GeV Nbins=160 100 eV - 10 GeV
0226 StepLim=DistanceToBoundary Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0227
0228 eIoni: for e- XStype:1 SubType=2
0229 dE/dx and range tables from 10 eV to 10 GeV in 180 bins
0230 Lambda tables from threshold to 10 GeV, 20 bins/decade, spline: 1
0231 StepFunction=(0.2, 0.01 mm), integ: 1, fluct: 1, linLossLim= 0.01
0232 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0233 myMollerBhabha : Emin= 0 eV Emax= 10 GeV
0234 CSDA range table up to 10 GeV in 180 bins
0235
0236 msc: for e+ SubType= 10
0237 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0238 UrbanMsc : Emin= 0 eV Emax= 10 GeV Nbins=160 100 eV - 10 GeV
0239 StepLim=DistanceToBoundary Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0240
0241 eIoni: for e+ XStype:1 SubType=2
0242 dE/dx and range tables from 10 eV to 10 GeV in 180 bins
0243 Lambda tables from threshold to 10 GeV, 20 bins/decade, spline: 1
0244 StepFunction=(0.2, 0.01 mm), integ: 1, fluct: 1, linLossLim= 0.01
0245 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0246 myMollerBhabha : Emin= 0 eV Emax= 10 GeV
0247 CSDA range table up to 10 GeV in 180 bins
0248
0249 annihil: for e+ XStype:2 SubType=5 AtRestModel:Simple BuildTable=0
0250 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0251 eplus2gg : Emin= 0 eV Emax= 10 GeV
0252
0253 msc: for proton SubType= 10
0254 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0255 UrbanMsc : Emin= 0 eV Emax= 10 GeV Nbins=160 100 eV - 10 GeV
0256 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0257
0258 hIoni: for proton XStype:1 SubType=2
0259 dE/dx and range tables from 10 eV to 10 GeV in 180 bins
0260 Lambda tables from threshold to 10 GeV, 20 bins/decade, spline: 1
0261 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01
0262 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0263 Bragg : Emin= 0 eV Emax= 2 MeV
0264 BetheBloch : Emin= 2 MeV Emax= 10 GeV
0265 CSDA range table up to 10 GeV in 180 bins
0266
0267 ========= Table of registered couples ============================
0268
0269 Index : 0 used in the geometry : Yes
0270 Material : G4_WATER
0271 Range cuts : gamma 10 km e- 10 km e+ 10 km proton 10 km
0272 Energy thresholds : gamma 9.54965 GeV e- 9.54965 GeV e+ 9.54965 GeV proton 1 TeV
0273 Region(s) which use this couple :
0274 DefaultRegionForTheWorld
0275
0276 Index : 1 used in the geometry : Yes
0277 Material : g4Water_gas
0278 Range cuts : gamma 10 km e- 10 km e+ 10 km proton 10 km
0279 Energy thresholds : gamma 9.54965 GeV e- 9.54965 GeV e+ 9.54965 GeV proton 1 TeV
0280 Region(s) which use this couple :
0281 DefaultRegionForTheWorld
0282
0283 ==================================================================
0284
0285 ### Run 0 starts.
0286 ### Run 0 start.
0287 --> Event 0 starts.
0288 --> Event 1000 starts.
0289 ---> NbofEvents= 1000 NbOfelectr= 70472 Tkin= 587.833 keV (0.000 %) NbOfelec in cav= 4882 Dose/EnFluence= 29.487 cm2/g (0.000 %)
0290
0291 --> Event 2000 starts.
0292 ---> NbofEvents= 2000 NbOfelectr= 140752 Tkin= 587.532 keV (-0.051 %) NbOfelec in cav= 9647 Dose/EnFluence= 29.645 cm2/g (0.538 %)
0293
0294 --> Event 3000 starts.
0295 ---> NbofEvents= 3000 NbOfelectr= 211683 Tkin= 587.487 keV (-0.008 %) NbOfelec in cav= 14584 Dose/EnFluence= 28.703 cm2/g (-3.180 %)
0296
0297 Run terminated.
0298 Run Summary
0299 Number of events processed : 4000
0300 User=8.790000s Real=9.084049s Sys=0.000000s
0301
0302 ======================== run summary ======================
0303
0304 The run consists of 4000 gamma of 1.250 MeV through 2*5.000 mm of G4_WATER (density: 1.000 g/cm3 )
0305
0306 the cavity is 2.000 mm of g4Water_gas (density: 1.000 mg/cm3) Mass = 0.628 mg
0307
0308 ============================================================
0309
0310 Process calls frequency ---> compt= 282020 Transportation= 78385 msc= 1049116 eIoni= 2098867
0311
0312 Gamma crossSections in wall material : compt= 63.464 cm2/g --> total= 63.464 cm2/g
0313
0314 Mean energy of secondary e- = 587.9031 keV +- 0.1030 % (--> range in wall material = 2.2194 mm )
0315 Mass_energy_transfer coef: 29.8484 cm2/g
0316
0317 StoppingPower in wall = 1.9646 MeV*cm2/g
0318 in cavity = 1.9646 MeV*cm2/g
0319
0320 Charged particle flow in cavity :
0321 Enter --> nbParticles = 19379 Energy = 9.8581 GeV
0322 Exit --> nbParticles = 19382 Energy = 9.8562 GeV
0323
0324 Total edep in cavity = 12.9357 MeV +- 3.2874 % Total charged trackLength = 53.5865 m (mean value = 2.7652 mm )
0325 Total dose in cavity = 20.5879 MeV/mg
0326 Dose/EnergyFluence = 29.1054 cm2/g
0327
0328 (Dose/EnergyFluence)/Mass_energy_transfer = 0.97511 +- 0.03207
0329
0330 StepSize of ch. tracks in wall = 148.8790 um +- 150.7159 um (nbSteps/track = 11.2614)
0331 StepSize of ch. tracks in cavity = 1.2519 mm +- 1.6595 mm (nbSteps/track = 2.2087)
0332 ================== Deleting memory pools ===================
0333 Number of memory pools allocated: 9 of which, static: 0
0334 Dynamic pools deleted: 9 / Total memory freed: 0.074 MB
0335 ============================================================
