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-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 = 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 Positron annihilation at rest model SimplePositronium
0158 Enable 3 gamma annihilation on fly 0
0159 Lowest triplet kinetic energy 1 MeV
0160 Enable sampling of gamma linear polarisation 0
0161 5D gamma conversion model type 0
0162 5D gamma conversion model on isolated ion 0
0163 Use RiGe 5D e+e- pair production model by muons 0
0164 Livermore data directory epics_2017
0165 =======================================================================
0166 ====== Ionisation Parameters ========
0167 =======================================================================
0168 Step function for e+- (0.2, 0.01 mm)
0169 Step function for muons/hadrons (0.2, 0.1 mm)
0170 Step function for light ions (0.2, 0.1 mm)
0171 Step function for general ions (0.2, 0.1 mm)
0172 Lowest e+e- kinetic energy 100 eV
0173 Lowest muon/hadron kinetic energy 1 keV
0174 Use ICRU90 data 0
0175 Fluctuations of dE/dx are enabled 1
0176 Type of fluctuation model for leptons and hadrons Urban
0177 Use built-in Birks satuaration 0
0178 Build CSDA range enabled 1
0179 Use cut as a final range enabled 0
0180 Enable angular generator interface 0
0181 Max kinetic energy for CSDA tables 10 GeV
0182 Max kinetic energy for NIEL computation 0 eV
0183 Linear loss limit 0.01
0184 Read data from file for e+e- pair production by mu 0
0185 =======================================================================
0186 ====== Multiple Scattering Parameters ========
0187 =======================================================================
0188 Type of msc step limit algorithm for e+- 3
0189 Type of msc step limit algorithm for muons/hadrons 0
0190 Msc lateral displacement for e+- enabled 1
0191 Msc lateral displacement for muons and hadrons 0
0192 Urban msc model lateral displacement alg96 1
0193 Range factor for msc step limit for e+- 0.04
0194 Range factor for msc step limit for muons/hadrons 0.2
0195 Geometry factor for msc step limitation of e+- 2.5
0196 Safety factor for msc step limit for e+- 0.6
0197 Skin parameter for msc step limitation of e+- 1
0198 Lambda limit for msc step limit for e+- 1 mm
0199 Use Mott correction for e- scattering 0
0200 Factor used for dynamic computation of angular
0201 limit between single and multiple scattering 1
0202 Fixed angular limit between single
0203 and multiple scattering 3.1416 rad
0204 Upper energy limit for e+- multiple scattering 100 MeV
0205 Type of electron single scattering model 0
0206 Type of nuclear form-factor 1
0207 Screening factor 1
0208 =======================================================================
0209
0210 phot: for gamma SubType=12 BuildTable=0
0211 LambdaPrime table from 200 keV to 10 GeV in 94 bins
0212 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0213 PhotoElectric : Emin= 0 eV Emax= 10 GeV SauterGavrila
0214
0215 compt: for gamma SubType=13 BuildTable=1
0216 Lambda table from 10 eV to 1 MeV, 20 bins/decade, spline: 1
0217 LambdaPrime table from 1 MeV to 10 GeV in 80 bins
0218 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0219 myKlein-Nishina : Emin= 0 eV Emax= 10 GeV
0220
0221 conv: for gamma SubType=14 BuildTable=1
0222 Lambda table from 1.022 MeV to 10 GeV, 25 bins/decade, spline: 1
0223 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0224 BetheHeitlerLPM : Emin= 0 eV Emax= 10 GeV ModifiedTsai
0225
0226 msc: for e- SubType= 10
0227 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0228 UrbanMsc : Emin= 0 eV Emax= 10 GeV Nbins=160 100 eV - 10 GeV
0229 StepLim=DistanceToBoundary Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0230
0231 eIoni: for e- XStype:1 SubType=2
0232 dE/dx and range tables from 10 eV to 10 GeV in 180 bins
0233 Lambda tables from threshold to 10 GeV, 20 bins/decade, spline: 1
0234 StepFunction=(0.2, 0.01 mm), integ: 1, fluct: 1, linLossLim= 0.01
0235 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0236 myMollerBhabha : Emin= 0 eV Emax= 10 GeV
0237 CSDA range table up to 10 GeV in 180 bins
0238
0239 msc: for e+ SubType= 10
0240 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0241 UrbanMsc : Emin= 0 eV Emax= 10 GeV Nbins=160 100 eV - 10 GeV
0242 StepLim=DistanceToBoundary Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0243
0244 eIoni: for e+ XStype:1 SubType=2
0245 dE/dx and range tables from 10 eV to 10 GeV in 180 bins
0246 Lambda tables from threshold to 10 GeV, 20 bins/decade, spline: 1
0247 StepFunction=(0.2, 0.01 mm), integ: 1, fluct: 1, linLossLim= 0.01
0248 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0249 myMollerBhabha : Emin= 0 eV Emax= 10 GeV
0250 CSDA range table up to 10 GeV in 180 bins
0251
0252 annihil: for e+ XStype:2 SubType=5 AtRestModel:Simple BuildTable=0
0253 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0254 eplus2gg : Emin= 0 eV Emax= 10 GeV
0255
0256 msc: for proton SubType= 10
0257 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0258 UrbanMsc : Emin= 0 eV Emax= 10 GeV Nbins=160 100 eV - 10 GeV
0259 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0260
0261 hIoni: for proton XStype:1 SubType=2
0262 dE/dx and range tables from 10 eV to 10 GeV in 180 bins
0263 Lambda tables from threshold to 10 GeV, 20 bins/decade, spline: 1
0264 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01
0265 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0266 Bragg : Emin= 0 eV Emax= 2 MeV
0267 BetheBloch : Emin= 2 MeV Emax= 10 GeV
0268 CSDA range table up to 10 GeV in 180 bins
0269
0270 ========= Table of registered couples ============================
0271
0272 Index : 0 used in the geometry : Yes
0273 Material : G4_WATER
0274 Range cuts : gamma 10 km e- 10 km e+ 10 km proton 10 km
0275 Energy thresholds : gamma 9.54965 GeV e- 9.54965 GeV e+ 9.54965 GeV proton 1 TeV
0276 Region(s) which use this couple :
0277 DefaultRegionForTheWorld
0278
0279 Index : 1 used in the geometry : Yes
0280 Material : g4Water_gas
0281 Range cuts : gamma 10 km e- 10 km e+ 10 km proton 10 km
0282 Energy thresholds : gamma 9.54965 GeV e- 9.54965 GeV e+ 9.54965 GeV proton 1 TeV
0283 Region(s) which use this couple :
0284 DefaultRegionForTheWorld
0285
0286 ==================================================================
0287
0288 ### Run 0 starts.
0289 ### Run 0 start.
0290 --> Event 0 starts.
0291 --> Event 1000 starts.
0292 ---> NbofEvents= 1000 NbOfelectr= 70472 Tkin= 587.833 keV (0.000 %) NbOfelec in cav= 4882 Dose/EnFluence= 29.487 cm2/g (0.000 %)
0293
0294 --> Event 2000 starts.
0295 ---> NbofEvents= 2000 NbOfelectr= 140752 Tkin= 587.532 keV (-0.051 %) NbOfelec in cav= 9647 Dose/EnFluence= 29.645 cm2/g (0.538 %)
0296
0297 --> Event 3000 starts.
0298 ---> NbofEvents= 3000 NbOfelectr= 211683 Tkin= 587.487 keV (-0.008 %) NbOfelec in cav= 14584 Dose/EnFluence= 28.703 cm2/g (-3.180 %)
0299
0300 Run terminated.
0301 Run Summary
0302 Number of events processed : 4000
0303 User=4.920000s Real=4.916582s Sys=0.000000s
0304
0305 ======================== run summary ======================
0306
0307 The run consists of 4000 gamma of 1.250 MeV through 2*5.000 mm of G4_WATER (density: 1.000 g/cm3 )
0308
0309 the cavity is 2.000 mm of g4Water_gas (density: 1.000 mg/cm3) Mass = 0.628 mg
0310
0311 ============================================================
0312
0313 Process calls frequency ---> compt= 282020 Transportation= 78385 msc= 1049116 eIoni= 2098867
0314
0315 Gamma crossSections in wall material : compt= 63.464 cm2/g --> total= 63.464 cm2/g
0316
0317 Mean energy of secondary e- = 587.9031 keV +- 0.1030 % (--> range in wall material = 2.2194 mm )
0318 Mass_energy_transfer coef: 29.8484 cm2/g
0319
0320 StoppingPower in wall = 1.9646 MeV*cm2/g
0321 in cavity = 1.9646 MeV*cm2/g
0322
0323 Charged particle flow in cavity :
0324 Enter --> nbParticles = 19379 Energy = 9.8581 GeV
0325 Exit --> nbParticles = 19382 Energy = 9.8562 GeV
0326
0327 Total edep in cavity = 12.9357 MeV +- 3.2874 % Total charged trackLength = 53.5865 m (mean value = 2.7652 mm )
0328 Total dose in cavity = 20.5879 MeV/mg
0329 Dose/EnergyFluence = 29.1054 cm2/g
0330
0331 (Dose/EnergyFluence)/Mass_energy_transfer = 0.97511 +- 0.03207
0332
0333 StepSize of ch. tracks in wall = 148.8790 um +- 150.7159 um (nbSteps/track = 11.2614)
0334 StepSize of ch. tracks in cavity = 1.2519 mm +- 1.6595 mm (nbSteps/track = 2.2087)
0335 ================== Deleting memory pools ===================
0336 Number of memory pools allocated: 9 of which, static: 0
0337 Dynamic pools deleted: 9 / Total memory freed: 0.074 MB
0338 ============================================================
