Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/hadrontherapy/hadrontherapy.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 Activating geometry default
 23 Going to register Parallel world...... done
 24 Going to register G4ParallelWorldPhysics
 25 Accumulable registered as "RBE"
 26 Visualization Manager instantiating with verbosity "warnings (3)"...
 27 Visualization Manager initialising...
 28 Registering graphics systems...
 29 
 30 You have successfully registered the following graphics systems.
 31 Registered graphics systems are:
 32   ASCIITree (ATree)
 33   DAWNFILE (DAWNFILE)
 34   G4HepRepFile (HepRepFile)
 35   RayTracer (RayTracer)
 36   VRML2FILE (VRML2FILE)
 37   gMocrenFile (gMocrenFile)
 38   TOOLSSG_OFFSCREEN (TSG_OFFSCREEN, TSG_FILE)
 39   OpenGLImmediateQt (OGLIQt, OGLI)
 40   OpenGLStoredQt (OGLSQt, OGL, OGLS)
 41   OpenGLImmediateXm (OGLIXm, OGLIQt_FALLBACK)
 42   OpenGLStoredXm (OGLSXm, OGLSQt_FALLBACK)
 43   OpenGLImmediateX (OGLIX, OGLIQt_FALLBACK, OGLIXm_FALLBACK)
 44   OpenGLStoredX (OGLSX, OGLSQt_FALLBACK, OGLSXm_FALLBACK)
 45   RayTracerX (RayTracerX)
 46   Qt3D (Qt3D)
 47   TOOLSSG_X11_GLES (TSG_X11_GLES, TSGX11, TSG_XT_GLES_FALLBACK)
 48   TOOLSSG_X11_ZB (TSG_X11_ZB, TSGX11ZB)
 49   TOOLSSG_XT_GLES (TSG_XT_GLES, TSGXt, TSG_QT_GLES_FALLBACK)
 50   TOOLSSG_XT_ZB (TSG_XT_ZB, TSGXtZB)
 51   TOOLSSG_QT_GLES (TSG_QT_GLES, TSGQt, TSG)
 52   TOOLSSG_QT_ZB (TSG_QT_ZB, TSGQtZB)
 53 You may choose a graphics system (driver) with a parameter of
 54 the command "/vis/open" or "/vis/sceneHandler/create",
 55 or you may omit the driver parameter and choose at run time:
 56 - by argument in the construction of G4VisExecutive
 57 - by environment variable "G4VIS_DEFAULT_DRIVER"
 58 - by entry in "~/.g4session"
 59 - by build flags.
 60 - Note: This feature is not allowed in batch mode.
 61 For further information see "examples/basic/B1/exampleB1.cc"
 62 and "vis.mac".
 63 
 64 Registering model factories...
 65 
 66 You have successfully registered the following model factories.
 67 Registered model factories:
 68   generic
 69   drawByAttribute
 70   drawByCharge
 71   drawByOriginVolume
 72   drawByParticleID
 73   drawByEncounteredVolume
 74 
 75 Registered models:
 76   None
 77 
 78 Registered filter factories:
 79   attributeFilter
 80   chargeFilter
 81   originVolumeFilter
 82   particleFilter
 83   encounteredVolumeFilter
 84 
 85 Registered filters:
 86   None
 87 
 88 You have successfully registered the following user vis actions.
 89 Run Duration User Vis Actions: none
 90 End of Event User Vis Actions: none
 91 End of Run User Vis Actions: none
 92 
 93 Some /vis commands (optionally) take a string to specify colour.
 94 "/vis/list" to see available colours.
 95 /tracking/verbose 0
 96 /run/verbose 1
 97 /event/verbose 0
 98 /Physics/addPhysics standard_opt4
 99 THE FOLLOWING ELECTROMAGNETIC PHYSICS LIST HAS BEEN ACTIVATED: G4EmStandardPhysics_option4
100 /run/initialize
101 /run/geometryModified
102 HadrontherapyMatrix: Memory space to store physical dose into 200 voxels has been allocated 
103 /run/geometryModified
104 The (X,Y,Z) dimensions of the phantom are : (40 cm ,40 cm ,40 cm )
105 The (X,Y,Z) dimensions of the detector are : (40 cm ,40 cm ,40 cm )
106 Displacement between Phantom and World is: DX= 20 cm DY= 0 fm DZ= 0 fm 
107 The (X,Y,Z) sizes of the Voxels are: (200 um ,4 cm ,4 cm )
108 The number of Voxels along (X,Y,Z) is: (200,1,1)
109 --- G4CoupledTransportation is used 
110 Warning : Region <DetectorLog> does not have specific production cuts,
111 even though it appears in the current tracking world.
112 Default cuts are used for this region.
113 /gps/pos/shape Circle
114 /gps/pos/centre -310. 0. 0. cm
115 /gps/pos/radius 0. mm
116 /gps/pos/sigma_r 2. mm
117 /gps/particle proton
118 /gps/pos/type Beam
119 /gps/pos/rot1 0 1 0
120 /gps/pos/rot2 0 0 1
121 /gps/ang/rot1 0 0 1
122 /gps/ang/rot2 0 1 0
123 /gps/ang/type beam1d
124 /gps/ang/sigma_r 0. deg
125 /gps/ene/type Gauss
126 /gps/ene/mono 62 MeV
127 /gps/ene/sigma 0.3 MeV
128 /Step/waterPhantomStepMax 1 mm
129 /changePhantom/size 40 40 40 cm
130 /changePhantom/position 20 0 0 cm
131 /changeDetector/size 4 4 4 cm
132 /changeDetector/voxelSize 1 40 40 mm
133 /changeDetector/displacement 0 18 18 cm
134 /changePhantom/update
135 HadrontherapyMatrix: Memory space to store physical dose into 40 voxels has been allocated 
136 /run/geometryModified
137 The (X,Y,Z) dimensions of the phantom are : (40 cm ,40 cm ,40 cm )
138 The (X,Y,Z) dimensions of the detector are : (4 cm ,4 cm ,4 cm )
139 Displacement between Phantom and World is: DX= 20 cm DY= 0 fm DZ= 0 fm 
140 The (X,Y,Z) sizes of the Voxels are: (1 mm ,4 cm ,4 cm )
141 The number of Voxels along (X,Y,Z) is: (40,1,1)
142 /event/printEventNumber 100
143 /run/beamOn 500
144 =======================================================================
145 ======                 Electromagnetic Physics Parameters      ========
146 =======================================================================
147 LPM effect enabled                                 1
148 Enable creation and use of sampling tables         0
149 Apply cuts on all EM processes                     0
150 Use combined TransportationWithMsc                 Disabled
151 Use general process                                1
152 Enable linear polarisation for gamma               0
153 Enable photoeffect sampling below K-shell          1
154 Enable sampling of quantum entanglement            0
155 X-section factor for integral approach             0.8
156 Min kinetic energy for tables                      100 eV 
157 Max kinetic energy for tables                      100 TeV
158 Number of bins per decade of a table               20
159 Verbose level                                      1
160 Verbose level for worker thread                    0
161 Bremsstrahlung energy threshold above which 
162   primary e+- is added to the list of secondary    100 TeV
163 Bremsstrahlung energy threshold above which primary
164   muon/hadron is added to the list of secondary    100 TeV
165 Positron annihilation at rest model                AllisonPositronium
166 Enable 3 gamma annihilation on fly                 1
167 Lowest triplet kinetic energy                      1 MeV
168 Enable sampling of gamma linear polarisation       0
169 5D gamma conversion model type                     0
170 5D gamma conversion model on isolated ion          0
171 Use Ricardo-Gerardo pair production model          0
172 Livermore data directory                           epics_2017
173 =======================================================================
174 ======                 Ionisation Parameters                   ========
175 =======================================================================
176 Step function for e+-                              (0.2, 0.01 mm)
177 Step function for muons/hadrons                    (0.1, 0.05 mm)
178 Step function for light ions                       (0.1, 0.02 mm)
179 Step function for general ions                     (0.1, 0.001 mm)
180 Lowest e+e- kinetic energy                         100 eV 
181 Lowest muon/hadron kinetic energy                  1 keV
182 Use ICRU90 data                                    1
183 Fluctuations of dE/dx are enabled                  1
184 Type of fluctuation model for leptons and hadrons  Urban
185 Use built-in Birks satuaration                     0
186 Build CSDA range enabled                           0
187 Use cut as a final range enabled                   0
188 Enable angular generator interface                 1
189 Max kinetic energy for CSDA tables                 1 GeV
190 Max kinetic energy for NIEL computation            1 MeV
191 Linear loss limit                                  0.01
192 Read data from file for e+e- pair production by mu 0
193 =======================================================================
194 ======                 Multiple Scattering Parameters          ========
195 =======================================================================
196 Type of msc step limit algorithm for e+-           2
197 Type of msc step limit algorithm for muons/hadrons 0
198 Msc lateral displacement for e+- enabled           1
199 Msc lateral displacement for muons and hadrons     1
200 Urban msc model lateral displacement alg96         1
201 Range factor for msc step limit for e+-            0.08
202 Range factor for msc step limit for muons/hadrons  0.2
203 Geometry factor for msc step limitation of e+-     2.5
204 Safety factor for msc step limit for e+-           0.6
205 Skin parameter for msc step limitation of e+-      3
206 Lambda limit for msc step limit for e+-            1 mm
207 Use Mott correction for e- scattering              1
208 Factor used for dynamic computation of angular 
209   limit between single and multiple scattering     1
210 Fixed angular limit between single 
211   and multiple scattering                          3.1416 rad
212 Upper energy limit for e+- multiple scattering     100 MeV
213 Type of electron single scattering model           0
214 Type of nuclear form-factor                        1
215 Screening factor                                   1
216 =======================================================================
217 ======                 Atomic Deexcitation Parameters          ========
218 =======================================================================
219 Fluorescence enabled                               1
220 Directory in G4LEDATA for fluorescence data files  fluor
221 Auger electron cascade enabled                     0
222 PIXE atomic de-excitation enabled                  0
223 De-excitation module ignores cuts                  0
224 Type of PIXE cross section for hadrons             Empirical
225 Type of PIXE cross section for e+-                 Livermore
226 =======================================================================
227 
228 ### ===  Deexcitation model UAtomDeexcitation is activated for 2 regions:
229           DefaultRegionForTheWorld  1  0  0
230           DetectorLog  1  0  0
231 ### ===  Ignore cuts flag:   0
232 
233 phot:  for gamma SubType=12 BuildTable=0
234       LambdaPrime table from 200 keV to 100 TeV in 174 bins 
235       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
236  LivermorePhElectric : Emin=    0 eV  Emax=  100 TeV  SauterGavrila Fluo
237 
238 compt:  for gamma SubType=13 BuildTable=1
239       Lambda table from 100 eV  to 1 MeV, 20 bins/decade, spline: 1
240       LambdaPrime table from 1 MeV to 100 TeV in 160 bins 
241       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
242    LowEPComptonModel : Emin=    0 eV  Emax=   20 MeV Fluo
243         KleinNishina : Emin=   20 MeV Emax=  100 TeV Fluo
244 
245 conv:  for gamma SubType=14 BuildTable=1
246       Lambda table from 1.022 MeV to 100 TeV, 20 bins/decade, spline: 1
247       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
248       BetheHeitler5D : Emin=    0 eV  Emax=  100 TeV  ModifiedTsai
249 
250 Rayl:  for gamma SubType=11 BuildTable=1
251       Lambda table from 100 eV  to 150 keV, 20 bins/decade, spline: 0
252       LambdaPrime table from 150 keV to 100 TeV in 176 bins 
253       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
254    LivermoreRayleigh : Emin=    0 eV  Emax=  100 TeV  CullenGenerator
255 
256 msc:  for e-  SubType= 10
257       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
258   GoudsmitSaunderson : Emin=    0 eV  Emax=  100 MeV Nbins=120 100 eV  - 100 MeV
259           StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
260         WentzelVIUni : Emin=  100 MeV Emax=  100 TeV Nbins=120 100 MeV - 100 TeV
261           StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
262 
263 eIoni:  for e-  XStype:3  SubType=2
264       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
265       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
266       StepFunction=(0.2, 0.01 mm), integ: 3, fluct: 1, linLossLim= 0.01
267       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
268              PenIoni : Emin=    0 eV  Emax=  100 keV
269         MollerBhabha : Emin=  100 keV Emax=  100 TeV  deltaVI
270 
271 eBrem:  for e-  XStype:4  SubType=3
272       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
273       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
274       LPM flag: 1 for E > 1 GeV,  VertexHighEnergyTh(GeV)= 100000
275       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
276              eBremSB : Emin=    0 eV  Emax=    1 GeV  AngularGen2BS
277             eBremLPM : Emin=    1 GeV Emax=  100 TeV  AngularGen2BS
278 
279 ePairProd:  for e-  XStype:1  SubType=4
280       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
281       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 0
282       Sampling table 25x1001 from 0.1 GeV to 100 TeV 
283       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
284            ePairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
285 
286 CoulombScat:  for e- XStype:1 SubType=1 BuildTable=1
287       Lambda table from 100 MeV to 100 TeV, 20 bins/decade, spline: 0
288       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
289       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
290   eCoulombScattering : Emin=  100 MeV Emax=  100 TeV
291 
292 msc:  for e+  SubType= 10
293       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
294   GoudsmitSaunderson : Emin=    0 eV  Emax=  100 MeV Nbins=120 100 eV  - 100 MeV
295           StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
296         WentzelVIUni : Emin=  100 MeV Emax=  100 TeV Nbins=120 100 MeV - 100 TeV
297           StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
298 
299 eIoni:  for e+  XStype:3  SubType=2
300       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
301       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
302       StepFunction=(0.2, 0.01 mm), integ: 3, fluct: 1, linLossLim= 0.01
303       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
304              PenIoni : Emin=    0 eV  Emax=  100 keV
305         MollerBhabha : Emin=  100 keV Emax=  100 TeV  deltaVI
306 
307 eBrem:  for e+  XStype:4  SubType=3
308       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
309       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
310       LPM flag: 1 for E > 1 GeV,  VertexHighEnergyTh(GeV)= 100000
311       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
312              eBremSB : Emin=    0 eV  Emax=    1 GeV  AngularGen2BS
313             eBremLPM : Emin=    1 GeV Emax=  100 TeV  AngularGen2BS
314 
315 ePairProd:  for e+  XStype:1  SubType=4
316       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
317       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 0
318       Sampling table 25x1001 from 0.1 GeV to 100 TeV 
319       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
320            ePairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
321 
322 annihil:  for e+ XStype:2 SubType=5 AtRestModel:Allison BuildTable=0
323       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
324     eplusTo2or3gamma : Emin=    0 eV  Emax=  100 TeV
325 
326 CoulombScat:  for e+ XStype:1 SubType=1 BuildTable=1
327       Lambda table from 100 MeV to 100 TeV, 20 bins/decade, spline: 0
328       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
329       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
330   eCoulombScattering : Emin=  100 MeV Emax=  100 TeV
331 
332 msc:  for proton  SubType= 10
333       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
334         WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=240 100 eV  - 100 TeV
335           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
336 
337 hIoni:  for proton  XStype:3  SubType=2
338       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
339       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
340       StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
341       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
342                Bragg : Emin=    0 eV  Emax=    2 MeV  deltaVI
343           BetheBloch : Emin=    2 MeV Emax=  100 TeV  deltaVI
344 
345 hBrems:  for proton  XStype:1  SubType=3
346       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
347       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
348       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
349                hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
350 
351 hPairProd:  for proton  XStype:1  SubType=4
352       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
353       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
354       Sampling table 17x1001 from 7.50618 GeV to 100 TeV 
355       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
356            hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
357 
358 CoulombScat:  for proton XStype:1 SubType=1 BuildTable=1
359       Lambda table from threshold  to 100 TeV, 20 bins/decade, spline: 0
360       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
361       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
362   eCoulombScattering : Emin=    0 eV  Emax=  100 TeV
363 
364 nuclearStopping:  for proton SubType=8 BuildTable=0
365       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
366    ICRU49NucStopping : Emin=    0 eV  Emax=    1 MeV
367 
368 msc:  for GenericIon  SubType= 10
369       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
370             UrbanMsc : Emin=    0 eV  Emax=  100 TeV
371           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
372 
373 ionIoni:  for GenericIon  XStype:3  SubType=2
374       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
375       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
376       StepFunction=(0.1, 0.001 mm), integ: 3, fluct: 1, linLossLim= 0.02
377       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
378     LindhardSorensen : Emin=    0 eV  Emax=  100 TeV  deltaVI
379 
380 nuclearStopping:  for GenericIon SubType=8 BuildTable=0
381       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
382    ICRU49NucStopping : Emin=    0 eV  Emax=    1 MeV
383 
384 msc:  for alpha  SubType= 10
385       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
386             UrbanMsc : Emin=    0 eV  Emax=  100 TeV
387           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
388 
389 ionIoni:  for alpha  XStype:3  SubType=2
390       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
391       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
392       StepFunction=(0.1, 0.02 mm), integ: 3, fluct: 1, linLossLim= 0.02
393       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
394             BraggIon : Emin=    0 eV  Emax=7.9452 MeV  deltaVI
395           BetheBloch : Emin=7.9452 MeV Emax=  100 TeV  deltaVI
396 
397 nuclearStopping:  for alpha SubType=8 BuildTable=0
398       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
399    ICRU49NucStopping : Emin=    0 eV  Emax=    1 MeV
400 
401 msc:  for anti_proton  SubType= 10
402       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
403         WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=240 100 eV  - 100 TeV
404           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
405 
406 hIoni:  for anti_proton  XStype:3  SubType=2
407       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
408       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
409       StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
410       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
411             ICRU73QO : Emin=    0 eV  Emax=    2 MeV  deltaVI
412           BetheBloch : Emin=    2 MeV Emax=  100 TeV  deltaVI
413 
414 hBrems:  for anti_proton  XStype:1  SubType=3
415       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
416       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
417       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
418                hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
419 
420 hPairProd:  for anti_proton  XStype:1  SubType=4
421       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
422       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
423       Sampling table 17x1001 from 7.50618 GeV to 100 TeV 
424       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
425            hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
426 
427 CoulombScat:  for anti_proton XStype:1 SubType=1 BuildTable=1
428       Lambda table from threshold  to 100 TeV, 20 bins/decade, spline: 0
429       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
430       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
431   eCoulombScattering : Emin=    0 eV  Emax=  100 TeV
432 
433 msc:  for kaon+  SubType= 10
434       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
435         WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=240 100 eV  - 100 TeV
436           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
437 
438 hIoni:  for kaon+  XStype:3  SubType=2
439       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
440       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
441       StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
442       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
443                Bragg : Emin=    0 eV  Emax=1.05231 MeV  deltaVI
444           BetheBloch : Emin=1.05231 MeV Emax=  100 TeV  deltaVI
445 
446 hBrems:  for kaon+  XStype:1  SubType=3
447       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
448       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
449       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
450                hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
451 
452 hPairProd:  for kaon+  XStype:1  SubType=4
453       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
454       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
455       Sampling table 18x1001 from 3.94942 GeV to 100 TeV 
456       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
457            hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
458 
459 CoulombScat:  for kaon+ XStype:1 SubType=1 BuildTable=1
460       Lambda table from threshold  to 100 TeV, 20 bins/decade, spline: 0
461       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
462       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
463   eCoulombScattering : Emin=    0 eV  Emax=  100 TeV
464 
465 msc:  for kaon-  SubType= 10
466       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
467         WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=240 100 eV  - 100 TeV
468           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
469 
470 hIoni:  for kaon-  XStype:3  SubType=2
471       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
472       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
473       StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
474       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
475             ICRU73QO : Emin=    0 eV  Emax=1.05231 MeV  deltaVI
476           BetheBloch : Emin=1.05231 MeV Emax=  100 TeV  deltaVI
477 
478 hBrems:  for kaon-  XStype:1  SubType=3
479       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
480       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
481       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
482                hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
483 
484 hPairProd:  for kaon-  XStype:1  SubType=4
485       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
486       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
487       Sampling table 18x1001 from 3.94942 GeV to 100 TeV 
488       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
489            hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
490 
491 CoulombScat:  for kaon- XStype:1 SubType=1 BuildTable=1
492       Used Lambda table of kaon+
493       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
494       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
495   eCoulombScattering : Emin=    0 eV  Emax=  100 TeV
496 
497 msc:  for mu+  SubType= 10
498       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
499         WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=240 100 eV  - 100 TeV
500           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
501 
502 muIoni:  for mu+  XStype:3  SubType=2
503       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
504       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
505       StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
506       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
507                Bragg : Emin=    0 eV  Emax=  200 keV  deltaVI
508         MuBetheBloch : Emin=  200 keV Emax=  100 TeV  deltaVI
509 
510 muBrems:  for mu+  XStype:1  SubType=3
511       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
512       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
513       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
514               MuBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
515 
516 muPairProd:  for mu+  XStype:1  SubType=4
517       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
518       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
519       Sampling table 21x1001 from 0.85 GeV to 100 TeV 
520       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
521           muPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
522 
523 CoulombScat:  for mu+ XStype:1 SubType=1 BuildTable=1
524       Lambda table from threshold  to 100 TeV, 20 bins/decade, spline: 0
525       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
526       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
527   eCoulombScattering : Emin=    0 eV  Emax=  100 TeV
528 
529 msc:  for mu-  SubType= 10
530       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
531         WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=240 100 eV  - 100 TeV
532           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
533 
534 muIoni:  for mu-  XStype:3  SubType=2
535       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
536       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
537       StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
538       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
539             ICRU73QO : Emin=    0 eV  Emax=  200 keV  deltaVI
540         MuBetheBloch : Emin=  200 keV Emax=  100 TeV  deltaVI
541 
542 muBrems:  for mu-  XStype:1  SubType=3
543       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
544       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
545       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
546               MuBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
547 
548 muPairProd:  for mu-  XStype:1  SubType=4
549       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
550       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
551       Sampling table 21x1001 from 0.85 GeV to 100 TeV 
552       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
553           muPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
554 
555 CoulombScat:  for mu- XStype:1 SubType=1 BuildTable=1
556       Used Lambda table of mu+
557       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
558       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
559   eCoulombScattering : Emin=    0 eV  Emax=  100 TeV
560 
561 msc:  for pi+  SubType= 10
562       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
563         WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=240 100 eV  - 100 TeV
564           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
565 
566 hIoni:  for pi+  XStype:3  SubType=2
567       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
568       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
569       StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
570       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
571                Bragg : Emin=    0 eV  Emax=297.505 keV  deltaVI
572           BetheBloch : Emin=297.505 keV Emax=  100 TeV  deltaVI
573 
574 hBrems:  for pi+  XStype:1  SubType=3
575       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
576       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
577       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
578                hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
579 
580 hPairProd:  for pi+  XStype:1  SubType=4
581       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
582       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
583       Sampling table 20x1001 from 1.11656 GeV to 100 TeV 
584       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
585            hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
586 
587 CoulombScat:  for pi+ XStype:1 SubType=1 BuildTable=1
588       Lambda table from threshold  to 100 TeV, 20 bins/decade, spline: 0
589       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
590       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
591   eCoulombScattering : Emin=    0 eV  Emax=  100 TeV
592 
593 msc:  for pi-  SubType= 10
594       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
595         WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=240 100 eV  - 100 TeV
596           StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
597 
598 hIoni:  for pi-  XStype:3  SubType=2
599       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
600       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
601       StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
602       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
603             ICRU73QO : Emin=    0 eV  Emax=297.505 keV  deltaVI
604           BetheBloch : Emin=297.505 keV Emax=  100 TeV  deltaVI
605 
606 hBrems:  for pi-  XStype:1  SubType=3
607       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
608       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
609       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
610                hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
611 
612 hPairProd:  for pi-  XStype:1  SubType=4
613       dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
614       Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
615       Sampling table 20x1001 from 1.11656 GeV to 100 TeV 
616       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
617            hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
618 
619 CoulombScat:  for pi- XStype:1 SubType=1 BuildTable=1
620       Used Lambda table of pi+
621       ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
622       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
623   eCoulombScattering : Emin=    0 eV  Emax=  100 TeV
624 
625 ========= Table of registered couples ============================
626 
627 Index : 0     used in the geometry : Yes
628  Material : G4_AIR
629  Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm 
630  Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 100 keV
631  Region(s) which use this couple : 
632     DefaultRegionForTheWorld
633 
634 Index : 1     used in the geometry : Yes
635  Material : G4_Al
636  Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm 
637  Energy thresholds :  gamma  6.91572 keV    e-  597.438 keV    e+  574.135 keV proton 100 keV
638  Region(s) which use this couple : 
639     DefaultRegionForTheWorld
640 
641 Index : 2     used in the geometry : Yes
642  Material : G4_Galactic
643  Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm 
644  Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 100 keV
645  Region(s) which use this couple : 
646     DefaultRegionForTheWorld
647 
648 Index : 3     used in the geometry : Yes
649  Material : G4_Ta
650  Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm 
651  Energy thresholds :  gamma  101.483 keV    e-  2.01234 MeV    e+  1.87285 MeV proton 100 keV
652  Region(s) which use this couple : 
653     DefaultRegionForTheWorld
654 
655 Index : 4     used in the geometry : Yes
656  Material : G4_KAPTON
657  Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm 
658  Energy thresholds :  gamma  2.97329 keV    e-  419.475 keV    e+  406.697 keV proton 100 keV
659  Region(s) which use this couple : 
660     DefaultRegionForTheWorld
661 
662 Index : 5     used in the geometry : Yes
663  Material : Brass
664  Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm 
665  Energy thresholds :  gamma  24.3022 keV    e-  1.32652 MeV    e+  1.25172 MeV proton 100 keV
666  Region(s) which use this couple : 
667     DefaultRegionForTheWorld
668 
669 Index : 6     used in the geometry : Yes
670  Material : G4_PLEXIGLASS
671  Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm 
672  Energy thresholds :  gamma  2.77786 keV    e-  389.158 keV    e+  377.95 keV proton 100 keV
673  Region(s) which use this couple : 
674     DefaultRegionForTheWorld
675 
676 Index : 7     used in the geometry : Yes
677  Material : G4_Cu
678  Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm 
679  Energy thresholds :  gamma  24.8228 keV    e-  1.40308 MeV    e+  1.3229 MeV proton 100 keV
680  Region(s) which use this couple : 
681     DefaultRegionForTheWorld
682 
683 Index : 8     used in the geometry : Yes
684  Material : G4_MYLAR
685  Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm 
686  Energy thresholds :  gamma  3.01225 keV    e-  420.616 keV    e+  407.833 keV proton 100 keV
687  Region(s) which use this couple : 
688     DefaultRegionForTheWorld
689 
690 Index : 9     used in the geometry : Yes
691  Material : G4_WATER
692  Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm 
693  Energy thresholds :  gamma  2.93547 keV    e-  351.353 keV    e+  341.657 keV proton 100 keV
694  Region(s) which use this couple : 
695     DefaultRegionForTheWorld
696     DetectorLog
697 
698 ==================================================================
699 
700 ### Run 0 starts.
701 Run 0 starts ...
702  Run terminated.
703 Run Summary
704   Number of events processed : 500
705   User=5.160000s Real=5.450592s Sys=0.080000s
706 Graphics systems deleted.
707 Visualization Manager deleting...
708 The simulation took: 8.16199 s to run (real time)
709 Dose is being written to Dose.out
710 i j k Dose(Gy)================== Deleting memory pools ===================
711 Number of memory pools allocated: 12 of which, static: 0
712 Dynamic pools deleted: 12 / Total memory freed: 0.32 MB
713 ============================================================