Warning, /geant4/examples/extended/medical/GammaTherapy/GammaTherapy.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 ===== GammaTherapy is started with 1 threads =====
0023
0024 ***** Table : Nb of materials = 6 *****
0025
0026 Material: G4_Be density: 1.848 g/cm3 RadL: 35.276 cm Nucl.Int.Length: 39.413 cm
0027 Imean: 63.700 eV temperature: 293.15 K pressure: 1.00 atm
0028
0029 ---> Element: Be (Be) Z = 4.0 N = 9 A = 9.012 g/mole
0030 ---> Isotope: Be9 Z = 4 N = 9 A = 9.01 g/mole abundance: 100.000 %
0031 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0032
0033
0034 Material: G4_W density: 19.300 g/cm3 RadL: 3.504 mm Nucl.Int.Length: 10.312 cm
0035 Imean: 727.000 eV temperature: 293.15 K pressure: 1.00 atm
0036
0037 ---> Element: W (W) Z = 74.0 N = 184 A = 183.842 g/mole
0038 ---> Isotope: W180 Z = 74 N = 180 A = 179.95 g/mole abundance: 0.120 %
0039 ---> Isotope: W182 Z = 74 N = 182 A = 181.95 g/mole abundance: 26.500 %
0040 ---> Isotope: W183 Z = 74 N = 183 A = 182.95 g/mole abundance: 14.310 %
0041 ---> Isotope: W184 Z = 74 N = 184 A = 183.95 g/mole abundance: 30.640 %
0042 ---> Isotope: W186 Z = 74 N = 186 A = 185.95 g/mole abundance: 28.430 %
0043 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0044
0045
0046 Material: G4_He density: 0.166 mg/cm3 RadL: 5.671 km Nucl.Int.Length: 3.341 km
0047 Imean: 41.800 eV temperature: 293.15 K pressure: 1.00 atm
0048
0049 ---> Element: He (He) Z = 2.0 N = 4 A = 4.003 g/mole
0050 ---> Isotope: He3 Z = 2 N = 3 A = 3.02 g/mole abundance: 0.000 %
0051 ---> Isotope: He4 Z = 2 N = 4 A = 4.00 g/mole abundance: 100.000 %
0052 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0053
0054
0055 Material: G4_WATER H_2O density: 1.000 g/cm3 RadL: 36.083 cm Nucl.Int.Length: 75.375 cm
0056 Imean: 78.000 eV temperature: 293.15 K pressure: 1.00 atm
0057
0058 ---> Element: H (H) Z = 1.0 N = 1 A = 1.008 g/mole
0059 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 %
0060 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 %
0061 ElmMassFraction: 11.19 % ElmAbundance 66.67 %
0062
0063 ---> Element: O (O) Z = 8.0 N = 16 A = 15.999 g/mole
0064 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 %
0065 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 %
0066 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 %
0067 ElmMassFraction: 88.81 % ElmAbundance 33.33 %
0068
0069
0070 Material: G4_AIR density: 1.205 mg/cm3 RadL: 303.921 m Nucl.Int.Length: 710.095 m
0071 Imean: 85.700 eV temperature: 293.15 K pressure: 1.00 atm
0072
0073 ---> Element: C (C) Z = 6.0 N = 12 A = 12.011 g/mole
0074 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.930 %
0075 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.070 %
0076 ElmMassFraction: 0.01 % ElmAbundance 0.02 %
0077
0078 ---> Element: N (N) Z = 7.0 N = 14 A = 14.007 g/mole
0079 ---> Isotope: N14 Z = 7 N = 14 A = 14.00 g/mole abundance: 99.632 %
0080 ---> Isotope: N15 Z = 7 N = 15 A = 15.00 g/mole abundance: 0.368 %
0081 ElmMassFraction: 75.53 % ElmAbundance 78.44 %
0082
0083 ---> Element: O (O) Z = 8.0 N = 16 A = 15.999 g/mole
0084 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 %
0085 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 %
0086 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 %
0087 ElmMassFraction: 23.18 % ElmAbundance 21.07 %
0088
0089 ---> Element: Ar (Ar) Z = 18.0 N = 40 A = 39.948 g/mole
0090 ---> Isotope: Ar36 Z = 18 N = 36 A = 35.97 g/mole abundance: 0.337 %
0091 ---> Isotope: Ar38 Z = 18 N = 38 A = 37.96 g/mole abundance: 0.063 %
0092 ---> Isotope: Ar40 Z = 18 N = 40 A = 39.96 g/mole abundance: 99.600 %
0093 ElmMassFraction: 1.28 % ElmAbundance 0.47 %
0094
0095
0096 Material: G4_MYLAR density: 1.400 g/cm3 RadL: 28.536 cm Nucl.Int.Length: 55.923 cm
0097 Imean: 78.700 eV temperature: 293.15 K pressure: 1.00 atm
0098
0099 ---> Element: C (C) Z = 6.0 N = 12 A = 12.011 g/mole
0100 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.930 %
0101 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.070 %
0102 ElmMassFraction: 62.50 % ElmAbundance 45.45 %
0103
0104 ---> Element: H (H) Z = 1.0 N = 1 A = 1.008 g/mole
0105 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 %
0106 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 %
0107 ElmMassFraction: 4.20 % ElmAbundance 36.36 %
0108
0109 ---> Element: O (O) Z = 8.0 N = 16 A = 15.999 g/mole
0110 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 %
0111 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 %
0112 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 %
0113 ElmMassFraction: 33.30 % ElmAbundance 18.18 %
0114
0115
0116
0117 ### PhysicsList Construte Particles
0118 /control/cout/ignoreThreadsExcept 0
0119 /run/verbose 1
0120 /tracking/verbose 0
0121 /event/verbose 0
0122 /testem/verbose 0
0123 #
0124 /testem/target1Material G4_Be
0125 /testem/target2Material G4_W
0126 /testem/target1Z 9 mm
0127 /testem/target2Z 6 mm
0128 #
0129 /testem/phys/addPhysics emstandard
0130 ### PhysicsList Add Physics <emstandard>
0131 #/testem/phys/addPhysics emstandard_opt2
0132 #/testem/phys/addPhysics empenelope
0133 #/testem/phys/addPhysics emlivermore
0134 #/testem/phys/addPhysics step_limit
0135 #
0136 #/testem/histoName be_w_st
0137 #
0138 #/testem/stepMax 1.0 mm
0139 #
0140 /run/setCut 0.1 mm
0141 #
0142 /run/initialize
0143 ===================================================
0144 # GammaTherapy Geometry #
0145 ===================================================
0146 World width= 747 mm
0147 Window width= 0.05 mm position = -622.525 mm:
0148 TargetV width= 45.001 mm position = -600 mm:
0149 Target1 width= 9 mm position = 11.9995 mm:
0150 Target2 width= 6 mm position = 19.4995 mm:
0151 Gas width= 705 mm position = -224.999 mm:
0152 Mylar width= 0.02 mm position = 127.512 mm:
0153 Check width= 0.1 mm position = 327.551 mm:
0154 Air width= 210 mm
0155 Phantom width= 300 mm position = 487.501 mm:
0156 Absorb width= 5 mm position = -77.5 mm:
0157 Absorb shift= 337.501 mm
0158 Target1 G4_Be
0159 Target2 G4_W
0160 Phantom G4_WATER
0161 ===================================================
0162 Phantom R= 200 dz= 150
0163 Absorber R= 200 dz= 2.5 posZ= -77.5
0164 ### PhysicsList Construte Processes
0165 #
0166 /gun/particle e-
0167 /gun/energy 50.0 MeV
0168 /run/beamOn 10000
0169 =======================================================================
0170 ====== Electromagnetic Physics Parameters ========
0171 =======================================================================
0172 LPM effect enabled 1
0173 Enable creation and use of sampling tables 0
0174 Apply cuts on all EM processes 0
0175 Use combined TransportationWithMsc Disabled
0176 Use general process 1
0177 Enable linear polarisation for gamma 0
0178 Enable photoeffect sampling below K-shell 1
0179 Enable sampling of quantum entanglement 0
0180 X-section factor for integral approach 0.8
0181 Min kinetic energy for tables 10 eV
0182 Max kinetic energy for tables 10 GeV
0183 Number of bins per decade of a table 10
0184 Verbose level 1
0185 Verbose level for worker thread 0
0186 Bremsstrahlung energy threshold above which
0187 primary e+- is added to the list of secondary 100 TeV
0188 Bremsstrahlung energy threshold above which primary
0189 muon/hadron is added to the list of secondary 100 TeV
0190 Positron annihilation at rest model SimplePositronium
0191 Enable 3 gamma annihilation on fly 0
0192 Lowest triplet kinetic energy 1 MeV
0193 Enable sampling of gamma linear polarisation 0
0194 5D gamma conversion model type 0
0195 5D gamma conversion model on isolated ion 0
0196 Use RiGe 5D e+e- pair production model by muons 0
0197 Livermore data directory epics_2017
0198 =======================================================================
0199 ====== Ionisation Parameters ========
0200 =======================================================================
0201 Step function for e+- (0.2, 1 mm)
0202 Step function for muons/hadrons (0.2, 0.1 mm)
0203 Step function for light ions (0.2, 0.1 mm)
0204 Step function for general ions (0.2, 0.1 mm)
0205 Lowest e+e- kinetic energy 1 keV
0206 Lowest muon/hadron kinetic energy 1 keV
0207 Use ICRU90 data 0
0208 Fluctuations of dE/dx are enabled 1
0209 Type of fluctuation model for leptons and hadrons Urban
0210 Use built-in Birks satuaration 0
0211 Build CSDA range enabled 0
0212 Use cut as a final range enabled 0
0213 Enable angular generator interface 0
0214 Max kinetic energy for CSDA tables 1 GeV
0215 Max kinetic energy for NIEL computation 0 eV
0216 Linear loss limit 0.01
0217 Read data from file for e+e- pair production by mu 0
0218 =======================================================================
0219 ====== Multiple Scattering Parameters ========
0220 =======================================================================
0221 Type of msc step limit algorithm for e+- 1
0222 Type of msc step limit algorithm for muons/hadrons 0
0223 Msc lateral displacement for e+- enabled 1
0224 Msc lateral displacement for muons and hadrons 0
0225 Urban msc model lateral displacement alg96 1
0226 Range factor for msc step limit for e+- 0.04
0227 Range factor for msc step limit for muons/hadrons 0.2
0228 Geometry factor for msc step limitation of e+- 2.5
0229 Safety factor for msc step limit for e+- 0.6
0230 Skin parameter for msc step limitation of e+- 1
0231 Lambda limit for msc step limit for e+- 1 mm
0232 Use Mott correction for e- scattering 0
0233 Factor used for dynamic computation of angular
0234 limit between single and multiple scattering 1
0235 Fixed angular limit between single
0236 and multiple scattering 3.1416 rad
0237 Upper energy limit for e+- multiple scattering 100 MeV
0238 Type of electron single scattering model 0
0239 Type of nuclear form-factor 1
0240 Screening factor 1
0241 =======================================================================
0242
0243 phot: for gamma SubType=12 BuildTable=0
0244 LambdaPrime table from 200 keV to 10 GeV in 47 bins
0245 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0246 LivermorePhElectric : Emin= 0 eV Emax= 10 GeV SauterGavrila Fluo
0247
0248 compt: for gamma SubType=13 BuildTable=1
0249 Lambda table from 10 eV to 1 MeV, 10 bins/decade, spline: 1
0250 LambdaPrime table from 1 MeV to 10 GeV in 40 bins
0251 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0252 Klein-Nishina : Emin= 0 eV Emax= 10 GeV
0253
0254 conv: for gamma SubType=14 BuildTable=1
0255 Lambda table from 1.022 MeV to 10 GeV, 25 bins/decade, spline: 1
0256 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0257 BetheHeitlerLPM : Emin= 0 eV Emax= 10 GeV ModifiedTsai
0258
0259 Rayl: for gamma SubType=11 BuildTable=1
0260 Lambda table from 10 eV to 150 keV, 10 bins/decade, spline: 0
0261 LambdaPrime table from 150 keV to 10 GeV in 48 bins
0262 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0263 LivermoreRayleigh : Emin= 0 eV Emax= 10 GeV CullenGenerator
0264
0265 msc: for e- SubType= 10
0266 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0267 UrbanMsc : Emin= 0 eV Emax= 100 MeV Nbins=60 100 eV - 100 MeV
0268 StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0269 WentzelVIUni : Emin= 100 MeV Emax= 10 GeV Nbins=20 100 MeV - 10 GeV
0270 StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0271
0272 eIoni: for e- XStype:3 SubType=2
0273 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0274 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0275 StepFunction=(0.2, 1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0276 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0277 MollerBhabha : Emin= 0 eV Emax= 10 GeV
0278
0279 eBrem: for e- XStype:4 SubType=3
0280 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0281 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0282 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000
0283 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0284 eBremSB : Emin= 0 eV Emax= 1 GeV ModifiedTsai
0285 eBremLPM : Emin= 1 GeV Emax= 10 GeV ModifiedTsai
0286
0287 CoulombScat: for e- XStype:1 SubType=1 BuildTable=1
0288 Lambda table from 100 MeV to 10 GeV, 10 bins/decade, spline: 0
0289 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0290 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0291 eCoulombScattering : Emin= 100 MeV Emax= 10 GeV
0292
0293 msc: for e+ SubType= 10
0294 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0295 UrbanMsc : Emin= 0 eV Emax= 100 MeV Nbins=60 100 eV - 100 MeV
0296 StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0297 WentzelVIUni : Emin= 100 MeV Emax= 10 GeV Nbins=20 100 MeV - 10 GeV
0298 StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0299
0300 eIoni: for e+ XStype:3 SubType=2
0301 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0302 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0303 StepFunction=(0.2, 1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0304 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0305 MollerBhabha : Emin= 0 eV Emax= 10 GeV
0306
0307 eBrem: for e+ XStype:4 SubType=3
0308 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0309 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0310 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000
0311 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0312 eBremSB : Emin= 0 eV Emax= 1 GeV ModifiedTsai
0313 eBremLPM : Emin= 1 GeV Emax= 10 GeV ModifiedTsai
0314
0315 annihil: for e+ XStype:2 SubType=5 AtRestModel:Simple BuildTable=0
0316 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0317 eplus2gg : Emin= 0 eV Emax= 10 GeV
0318
0319 CoulombScat: for e+ XStype:1 SubType=1 BuildTable=1
0320 Lambda table from 100 MeV to 10 GeV, 10 bins/decade, spline: 0
0321 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0322 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0323 eCoulombScattering : Emin= 100 MeV Emax= 10 GeV
0324
0325 msc: for proton SubType= 10
0326 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0327 WentzelVIUni : Emin= 0 eV Emax= 10 GeV Nbins=80 100 eV - 10 GeV
0328 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0329
0330 hIoni: for proton XStype:3 SubType=2
0331 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0332 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0333 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0334 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0335 Bragg : Emin= 0 eV Emax= 2 MeV
0336 BetheBloch : Emin= 2 MeV Emax= 10 GeV
0337
0338 hBrems: for proton XStype:1 SubType=3
0339 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0340 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0341 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0342 hBrem : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0343
0344 hPairProd: for proton XStype:1 SubType=4
0345 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0346 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0347 Sampling table 4x1001 from 7.50618 GeV to 0.0150124 TeV
0348 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0349 hPairProd : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0350
0351 CoulombScat: for proton XStype:1 SubType=1 BuildTable=1
0352 Lambda table from threshold to 10 GeV, 10 bins/decade, spline: 0
0353 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0354 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0355 eCoulombScattering : Emin= 0 eV Emax= 10 GeV
0356
0357 msc: for GenericIon SubType= 10
0358 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0359 UrbanMsc : Emin= 0 eV Emax= 10 GeV
0360 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0361
0362 ionIoni: for GenericIon XStype:3 SubType=2
0363 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0364 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0365 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.02
0366 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0367 Bragg : Emin= 0 eV Emax= 2 MeV
0368 BetheBloch : Emin= 2 MeV Emax= 10 GeV
0369
0370 msc: for alpha SubType= 10
0371 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0372 UrbanMsc : Emin= 0 eV Emax= 10 GeV
0373 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0374
0375 ionIoni: for alpha XStype:3 SubType=2
0376 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0377 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0378 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.02
0379 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0380 BraggIon : Emin= 0 eV Emax=7.9452 MeV
0381 BetheBloch : Emin=7.9452 MeV Emax= 10 GeV
0382
0383 msc: for anti_proton SubType= 10
0384 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0385 WentzelVIUni : Emin= 0 eV Emax= 10 GeV Nbins=80 100 eV - 10 GeV
0386 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0387
0388 hIoni: for anti_proton XStype:3 SubType=2
0389 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0390 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0391 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0392 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0393 ICRU73QO : Emin= 0 eV Emax= 2 MeV
0394 BetheBloch : Emin= 2 MeV Emax= 10 GeV
0395
0396 hBrems: for anti_proton XStype:1 SubType=3
0397 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0398 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0399 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0400 hBrem : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0401
0402 hPairProd: for anti_proton XStype:1 SubType=4
0403 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0404 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0405 Sampling table 4x1001 from 7.50618 GeV to 0.0150124 TeV
0406 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0407 hPairProd : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0408
0409 CoulombScat: for anti_proton XStype:1 SubType=1 BuildTable=1
0410 Lambda table from threshold to 10 GeV, 10 bins/decade, spline: 0
0411 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0412 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0413 eCoulombScattering : Emin= 0 eV Emax= 10 GeV
0414
0415 msc: for kaon+ SubType= 10
0416 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0417 WentzelVIUni : Emin= 0 eV Emax= 10 GeV Nbins=80 100 eV - 10 GeV
0418 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0419
0420 hIoni: for kaon+ XStype:3 SubType=2
0421 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0422 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0423 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0424 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0425 Bragg : Emin= 0 eV Emax=1.05231 MeV
0426 BetheBloch : Emin=1.05231 MeV Emax= 10 GeV
0427
0428 hBrems: for kaon+ XStype:1 SubType=3
0429 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0430 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0431 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0432 hBrem : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0433
0434 hPairProd: for kaon+ XStype:1 SubType=4
0435 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0436 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0437 Sampling table 4x1001 from 3.94942 GeV to 0.01 TeV
0438 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0439 hPairProd : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0440
0441 CoulombScat: for kaon+ XStype:1 SubType=1 BuildTable=1
0442 Lambda table from threshold to 10 GeV, 10 bins/decade, spline: 0
0443 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0444 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0445 eCoulombScattering : Emin= 0 eV Emax= 10 GeV
0446
0447 msc: for kaon- SubType= 10
0448 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0449 WentzelVIUni : Emin= 0 eV Emax= 10 GeV Nbins=80 100 eV - 10 GeV
0450 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0451
0452 hIoni: for kaon- XStype:3 SubType=2
0453 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0454 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0455 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0456 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0457 ICRU73QO : Emin= 0 eV Emax=1.05231 MeV
0458 BetheBloch : Emin=1.05231 MeV Emax= 10 GeV
0459
0460 hBrems: for kaon- XStype:1 SubType=3
0461 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0462 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0463 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0464 hBrem : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0465
0466 hPairProd: for kaon- XStype:1 SubType=4
0467 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0468 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0469 Sampling table 4x1001 from 3.94942 GeV to 0.01 TeV
0470 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0471 hPairProd : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0472
0473 CoulombScat: for kaon- XStype:1 SubType=1 BuildTable=1
0474 Used Lambda table of kaon+
0475 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0476 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0477 eCoulombScattering : Emin= 0 eV Emax= 10 GeV
0478
0479 msc: for mu+ SubType= 10
0480 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0481 WentzelVIUni : Emin= 0 eV Emax= 10 GeV Nbins=80 100 eV - 10 GeV
0482 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0483
0484 muIoni: for mu+ XStype:3 SubType=2
0485 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0486 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0487 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0488 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0489 Bragg : Emin= 0 eV Emax= 200 keV
0490 MuBetheBloch : Emin= 200 keV Emax= 10 GeV
0491
0492 muBrems: for mu+ XStype:1 SubType=3
0493 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0494 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0495 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0496 MuBrem : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0497
0498 muPairProd: for mu+ XStype:1 SubType=4
0499 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0500 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0501 Sampling table 5x1001 from 0.85 GeV to 0.01 TeV
0502 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0503 muPairProd : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0504
0505 CoulombScat: for mu+ XStype:1 SubType=1 BuildTable=1
0506 Lambda table from threshold to 10 GeV, 10 bins/decade, spline: 0
0507 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0508 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0509 eCoulombScattering : Emin= 0 eV Emax= 10 GeV
0510
0511 msc: for mu- SubType= 10
0512 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0513 WentzelVIUni : Emin= 0 eV Emax= 10 GeV Nbins=80 100 eV - 10 GeV
0514 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0515
0516 muIoni: for mu- XStype:3 SubType=2
0517 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0518 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0519 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0520 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0521 ICRU73QO : Emin= 0 eV Emax= 200 keV
0522 MuBetheBloch : Emin= 200 keV Emax= 10 GeV
0523
0524 muBrems: for mu- XStype:1 SubType=3
0525 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0526 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0527 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0528 MuBrem : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0529
0530 muPairProd: for mu- XStype:1 SubType=4
0531 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0532 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0533 Sampling table 5x1001 from 0.85 GeV to 0.01 TeV
0534 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0535 muPairProd : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0536
0537 CoulombScat: for mu- XStype:1 SubType=1 BuildTable=1
0538 Used Lambda table of mu+
0539 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0540 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0541 eCoulombScattering : Emin= 0 eV Emax= 10 GeV
0542
0543 msc: for pi+ SubType= 10
0544 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0545 WentzelVIUni : Emin= 0 eV Emax= 10 GeV Nbins=80 100 eV - 10 GeV
0546 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0547
0548 hIoni: for pi+ XStype:3 SubType=2
0549 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0550 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0551 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0552 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0553 Bragg : Emin= 0 eV Emax=297.505 keV
0554 BetheBloch : Emin=297.505 keV Emax= 10 GeV
0555
0556 hBrems: for pi+ XStype:1 SubType=3
0557 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0558 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0559 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0560 hBrem : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0561
0562 hPairProd: for pi+ XStype:1 SubType=4
0563 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0564 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0565 Sampling table 4x1001 from 1.11656 GeV to 0.01 TeV
0566 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0567 hPairProd : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0568
0569 CoulombScat: for pi+ XStype:1 SubType=1 BuildTable=1
0570 Lambda table from threshold to 10 GeV, 10 bins/decade, spline: 0
0571 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0572 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0573 eCoulombScattering : Emin= 0 eV Emax= 10 GeV
0574
0575 msc: for pi- SubType= 10
0576 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0577 WentzelVIUni : Emin= 0 eV Emax= 10 GeV Nbins=80 100 eV - 10 GeV
0578 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0579
0580 hIoni: for pi- XStype:3 SubType=2
0581 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0582 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0583 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0584 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0585 ICRU73QO : Emin= 0 eV Emax=297.505 keV
0586 BetheBloch : Emin=297.505 keV Emax= 10 GeV
0587
0588 hBrems: for pi- XStype:1 SubType=3
0589 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0590 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0591 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0592 hBrem : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0593
0594 hPairProd: for pi- XStype:1 SubType=4
0595 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0596 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0597 Sampling table 4x1001 from 1.11656 GeV to 0.01 TeV
0598 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0599 hPairProd : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0600
0601 CoulombScat: for pi- XStype:1 SubType=1 BuildTable=1
0602 Used Lambda table of pi+
0603 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0604 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0605 eCoulombScattering : Emin= 0 eV Emax= 10 GeV
0606
0607 ========= Table of registered couples ============================
0608
0609 Index : 0 used in the geometry : Yes
0610 Material : G4_AIR
0611 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 100 um
0612 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 10 keV
0613 Region(s) which use this couple :
0614 DefaultRegionForTheWorld
0615
0616 Index : 1 used in the geometry : Yes
0617 Material : G4_Be
0618 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 100 um
0619 Energy thresholds : gamma 1.0488 keV e- 108.985 keV e+ 107.379 keV proton 10 keV
0620 Region(s) which use this couple :
0621 DefaultRegionForTheWorld
0622
0623 Index : 2 used in the geometry : Yes
0624 Material : G4_He
0625 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 100 um
0626 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 10 keV
0627 Region(s) which use this couple :
0628 DefaultRegionForTheWorld
0629
0630 Index : 3 used in the geometry : Yes
0631 Material : G4_W
0632 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 100 um
0633 Energy thresholds : gamma 36.2784 keV e- 351.887 keV e+ 338.074 keV proton 10 keV
0634 Region(s) which use this couple :
0635 DefaultRegionForTheWorld
0636
0637 Index : 4 used in the geometry : Yes
0638 Material : G4_MYLAR
0639 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 100 um
0640 Energy thresholds : gamma 1.20959 keV e- 99.1194 keV e+ 97.66 keV proton 10 keV
0641 Region(s) which use this couple :
0642 DefaultRegionForTheWorld
0643
0644 Index : 5 used in the geometry : Yes
0645 Material : G4_WATER
0646 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 100 um
0647 Energy thresholds : gamma 1.11075 keV e- 85.0445 keV e+ 83.9021 keV proton 10 keV
0648 Region(s) which use this couple :
0649 DefaultRegionForTheWorld
0650
0651 ==================================================================
0652
0653 80 bins R stepR= 2.5 mm
0654 60 bins Z stepZ= 5 mm
0655 200 bins E stepE= 0.25 MeV
0656 14-th bin in Z is used for R distribution
0657 ### Run 0 starts.
0658
0659 ------- MixMaxRng engine status -------
0660 Current state vector is:
0661 mixmax state, file version 1.0
0662 N=17 V[N]={906770732717044781, 629165745432651234, 1235682547346241386, 68420008233404568, 2262190991329242458, 2266470399991071809, 1976726662926872232, 245458862506414172, 1955974201201518530, 2155248512522080758, 604170912935414061, 1116171330120743511, 1861018313684488333, 1296715403254578286, 1549011045957234151, 370819759640195970, 2230139271784837643} counter= 17sumtot= 1977567618660788324
0663 ---------------------------------------
0664 Run terminated.
0665 Run Summary
0666 Number of events processed : 10000
0667 User=1.600000s Real=1.628730s Sys=0.000000s
0668 Histo: End of run actions are started
0669 ========================================================
0670 Number of events 10000
0671 Average number of e- 33.13
0672 Average number of gamma 19.58
0673 Average number of e+ 0.7605
0674 Average number of steps in the phantom 11.51
0675 Average number of e- steps in the target 17.98
0676 Average number of g produced in the target 18.72
0677 Average number of e- produced in the target 20.83
0678 Average number of g produced in the phantom 0.4933
0679 Average number of e- produced in the phantom 5.953
0680 Total fGamma fluence in front of the phantom 7.954 MeV
0681 ========================================================
0682
0683
0684
0685 ------- MixMaxRng engine status -------
0686 Current state vector is:
0687 mixmax state, file version 1.0
0688 N=17 V[N]={961145400483405076, 699745897242935866, 1057413413739364408, 662555405695072145, 1718128862343597050, 1055151940051380955, 2185583859974491267, 673048765014179189, 467115907169885361, 358922768063147600, 1116951267450062581, 2116545342320611943, 2097059162164451181, 93621624613885263, 529728786058221144, 735650886969976136, 590625775673630549} counter= 4sumtot= 978094000532440057
0689 ---------------------------------------
0690 #
0691 ================== Deleting memory pools ===================
0692 Number of memory pools allocated: 10 of which, static: 0
0693 Dynamic pools deleted: 10 / Total memory freed: 0.058 MB
0694 ============================================================
