Geant4 Cross Reference |
1 Environment variable "G4FORCE_RUN_MANAGER_TYPE" enabled with value == Serial. Forcing G4RunManager type... 2 3 ############################################ 4 !!! WARNING - FPE detection is activated !!! 5 ############################################ 6 7 8 ################################ 9 !!! G4Backtrace is activated !!! 10 ################################ 11 12 13 ************************************************************** 14 Geant4 version Name: geant4-11-03-ref-00 (6-December-2024) 15 Copyright : Geant4 Collaboration 16 References : NIM A 506 (2003), 250-303 17 : IEEE-TNS 53 (2006), 270-278 18 : NIM A 835 (2016), 186-225 19 WWW : http://geant4.org/ 20 ************************************************************** 21 22 Accumulable registered as "accumulable_0" 23 Accumulable registered as "accumulable_1" 24 Visualization Manager instantiating with verbosity "warnings (3)"... 25 Visualization Manager initialising... 26 Registering graphics systems... 27 28 You have successfully registered the following graphics systems. 29 Registered graphics systems are: 30 ASCIITree (ATree) 31 DAWNFILE (DAWNFILE) 32 G4HepRepFile (HepRepFile) 33 RayTracer (RayTracer) 34 VRML2FILE (VRML2FILE) 35 gMocrenFile (gMocrenFile) 36 TOOLSSG_OFFSCREEN (TSG_OFFSCREEN, TSG_FILE) 37 OpenGLImmediateQt (OGLIQt, OGLI) 38 OpenGLStoredQt (OGLSQt, OGL, OGLS) 39 OpenGLImmediateXm (OGLIXm, OGLIQt_FALLBACK) 40 OpenGLStoredXm (OGLSXm, OGLSQt_FALLBACK) 41 OpenGLImmediateX (OGLIX, OGLIQt_FALLBACK, OGLIXm_FALLBACK) 42 OpenGLStoredX (OGLSX, OGLSQt_FALLBACK, OGLSXm_FALLBACK) 43 RayTracerX (RayTracerX) 44 Qt3D (Qt3D) 45 TOOLSSG_X11_GLES (TSG_X11_GLES, TSGX11, TSG_XT_GLES_FALLBACK) 46 TOOLSSG_X11_ZB (TSG_X11_ZB, TSGX11ZB) 47 TOOLSSG_XT_GLES (TSG_XT_GLES, TSGXt, TSG_QT_GLES_FALLBACK) 48 TOOLSSG_XT_ZB (TSG_XT_ZB, TSGXtZB) 49 TOOLSSG_QT_GLES (TSG_QT_GLES, TSGQt, TSG) 50 TOOLSSG_QT_ZB (TSG_QT_ZB, TSGQtZB) 51 You may choose a graphics system (driver) with a parameter of 52 the command "/vis/open" or "/vis/sceneHandler/create", 53 or you may omit the driver parameter and choose at run time: 54 - by argument in the construction of G4VisExecutive 55 - by environment variable "G4VIS_DEFAULT_DRIVER" 56 - by entry in "~/.g4session" 57 - by build flags. 58 - Note: This feature is not allowed in batch mode. 59 For further information see "examples/basic/B1/exampleB1.cc" 60 and "vis.mac". 61 62 Registering model factories... 63 64 You have successfully registered the following model factories. 65 Registered model factories: 66 generic 67 drawByAttribute 68 drawByCharge 69 drawByOriginVolume 70 drawByParticleID 71 drawByEncounteredVolume 72 73 Registered models: 74 None 75 76 Registered filter factories: 77 attributeFilter 78 chargeFilter 79 originVolumeFilter 80 particleFilter 81 encounteredVolumeFilter 82 83 Registered filters: 84 None 85 86 You have successfully registered the following user vis actions. 87 Run Duration User Vis Actions: none 88 End of Event User Vis Actions: none 89 End of Run User Vis Actions: none 90 91 Some /vis commands (optionally) take a string to specify colour. 92 "/vis/list" to see available colours. 93 Checking overlaps for volume crystal:0 (G4Box) ... OK! 94 Checking overlaps for volume crystal:1 (G4Box) ... OK! 95 Checking overlaps for volume crystal:2 (G4Box) ... OK! 96 Checking overlaps for volume crystal:3 (G4Box) ... OK! 97 Checking overlaps for volume crystal:4 (G4Box) ... OK! 98 Checking overlaps for volume crystal:5 (G4Box) ... OK! 99 Checking overlaps for volume crystal:6 (G4Box) ... OK! 100 Checking overlaps for volume crystal:7 (G4Box) ... OK! 101 Checking overlaps for volume crystal:8 (G4Box) ... OK! 102 Checking overlaps for volume crystal:9 (G4Box) ... OK! 103 Checking overlaps for volume crystal:10 (G4Box) ... OK! 104 Checking overlaps for volume crystal:11 (G4Box) ... OK! 105 Checking overlaps for volume crystal:12 (G4Box) ... OK! 106 Checking overlaps for volume crystal:13 (G4Box) ... OK! 107 Checking overlaps for volume crystal:14 (G4Box) ... OK! 108 Checking overlaps for volume crystal:15 (G4Box) ... OK! 109 Checking overlaps for volume crystal:16 (G4Box) ... OK! 110 Checking overlaps for volume crystal:17 (G4Box) ... OK! 111 Checking overlaps for volume crystal:18 (G4Box) ... OK! 112 Checking overlaps for volume crystal:19 (G4Box) ... OK! 113 Checking overlaps for volume crystal:20 (G4Box) ... OK! 114 Checking overlaps for volume crystal:21 (G4Box) ... OK! 115 Checking overlaps for volume crystal:22 (G4Box) ... OK! 116 Checking overlaps for volume crystal:23 (G4Box) ... OK! 117 Checking overlaps for volume crystal:24 (G4Box) ... OK! 118 Checking overlaps for volume crystal:25 (G4Box) ... OK! 119 Checking overlaps for volume crystal:26 (G4Box) ... OK! 120 Checking overlaps for volume crystal:27 (G4Box) ... OK! 121 Checking overlaps for volume crystal:28 (G4Box) ... OK! 122 Checking overlaps for volume crystal:29 (G4Box) ... OK! 123 Checking overlaps for volume crystal:30 (G4Box) ... OK! 124 Checking overlaps for volume crystal:31 (G4Box) ... OK! 125 Checking overlaps for volume ring:0 (G4Tubs) ... OK! 126 Checking overlaps for volume ring:1 (G4Tubs) ... OK! 127 Checking overlaps for volume ring:2 (G4Tubs) ... OK! 128 Checking overlaps for volume ring:3 (G4Tubs) ... OK! 129 Checking overlaps for volume ring:4 (G4Tubs) ... OK! 130 Checking overlaps for volume ring:5 (G4Tubs) ... OK! 131 Checking overlaps for volume ring:6 (G4Tubs) ... OK! 132 Checking overlaps for volume ring:7 (G4Tubs) ... OK! 133 Checking overlaps for volume ring:8 (G4Tubs) ... OK! 134 Checking overlaps for volume Detector:0 (G4Tubs) ... OK! 135 Checking overlaps for volume Patient:0 (G4Tubs) ... OK! 136 137 ***** Table : Nb of materials = 3 ***** 138 139 Material: Lu2SiO5 density: 7.400 g/cm3 RadL: 1.143 cm Nucl.Int.Length: 20.920 cm 140 Imean: 418.807 eV temperature: 293.15 K pressure: 1.00 atm 141 142 ---> Element: Lu (Lu) Z = 71.0 N = 175 A = 174.967 g/mole 143 ---> Isotope: Lu175 Z = 71 N = 175 A = 174.94 g/mole abundance: 97.410 % 144 ---> Isotope: Lu176 Z = 71 N = 176 A = 175.94 g/mole abundance: 2.590 % 145 ElmMassFraction: 76.40 % ElmAbundance 25.00 % 146 147 ---> Element: Si (Si) Z = 14.0 N = 28 A = 28.085 g/mole 148 ---> Isotope: Si28 Z = 14 N = 28 A = 27.98 g/mole abundance: 92.230 % 149 ---> Isotope: Si29 Z = 14 N = 29 A = 28.98 g/mole abundance: 4.683 % 150 ---> Isotope: Si30 Z = 14 N = 30 A = 29.97 g/mole abundance: 3.087 % 151 ElmMassFraction: 6.13 % ElmAbundance 12.50 % 152 153 ---> Element: O (O) Z = 8.0 N = 16 A = 15.999 g/mole 154 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 % 155 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 % 156 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 % 157 ElmMassFraction: 17.47 % ElmAbundance 62.50 % 158 159 160 Material: G4_AIR density: 1.205 mg/cm3 RadL: 303.921 m Nucl.Int.Length: 710.095 m 161 Imean: 85.700 eV temperature: 293.15 K pressure: 1.00 atm 162 163 ---> Element: C (C) Z = 6.0 N = 12 A = 12.011 g/mole 164 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.930 % 165 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.070 % 166 ElmMassFraction: 0.01 % ElmAbundance 0.02 % 167 168 ---> Element: N (N) Z = 7.0 N = 14 A = 14.007 g/mole 169 ---> Isotope: N14 Z = 7 N = 14 A = 14.00 g/mole abundance: 99.632 % 170 ---> Isotope: N15 Z = 7 N = 15 A = 15.00 g/mole abundance: 0.368 % 171 ElmMassFraction: 75.53 % ElmAbundance 78.44 % 172 173 ---> Element: O (O) Z = 8.0 N = 16 A = 15.999 g/mole 174 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 % 175 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 % 176 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 % 177 ElmMassFraction: 23.18 % ElmAbundance 21.07 % 178 179 ---> Element: Ar (Ar) Z = 18.0 N = 40 A = 39.948 g/mole 180 ---> Isotope: Ar36 Z = 18 N = 36 A = 35.97 g/mole abundance: 0.337 % 181 ---> Isotope: Ar38 Z = 18 N = 38 A = 37.96 g/mole abundance: 0.063 % 182 ---> Isotope: Ar40 Z = 18 N = 40 A = 39.96 g/mole abundance: 99.600 % 183 ElmMassFraction: 1.28 % ElmAbundance 0.47 % 184 185 186 Material: G4_BRAIN_ICRP density: 1.040 g/cm3 RadL: 35.403 cm Nucl.Int.Length: 72.156 cm 187 Imean: 73.300 eV temperature: 293.15 K pressure: 1.00 atm 188 189 ---> Element: H (H) Z = 1.0 N = 1 A = 1.008 g/mole 190 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 % 191 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 % 192 ElmMassFraction: 10.70 % ElmAbundance 64.44 % 193 194 ---> Element: C (C) Z = 6.0 N = 12 A = 12.011 g/mole 195 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.930 % 196 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.070 % 197 ElmMassFraction: 14.50 % ElmAbundance 7.33 % 198 199 ---> Element: N (N) Z = 7.0 N = 14 A = 14.007 g/mole 200 ---> Isotope: N14 Z = 7 N = 14 A = 14.00 g/mole abundance: 99.632 % 201 ---> Isotope: N15 Z = 7 N = 15 A = 15.00 g/mole abundance: 0.368 % 202 ElmMassFraction: 2.20 % ElmAbundance 0.95 % 203 204 ---> Element: O (O) Z = 8.0 N = 16 A = 15.999 g/mole 205 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 % 206 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 % 207 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 % 208 ElmMassFraction: 71.20 % ElmAbundance 27.01 % 209 210 ---> Element: Na (Na) Z = 11.0 N = 23 A = 22.990 g/mole 211 ---> Isotope: Na23 Z = 11 N = 23 A = 22.99 g/mole abundance: 100.000 % 212 ElmMassFraction: 0.20 % ElmAbundance 0.05 % 213 214 ---> Element: P (P) Z = 15.0 N = 31 A = 30.974 g/mole 215 ---> Isotope: P31 Z = 15 N = 31 A = 30.97 g/mole abundance: 100.000 % 216 ElmMassFraction: 0.40 % ElmAbundance 0.08 % 217 218 ---> Element: S (S) Z = 16.0 N = 32 A = 32.066 g/mole 219 ---> Isotope: S32 Z = 16 N = 32 A = 31.97 g/mole abundance: 94.930 % 220 ---> Isotope: S33 Z = 16 N = 33 A = 32.97 g/mole abundance: 0.760 % 221 ---> Isotope: S34 Z = 16 N = 34 A = 33.97 g/mole abundance: 4.290 % 222 ---> Isotope: S36 Z = 16 N = 36 A = 35.97 g/mole abundance: 0.020 % 223 ElmMassFraction: 0.20 % ElmAbundance 0.04 % 224 225 ---> Element: Cl (Cl) Z = 17.0 N = 35 A = 35.453 g/mole 226 ---> Isotope: Cl35 Z = 17 N = 35 A = 34.97 g/mole abundance: 75.780 % 227 ---> Isotope: Cl37 Z = 17 N = 37 A = 36.97 g/mole abundance: 24.220 % 228 ElmMassFraction: 0.30 % ElmAbundance 0.05 % 229 230 ---> Element: K (K) Z = 19.0 N = 39 A = 39.098 g/mole 231 ---> Isotope: K39 Z = 19 N = 39 A = 38.96 g/mole abundance: 93.258 % 232 ---> Isotope: K40 Z = 19 N = 40 A = 39.96 g/mole abundance: 0.012 % 233 ---> Isotope: K41 Z = 19 N = 41 A = 40.96 g/mole abundance: 6.730 % 234 ElmMassFraction: 0.30 % ElmAbundance 0.05 % 235 236 237 238 G4SDManager::AddNewCollection : the collection <crystal/edep> is registered at 1 239 G4SDManager::AddNewCollection : the collection <patient/dose> is registered at 2 240 # 241 /run/printProgress 1000 242 /run/beamOn 10000 243 ======================================================================= 244 ====== Electromagnetic Physics Parameters ======== 245 ======================================================================= 246 LPM effect enabled 1 247 Enable creation and use of sampling tables 0 248 Apply cuts on all EM processes 0 249 Use combined TransportationWithMsc Disabled 250 Use general process 1 251 Enable linear polarisation for gamma 0 252 Enable photoeffect sampling below K-shell 1 253 Enable sampling of quantum entanglement 0 254 X-section factor for integral approach 0.8 255 Min kinetic energy for tables 100 eV 256 Max kinetic energy for tables 100 TeV 257 Number of bins per decade of a table 7 258 Verbose level 1 259 Verbose level for worker thread 0 260 Bremsstrahlung energy threshold above which 261 primary e+- is added to the list of secondary 100 TeV 262 Bremsstrahlung energy threshold above which primary 263 muon/hadron is added to the list of secondary 100 TeV 264 Positron annihilation at rest model SimplePositronium 265 Enable 3 gamma annihilation on fly 0 266 Lowest triplet kinetic energy 1 MeV 267 Enable sampling of gamma linear polarisation 0 268 5D gamma conversion model type 0 269 5D gamma conversion model on isolated ion 0 270 Use Ricardo-Gerardo pair production model 0 271 Livermore data directory epics_2017 272 ======================================================================= 273 ====== Ionisation Parameters ======== 274 ======================================================================= 275 Step function for e+- (0.2, 1 mm) 276 Step function for muons/hadrons (0.2, 0.1 mm) 277 Step function for light ions (0.2, 0.1 mm) 278 Step function for general ions (0.2, 0.1 mm) 279 Lowest e+e- kinetic energy 1 keV 280 Lowest muon/hadron kinetic energy 1 keV 281 Use ICRU90 data 0 282 Fluctuations of dE/dx are enabled 1 283 Type of fluctuation model for leptons and hadrons Urban 284 Use built-in Birks satuaration 0 285 Build CSDA range enabled 0 286 Use cut as a final range enabled 0 287 Enable angular generator interface 0 288 Max kinetic energy for CSDA tables 1 GeV 289 Max kinetic energy for NIEL computation 0 eV 290 Linear loss limit 0.01 291 Read data from file for e+e- pair production by mu 0 292 ======================================================================= 293 ====== Multiple Scattering Parameters ======== 294 ======================================================================= 295 Type of msc step limit algorithm for e+- 1 296 Type of msc step limit algorithm for muons/hadrons 0 297 Msc lateral displacement for e+- enabled 1 298 Msc lateral displacement for muons and hadrons 0 299 Urban msc model lateral displacement alg96 1 300 Range factor for msc step limit for e+- 0.04 301 Range factor for msc step limit for muons/hadrons 0.2 302 Geometry factor for msc step limitation of e+- 2.5 303 Safety factor for msc step limit for e+- 0.6 304 Skin parameter for msc step limitation of e+- 1 305 Lambda limit for msc step limit for e+- 1 mm 306 Use Mott correction for e- scattering 0 307 Factor used for dynamic computation of angular 308 limit between single and multiple scattering 1 309 Fixed angular limit between single 310 and multiple scattering 3.1416 rad 311 Upper energy limit for e+- multiple scattering 100 MeV 312 Type of electron single scattering model 0 313 Type of nuclear form-factor 1 314 Screening factor 1 315 ======================================================================= 316 ====== Atomic Deexcitation Parameters ======== 317 ======================================================================= 318 Fluorescence enabled 1 319 Directory in G4LEDATA for fluorescence data files fluor 320 Auger electron cascade enabled 1 321 PIXE atomic de-excitation enabled 0 322 De-excitation module ignores cuts 1 323 Type of PIXE cross section for hadrons Empirical 324 Type of PIXE cross section for e+- Livermore 325 ======================================================================= 326 327 ### === Deexcitation model UAtomDeexcitation is activated for 1 region: 328 DefaultRegionForTheWorld 1 1 0 329 ### === Auger flag: 1 330 ### === Ignore cuts flag: 1 331 332 phot: for gamma SubType=12 BuildTable=0 333 LambdaPrime table from 200 keV to 100 TeV in 61 bins 334 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 335 LivermorePhElectric : Emin= 0 eV Emax= 100 TeV SauterGavrila Fluo 336 337 compt: for gamma SubType=13 BuildTable=1 338 Lambda table from 100 eV to 1 MeV, 7 bins/decade, spline: 1 339 LambdaPrime table from 1 MeV to 100 TeV in 56 bins 340 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 341 Klein-Nishina : Emin= 0 eV Emax= 100 TeV 342 343 conv: for gamma SubType=14 BuildTable=1 344 Lambda table from 1.022 MeV to 100 TeV, 18 bins/decade, spline: 1 345 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 346 BetheHeitlerLPM : Emin= 0 eV Emax= 100 TeV ModifiedTsai 347 348 Rayl: for gamma SubType=11 BuildTable=1 349 Lambda table from 100 eV to 150 keV, 7 bins/decade, spline: 0 350 LambdaPrime table from 150 keV to 100 TeV in 62 bins 351 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 352 LivermoreRayleigh : Emin= 0 eV Emax= 100 TeV CullenGenerator 353 354 msc: for e- SubType= 10 355 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 356 UrbanMsc : Emin= 0 eV Emax= 100 MeV Nbins=42 100 eV - 100 MeV 357 StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm 358 WentzelVIUni : Emin= 100 MeV Emax= 100 TeV Nbins=42 100 MeV - 100 TeV 359 StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm 360 361 eIoni: for e- XStype:3 SubType=2 362 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 363 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 364 StepFunction=(0.2, 1 mm), integ: 3, fluct: 1, linLossLim= 0.01 365 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 366 MollerBhabha : Emin= 0 eV Emax= 100 TeV 367 368 eBrem: for e- XStype:4 SubType=3 369 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 370 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 371 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000 372 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 373 eBremSB : Emin= 0 eV Emax= 1 GeV ModifiedTsai 374 eBremLPM : Emin= 1 GeV Emax= 100 TeV ModifiedTsai 375 376 CoulombScat: for e- XStype:1 SubType=1 BuildTable=1 377 Lambda table from 100 MeV to 100 TeV, 7 bins/decade, spline: 0 378 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 379 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 380 eCoulombScattering : Emin= 100 MeV Emax= 100 TeV 381 382 msc: for e+ SubType= 10 383 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 384 UrbanMsc : Emin= 0 eV Emax= 100 MeV Nbins=42 100 eV - 100 MeV 385 StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm 386 WentzelVIUni : Emin= 100 MeV Emax= 100 TeV Nbins=42 100 MeV - 100 TeV 387 StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm 388 389 eIoni: for e+ XStype:3 SubType=2 390 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 391 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 392 StepFunction=(0.2, 1 mm), integ: 3, fluct: 1, linLossLim= 0.01 393 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 394 MollerBhabha : Emin= 0 eV Emax= 100 TeV 395 396 eBrem: for e+ XStype:4 SubType=3 397 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 398 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 399 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000 400 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 401 eBremSB : Emin= 0 eV Emax= 1 GeV ModifiedTsai 402 eBremLPM : Emin= 1 GeV Emax= 100 TeV ModifiedTsai 403 404 annihil: for e+ XStype:2 SubType=5 AtRestModel:Simple BuildTable=0 405 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 406 eplus2gg : Emin= 0 eV Emax= 100 TeV 407 408 CoulombScat: for e+ XStype:1 SubType=1 BuildTable=1 409 Lambda table from 100 MeV to 100 TeV, 7 bins/decade, spline: 0 410 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 411 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 412 eCoulombScattering : Emin= 100 MeV Emax= 100 TeV 413 414 msc: for proton SubType= 10 415 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 416 WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV 417 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm 418 419 hIoni: for proton XStype:3 SubType=2 420 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 421 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 422 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01 423 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 424 Bragg : Emin= 0 eV Emax= 2 MeV 425 BetheBloch : Emin= 2 MeV Emax= 100 TeV 426 427 hBrems: for proton XStype:1 SubType=3 428 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 429 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 430 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 431 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi 432 433 hPairProd: for proton XStype:1 SubType=4 434 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 435 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 436 Sampling table 17x1001 from 7.50618 GeV to 100 TeV 437 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 438 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi 439 440 CoulombScat: for proton XStype:1 SubType=1 BuildTable=1 441 Lambda table from threshold to 100 TeV, 7 bins/decade, spline: 0 442 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 443 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 444 eCoulombScattering : Emin= 0 eV Emax= 100 TeV 445 446 msc: for GenericIon SubType= 10 447 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 448 UrbanMsc : Emin= 0 eV Emax= 100 TeV 449 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm 450 451 ionIoni: for GenericIon XStype:3 SubType=2 452 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 453 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 454 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.02 455 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 456 Bragg : Emin= 0 eV Emax= 2 MeV 457 BetheBloch : Emin= 2 MeV Emax= 100 TeV 458 ====================================================================== 459 ====== Radioactive Decay Physics Parameters ======= 460 ====================================================================== 461 min MeanLife (from G4NuclideTable) 1 ns 462 Max life time (from G4DeexPrecoParameters) 1000 ps 463 Internal e- conversion flag 1 464 Stored internal conversion coefficients 1 465 Enabled atomic relaxation mode 1 466 Enable correlated gamma emission 0 467 Max 2J for sampling of angular correlations 10 468 Atomic de-excitation enabled 1 469 Auger electron emission enabled 1 470 Check EM cuts disabled for atomic de-excitation 1 471 Use Bearden atomic level energies 0 472 Use ANSTO fluorescence model 0 473 Threshold for very long decay time at rest 1 y 474 ====================================================================== 475 476 ==================================================================== 477 HADRONIC PROCESSES SUMMARY (verbose level 1) 478 ------------------------------------------------------------------------- 479 Hadronic Processes for GenericIon 480 Process: RadioactiveDecay 481 482 msc: for alpha SubType= 10 483 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 484 UrbanMsc : Emin= 0 eV Emax= 100 TeV 485 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm 486 487 ionIoni: for alpha XStype:3 SubType=2 488 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 489 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 490 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.02 491 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 492 BraggIon : Emin= 0 eV Emax=7.9452 MeV 493 BetheBloch : Emin=7.9452 MeV Emax= 100 TeV 494 495 msc: for anti_proton SubType= 10 496 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 497 WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV 498 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm 499 500 hIoni: for anti_proton XStype:3 SubType=2 501 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 502 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 503 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01 504 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 505 ICRU73QO : Emin= 0 eV Emax= 2 MeV 506 BetheBloch : Emin= 2 MeV Emax= 100 TeV 507 508 hBrems: for anti_proton XStype:1 SubType=3 509 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 510 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 511 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 512 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi 513 514 hPairProd: for anti_proton XStype:1 SubType=4 515 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 516 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 517 Sampling table 17x1001 from 7.50618 GeV to 100 TeV 518 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 519 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi 520 521 CoulombScat: for anti_proton XStype:1 SubType=1 BuildTable=1 522 Lambda table from threshold to 100 TeV, 7 bins/decade, spline: 0 523 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 524 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 525 eCoulombScattering : Emin= 0 eV Emax= 100 TeV 526 527 msc: for kaon+ SubType= 10 528 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 529 WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV 530 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm 531 532 hIoni: for kaon+ XStype:3 SubType=2 533 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 534 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 535 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01 536 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 537 Bragg : Emin= 0 eV Emax=1.05231 MeV 538 BetheBloch : Emin=1.05231 MeV Emax= 100 TeV 539 540 hBrems: for kaon+ XStype:1 SubType=3 541 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 542 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 543 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 544 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi 545 546 hPairProd: for kaon+ XStype:1 SubType=4 547 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 548 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 549 Sampling table 18x1001 from 3.94942 GeV to 100 TeV 550 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 551 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi 552 553 CoulombScat: for kaon+ XStype:1 SubType=1 BuildTable=1 554 Lambda table from threshold to 100 TeV, 7 bins/decade, spline: 0 555 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 556 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 557 eCoulombScattering : Emin= 0 eV Emax= 100 TeV 558 559 msc: for kaon- SubType= 10 560 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 561 WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV 562 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm 563 564 hIoni: for kaon- XStype:3 SubType=2 565 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 566 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 567 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01 568 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 569 ICRU73QO : Emin= 0 eV Emax=1.05231 MeV 570 BetheBloch : Emin=1.05231 MeV Emax= 100 TeV 571 572 hBrems: for kaon- XStype:1 SubType=3 573 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 574 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 575 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 576 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi 577 578 hPairProd: for kaon- XStype:1 SubType=4 579 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 580 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 581 Sampling table 18x1001 from 3.94942 GeV to 100 TeV 582 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 583 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi 584 585 CoulombScat: for kaon- XStype:1 SubType=1 BuildTable=1 586 Used Lambda table of kaon+ 587 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 588 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 589 eCoulombScattering : Emin= 0 eV Emax= 100 TeV 590 591 msc: for mu+ SubType= 10 592 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 593 WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV 594 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm 595 596 muIoni: for mu+ XStype:3 SubType=2 597 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 598 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 599 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01 600 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 601 Bragg : Emin= 0 eV Emax= 200 keV 602 MuBetheBloch : Emin= 200 keV Emax= 100 TeV 603 604 muBrems: for mu+ XStype:1 SubType=3 605 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 606 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 607 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 608 MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi 609 610 muPairProd: for mu+ XStype:1 SubType=4 611 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 612 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 613 Sampling table 21x1001 from 0.85 GeV to 100 TeV 614 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 615 muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi 616 617 CoulombScat: for mu+ XStype:1 SubType=1 BuildTable=1 618 Lambda table from threshold to 100 TeV, 7 bins/decade, spline: 0 619 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 620 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 621 eCoulombScattering : Emin= 0 eV Emax= 100 TeV 622 623 msc: for mu- SubType= 10 624 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 625 WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV 626 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm 627 628 muIoni: for mu- XStype:3 SubType=2 629 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 630 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 631 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01 632 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 633 ICRU73QO : Emin= 0 eV Emax= 200 keV 634 MuBetheBloch : Emin= 200 keV Emax= 100 TeV 635 636 muBrems: for mu- XStype:1 SubType=3 637 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 638 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 639 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 640 MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi 641 642 muPairProd: for mu- XStype:1 SubType=4 643 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 644 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 645 Sampling table 21x1001 from 0.85 GeV to 100 TeV 646 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 647 muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi 648 649 CoulombScat: for mu- XStype:1 SubType=1 BuildTable=1 650 Used Lambda table of mu+ 651 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 652 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 653 eCoulombScattering : Emin= 0 eV Emax= 100 TeV 654 655 msc: for pi+ SubType= 10 656 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 657 WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV 658 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm 659 660 hIoni: for pi+ XStype:3 SubType=2 661 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 662 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 663 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01 664 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 665 Bragg : Emin= 0 eV Emax=297.505 keV 666 BetheBloch : Emin=297.505 keV Emax= 100 TeV 667 668 hBrems: for pi+ XStype:1 SubType=3 669 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 670 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 671 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 672 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi 673 674 hPairProd: for pi+ XStype:1 SubType=4 675 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 676 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 677 Sampling table 20x1001 from 1.11656 GeV to 100 TeV 678 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 679 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi 680 681 CoulombScat: for pi+ XStype:1 SubType=1 BuildTable=1 682 Lambda table from threshold to 100 TeV, 7 bins/decade, spline: 0 683 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 684 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 685 eCoulombScattering : Emin= 0 eV Emax= 100 TeV 686 687 msc: for pi- SubType= 10 688 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 689 WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV 690 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm 691 692 hIoni: for pi- XStype:3 SubType=2 693 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 694 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 695 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01 696 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 697 ICRU73QO : Emin= 0 eV Emax=297.505 keV 698 BetheBloch : Emin=297.505 keV Emax= 100 TeV 699 700 hBrems: for pi- XStype:1 SubType=3 701 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 702 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 703 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 704 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi 705 706 hPairProd: for pi- XStype:1 SubType=4 707 dE/dx and range tables from 100 eV to 100 TeV in 84 bins 708 Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 709 Sampling table 20x1001 from 1.11656 GeV to 100 TeV 710 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 711 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi 712 713 CoulombScat: for pi- XStype:1 SubType=1 BuildTable=1 714 Used Lambda table of pi+ 715 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 716 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 717 eCoulombScattering : Emin= 0 eV Emax= 100 TeV 718 ### Run 0 starts. 719 ### Run 0 start. 720 721 -------- WWWW ------- G4Exception-START -------- WWWW ------- 722 *** G4Exception : Analysis_W001 723 issued by : G4RootNtupleFileManager::SetNtupleMergingMode 724 Merging ntuples is not applicable in sequential application. 725 Setting was ignored. 726 *** This is just a warning message. *** 727 -------- WWWW -------- G4Exception-END --------- WWWW ------- 728 729 ... set ntuple merging row mode : row-wise - done 730 ... create file : scoring.root - done 731 ... open analysis file : scoring.root - done 732 ... open analysis file : scoring.root - done 733 --> Event 0 starts. 734 --> Event 1000 starts. 735 --> Event 2000 starts. 736 --> Event 3000 starts. 737 --> Event 4000 starts. 738 --> Event 5000 starts. 739 --> Event 6000 starts. 740 --> Event 7000 starts. 741 --> Event 8000 starts. 742 --> Event 9000 starts. 743 744 --------------------End of Global Run----------------------- 745 The run was 10000 events Nb of 'good' e+ annihilations: 1306 746 Total dose in patient : 306.658 picoGy 747 ------------------------------------------------------------ 748 749 ... write file : scoring.root - done 750 ... close file : scoring.root - done 751 Graphics systems deleted. 752 Visualization Manager deleting...