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-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 ===== 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 Lowest triplet kinetic energy 1 MeV
0191 Enable sampling of gamma linear polarisation 0
0192 5D gamma conversion model type 0
0193 5D gamma conversion model on isolated ion 0
0194 Livermore data directory epics_2017
0195 =======================================================================
0196 ====== Ionisation Parameters ========
0197 =======================================================================
0198 Step function for e+- (0.2, 1 mm)
0199 Step function for muons/hadrons (0.2, 0.1 mm)
0200 Step function for light ions (0.2, 0.1 mm)
0201 Step function for general ions (0.2, 0.1 mm)
0202 Lowest e+e- kinetic energy 1 keV
0203 Lowest muon/hadron kinetic energy 1 keV
0204 Use ICRU90 data 0
0205 Fluctuations of dE/dx are enabled 1
0206 Type of fluctuation model for leptons and hadrons Urban
0207 Use built-in Birks satuaration 0
0208 Build CSDA range enabled 0
0209 Use cut as a final range enabled 0
0210 Enable angular generator interface 0
0211 Max kinetic energy for CSDA tables 1 GeV
0212 Max kinetic energy for NIEL computation 0 eV
0213 Linear loss limit 0.01
0214 Read data from file for e+e- pair production by mu 0
0215 =======================================================================
0216 ====== Multiple Scattering Parameters ========
0217 =======================================================================
0218 Type of msc step limit algorithm for e+- 1
0219 Type of msc step limit algorithm for muons/hadrons 0
0220 Msc lateral displacement for e+- enabled 1
0221 Msc lateral displacement for muons and hadrons 0
0222 Urban msc model lateral displacement alg96 1
0223 Range factor for msc step limit for e+- 0.04
0224 Range factor for msc step limit for muons/hadrons 0.2
0225 Geometry factor for msc step limitation of e+- 2.5
0226 Safety factor for msc step limit for e+- 0.6
0227 Skin parameter for msc step limitation of e+- 1
0228 Lambda limit for msc step limit for e+- 1 mm
0229 Use Mott correction for e- scattering 0
0230 Factor used for dynamic computation of angular
0231 limit between single and multiple scattering 1
0232 Fixed angular limit between single
0233 and multiple scattering 3.1416 rad
0234 Upper energy limit for e+- multiple scattering 100 MeV
0235 Type of electron single scattering model 0
0236 Type of nuclear form-factor 1
0237 Screening factor 1
0238 =======================================================================
0239
0240 phot: for gamma SubType=12 BuildTable=0
0241 LambdaPrime table from 200 keV to 10 GeV in 47 bins
0242 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0243 LivermorePhElectric : Emin= 0 eV Emax= 10 GeV SauterGavrila Fluo
0244
0245 compt: for gamma SubType=13 BuildTable=1
0246 Lambda table from 10 eV to 1 MeV, 10 bins/decade, spline: 1
0247 LambdaPrime table from 1 MeV to 10 GeV in 40 bins
0248 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0249 Klein-Nishina : Emin= 0 eV Emax= 10 GeV
0250
0251 conv: for gamma SubType=14 BuildTable=1
0252 Lambda table from 1.022 MeV to 10 GeV, 25 bins/decade, spline: 1
0253 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0254 BetheHeitlerLPM : Emin= 0 eV Emax= 10 GeV ModifiedTsai
0255
0256 Rayl: for gamma SubType=11 BuildTable=1
0257 Lambda table from 10 eV to 150 keV, 10 bins/decade, spline: 0
0258 LambdaPrime table from 150 keV to 10 GeV in 48 bins
0259 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0260 LivermoreRayleigh : Emin= 0 eV Emax= 10 GeV CullenGenerator
0261
0262 msc: for e- SubType= 10
0263 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0264 UrbanMsc : Emin= 0 eV Emax= 100 MeV Nbins=60 100 eV - 100 MeV
0265 StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0266 WentzelVIUni : Emin= 100 MeV Emax= 10 GeV Nbins=20 100 MeV - 10 GeV
0267 StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0268
0269 eIoni: for e- XStype:3 SubType=2
0270 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0271 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0272 StepFunction=(0.2, 1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0273 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0274 MollerBhabha : Emin= 0 eV Emax= 10 GeV
0275
0276 eBrem: for e- XStype:4 SubType=3
0277 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0278 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0279 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000
0280 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0281 eBremSB : Emin= 0 eV Emax= 1 GeV ModifiedTsai
0282 eBremLPM : Emin= 1 GeV Emax= 10 GeV ModifiedTsai
0283
0284 CoulombScat: for e- XStype:1 SubType=1 BuildTable=1
0285 Lambda table from 100 MeV to 10 GeV, 10 bins/decade, spline: 0
0286 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0287 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0288 eCoulombScattering : Emin= 100 MeV Emax= 10 GeV
0289
0290 msc: for e+ SubType= 10
0291 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0292 UrbanMsc : Emin= 0 eV Emax= 100 MeV Nbins=60 100 eV - 100 MeV
0293 StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0294 WentzelVIUni : Emin= 100 MeV Emax= 10 GeV Nbins=20 100 MeV - 10 GeV
0295 StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0296
0297 eIoni: for e+ XStype:3 SubType=2
0298 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0299 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0300 StepFunction=(0.2, 1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0301 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0302 MollerBhabha : Emin= 0 eV Emax= 10 GeV
0303
0304 eBrem: for e+ XStype:4 SubType=3
0305 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0306 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0307 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000
0308 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0309 eBremSB : Emin= 0 eV Emax= 1 GeV ModifiedTsai
0310 eBremLPM : Emin= 1 GeV Emax= 10 GeV ModifiedTsai
0311
0312 annihil: for e+ XStype:2 SubType=5 AtRestModel:Simple BuildTable=0
0313 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0314 eplus2gg : Emin= 0 eV Emax= 10 GeV
0315
0316 CoulombScat: for e+ XStype:1 SubType=1 BuildTable=1
0317 Lambda table from 100 MeV to 10 GeV, 10 bins/decade, spline: 0
0318 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0319 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0320 eCoulombScattering : Emin= 100 MeV Emax= 10 GeV
0321
0322 msc: for proton SubType= 10
0323 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0324 WentzelVIUni : Emin= 0 eV Emax= 10 GeV Nbins=80 100 eV - 10 GeV
0325 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0326
0327 hIoni: for proton XStype:3 SubType=2
0328 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0329 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0330 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0331 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0332 Bragg : Emin= 0 eV Emax= 2 MeV
0333 BetheBloch : Emin= 2 MeV Emax= 10 GeV
0334
0335 hBrems: for proton XStype:1 SubType=3
0336 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0337 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0338 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0339 hBrem : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0340
0341 hPairProd: for proton XStype:1 SubType=4
0342 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0343 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0344 Sampling table 4x1001 from 7.50618 GeV to 0.0150124 TeV
0345 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0346 hPairProd : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0347
0348 CoulombScat: for proton XStype:1 SubType=1 BuildTable=1
0349 Lambda table from threshold to 10 GeV, 10 bins/decade, spline: 0
0350 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0351 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0352 eCoulombScattering : Emin= 0 eV Emax= 10 GeV
0353
0354 msc: for GenericIon SubType= 10
0355 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0356 UrbanMsc : Emin= 0 eV Emax= 10 GeV
0357 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0358
0359 ionIoni: for GenericIon XStype:3 SubType=2
0360 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0361 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0362 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.02
0363 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0364 BraggIon : Emin= 0 eV Emax= 2 MeV
0365 BetheBloch : Emin= 2 MeV Emax= 10 GeV
0366
0367 msc: for alpha SubType= 10
0368 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0369 UrbanMsc : Emin= 0 eV Emax= 10 GeV
0370 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0371
0372 ionIoni: for alpha XStype:3 SubType=2
0373 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0374 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0375 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.02
0376 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0377 BraggIon : Emin= 0 eV Emax=7.9452 MeV
0378 BetheBloch : Emin=7.9452 MeV Emax= 10 GeV
0379
0380 msc: for anti_proton SubType= 10
0381 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0382 WentzelVIUni : Emin= 0 eV Emax= 10 GeV Nbins=80 100 eV - 10 GeV
0383 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0384
0385 hIoni: for anti_proton XStype:3 SubType=2
0386 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0387 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0388 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0389 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0390 ICRU73QO : Emin= 0 eV Emax= 2 MeV
0391 BetheBloch : Emin= 2 MeV Emax= 10 GeV
0392
0393 hBrems: for anti_proton XStype:1 SubType=3
0394 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0395 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0396 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0397 hBrem : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0398
0399 hPairProd: for anti_proton XStype:1 SubType=4
0400 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0401 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0402 Sampling table 4x1001 from 7.50618 GeV to 0.0150124 TeV
0403 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0404 hPairProd : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0405
0406 CoulombScat: for anti_proton XStype:1 SubType=1 BuildTable=1
0407 Lambda table from threshold to 10 GeV, 10 bins/decade, spline: 0
0408 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0409 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0410 eCoulombScattering : Emin= 0 eV Emax= 10 GeV
0411
0412 msc: for kaon+ SubType= 10
0413 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0414 WentzelVIUni : Emin= 0 eV Emax= 10 GeV Nbins=80 100 eV - 10 GeV
0415 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0416
0417 hIoni: for kaon+ XStype:3 SubType=2
0418 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0419 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0420 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0421 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0422 Bragg : Emin= 0 eV Emax=1.05231 MeV
0423 BetheBloch : Emin=1.05231 MeV Emax= 10 GeV
0424
0425 hBrems: for kaon+ XStype:1 SubType=3
0426 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0427 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0428 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0429 hBrem : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0430
0431 hPairProd: for kaon+ XStype:1 SubType=4
0432 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0433 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0434 Sampling table 4x1001 from 3.94942 GeV to 0.01 TeV
0435 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0436 hPairProd : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0437
0438 CoulombScat: for kaon+ XStype:1 SubType=1 BuildTable=1
0439 Lambda table from threshold to 10 GeV, 10 bins/decade, spline: 0
0440 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0441 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0442 eCoulombScattering : Emin= 0 eV Emax= 10 GeV
0443
0444 msc: for kaon- SubType= 10
0445 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0446 WentzelVIUni : Emin= 0 eV Emax= 10 GeV Nbins=80 100 eV - 10 GeV
0447 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0448
0449 hIoni: for kaon- XStype:3 SubType=2
0450 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0451 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0452 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0453 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0454 ICRU73QO : Emin= 0 eV Emax=1.05231 MeV
0455 BetheBloch : Emin=1.05231 MeV Emax= 10 GeV
0456
0457 hBrems: for kaon- XStype:1 SubType=3
0458 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0459 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0460 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0461 hBrem : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0462
0463 hPairProd: for kaon- XStype:1 SubType=4
0464 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0465 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0466 Sampling table 4x1001 from 3.94942 GeV to 0.01 TeV
0467 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0468 hPairProd : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0469
0470 CoulombScat: for kaon- XStype:1 SubType=1 BuildTable=1
0471 Used Lambda table of kaon+
0472 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0473 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0474 eCoulombScattering : Emin= 0 eV Emax= 10 GeV
0475
0476 msc: for mu+ SubType= 10
0477 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0478 WentzelVIUni : Emin= 0 eV Emax= 10 GeV Nbins=80 100 eV - 10 GeV
0479 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0480
0481 muIoni: for mu+ XStype:3 SubType=2
0482 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0483 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0484 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0485 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0486 Bragg : Emin= 0 eV Emax= 200 keV
0487 MuBetheBloch : Emin= 200 keV Emax= 10 GeV
0488
0489 muBrems: for mu+ XStype:1 SubType=3
0490 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0491 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0492 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0493 MuBrem : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0494
0495 muPairProd: for mu+ XStype:1 SubType=4
0496 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0497 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0498 Sampling table 5x1001 from 0.85 GeV to 0.01 TeV
0499 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0500 muPairProd : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0501
0502 CoulombScat: for mu+ XStype:1 SubType=1 BuildTable=1
0503 Lambda table from threshold to 10 GeV, 10 bins/decade, spline: 0
0504 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0505 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0506 eCoulombScattering : Emin= 0 eV Emax= 10 GeV
0507
0508 msc: for mu- SubType= 10
0509 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0510 WentzelVIUni : Emin= 0 eV Emax= 10 GeV Nbins=80 100 eV - 10 GeV
0511 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0512
0513 muIoni: for mu- XStype:3 SubType=2
0514 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0515 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0516 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0517 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0518 ICRU73QO : Emin= 0 eV Emax= 200 keV
0519 MuBetheBloch : Emin= 200 keV Emax= 10 GeV
0520
0521 muBrems: for mu- XStype:1 SubType=3
0522 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0523 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0524 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0525 MuBrem : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0526
0527 muPairProd: for mu- XStype:1 SubType=4
0528 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0529 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0530 Sampling table 5x1001 from 0.85 GeV to 0.01 TeV
0531 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0532 muPairProd : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0533
0534 CoulombScat: for mu- XStype:1 SubType=1 BuildTable=1
0535 Used Lambda table of mu+
0536 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0537 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0538 eCoulombScattering : Emin= 0 eV Emax= 10 GeV
0539
0540 msc: for pi+ SubType= 10
0541 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0542 WentzelVIUni : Emin= 0 eV Emax= 10 GeV Nbins=80 100 eV - 10 GeV
0543 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0544
0545 hIoni: for pi+ XStype:3 SubType=2
0546 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0547 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0548 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0549 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0550 Bragg : Emin= 0 eV Emax=297.505 keV
0551 BetheBloch : Emin=297.505 keV Emax= 10 GeV
0552
0553 hBrems: for pi+ XStype:1 SubType=3
0554 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0555 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0556 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0557 hBrem : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0558
0559 hPairProd: for pi+ XStype:1 SubType=4
0560 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0561 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0562 Sampling table 4x1001 from 1.11656 GeV to 0.01 TeV
0563 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0564 hPairProd : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0565
0566 CoulombScat: for pi+ XStype:1 SubType=1 BuildTable=1
0567 Lambda table from threshold to 10 GeV, 10 bins/decade, spline: 0
0568 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0569 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0570 eCoulombScattering : Emin= 0 eV Emax= 10 GeV
0571
0572 msc: for pi- SubType= 10
0573 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0574 WentzelVIUni : Emin= 0 eV Emax= 10 GeV Nbins=80 100 eV - 10 GeV
0575 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
0576
0577 hIoni: for pi- XStype:3 SubType=2
0578 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0579 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0580 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0581 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0582 ICRU73QO : Emin= 0 eV Emax=297.505 keV
0583 BetheBloch : Emin=297.505 keV Emax= 10 GeV
0584
0585 hBrems: for pi- XStype:1 SubType=3
0586 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0587 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0588 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0589 hBrem : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0590
0591 hPairProd: for pi- XStype:1 SubType=4
0592 dE/dx and range tables from 10 eV to 10 GeV in 90 bins
0593 Lambda tables from threshold to 10 GeV, 10 bins/decade, spline: 1
0594 Sampling table 4x1001 from 1.11656 GeV to 0.01 TeV
0595 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0596 hPairProd : Emin= 0 eV Emax= 10 GeV ModifiedMephi
0597
0598 CoulombScat: for pi- XStype:1 SubType=1 BuildTable=1
0599 Used Lambda table of pi+
0600 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0601 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0602 eCoulombScattering : Emin= 0 eV Emax= 10 GeV
0603
0604 ========= Table of registered couples ============================
0605
0606 Index : 0 used in the geometry : Yes
0607 Material : G4_AIR
0608 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 100 um
0609 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 10 keV
0610 Region(s) which use this couple :
0611 DefaultRegionForTheWorld
0612
0613 Index : 1 used in the geometry : Yes
0614 Material : G4_Be
0615 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 100 um
0616 Energy thresholds : gamma 1.0488 keV e- 108.985 keV e+ 107.379 keV proton 10 keV
0617 Region(s) which use this couple :
0618 DefaultRegionForTheWorld
0619
0620 Index : 2 used in the geometry : Yes
0621 Material : G4_He
0622 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 100 um
0623 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 10 keV
0624 Region(s) which use this couple :
0625 DefaultRegionForTheWorld
0626
0627 Index : 3 used in the geometry : Yes
0628 Material : G4_W
0629 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 100 um
0630 Energy thresholds : gamma 36.2784 keV e- 351.887 keV e+ 338.074 keV proton 10 keV
0631 Region(s) which use this couple :
0632 DefaultRegionForTheWorld
0633
0634 Index : 4 used in the geometry : Yes
0635 Material : G4_MYLAR
0636 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 100 um
0637 Energy thresholds : gamma 1.20959 keV e- 99.1194 keV e+ 97.66 keV proton 10 keV
0638 Region(s) which use this couple :
0639 DefaultRegionForTheWorld
0640
0641 Index : 5 used in the geometry : Yes
0642 Material : G4_WATER
0643 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 100 um
0644 Energy thresholds : gamma 1.11075 keV e- 85.0445 keV e+ 83.9021 keV proton 10 keV
0645 Region(s) which use this couple :
0646 DefaultRegionForTheWorld
0647
0648 ==================================================================
0649
0650 80 bins R stepR= 2.5 mm
0651 60 bins Z stepZ= 5 mm
0652 200 bins E stepE= 0.25 MeV
0653 14-th bin in Z is used for R distribution
0654 ### Run 0 starts.
0655
0656 ------- MixMaxRng engine status -------
0657 Current state vector is:
0658 mixmax state, file version 1.0
0659 N=17 V[N]={906770732717044781, 629165745432651234, 1235682547346241386, 68420008233404568, 2262190991329242458, 2266470399991071809, 1976726662926872232, 245458862506414172, 1955974201201518530, 2155248512522080758, 604170912935414061, 1116171330120743511, 1861018313684488333, 1296715403254578286, 1549011045957234151, 370819759640195970, 2230139271784837643} counter= 17sumtot= 1977567618660788324
0660 ---------------------------------------
0661 Run terminated.
0662 Run Summary
0663 Number of events processed : 10000
0664 User=3.370000s Real=3.442387s Sys=0.000000s
0665 Histo: End of run actions are started
0666 ========================================================
0667 Number of events 10000
0668 Average number of e- 33.13
0669 Average number of gamma 19.58
0670 Average number of e+ 0.7605
0671 Average number of steps in the phantom 11.51
0672 Average number of e- steps in the target 17.98
0673 Average number of g produced in the target 18.72
0674 Average number of e- produced in the target 20.83
0675 Average number of g produced in the phantom 0.4933
0676 Average number of e- produced in the phantom 5.953
0677 Total fGamma fluence in front of the phantom 7.954 MeV
0678 ========================================================
0679
0680
0681
0682 ------- MixMaxRng engine status -------
0683 Current state vector is:
0684 mixmax state, file version 1.0
0685 N=17 V[N]={961145400483405076, 699745897242935866, 1057413413739364408, 662555405695072145, 1718128862343597050, 1055151940051380955, 2185583859974491267, 673048765014179189, 467115907169885361, 358922768063147600, 1116951267450062581, 2116545342320611943, 2097059162164451181, 93621624613885263, 529728786058221144, 735650886969976136, 590625775673630549} counter= 4sumtot= 978094000532440057
0686 ---------------------------------------
0687 #
0688 ================== Deleting memory pools ===================
0689 Number of memory pools allocated: 10 of which, static: 0
0690 Dynamic pools deleted: 10 / Total memory freed: 0.058 MB
0691 ============================================================
