Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/human_phantom/human_phantom.out

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  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 Visualization Manager instantiating with verbosity "warnings (3)"...
 23 Visualization Manager initialising...
 24 Registering graphics systems...
 25 
 26 You have successfully registered the following graphics systems.
 27 Registered graphics systems are:
 28   ASCIITree (ATree)
 29   DAWNFILE (DAWNFILE)
 30   G4HepRepFile (HepRepFile)
 31   RayTracer (RayTracer)
 32   VRML2FILE (VRML2FILE)
 33   gMocrenFile (gMocrenFile)
 34   TOOLSSG_OFFSCREEN (TSG_OFFSCREEN, TSG_FILE)
 35   OpenGLImmediateQt (OGLIQt, OGLI)
 36   OpenGLStoredQt (OGLSQt, OGL, OGLS)
 37   OpenGLImmediateXm (OGLIXm, OGLIQt_FALLBACK)
 38   OpenGLStoredXm (OGLSXm, OGLSQt_FALLBACK)
 39   OpenGLImmediateX (OGLIX, OGLIQt_FALLBACK, OGLIXm_FALLBACK)
 40   OpenGLStoredX (OGLSX, OGLSQt_FALLBACK, OGLSXm_FALLBACK)
 41   RayTracerX (RayTracerX)
 42   Qt3D (Qt3D)
 43   TOOLSSG_X11_GLES (TSG_X11_GLES, TSGX11, TSG_XT_GLES_FALLBACK)
 44   TOOLSSG_X11_ZB (TSG_X11_ZB, TSGX11ZB)
 45   TOOLSSG_XT_GLES (TSG_XT_GLES, TSGXt, TSG_QT_GLES_FALLBACK)
 46   TOOLSSG_XT_ZB (TSG_XT_ZB, TSGXtZB)
 47   TOOLSSG_QT_GLES (TSG_QT_GLES, TSGQt, TSG)
 48   TOOLSSG_QT_ZB (TSG_QT_ZB, TSGQtZB)
 49 You may choose a graphics system (driver) with a parameter of
 50 the command "/vis/open" or "/vis/sceneHandler/create",
 51 or you may omit the driver parameter and choose at run time:
 52 - by argument in the construction of G4VisExecutive
 53 - by environment variable "G4VIS_DEFAULT_DRIVER"
 54 - by entry in "~/.g4session"
 55 - by build flags.
 56 - Note: This feature is not allowed in batch mode.
 57 For further information see "examples/basic/B1/exampleB1.cc"
 58 and "vis.mac".
 59 
 60 Registering model factories...
 61 
 62 You have successfully registered the following model factories.
 63 Registered model factories:
 64   generic
 65   drawByAttribute
 66   drawByCharge
 67   drawByOriginVolume
 68   drawByParticleID
 69   drawByEncounteredVolume
 70 
 71 Registered models:
 72   None
 73 
 74 Registered filter factories:
 75   attributeFilter
 76   chargeFilter
 77   originVolumeFilter
 78   particleFilter
 79   encounteredVolumeFilter
 80 
 81 Registered filters:
 82   None
 83 
 84 You have successfully registered the following user vis actions.
 85 Run Duration User Vis Actions: none
 86 End of Event User Vis Actions: none
 87 End of Run User Vis Actions: none
 88 
 89 Some /vis commands (optionally) take a string to specify colour.
 90 "/vis/list" to see available colours.
 91 #
 92 #
 93 #
 94 # Initialize New Phantom
 95 #
 96 # Choose model: ORNLFemale, ORNLMale, MIRD, MIRDHead, ORNLHead
 97 #
 98 /phantom/setPhantomModel MIRD
 99  >> Phantom MIRD will be built.
100 #
101 # Choose Sex of Phantom: Male or Female
102 #
103 /phantom/setPhantomSex Male
104 >> Male Phantom will be built.
105 # Construct your Phantom
106 #
107 /phantom/buildNewPhantom
108  ****************>>>> NEW PHANTOM CONSTRUCTION <<<<***************** 
109 /run/initialize
110 MotherVolume: physicalWorld
111 sensitivity : 0
112 Construct Head with mother physicalWorld
113 Checking overlaps for volume physicalHead:0 (G4UnionSolid) ... OK! 
114 Volume of Head = 4654.79 cm^3
115 Material of Head = soft_tissue
116 Density of Material = 0.9869 g/cm^3
117 Mass of Head = 4593.82 g
118 MotherVolume: physicalHead
119 sensitivity : 0
120 Construct Skull with mother volume physicalHead
121 Checking overlaps for volume physicalSkull:0 (G4SubtractionSolid) ... OK! 
122 Skull created !!!!!!
123 Volume of Skull = 846.921 cm^3
124 Material of Skull = skeleton
125 Density of Material = 1.4862 g/cm^3
126 Mass of Skull = 1258.69 g
127 Construct Brain with mother physicalHead
128 Checking overlaps for volume physicalBrain:0 (G4Ellipsoid) ... OK! 
129 Volume of Brain = 1470.27 cm^3
130 Material of Brain = soft_tissue
131 Density of Material = 0.9869 g/cm^3
132 Mass of Brain = 1451 g
133 MotherVolume: physicalWorld
134 sensitivity : 0
135 Construct Trunk with mother volume physicalWorld
136 Checking overlaps for volume physicalTrunk:0 (G4EllipticalTube) ... OK! 
137 Trunk created !!!!!!
138 Volume of Trunk = 43982.3 cm^3
139 Material of Trunk = soft_tissue
140 Density of Material = 0.9869 g/cm^3
141 Mass of Trunk = 43406.1 g
142 MotherVolume: physicalWorld
143 sensitivity : 0
144 Construct LeftLeg with mother volume physicalWorld
145 Checking overlaps for volume physicalLeftLeg:0 (G4Cons) ... OK! 
146 LeftLeg created !!!!!!
147 Volume of LeftLeg = 10388.2 cm^3
148 Material of LeftLeg = soft_tissue
149 Density of Material = 0.9869 g/cm^3
150 Mass of LeftLeg = 10252.1 g
151 MotherVolume: physicalWorld
152 sensitivity : 0
153 Construct RightLeg with mother volume physicalWorld
154 Checking overlaps for volume physicalRightLeg:0 (G4Cons) ... OK! 
155 RightLeg created !!!!!!
156 Volume of RightLeg = 10388.2 cm^3
157 Material of RightLeg = soft_tissue
158 Density of Material = 0.9869 g/cm^3
159 Mass of RightLeg = 10252.1 g
160 MotherVolume: physicalTrunk
161 sensitivity : 0
162 Construct LeftArmBone with mother physicalTrunk
163 Checking overlaps for volume physicalLeftArmBone:0 (G4EllipticalTube) ... OK! 
164 LeftArmBone created !!!!!!
165 Volume of LeftArmBone = 819.39 cm^3
166 Material of LeftArmBone = skeleton
167 Density of Material = 1.4862 g/cm^3
168 Mass of LeftArmBone = 1217.78 g
169 MotherVolume: physicalTrunk
170 sensitivity : 0
171 Construct RightArmBone with mother volume physicalTrunk
172 Checking overlaps for volume physicalRightArmBone:0 (G4EllipticalTube) ... OK! 
173 RightArmBone created !!!!!!
174 Volume of RightArmBone = 819.39 cm^3
175 Material of RightArmBone = skeleton
176 Density of Material = 1.4862 g/cm^3
177 Mass of RightArmBone = 1217.78 g
178 MotherVolume: physicalLeftLeg
179 sensitivity : 0
180 Construct LeftLegBone with mother volume physicalLeftLeg
181 Checking overlaps for volume physicalLeftLegBone:0 (G4Cons) ... OK! 
182 LeftLegBone created !!!!!!
183 Volume of LeftLegBone = 1399.74 cm^3
184 Material of LeftLegBone = skeleton
185 Density of Material = 1.4862 g/cm^3
186 Mass of LeftLegBone = 2080.29 g
187 MotherVolume: physicalRightLeg
188 sensitivity : 0
189 Construct RightLegBone with mother volume physicalRightLeg
190 Checking overlaps for volume physicalRightLegBone:0 (G4Cons) ... OK! 
191 RightLegBone created !!!!!!
192 Volume of RightLegBone = 1399.74 cm^3
193 Material of RightLegBone = skeleton
194 Density of Material = 1.4862 g/cm^3
195 Mass of RightLegBone = 2080.29 g
196 MotherVolume: physicalHead
197 sensitivity : 0
198 Construct UpperSpine with mother volume physicalHead
199 Checking overlaps for volume physicalUpperSpine:0 (G4SubtractionSolid) ... OK! 
200 UpperSpine created !!!!!!
201 Volume of UpperSpine = 126.45 cm^3
202 Material of UpperSpine = skeleton
203 Density of Material = 1.4862 g/cm^3
204 Mass of UpperSpine = 187.93 g
205 MotherVolume: physicalTrunk
206 sensitivity : 0
207 Construct LeftScapula with mother volume physicalTrunk
208 Checking overlaps for volume physicalLeftScapula:0 (G4SubtractionSolid) ... OK! 
209 LeftScapula created !!!!!!
210 Volume of LeftScapula = 91.848 cm^3
211 Material of LeftScapula = skeleton
212 Density of Material = 1.4862 g/cm^3
213 Mass of LeftScapula = 136.505 g
214 MotherVolume: physicalTrunk
215 sensitivity : 0
216 Construct RightScapula with mother volume physicalTrunk
217 Checking overlaps for volume physicalRightScapula:0 (G4SubtractionSolid) ... OK! 
218 RightScapula created !!!!!!
219 Volume of RightScapula = 106.747 cm^3
220 Material of RightScapula = skeleton
221 Density of Material = 1.4862 g/cm^3
222 Mass of RightScapula = 158.647 g
223 Construct LeftAdrenal with mother physicalTrunk
224 Checking overlaps for volume physicalLeftAdrenal:0 (G4Ellipsoid) ... OK! 
225 Left LeftAdrenal created !!!!!!
226 Volume of LeftAdrenal = 7.85398 cm^3
227 Material of LeftAdrenal = soft_tissue
228 Density of Material = 0.9869 g/cm^3
229 Mass of LeftAdrenal = 7.75109 g
230 Construct RightAdrenal with mother volume physicalTrunk
231 Checking overlaps for volume physicalRightAdrenal:0 (G4Ellipsoid) ... OK! 
232 Right RightAdrenal created !!!!!!
233 Volume of RightAdrenal = 7.85398 cm^3
234 Material of RightAdrenal = soft_tissue
235 Density of Material = 0.9869 g/cm^3
236 Mass of RightAdrenal = 7.75109 g
237 Construct Thymus with mother volume physicalTrunk
238 Checking overlaps for volume physicalThymus:0 (G4Ellipsoid) ... OK! 
239 Thymus created !!!!!!
240 Volume of Thymus = 25.1327 cm^3
241 Material of Thymus = soft_tissue
242 Density of Material = 0.9869 g/cm^3
243 Mass of Thymus = 24.8035 g
244 MotherVolume: physicalTrunk
245 sensitivity : 0
246 Construct LeftClavicle with mother volume physicalTrunk
247 Checking overlaps for volume physicalLeftClavicle:0 (G4Torus) ... OK! 
248 LeftClavicle created !!!!!!
249 Volume of LeftClavicle = 13.6657 cm^3
250 Material of LeftClavicle = skeleton
251 Density of Material = 1.4862 g/cm^3
252 Mass of LeftClavicle = 20.3099 g
253 MotherVolume: physicalTrunk
254 sensitivity : 0
255 Construct RightClavicle with mother volume physicalTrunk
256 Checking overlaps for volume physicalRightClavicle:0 (G4Torus) ... OK! 
257 RightClavicle created !!!!!!
258 Volume of RightClavicle = 13.6657 cm^3
259 Material of RightClavicle = skeleton
260 Density of Material = 1.4862 g/cm^3
261 Mass of RightClavicle = 20.3099 g
262 Construct SmallIntestine with mother volume physicalTrunk
263 Checking overlaps for volume physicalSmallIntestine:0 (G4SubtractionSolid) ... OK! 
264 SmallIntestine created !!!!!!
265 Volume of SmallIntestine = 1020.2 cm^3
266 Material of SmallIntestine = soft_tissue
267 Density of Material = 0.9869 g/cm^3
268 Mass of SmallIntestine = 1006.83 g
269 Construct RibCage with mother volume physicalTrunk
270 Checking overlaps for volume physicalRibCage:0 (G4SubtractionSolid) ... OK! 
271 Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! 
272 Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! 
273 Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! 
274 Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! 
275 Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! 
276 Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! 
277 Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! 
278 Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! 
279 Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! 
280 Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! 
281 Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! 
282 Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! 
283 RibCage created !!!!!!
284 Volume of RibCage = 692.88 cm^3
285 Material of RibCage = skeleton
286 Density of Material = 1.4862 g/cm^3
287 Mass of RibCage = 1029.76 g
288 MotherVolume: physicalTrunk
289 sensitivity : 0
290 Construct MiddleLowerSpine with mother volume physicalTrunk
291 Checking overlaps for volume physicalMiddleLowerSpine:0 (G4EllipticalTube) ... OK! 
292 MiddleLowerSpine created !!!!!!
293 Volume of MiddleLowerSpine = 753.982 cm^3
294 Material of MiddleLowerSpine = skeleton
295 Density of Material = 1.4862 g/cm^3
296 Mass of MiddleLowerSpine = 1120.57 g
297 Construct Pelvis with mother volume physicalTrunk
298 Checking overlaps for volume physicalPelvis:0 (G4SubtractionSolid) ... OK! 
299 Pelvis created !!!!!!
300 Volume of Pelvis = 606.07 cm^3
301 Material of Pelvis = skeleton
302 Density of Material = 1.4862 g/cm^3
303 Mass of Pelvis = 900.741 g
304 Construct Stomach with mother volume physicalTrunk
305 Checking overlaps for volume physicalStomach:0 (G4Ellipsoid) ... OK! 
306 Stomach created !!!!!!
307 Volume of Stomach = 402.124 cm^3
308 Material of Stomach = soft_tissue
309 Density of Material = 0.9869 g/cm^3
310 Mass of Stomach = 396.856 g
311 Construct UpperLargeIntestine with mother volume physicalTrunk
312 Checking overlaps for volume physicalUpperLargeIntestine:0 (G4UnionSolid) ... OK! 
313 UpperLargeIntestine created !!!!!!
314 Volume of UpperLargeIntestine = 434.914 cm^3
315 Material of UpperLargeIntestine = soft_tissue
316 Density of Material = 0.9869 g/cm^3
317 Mass of UpperLargeIntestine = 429.217 g
318 Construct LowerLargeIntestine with mother volume physicalTrunk
319 Checking overlaps for volume physicalLowerLargeIntestine:0 (G4UnionSolid) ... OK! 
320 LowerLargeIntestine created !!!!!!
321 Volume of LowerLargeIntestine = 344.316 cm^3
322 Material of LowerLargeIntestine = soft_tissue
323 Density of Material = 0.9869 g/cm^3
324 Mass of LowerLargeIntestine = 339.805 g
325 Construct Spleen with mother volume physicalTrunk
326 Checking overlaps for volume physicalSpleen:0 (G4Ellipsoid) ... OK! 
327 Spleen created !!!!!!
328 Volume of Spleen = 175.929 cm^3
329 Material of Spleen = soft_tissue
330 Density of Material = 0.9869 g/cm^3
331 Mass of Spleen = 173.625 g
332 Construct Pancreas with mother volume physicalTrunk
333 Checking overlaps for volume physicalPancreas:0 (G4SubtractionSolid) ... OK! 
334 Pancreas created !!!!!!
335 Volume of Pancreas = 61.101 cm^3
336 Material of Pancreas = soft_tissue
337 Density of Material = 0.9869 g/cm^3
338 Mass of Pancreas = 60.3006 g
339 Construct LeftKidney with mother volume physicalTrunk
340 Checking overlaps for volume physicalLeftKidney:0 (G4SubtractionSolid) ... OK! 
341 Left LeftKidney created !!!!!!
342 Volume of LeftKidney = 144.001 cm^3
343 Material of LeftKidney = soft_tissue
344 Density of Material = 0.9869 g/cm^3
345 Mass of LeftKidney = 142.115 g
346 Construct RightKidney with mother volume physicalTrunk
347 Checking overlaps for volume physicalRightKidney:0 (G4SubtractionSolid) ... OK! 
348 RightKidney created !!!!!!
349 Volume of RightKidney = 143.986 cm^3
350 Material of RightKidney = soft_tissue
351 Density of Material = 0.9869 g/cm^3
352 Mass of RightKidney = 142.1 g
353 Construct UrinaryBladder with mother volume physicalTrunk
354 Checking overlaps for volume physicalUrinaryBladder:0 (G4SubtractionSolid) ... OK! 
355 UrinaryBladder created !!!!!!
356 Volume of UrinaryBladder = 45.7615 cm^3
357 Material of UrinaryBladder = soft_tissue
358 Density of Material = 0.9869 g/cm^3
359 Mass of UrinaryBladder = 45.1621 g
360 MotherVolume: physicalWorld
361 sensitivity : 0
362 Construct MaleGenitalia with mother volume physicalWorld
363 Checking overlaps for volume physicalMaleGenitalia:0 (G4SubtractionSolid) ... OK! 
364 MaleGenitalia created !!!!!!
365 Volume of MaleGenitalia = 229.106 cm^3
366 Material of MaleGenitalia = soft_tissue
367 Density of Material = 0.9869 g/cm^3
368 Mass of MaleGenitalia = 226.105 g
369 MotherVolume: physicalWorld
370 sensitivity : 0
371 Construct LeftTeste with mother volume physicalMaleGenitalia
372 Checking overlaps for volume physicalLeftTeste:0 (G4Ellipsoid) ... OK! 
373 LeftTeste created !!!!!!
374 Volume of LeftTeste = 18.7867 cm^3
375 Material of LeftTeste = soft_tissue
376 Density of Material = 0.9869 g/cm^3
377 Mass of LeftTeste = 18.5406 g
378 MotherVolume: physicalWorld
379 sensitivity : 0
380 Construct RightTeste with mother volume physicalMaleGenitalia
381 Checking overlaps for volume physicalRightTeste:0 (G4Ellipsoid) ... OK! 
382 RightTeste created !!!!!!
383 Volume of RightTeste = 18.7867 cm^3
384 Material of RightTeste = soft_tissue
385 Density of Material = 0.9869 g/cm^3
386 Mass of RightTeste = 18.5406 g
387 # or ...
388 #/control/execute adultMIRDFemale.mac
389 # or ...
390 #/control/execute adultORNLMale.mac
391 # or...
392 #/control/execute adultORNLFemale.mac
393 #
394 # Primary radiation field
395 /control/execute primary.mac
396 #
397 ### Radiation field defined
398 ### by means of GPS 
399 /gps/verbose 0
400 /gps/particle e-
401 /gps/energy 1 MeV
402 #
403 /gps/pos/type Volume
404 #
405 /gps/pos/shape Para
406 /gps/pos/halfx 20 cm
407 /gps/pos/halfy 1 mm
408 /gps/pos/halfz 100 cm
409 /gps/pos/centre 0. 30. 0. cm
410 /gps/direction 0 -1 0
411 #
412 /run/beamOn 100
413 =======================================================================
414 ======                 Electromagnetic Physics Parameters      ========
415 =======================================================================
416 LPM effect enabled                                 1
417 Enable creation and use of sampling tables         0
418 Apply cuts on all EM processes                     0
419 Use combined TransportationWithMsc                 Disabled
420 Use general process                                1
421 Enable linear polarisation for gamma               0
422 Enable photoeffect sampling below K-shell          1
423 Enable sampling of quantum entanglement            0
424 X-section factor for integral approach             0.8
425 Min kinetic energy for tables                      100 eV 
426 Max kinetic energy for tables                      100 TeV
427 Number of bins per decade of a table               20
428 Verbose level                                      1
429 Verbose level for worker thread                    0
430 Bremsstrahlung energy threshold above which 
431   primary e+- is added to the list of secondary    100 TeV
432 Bremsstrahlung energy threshold above which primary
433   muon/hadron is added to the list of secondary    100 TeV
434 Positron annihilation at rest model                AllisonPositronium
435 Enable 3 gamma annihilation on fly                 1
436 Lowest triplet kinetic energy                      1 MeV
437 Enable sampling of gamma linear polarisation       0
438 5D gamma conversion model type                     0
439 5D gamma conversion model on isolated ion          0
440 Use Ricardo-Gerardo pair production model          0
441 Livermore data directory                           epics_2017
442 =======================================================================
443 ======                 Ionisation Parameters                   ========
444 =======================================================================
445 Step function for e+-                              (0.2, 0.01 mm)
446 Step function for muons/hadrons                    (0.1, 0.05 mm)
447 Step function for light ions                       (0.1, 0.02 mm)
448 Step function for general ions                     (0.1, 0.001 mm)
449 Lowest e+e- kinetic energy                         100 eV 
450 Lowest muon/hadron kinetic energy                  1 keV
451 Use ICRU90 data                                    1
452 Fluctuations of dE/dx are enabled                  1
453 Type of fluctuation model for leptons and hadrons  Urban
454 Use built-in Birks satuaration                     0
455 Build CSDA range enabled                           0
456 Use cut as a final range enabled                   0
457 Enable angular generator interface                 1
458 Max kinetic energy for CSDA tables                 1 GeV
459 Max kinetic energy for NIEL computation            1 MeV
460 Linear loss limit                                  0.01
461 Read data from file for e+e- pair production by mu 0
462 =======================================================================
463 ======                 Multiple Scattering Parameters          ========
464 =======================================================================
465 Type of msc step limit algorithm for e+-           2
466 Type of msc step limit algorithm for muons/hadrons 0
467 Msc lateral displacement for e+- enabled           1
468 Msc lateral displacement for muons and hadrons     1
469 Urban msc model lateral displacement alg96         1
470 Range factor for msc step limit for e+-            0.08
471 Range factor for msc step limit for muons/hadrons  0.2
472 Geometry factor for msc step limitation of e+-     2.5
473 Safety factor for msc step limit for e+-           0.6
474 Skin parameter for msc step limitation of e+-      3
475 Lambda limit for msc step limit for e+-            1 mm
476 Use Mott correction for e- scattering              1
477 Factor used for dynamic computation of angular 
478   limit between single and multiple scattering     1
479 Fixed angular limit between single 
480   and multiple scattering                          3.1416 rad
481 Upper energy limit for e+- multiple scattering     100 MeV
482 Type of electron single scattering model           0
483 Type of nuclear form-factor                        1
484 Screening factor                                   1
485 =======================================================================
486 ======                 Atomic Deexcitation Parameters          ========
487 =======================================================================
488 Fluorescence enabled                               1
489 Directory in G4LEDATA for fluorescence data files  fluor
490 Auger electron cascade enabled                     0
491 PIXE atomic de-excitation enabled                  0
492 De-excitation module ignores cuts                  0
493 Type of PIXE cross section for hadrons             Empirical
494 Type of PIXE cross section for e+-                 Livermore
495 =======================================================================
496 
497 ### ===  Deexcitation model UAtomDeexcitation is activated for 1 region:
498           DefaultRegionForTheWorld  1  0  0
499 ### ===  Ignore cuts flag:   0
500 
501 phot:  for gamma SubType=12 BuildTable=0
502       LambdaPrime table from 200 keV to 100 TeV in 174 bins 
503       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
504  LivermorePhElectric : Emin=    0 eV  Emax=  100 TeV  SauterGavrila Fluo
505 
506 compt:  for gamma SubType=13 BuildTable=1
507       Lambda table from 100 eV  to 1 MeV, 20 bins/decade, spline: 1
508       LambdaPrime table from 1 MeV to 100 TeV in 160 bins 
509       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
510    LowEPComptonModel : Emin=    0 eV  Emax=   20 MeV Fluo
511         KleinNishina : Emin=   20 MeV Emax=  100 TeV Fluo
512 
513 conv:  for gamma SubType=14 BuildTable=1
514       Lambda table from 1.022 MeV to 100 TeV, 20 bins/decade, spline: 1
515       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
516       BetheHeitler5D : Emin=    0 eV  Emax=  100 TeV  ModifiedTsai
517 
518 Rayl:  for gamma SubType=11 BuildTable=1
519       Lambda table from 100 eV  to 150 keV, 20 bins/decade, spline: 0
520       LambdaPrime table from 150 keV to 100 TeV in 176 bins 
521       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
522    LivermoreRayleigh : Emin=    0 eV  Emax=  100 TeV  CullenGenerator
523 
524 msc:  for e-  SubType= 10
525       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
526   GoudsmitSaunderson : Emin=    0 eV  Emax=  100 MeV Nbins=120 100 eV  - 100 MeV
527           StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
528         WentzelVIUni : Emin=  100 MeV Emax=  100 TeV Nbins=120 100 MeV - 100 TeV
529           StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
530 
531 eIoni:  for e-  XStype:3  SubType=2
532       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
533       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
534       StepFunction=(0.2, 0.01 mm), integ: 3, fluct: 1, linLossLim= 0.01
535       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
536              PenIoni : Emin=    0 eV  Emax=  100 keV
537         MollerBhabha : Emin=  100 keV Emax=  100 TeV  deltaVI
538 
539 eBrem:  for e-  XStype:4  SubType=3
540       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
541       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
542       LPM flag: 1 for E > 1 GeV,  VertexHighEnergyTh(GeV)= 100000
543       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
544              eBremSB : Emin=    0 eV  Emax=    1 GeV  AngularGen2BS
545             eBremLPM : Emin=    1 GeV Emax=  100 TeV  AngularGen2BS
546 
547 ePairProd:  for e-  XStype:1  SubType=4
548       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
549       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 0
550       Sampling table 25x1001 from 0.1 GeV to 100 TeV 
551       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
552            ePairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
553 
554 CoulombScat:  for e- XStype:1 SubType=1 BuildTable=1
555       Lambda table from 100 MeV to 100 TeV, 20 bins/decade, spline: 0
556       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
557       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
558   eCoulombScattering : Emin=  100 MeV Emax=  100 TeV
559 
560 msc:  for e+  SubType= 10
561       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
562   GoudsmitSaunderson : Emin=    0 eV  Emax=  100 MeV Nbins=120 100 eV  - 100 MeV
563           StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
564         WentzelVIUni : Emin=  100 MeV Emax=  100 TeV Nbins=120 100 MeV - 100 TeV
565           StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
566 
567 eIoni:  for e+  XStype:3  SubType=2
568       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
569       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
570       StepFunction=(0.2, 0.01 mm), integ: 3, fluct: 1, linLossLim= 0.01
571       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
572              PenIoni : Emin=    0 eV  Emax=  100 keV
573         MollerBhabha : Emin=  100 keV Emax=  100 TeV  deltaVI
574 
575 eBrem:  for e+  XStype:4  SubType=3
576       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
577       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
578       LPM flag: 1 for E > 1 GeV,  VertexHighEnergyTh(GeV)= 100000
579       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
580              eBremSB : Emin=    0 eV  Emax=    1 GeV  AngularGen2BS
581             eBremLPM : Emin=    1 GeV Emax=  100 TeV  AngularGen2BS
582 
583 ePairProd:  for e+  XStype:1  SubType=4
584       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
585       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 0
586       Sampling table 25x1001 from 0.1 GeV to 100 TeV 
587       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
588            ePairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
589 
590 annihil:  for e+ XStype:2 SubType=5 AtRestModel:Allison BuildTable=0
591       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
592     eplusTo2or3gamma : Emin=    0 eV  Emax=  100 TeV
593 
594 CoulombScat:  for e+ XStype:1 SubType=1 BuildTable=1
595       Lambda table from 100 MeV to 100 TeV, 20 bins/decade, spline: 0
596       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
597       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
598   eCoulombScattering : Emin=  100 MeV Emax=  100 TeV
599 
600 msc:  for proton  SubType= 10
601       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
602         WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=240 100 eV  - 100 TeV
603           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
604 
605 hIoni:  for proton  XStype:3  SubType=2
606       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
607       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
608       StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
609       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
610                Bragg : Emin=    0 eV  Emax=    2 MeV  deltaVI
611           BetheBloch : Emin=    2 MeV Emax=  100 TeV  deltaVI
612 
613 hBrems:  for proton  XStype:1  SubType=3
614       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
615       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
616       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
617                hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
618       ===== Limit on energy threshold has been applied 
619 
620 hPairProd:  for proton  XStype:1  SubType=4
621       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
622       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
623       Sampling table 17x1001 from 7.50618 GeV to 100 TeV 
624       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
625            hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
626 
627 CoulombScat:  for proton XStype:1 SubType=1 BuildTable=1
628       Lambda table from threshold  to 100 TeV, 20 bins/decade, spline: 0
629       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
630       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
631   eCoulombScattering : Emin=    0 eV  Emax=  100 TeV
632 
633 nuclearStopping:  for proton SubType=8 BuildTable=0
634       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
635    ICRU49NucStopping : Emin=    0 eV  Emax=    1 MeV
636 
637 msc:  for GenericIon  SubType= 10
638       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
639             UrbanMsc : Emin=    0 eV  Emax=  100 TeV
640           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
641 
642 ionIoni:  for GenericIon  XStype:3  SubType=2
643       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
644       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
645       StepFunction=(0.1, 0.001 mm), integ: 3, fluct: 1, linLossLim= 0.02
646       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
647     LindhardSorensen : Emin=    0 eV  Emax=  100 TeV  deltaVI
648 
649 nuclearStopping:  for GenericIon SubType=8 BuildTable=0
650       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
651    ICRU49NucStopping : Emin=    0 eV  Emax=    1 MeV
652 
653 msc:  for alpha  SubType= 10
654       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
655             UrbanMsc : Emin=    0 eV  Emax=  100 TeV
656           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
657 
658 ionIoni:  for alpha  XStype:3  SubType=2
659       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
660       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
661       StepFunction=(0.1, 0.02 mm), integ: 3, fluct: 1, linLossLim= 0.02
662       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
663             BraggIon : Emin=    0 eV  Emax=7.9452 MeV  deltaVI
664           BetheBloch : Emin=7.9452 MeV Emax=  100 TeV  deltaVI
665 
666 nuclearStopping:  for alpha SubType=8 BuildTable=0
667       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
668    ICRU49NucStopping : Emin=    0 eV  Emax=    1 MeV
669 
670 msc:  for anti_proton  SubType= 10
671       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
672         WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=240 100 eV  - 100 TeV
673           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
674 
675 hIoni:  for anti_proton  XStype:3  SubType=2
676       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
677       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
678       StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
679       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
680             ICRU73QO : Emin=    0 eV  Emax=    2 MeV  deltaVI
681           BetheBloch : Emin=    2 MeV Emax=  100 TeV  deltaVI
682 
683 hBrems:  for anti_proton  XStype:1  SubType=3
684       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
685       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
686       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
687                hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
688       ===== Limit on energy threshold has been applied 
689 
690 hPairProd:  for anti_proton  XStype:1  SubType=4
691       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
692       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
693       Sampling table 17x1001 from 7.50618 GeV to 100 TeV 
694       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
695            hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
696 
697 CoulombScat:  for anti_proton XStype:1 SubType=1 BuildTable=1
698       Lambda table from threshold  to 100 TeV, 20 bins/decade, spline: 0
699       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
700       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
701   eCoulombScattering : Emin=    0 eV  Emax=  100 TeV
702 
703 msc:  for kaon+  SubType= 10
704       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
705         WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=240 100 eV  - 100 TeV
706           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
707 
708 hIoni:  for kaon+  XStype:3  SubType=2
709       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
710       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
711       StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
712       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
713                Bragg : Emin=    0 eV  Emax=1.05231 MeV  deltaVI
714           BetheBloch : Emin=1.05231 MeV Emax=  100 TeV  deltaVI
715 
716 hBrems:  for kaon+  XStype:1  SubType=3
717       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
718       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
719       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
720                hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
721       ===== Limit on energy threshold has been applied 
722 
723 hPairProd:  for kaon+  XStype:1  SubType=4
724       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
725       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
726       Sampling table 18x1001 from 3.94942 GeV to 100 TeV 
727       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
728            hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
729 
730 CoulombScat:  for kaon+ XStype:1 SubType=1 BuildTable=1
731       Lambda table from threshold  to 100 TeV, 20 bins/decade, spline: 0
732       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
733       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
734   eCoulombScattering : Emin=    0 eV  Emax=  100 TeV
735 
736 msc:  for kaon-  SubType= 10
737       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
738         WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=240 100 eV  - 100 TeV
739           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
740 
741 hIoni:  for kaon-  XStype:3  SubType=2
742       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
743       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
744       StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
745       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
746             ICRU73QO : Emin=    0 eV  Emax=1.05231 MeV  deltaVI
747           BetheBloch : Emin=1.05231 MeV Emax=  100 TeV  deltaVI
748 
749 hBrems:  for kaon-  XStype:1  SubType=3
750       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
751       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
752       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
753                hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
754       ===== Limit on energy threshold has been applied 
755 
756 hPairProd:  for kaon-  XStype:1  SubType=4
757       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
758       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
759       Sampling table 18x1001 from 3.94942 GeV to 100 TeV 
760       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
761            hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
762 
763 CoulombScat:  for kaon- XStype:1 SubType=1 BuildTable=1
764       Used Lambda table of kaon+
765       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
766       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
767   eCoulombScattering : Emin=    0 eV  Emax=  100 TeV
768 
769 msc:  for mu+  SubType= 10
770       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
771         WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=240 100 eV  - 100 TeV
772           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
773 
774 muIoni:  for mu+  XStype:3  SubType=2
775       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
776       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
777       StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
778       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
779                Bragg : Emin=    0 eV  Emax=  200 keV  deltaVI
780         MuBetheBloch : Emin=  200 keV Emax=  100 TeV  deltaVI
781 
782 muBrems:  for mu+  XStype:1  SubType=3
783       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
784       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
785       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
786               MuBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
787       ===== Limit on energy threshold has been applied 
788 
789 muPairProd:  for mu+  XStype:1  SubType=4
790       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
791       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
792       Sampling table 21x1001 from 0.85 GeV to 100 TeV 
793       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
794           muPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
795 
796 CoulombScat:  for mu+ XStype:1 SubType=1 BuildTable=1
797       Lambda table from threshold  to 100 TeV, 20 bins/decade, spline: 0
798       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
799       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
800   eCoulombScattering : Emin=    0 eV  Emax=  100 TeV
801 
802 msc:  for mu-  SubType= 10
803       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
804         WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=240 100 eV  - 100 TeV
805           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
806 
807 muIoni:  for mu-  XStype:3  SubType=2
808       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
809       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
810       StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
811       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
812             ICRU73QO : Emin=    0 eV  Emax=  200 keV  deltaVI
813         MuBetheBloch : Emin=  200 keV Emax=  100 TeV  deltaVI
814 
815 muBrems:  for mu-  XStype:1  SubType=3
816       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
817       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
818       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
819               MuBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
820       ===== Limit on energy threshold has been applied 
821 
822 muPairProd:  for mu-  XStype:1  SubType=4
823       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
824       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
825       Sampling table 21x1001 from 0.85 GeV to 100 TeV 
826       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
827           muPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
828 
829 CoulombScat:  for mu- XStype:1 SubType=1 BuildTable=1
830       Used Lambda table of mu+
831       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
832       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
833   eCoulombScattering : Emin=    0 eV  Emax=  100 TeV
834 
835 msc:  for pi+  SubType= 10
836       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
837         WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=240 100 eV  - 100 TeV
838           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
839 
840 hIoni:  for pi+  XStype:3  SubType=2
841       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
842       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
843       StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
844       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
845                Bragg : Emin=    0 eV  Emax=297.505 keV  deltaVI
846           BetheBloch : Emin=297.505 keV Emax=  100 TeV  deltaVI
847 
848 hBrems:  for pi+  XStype:1  SubType=3
849       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
850       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
851       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
852                hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
853       ===== Limit on energy threshold has been applied 
854 
855 hPairProd:  for pi+  XStype:1  SubType=4
856       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
857       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
858       Sampling table 20x1001 from 1.11656 GeV to 100 TeV 
859       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
860            hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
861 
862 CoulombScat:  for pi+ XStype:1 SubType=1 BuildTable=1
863       Lambda table from threshold  to 100 TeV, 20 bins/decade, spline: 0
864       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
865       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
866   eCoulombScattering : Emin=    0 eV  Emax=  100 TeV
867 
868 msc:  for pi-  SubType= 10
869       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
870         WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=240 100 eV  - 100 TeV
871           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
872 
873 hIoni:  for pi-  XStype:3  SubType=2
874       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
875       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
876       StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
877       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
878             ICRU73QO : Emin=    0 eV  Emax=297.505 keV  deltaVI
879           BetheBloch : Emin=297.505 keV Emax=  100 TeV  deltaVI
880 
881 hBrems:  for pi-  XStype:1  SubType=3
882       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
883       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
884       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
885                hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
886       ===== Limit on energy threshold has been applied 
887 
888 hPairProd:  for pi-  XStype:1  SubType=4
889       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
890       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
891       Sampling table 20x1001 from 1.11656 GeV to 100 TeV 
892       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
893            hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
894 
895 CoulombScat:  for pi- XStype:1 SubType=1 BuildTable=1
896       Used Lambda table of pi+
897       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
898       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
899   eCoulombScattering : Emin=    0 eV  Emax=  100 TeV
900 
901 ========= Table of registered couples ============================
902 
903 Index : 0     used in the geometry : Yes
904  Material : Air
905  Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm 
906  Energy thresholds :  gamma  250 eV     e-  329.56 eV     e+  325.472 eV  proton 100 keV
907  Region(s) which use this couple : 
908     DefaultRegionForTheWorld
909 
910 Index : 1     used in the geometry : Yes
911  Material : soft_tissue
912  Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm 
913  Energy thresholds :  gamma  2.80001 keV    e-  347.402 keV    e+  337.896 keV proton 100 keV
914  Region(s) which use this couple : 
915     DefaultRegionForTheWorld
916 
917 Index : 2     used in the geometry : Yes
918  Material : skeleton
919  Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm 
920  Energy thresholds :  gamma  4.00973 keV    e-  441.639 keV    e+  427.712 keV proton 100 keV
921  Region(s) which use this couple : 
922     DefaultRegionForTheWorld
923 
924 ==================================================================
925 
926 ### Run 0 start.
927 ... set ntuple merging mode : G4NtupleMergeMode::kNone - done
928 ... set ntuple merging mode : G4NtupleMergeMode::kNone - done
929 ... set ntuple merging row mode : row-wise - done
930 ... create directory for ntuples : human_phantom_ntuple - done
931 ... create file : human_phantom.root - done
932 ... open analysis file : human_phantom.root - done
933 ... open analysis file : human_phantom.root - done
934 ... create ntuple booking : 1 ntupleId 1 - done
935 ... create ntuple T column : organID ntupleId 1 - done
936 ... create ntuple T column : edep ntupleId 1 - done
937 Number of events = 100
938 Energy Total in Run:logicalBrain, ID: 0, Energy Deposition (MeV): 0
939 Energy Total in Run:logicalHead, ID: 1, Energy Deposition (MeV): 4.63096
940 Analysis: organ 1 edep: 4.63096
941 Energy Total in Run:logicalHeart, ID: 2, Energy Deposition (MeV): 0
942 Energy Total in Run:logicalLeftAdrenal, ID: 3, Energy Deposition (MeV): 0
943 Energy Total in Run:logicalLeftArmBone, ID: 4, Energy Deposition (MeV): 0
944 Energy Total in Run:logicalLeftBreast, ID: 5, Energy Deposition (MeV): 0
945 Energy Total in Run:logicalLeftClavicle, ID: 6, Energy Deposition (MeV): 0
946 Energy Total in Run:logicalLeftKidney, ID: 7, Energy Deposition (MeV): 0
947 Energy Total in Run:logicalLeftLeg, ID: 8, Energy Deposition (MeV): 14.3441
948 Analysis: organ 8 edep: 14.3441
949 Energy Total in Run:logicalLeftLegBone, ID: 9, Energy Deposition (MeV): 0
950 Energy Total in Run:logicalLeftLung, ID: 10, Energy Deposition (MeV): 0
951 Energy Total in Run:logicalLeftOvary, ID: 11, Energy Deposition (MeV): 0
952 Energy Total in Run:logicalLeftScapula, ID: 12, Energy Deposition (MeV): 0
953 Energy Total in Run:logicalLeftTeste, ID: 13, Energy Deposition (MeV): 0
954 Energy Total in Run:logicalLowerLargeIntestine, ID: 14, Energy Deposition (MeV): 0
955 Energy Total in Run:logicalMaleGenitalia, ID: 15, Energy Deposition (MeV): 0
956 Energy Total in Run:logicalMiddleLowerSpine, ID: 16, Energy Deposition (MeV): 0
957 Energy Total in Run:logicalPancreas, ID: 17, Energy Deposition (MeV): 0
958 Energy Total in Run:logicalPelvis, ID: 18, Energy Deposition (MeV): 0
959 Energy Total in Run:logicalRibCage, ID: 19, Energy Deposition (MeV): 0.659836
960 Analysis: organ 19 edep: 0.659836
961 Energy Total in Run:logicalRightAdrenal, ID: 20, Energy Deposition (MeV): 0
962 Energy Total in Run:logicalRightArmBone, ID: 21, Energy Deposition (MeV): 0
963 Energy Total in Run:logicalRightBreast, ID: 22, Energy Deposition (MeV): 0
964 Energy Total in Run:logicalRightClavicle, ID: 23, Energy Deposition (MeV): 0
965 Energy Total in Run:logicalRightKidney, ID: 24, Energy Deposition (MeV): 0
966 Energy Total in Run:logicalRightLeg, ID: 25, Energy Deposition (MeV): 10.6217
967 Analysis: organ 25 edep: 10.6217
968 Energy Total in Run:logicalRightLegBone, ID: 26, Energy Deposition (MeV): 0
969 Energy Total in Run:logicalRightLung, ID: 27, Energy Deposition (MeV): 0
970 Energy Total in Run:logicalRightOvary, ID: 28, Energy Deposition (MeV): 0
971 Energy Total in Run:logicalRightScapula, ID: 29, Energy Deposition (MeV): 0
972 Energy Total in Run:logicalRightTeste, ID: 30, Energy Deposition (MeV): 0
973 Energy Total in Run:logicalSkull, ID: 31, Energy Deposition (MeV): 0
974 Energy Total in Run:logicalSmallIntestine, ID: 32, Energy Deposition (MeV): 0
975 Energy Total in Run:logicalSpleen, ID: 33, Energy Deposition (MeV): 0
976 Energy Total in Run:logicalStomach, ID: 34, Energy Deposition (MeV): 0
977 Energy Total in Run:logicalThymus, ID: 35, Energy Deposition (MeV): 0
978 Energy Total in Run:logicalThyroid, ID: 36, Energy Deposition (MeV): 0
979 Energy Total in Run:logicalTrunk, ID: 37, Energy Deposition (MeV): 23.6957
980 Analysis: organ 37 edep: 23.6957
981 Energy Total in Run:logicalUpperLargeIntestine, ID: 38, Energy Deposition (MeV): 0
982 Energy Total in Run:logicalUpperSpine, ID: 39, Energy Deposition (MeV): 0
983 Energy Total in Run:logicalUrinaryBladder, ID: 40, Energy Deposition (MeV): 0
984 Energy Total in Run:logicalUterus, ID: 41, Energy Deposition (MeV): 0
985 Total Energy deposit in the body is: 53.9523 MeV
986 ... write file : human_phantom.root - done
987 ... close file : human_phantom.root - done
988 ... clear files - done
989 ... clear h1 - done
990 ... clear h2 - done
991 ... clear h3 - done
992 ... clear p1 - done
993 ... clear p2 - done
994 ... clear ntupleBookings - done
995 ... clear ntuples - done
996 ... clear files - done
997 ... clear all data - done
998 Graphics systems deleted.
999 Visualization Manager deleting...