Warning, /geant4/examples/extended/medical/fanoCavity2/fanoCavity2.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 /run/initialize
0125
0126 ---------------------------------------------------------
0127 ---> The Wall is 5 cm of G4_WATER ( 1 g/cm3 )
0128 The Cavity is 2 mm of g4Water_gas ( 1 mg/cm3 )
0129 ---------------------------------------------------------
0130
0131 /gun/energy 1 MeV
0132 /run/printProgress 10000
0133 /run/beamOn 40000
0134 =======================================================================
0135 ====== Electromagnetic Physics Parameters ========
0136 =======================================================================
0137 LPM effect enabled 1
0138 Enable creation and use of sampling tables 0
0139 Apply cuts on all EM processes 0
0140 Use combined TransportationWithMsc Disabled
0141 Use general process 0
0142 Enable linear polarisation for gamma 0
0143 Enable photoeffect sampling below K-shell 1
0144 Enable sampling of quantum entanglement 0
0145 X-section factor for integral approach 0.8
0146 Min kinetic energy for tables 10 eV
0147 Max kinetic energy for tables 10 GeV
0148 Number of bins per decade of a table 20
0149 Verbose level 1
0150 Verbose level for worker thread 0
0151 Bremsstrahlung energy threshold above which
0152 primary e+- is added to the list of secondary 100 TeV
0153 Bremsstrahlung energy threshold above which primary
0154 muon/hadron is added to the list of secondary 100 TeV
0155 Positron annihilation at rest model SimplePositronium
0156 Enable 3 gamma annihilation on fly 0
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 Use RiGe 5D e+e- pair production model by muons 0
0162 Livermore data directory epics_2017
0163 =======================================================================
0164 ====== Ionisation Parameters ========
0165 =======================================================================
0166 Step function for e+- (0.2, 0.1 mm)
0167 Step function for muons/hadrons (0.2, 0.1 mm)
0168 Step function for light ions (0.2, 0.1 mm)
0169 Step function for general ions (0.2, 0.1 mm)
0170 Lowest e+e- kinetic energy 100 eV
0171 Lowest muon/hadron kinetic energy 1 keV
0172 Use ICRU90 data 0
0173 Fluctuations of dE/dx are enabled 1
0174 Type of fluctuation model for leptons and hadrons Urban
0175 Use built-in Birks satuaration 0
0176 Build CSDA range enabled 1
0177 Use cut as a final range enabled 0
0178 Enable angular generator interface 0
0179 Max kinetic energy for CSDA tables 10 GeV
0180 Max kinetic energy for NIEL computation 0 eV
0181 Linear loss limit 0.01
0182 Read data from file for e+e- pair production by mu 0
0183 =======================================================================
0184 ====== Multiple Scattering Parameters ========
0185 =======================================================================
0186 Type of msc step limit algorithm for e+- 3
0187 Type of msc step limit algorithm for muons/hadrons 0
0188 Msc lateral displacement for e+- enabled 1
0189 Msc lateral displacement for muons and hadrons 0
0190 Urban msc model lateral displacement alg96 1
0191 Range factor for msc step limit for e+- 0.04
0192 Range factor for msc step limit for muons/hadrons 0.2
0193 Geometry factor for msc step limitation of e+- 2.5
0194 Safety factor for msc step limit for e+- 0.6
0195 Skin parameter for msc step limitation of e+- 1
0196 Lambda limit for msc step limit for e+- 1 mm
0197 Use Mott correction for e- scattering 0
0198 Factor used for dynamic computation of angular
0199 limit between single and multiple scattering 1
0200 Fixed angular limit between single
0201 and multiple scattering 3.1416 rad
0202 Upper energy limit for e+- multiple scattering 100 MeV
0203 Type of electron single scattering model 0
0204 Type of nuclear form-factor 1
0205 Screening factor 1
0206 =======================================================================
0207
0208 phot: for gamma SubType=12 BuildTable=0
0209 LambdaPrime table from 200 keV to 10 GeV in 94 bins
0210 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0211 PhotoElectric : Emin= 0 eV Emax= 10 GeV SauterGavrila
0212
0213 compt: for gamma SubType=13 BuildTable=1
0214 Lambda table from 10 eV to 1 MeV, 20 bins/decade, spline: 1
0215 LambdaPrime table from 1 MeV to 10 GeV in 80 bins
0216 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0217 Klein-Nishina : Emin= 0 eV Emax= 10 GeV
0218
0219 conv: for gamma SubType=14 BuildTable=1
0220 Lambda table from 1.022 MeV to 10 GeV, 25 bins/decade, spline: 1
0221 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0222 BetheHeitlerLPM : Emin= 0 eV Emax= 10 GeV ModifiedTsai
0223
0224 msc: for e- SubType= 10
0225 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0226 UrbanMsc : Emin= 0 eV Emax= 10 GeV Nbins=160 100 eV - 10 GeV
0227 StepLim=DistanceToBoundary Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0228
0229 eIoni: for e- XStype:1 SubType=2
0230 dE/dx and range tables from 10 eV to 10 GeV in 180 bins
0231 Lambda tables from threshold to 10 GeV, 20 bins/decade, spline: 1
0232 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01
0233 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0234 myMollerBhabha : Emin= 0 eV Emax= 10 GeV
0235 CSDA range table up to 10 GeV in 180 bins
0236
0237 msc: for e+ SubType= 10
0238 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0239 UrbanMsc : Emin= 0 eV Emax= 10 GeV Nbins=160 100 eV - 10 GeV
0240 StepLim=DistanceToBoundary Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0241
0242 eIoni: for e+ XStype:1 SubType=2
0243 dE/dx and range tables from 10 eV to 10 GeV in 180 bins
0244 Lambda tables from threshold to 10 GeV, 20 bins/decade, spline: 1
0245 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01
0246 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0247 myMollerBhabha : Emin= 0 eV Emax= 10 GeV
0248 CSDA range table up to 10 GeV in 180 bins
0249
0250 annihil: for e+ XStype:2 SubType=5 AtRestModel:Simple BuildTable=0
0251 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0252 eplus2gg : Emin= 0 eV Emax= 10 GeV
0253
0254 msc: for proton SubType= 10
0255 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0256 UrbanMsc : Emin= 0 eV Emax= 10 GeV Nbins=160 100 eV - 10 GeV
0257 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0258
0259 hIoni: for proton XStype:1 SubType=2
0260 dE/dx and range tables from 10 eV to 10 GeV in 180 bins
0261 Lambda tables from threshold to 10 GeV, 20 bins/decade, spline: 1
0262 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01
0263 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0264 Bragg : Emin= 0 eV Emax= 2 MeV
0265 BetheBloch : Emin= 2 MeV Emax= 10 GeV
0266 CSDA range table up to 10 GeV in 180 bins
0267
0268 ========= Table of registered couples ============================
0269
0270 Index : 0 used in the geometry : Yes
0271 Material : G4_WATER
0272 Range cuts : gamma 10 km e- 10 km e+ 10 km proton 10 km
0273 Energy thresholds : gamma 9.54965 GeV e- 9.54965 GeV e+ 9.54965 GeV proton 1 TeV
0274 Region(s) which use this couple :
0275 DefaultRegionForTheWorld
0276
0277 Index : 1 used in the geometry : Yes
0278 Material : g4Water_gas
0279 Range cuts : gamma 10 km e- 10 km e+ 10 km proton 10 km
0280 Energy thresholds : gamma 9.54965 GeV e- 9.54965 GeV e+ 9.54965 GeV proton 1 TeV
0281 Region(s) which use this couple :
0282 DefaultRegionForTheWorld
0283
0284 ==================================================================
0285
0286
0287 ===================== run conditions =====================
0288
0289 The run will be 0 e- of 1.000 MeV through 2*5.268 mm of G4_WATER (density: 1.000 g/cm3 ) Mass/cm2 = 1.054 g
0290 csdaRange: 4.390 mm
0291
0292 the cavity is 2.000 mm of g4Water_gas (density: 1.000 mg/cm3) Mass/cm2 = 0.200 mg --> massRatio = 0.000190
0293 Wall radius: 10.000 m range in cavity: 4.390 m
0294
0295 ==========================================================
0296
0297 StoppingPower in wall = 1.8649 MeV*cm2/g
0298 in cavity = 1.8649 MeV*cm2/g
0299 ### Run 0 starts.
0300
0301 ------- MixMaxRng engine status -------
0302 Current state vector is:
0303 mixmax state, file version 1.0
0304 N=17 V[N]={906770732717044781, 629165745432651234, 1235682547346241386, 68420008233404568, 2262190991329242458, 2266470399991071809, 1976726662926872232, 245458862506414172, 1955974201201518530, 2155248512522080758, 604170912935414061, 1116171330120743511, 1861018313684488333, 1296715403254578286, 1549011045957234151, 370819759640195970, 2230139271784837643} counter= 17sumtot= 1977567618660788324
0305 ---------------------------------------
0306 ### Run 0 start.
0307 --> Event 0 starts.
0308 --> Event 10000 starts.
0309 --->evntNb= 10000 Nwall= 9997 Ncav= 4 Ic/Iw= 2.10768 Ne-_cav= 2958 doseCavity/Ebeam= 1.18242 (100*(ratio-1) = 18.24224 %)
0310
0311 --> Event 20000 starts.
0312 --->evntNb= 20000 Nwall= 19995 Ncav= 6 Ic/Iw= 1.58068 Ne-_cav= 5838 doseCavity/Ebeam= 1.15579 (100*(ratio-1) = 15.57930 %)
0313
0314 --> Event 30000 starts.
0315 --->evntNb= 30000 Nwall= 29992 Ncav= 9 Ic/Iw= 1.58071 Ne-_cav= 8657 doseCavity/Ebeam= 1.07851 (100*(ratio-1) = 7.85147 %)
0316
0317 Run terminated.
0318 Run Summary
0319 Number of events processed : 40000
0320 User=0.930000s Real=0.940515s Sys=0.000000s
0321
0322 Process calls frequency ---> eIoni= 315000 msc= 275995 Transportation= 22864
0323
0324 Charged particle flow in cavity :
0325 Enter --> nbParticles = 11430 Energy = 7.218 GeV
0326 Exit --> nbParticles = 11434 Energy = 7.223 GeV
0327
0328 beamFluence in wall = 39987 in cavity = 13 Icav/Iwall = 1.71253 energyFluence = 37.96066 MeV*cm2/mg
0329
0330 Total edep in cavity = 8.08834 MeV +- 4.39251 %
0331 Total dose in cavity = 40.44169 MeV*cm2/mg +- 4.39251 %
0332
0333 DoseCavity/EnergyFluence = 1.06536 +- 0.04680
0334
0335 Total charged trackLength in cavity = 36.6239 m (mean value = 3.2042 mm )
0336
0337 StepSize of ch. tracks in wall = 299.9595 um +- 255.1711 um (nbSteps/track = 14.7233)
0338 StepSize of ch. tracks in cavity = 1.4581 mm +- 2.6789 mm (nbSteps/track = 2.1976)
0339
0340 ------- MixMaxRng engine status -------
0341 Current state vector is:
0342 mixmax state, file version 1.0
0343 N=17 V[N]={1122560079547894776, 489338179114589148, 962953981966901183, 1371396244330764127, 1216257719796551649, 871351706115100433, 1331924270603534972, 1161597106298325732, 1677528997079743374, 561171356176391074, 1761454154557839398, 593086082936043028, 652088176841567984, 1436493037053116940, 1396378354312194542, 1141314676005422410, 747120547568034109} counter= 17sumtot= 47270596594463271
0344 ---------------------------------------
0345
0346 ------- MixMaxRng engine status -------
0347 Current state vector is:
0348 mixmax state, file version 1.0
0349 N=17 V[N]={1122560079547894776, 489338179114589148, 962953981966901183, 1371396244330764127, 1216257719796551649, 871351706115100433, 1331924270603534972, 1161597106298325732, 1677528997079743374, 561171356176391074, 1761454154557839398, 593086082936043028, 652088176841567984, 1436493037053116940, 1396378354312194542, 1141314676005422410, 747120547568034109} counter= 17sumtot= 47270596594463271
0350 ---------------------------------------
0351 ================== Deleting memory pools ===================
0352 Number of memory pools allocated: 9 of which, static: 0
0353 Dynamic pools deleted: 9 / Total memory freed: 0.014 MB
0354 ============================================================
