source: trunk/examples/extended/electromagnetic/TestEm5/TestEm5.out @ 1337

Last change on this file since 1337 was 1337, checked in by garnier, 14 years ago

tag geant4.9.4 beta 1 + modifs locales
File size: 94.3 KB
Line 
1
2        ############################################
3        !!! WARNING - FPE detection is activated !!!
4        ############################################
5
6*************************************************************
7 Geant4 version Name: geant4-09-03-ref-06    (25-June-2010)
8                      Copyright : Geant4 Collaboration
9                      Reference : NIM A 506 (2003), 250-303
10                            WWW : http://cern.ch/geant4
11*************************************************************
12
13/run/verbose 2
14#
15/testem/det/setAbsMat Silicon
16/testem/det/setAbsThick 20 um
17/testem/det/setAbsYZ 50 um
18#
19/testem/phys/addPhysics emstandard_opt3
20PhysicsList::AddPhysicsList: <emstandard_opt3>
21#
22/testem/phys/setCuts 7 um
23#
24/run/initialize
25userDetector->Construct() start.
26
27 Material: Galactic     density:  0.000 kg/m3   RadL: 204727512.315 pc   Nucl.Int.Length: 113804112.800 pc   Imean:  19.200 eV   temperature:   2.73 K  pressure:   0.00 atm
28   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
29
30 Material:  Silicon     density:  2.330 g/cm3   RadL:   9.368 cm   Nucl.Int.Length:  45.663 cm   Imean: 173.000 eV
31   --->  Element: Silicon ( )   Z = 14.0   N =  28.1   A =  28.09 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
32
33 The  WORLD   is made of 30 um  of Galactic. The transverse size (YZ) of the world is 60 um
34 The ABSORBER is made of 20 um  of Silicon. The transverse size (YZ) is 50 um
35 X position of the middle of the absorber 0 fm
36World is registered to the default region.
37physicsList->Construct() start.
38physicsList->Construct() start.
39physicsList->CheckParticleList() start.
40physicsList->setCut() start.
41PhysicsList::SetCuts:CutLength : 1 mm
42#
43/testem/gun/setDefault
44/gun/particle pi+
45/gun/energy 5 GeV
46#
47/testem/histo/setFileName testem5
48/testem/histo/setFileType root
49/testem/histo/setHisto 13 100 -0.25 0.25 mrad
50----> SetHisto 13: (transmit, charged) : projected angle at exit (mrad);  100 bins from -0.25 mrad to 0.25 mrad
51#
52/run/beamOn 100000
53
54phot:   for  gamma    SubType= 12
55      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
56       PhotoElectric :     Emin=          0 eV         Emax=   10 TeV
57
58compt:   for  gamma    SubType= 13
59      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
60      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
61       Klein-Nishina :     Emin=          0 eV         Emax=   10 TeV
62
63conv:   for  gamma    SubType= 14
64      Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1
65      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
66       Bethe-Heitler :     Emin=          0 eV         Emax=   10 TeV
67
68Rayl:   for  gamma    SubType= 11
69      Lambda tables from 100 eV  to 10 TeV in 200 bins, spline: 1
70      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
71   LivermoreRayleigh :     Emin=          0 eV         Emax=   100 GeV
72
73msc:   for e-    SubType= 10
74      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
75      RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
76      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
77          UrbanMsc93 :     Emin=          0 eV         Emax=   10 TeV
78
79eIoni:   for  e-    SubType= 2
80      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
81      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
82      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
83      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
84        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
85
86eBrem:   for  e-    SubType= 3
87      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
88      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
89      LPM flag: 1 for E > 1 GeV
90      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
91               eBrem :     Emin=          0 eV         Emax=   1 GeV
92            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
93
94eIoni:   for  e+    SubType= 2
95      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
96      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
97      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
98      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
99        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
100
101eBrem:   for  e+    SubType= 3
102      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
103      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
104      LPM flag: 1 for E > 1 GeV
105      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
106               eBrem :     Emin=          0 eV         Emax=   1 GeV
107            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
108
109annihil:   for  e+    SubType= 5
110      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
111      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
112            eplus2gg :     Emin=          0 eV         Emax=   10 TeV
113
114msc:   for proton    SubType= 10
115      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
116      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
117      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
118          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
119
120hIoni:   for  proton    SubType= 2
121      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
122      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
123      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
124      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
125               Bragg :     Emin=          0 eV         Emax=   2 MeV
126          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
127
128hBrems:   for  proton    SubType= 3
129      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
130      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
131      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
132               hBrem :     Emin=          0 eV         Emax=   10 TeV
133
134hPairProd:   for  proton    SubType= 4
135      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
136      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
137      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
138           hPairProd :     Emin=          0 eV         Emax=   10 TeV
139
140msc:   for GenericIon    SubType= 10
141      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
142      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
143          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
144
145ionIoni:   for  GenericIon    SubType= 2
146      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
147      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
148      finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
149      Stopping Power data for 17 ion/material pairs, nuclearStopping: 1
150      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
151         ParamICRU73 :     Emin=          0 eV         Emax=   10 TeV
152
153nuclearStopping:   for  GenericIon    SubType= 1
154      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
155   ICRU49NucStopping :     Emin=          0 eV         Emax=   10 TeV
156
157nuclearStopping:   for  alpha    SubType= 1
158      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
159   ICRU49NucStopping :     Emin=          0 eV         Emax=   10 TeV
160
161hIoni:   for  anti_proton    SubType= 2
162      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
163      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
164      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
165      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
166            ICRU73QO :     Emin=          0 eV         Emax=   2 MeV
167          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
168
169hIoni:   for  kaon+    SubType= 2
170      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
171      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
172      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
173      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
174               Bragg :     Emin=          0 eV         Emax=   1.05231 MeV
175          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
176
177hBrems:   for  kaon+    SubType= 3
178      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
179      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
180      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
181               hBrem :     Emin=          0 eV         Emax=   10 TeV
182
183hPairProd:   for  kaon+    SubType= 4
184      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
185      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
186      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
187           hPairProd :     Emin=          0 eV         Emax=   10 TeV
188
189hIoni:   for  kaon-    SubType= 2
190      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
191      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
192      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
193      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
194            ICRU73QO :     Emin=          0 eV         Emax=   1.05231 MeV
195          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
196
197hBrems:   for  kaon-    SubType= 3
198      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
199      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
200      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
201               hBrem :     Emin=          0 eV         Emax=   10 TeV
202
203hPairProd:   for  kaon-    SubType= 4
204      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
205      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
206      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
207           hPairProd :     Emin=          0 eV         Emax=   10 TeV
208
209muMsc:   for mu+    SubType= 10
210      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
211      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592
212      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
213        WentzelVIUni :     Emin=          0 eV         Emax=   10 TeV
214
215muIoni:   for  mu+    SubType= 2
216      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
217      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
218      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
219      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
220               Bragg :     Emin=          0 eV         Emax=   200 keV
221          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
222        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
223
224muBrems:   for  mu+    SubType= 3
225      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
226      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
227      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
228              MuBrem :     Emin=          0 eV         Emax=   10 TeV
229
230muPairProd:   for  mu+    SubType= 4
231      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
232      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
233      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
234          muPairProd :     Emin=          0 eV         Emax=   10 TeV
235
236CoulombScat:   for  mu+    SubType= 1
237      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
238      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
239      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
240  eCoulombScattering :     Emin=          0 eV         Emax=   10 TeV
241
242muIoni:   for  mu-    SubType= 2
243      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
244      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
245      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
246      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
247            ICRU73QO :     Emin=          0 eV         Emax=   200 keV
248          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
249        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
250
251muBrems:   for  mu-    SubType= 3
252      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
253      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
254      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
255              MuBrem :     Emin=          0 eV         Emax=   10 TeV
256
257muPairProd:   for  mu-    SubType= 4
258      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
259      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
260      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
261          muPairProd :     Emin=          0 eV         Emax=   10 TeV
262
263CoulombScat:   for  mu-    SubType= 1
264      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
265      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
266      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
267  eCoulombScattering :     Emin=          0 eV         Emax=   10 TeV
268
269hIoni:   for  pi+    SubType= 2
270      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
271      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
272      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
273      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
274               Bragg :     Emin=          0 eV         Emax=   297.504 keV
275          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
276
277hBrems:   for  pi+    SubType= 3
278      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
279      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
280      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
281               hBrem :     Emin=          0 eV         Emax=   10 TeV
282
283hPairProd:   for  pi+    SubType= 4
284      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
285      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
286      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
287           hPairProd :     Emin=          0 eV         Emax=   10 TeV
288
289msc:   for pi-    SubType= 10
290      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
291      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
292      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
293          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
294
295hIoni:   for  pi-    SubType= 2
296      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
297      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
298      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
299      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
300            ICRU73QO :     Emin=          0 eV         Emax=   297.504 keV
301          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
302
303hBrems:   for  pi-    SubType= 3
304      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
305      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
306      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
307               hBrem :     Emin=          0 eV         Emax=   10 TeV
308
309hPairProd:   for  pi-    SubType= 4
310      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
311      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
312      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
313           hPairProd :     Emin=          0 eV         Emax=   10 TeV
314
315Region <DefaultRegionForTheWorld> -- appears in <World> world volume
316 Root logical volume(s) : World
317 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
318 Materials : Galactic Silicon
319 Production cuts :   gamma 7 um      e- 7 um      e+ 7 um  proton 0 fm
320
321========= Table of registered couples ==============================
322
323Index : 0     used in the geometry : Yes     recalculation needed : No
324 Material : Galactic
325 Range cuts        :  gamma  7 um     e-  7 um     e+  7 um  proton 0 fm
326 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
327 Region(s) which use this couple :
328    DefaultRegionForTheWorld
329
330Index : 1     used in the geometry : Yes     recalculation needed : No
331 Material : Silicon
332 Range cuts        :  gamma  7 um     e-  7 um     e+  7 um  proton 0 fm
333 Energy thresholds :  gamma  990 eV     e-  21.4886 keV    e+  20.8954 keV proton 0 eV
334 Region(s) which use this couple :
335    DefaultRegionForTheWorld
336
337====================================================================
338
339Start closing geometry.
340G4GeometryManager::ReportVoxelStats -- Voxel Statistics
341
342    Total memory consumed for geometry optimisation:   0 kByte
343    Total CPU time elapsed for geometry optimisation: 0 seconds
344### Run 0 start.
345
346--------- Ranecu engine status ---------
347 Initial seed (index) = 0
348 Current couple of seeds = 9876, 54321
349----------------------------------------
350Start Run processing.
351
352---> Begin of Event: 0
353
354---> Begin of Event: 10000
355
356---> Begin of Event: 20000
357
358---> Begin of Event: 30000
359
360---> Begin of Event: 40000
361
362---> Begin of Event: 50000
363
364---> Begin of Event: 60000
365
366---> Begin of Event: 70000
367
368---> Begin of Event: 80000
369
370---> Begin of Event: 90000
371Run terminated.
372Run Summary
373  Number of events processed : 100000
374  User=4.22s Real=11.65s Sys=2.4s
375
376 ======================== run summary ======================
377
378 The run was 100000 pi+ of 5 GeV through 20 um  of Silicon (density: 2.33 g/cm3 )
379
380 Total energy deposit in absorber per event = 6.329 keV +- 15.51 eV
381
382 -----> Mean dE/dx = 3.165 MeV/cm       (1.358 MeV*cm2/g)
383
384 From formulas :
385   restricted dEdx = 2.963 MeV/cm       (1.272 MeV*cm2/g)
386   full dEdx       = 4.726 MeV/cm       (2.028 MeV*cm2/g)
387
388 Leakage :  primary = 5 GeV +- 294.7 eV    secondaries = 1.778 keV +- 293.3 eV
389 Energy balance :  edep + eleak = 5 GeV
390
391 Total track length (charged) in absorber per event = 20.26 um  +- 7.916 nm
392 Total track length (neutral) in absorber per event = 2.089 Ang +- 1.085 Ang
393
394 Number of steps (charged) in absorber per event = 1.156 +- 8.897e-06
395 Number of steps (neutral) in absorber per event = 6e-05 +- 8.897e-06
396
397 Number of secondaries per event : Gammas = 6e-05;   electrons = 0.01646;   positrons = 0
398
399 Number of events with the primary particle transmitted = 100 %
400 Number of events with at least  1 particle transmitted (same charge as primary) = 100 %
401
402 Number of events with the primary particle reflected = 0 %
403 Number of events with at least  1 particle reflected (same charge as primary) = 0 %
404
405 MultipleScattering:
406  rms proj angle of transmit primary particle = 0.02452 mrad (central part only)
407  computed theta0 (Highland formula)          = 0.02627 mrad
408  central part defined as +- 0.0788 mrad;   Tail ratio = 1.573 %
409
410--------- Ranecu engine status ---------
411 Initial seed (index) = 0
412 Current couple of seeds = 206812152, 1002169913
413----------------------------------------
414#
415/testem/det/setAbsMat Gold
416/testem/det/setAbsThick 9.658 um
417/testem/det/setAbsYZ 1 mm
418/testem/det/update
419
420 Material: Galactic     density:  0.000 kg/m3   RadL: 204727512.315 pc   Nucl.Int.Length: 113804112.800 pc   Imean:  19.200 eV   temperature:   2.73 K  pressure:   0.00 atm
421   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
422
423 Material:     Gold     density: 19.320 g/cm3   RadL:   3.344 mm   Nucl.Int.Length:  10.540 cm   Imean: 790.000 eV
424   --->  Element: Gold ( )   Z = 79.0   N = 197.0   A = 196.97 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
425
426 The  WORLD   is made of 14.487 um  of Galactic. The transverse size (YZ) of the world is 1.2 mm
427 The ABSORBER is made of 9.658 um  of Gold. The transverse size (YZ) is 1 mm
428 X position of the middle of the absorber 0 fm
429Obsolete world logical volume is removed from the default region.
430World is registered to the default region.
431/testem/gun/setDefault
432/gun/particle e-
433/gun/energy 15.7 MeV
434/run/beamOn 50000
435
436phot:   for  gamma    SubType= 12
437      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
438       PhotoElectric :     Emin=          0 eV         Emax=   10 TeV
439
440compt:   for  gamma    SubType= 13
441      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
442      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
443       Klein-Nishina :     Emin=          0 eV         Emax=   10 TeV
444
445conv:   for  gamma    SubType= 14
446      Lambda tables from 1.0219978 MeV to 10 TeV in 220 bins, spline: 1
447      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
448       Bethe-Heitler :     Emin=          0 eV         Emax=   10 TeV
449
450Rayl:   for  gamma    SubType= 11
451      Lambda tables from 100 eV  to 10 TeV in 200 bins, spline: 1
452      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
453   LivermoreRayleigh :     Emin=          0 eV         Emax=   100 GeV
454
455msc:   for e-    SubType= 10
456      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
457      RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
458      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
459          UrbanMsc93 :     Emin=          0 eV         Emax=   10 TeV
460
461eIoni:   for  e-    SubType= 2
462      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
463      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
464      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
465      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
466        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
467
468eBrem:   for  e-    SubType= 3
469      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
470      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
471      LPM flag: 1 for E > 1 GeV
472      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
473               eBrem :     Emin=          0 eV         Emax=   1 GeV
474            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
475
476eIoni:   for  e+    SubType= 2
477      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
478      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
479      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
480      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
481        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
482
483eBrem:   for  e+    SubType= 3
484      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
485      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
486      LPM flag: 1 for E > 1 GeV
487      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
488               eBrem :     Emin=          0 eV         Emax=   1 GeV
489            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
490
491annihil:   for  e+    SubType= 5
492      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
493      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
494            eplus2gg :     Emin=          0 eV         Emax=   10 TeV
495
496msc:   for proton    SubType= 10
497      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
498      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
499      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
500          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
501
502hIoni:   for  proton    SubType= 2
503      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
504      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
505      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
506      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
507               Bragg :     Emin=          0 eV         Emax=   2 MeV
508          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
509
510hBrems:   for  proton    SubType= 3
511      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
512      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
513      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
514               hBrem :     Emin=          0 eV         Emax=   10 TeV
515
516hPairProd:   for  proton    SubType= 4
517      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
518      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
519      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
520           hPairProd :     Emin=          0 eV         Emax=   10 TeV
521
522msc:   for GenericIon    SubType= 10
523      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
524      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
525          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
526
527ionIoni:   for  GenericIon    SubType= 2
528      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
529      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
530      finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
531      Stopping Power data for 17 ion/material pairs, nuclearStopping: 1
532      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
533         ParamICRU73 :     Emin=          0 eV         Emax=   10 TeV
534
535nuclearStopping:   for  GenericIon    SubType= 1
536      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
537   ICRU49NucStopping :     Emin=          0 eV         Emax=   10 TeV
538
539nuclearStopping:   for  alpha    SubType= 1
540      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
541   ICRU49NucStopping :     Emin=          0 eV         Emax=   10 TeV
542
543hIoni:   for  anti_proton    SubType= 2
544      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
545      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
546      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
547      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
548            ICRU73QO :     Emin=          0 eV         Emax=   2 MeV
549          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
550
551hIoni:   for  kaon+    SubType= 2
552      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
553      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
554      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
555      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
556               Bragg :     Emin=          0 eV         Emax=   1.052311 MeV
557          BetheBloch :     Emin=   1.052311 MeV        Emax=   10 TeV
558
559hBrems:   for  kaon+    SubType= 3
560      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
561      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
562      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
563               hBrem :     Emin=          0 eV         Emax=   10 TeV
564
565hPairProd:   for  kaon+    SubType= 4
566      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
567      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
568      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
569           hPairProd :     Emin=          0 eV         Emax=   10 TeV
570
571hIoni:   for  kaon-    SubType= 2
572      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
573      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
574      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
575      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
576            ICRU73QO :     Emin=          0 eV         Emax=   1.052311 MeV
577          BetheBloch :     Emin=   1.052311 MeV        Emax=   10 TeV
578
579hBrems:   for  kaon-    SubType= 3
580      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
581      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
582      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
583               hBrem :     Emin=          0 eV         Emax=   10 TeV
584
585hPairProd:   for  kaon-    SubType= 4
586      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
587      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
588      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
589           hPairProd :     Emin=          0 eV         Emax=   10 TeV
590
591muMsc:   for mu+    SubType= 10
592      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
593      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.459156
594      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
595        WentzelVIUni :     Emin=          0 eV         Emax=   10 TeV
596
597muIoni:   for  mu+    SubType= 2
598      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
599      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
600      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
601      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
602               Bragg :     Emin=          0 eV         Emax=   200 keV
603          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
604        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
605
606muBrems:   for  mu+    SubType= 3
607      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
608      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
609      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
610              MuBrem :     Emin=          0 eV         Emax=   10 TeV
611
612muPairProd:   for  mu+    SubType= 4
613      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
614      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
615      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
616          muPairProd :     Emin=          0 eV         Emax=   10 TeV
617
618CoulombScat:   for  mu+    SubType= 1
619      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
620      11.459156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123
621      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
622  eCoulombScattering :     Emin=          0 eV         Emax=   10 TeV
623
624muIoni:   for  mu-    SubType= 2
625      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
626      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
627      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
628      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
629            ICRU73QO :     Emin=          0 eV         Emax=   200 keV
630          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
631        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
632
633muBrems:   for  mu-    SubType= 3
634      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
635      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
636      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
637              MuBrem :     Emin=          0 eV         Emax=   10 TeV
638
639muPairProd:   for  mu-    SubType= 4
640      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
641      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
642      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
643          muPairProd :     Emin=          0 eV         Emax=   10 TeV
644
645CoulombScat:   for  mu-    SubType= 1
646      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
647      11.459156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123
648      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
649  eCoulombScattering :     Emin=          0 eV         Emax=   10 TeV
650
651hIoni:   for  pi+    SubType= 2
652      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
653      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
654      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
655      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
656               Bragg :     Emin=          0 eV         Emax=   297.50434 keV
657          BetheBloch :     Emin=  297.50434 keV        Emax=   10 TeV
658
659hBrems:   for  pi+    SubType= 3
660      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
661      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
662      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
663               hBrem :     Emin=          0 eV         Emax=   10 TeV
664
665hPairProd:   for  pi+    SubType= 4
666      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
667      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
668      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
669           hPairProd :     Emin=          0 eV         Emax=   10 TeV
670
671msc:   for pi-    SubType= 10
672      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
673      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
674      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
675          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
676
677hIoni:   for  pi-    SubType= 2
678      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
679      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
680      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
681      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
682            ICRU73QO :     Emin=          0 eV         Emax=   297.50434 keV
683          BetheBloch :     Emin=  297.50434 keV        Emax=   10 TeV
684
685hBrems:   for  pi-    SubType= 3
686      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
687      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
688      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
689               hBrem :     Emin=          0 eV         Emax=   10 TeV
690
691hPairProd:   for  pi-    SubType= 4
692      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
693      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
694      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
695           hPairProd :     Emin=          0 eV         Emax=   10 TeV
696
697Region <DefaultRegionForTheWorld> -- appears in <World> world volume
698 Root logical volume(s) : World
699 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
700 Materials : Galactic Gold
701 Production cuts :   gamma 7 um      e- 7 um      e+ 7 um  proton 0 fm
702
703========= Table of registered couples ==============================
704
705Index : 0     used in the geometry : Yes     recalculation needed : No
706 Material : Galactic
707 Range cuts        :  gamma  7 um     e-  7 um     e+  7 um  proton 0 fm
708 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
709 Region(s) which use this couple :
710    DefaultRegionForTheWorld
711
712Index : 1     used in the geometry : No      recalculation needed : No
713 Material : Silicon
714 Range cuts        :  gamma  7 um     e-  7 um     e+  7 um  proton 0 fm
715 Energy thresholds :  gamma  990 eV     e-  21.488579 keV    e+  20.895427 keV proton 0 eV
716
717Index : 2     used in the geometry : Yes     recalculation needed : No
718 Material : Gold
719 Range cuts        :  gamma  7 um     e-  7 um     e+  7 um  proton 0 fm
720 Energy thresholds :  gamma  6.5422867 keV    e-  65.140571 keV    e+  63.840474 keV proton 0 eV
721 Region(s) which use this couple :
722    DefaultRegionForTheWorld
723
724====================================================================
725
726Start closing geometry.
727G4GeometryManager::ReportVoxelStats -- Voxel Statistics
728
729    Total memory consumed for geometry optimisation:   0 kByte
730    Total CPU time elapsed for geometry optimisation: 0 seconds
731### Run 1 start.
732
733--------- Ranecu engine status ---------
734 Initial seed (index) = 0
735 Current couple of seeds = 206812152, 1002169913
736----------------------------------------
737Start Run processing.
738
739---> Begin of Event: 0
740
741---> Begin of Event: 10000
742
743---> Begin of Event: 20000
744
745---> Begin of Event: 30000
746
747---> Begin of Event: 40000
748Run terminated.
749Run Summary
750  Number of events processed : 50000
751  User=3.27s Real=4.5s Sys=1.22s
752
753 ======================== run summary ======================
754
755 The run was 50000 e- of 15.7 MeV through 9.66 um  of Gold (density: 19.3 g/cm3 )
756
757 Total energy deposit in absorber per event = 18.6 keV +- 64.68 eV
758
759 -----> Mean dE/dx = 19.26 MeV/cm       (0.9971 MeV*cm2/g)
760
761 From formulas :
762   restricted dEdx = 17.9 MeV/cm        (0.9266 MeV*cm2/g)
763   full dEdx       = 59.7 MeV/cm        (3.09 MeV*cm2/g)
764
765 Leakage :  primary = 15.64 MeV +- 2.154 keV   secondaries = 40.1 keV +- 2.145 keV
766 Energy balance :  edep + eleak = 15.7 MeV
767
768 Total track length (charged) in absorber per event = 9.95 um  +- 12.57 nm
769 Total track length (neutral) in absorber per event = 133.2 nm  +- 6.153 nm
770
771 Number of steps (charged) in absorber per event = 5.071 +- 1.585e-05
772 Number of steps (neutral) in absorber per event = 0.02956 +- 1.585e-05
773
774 Number of secondaries per event : Gammas = 0.02944;   electrons = 0.02198;   positrons = 0
775
776 Number of events with the primary particle transmitted = 99.97 %
777 Number of events with at least  1 particle transmitted (same charge as primary) = 99.97 %
778
779 Number of events with the primary particle reflected = 0.028 %
780 Number of events with at least  1 particle reflected (same charge as primary) = 0.276 %
781
782 MultipleScattering:
783  rms proj angle of transmit primary particle = 32.92 mrad (central part only)
784  computed theta0 (Highland formula)          = 35.1 mrad
785  central part defined as +- 105.3 mrad;   Tail ratio = 2.702 %
786
787--------- Ranecu engine status ---------
788 Initial seed (index) = 0
789 Current couple of seeds = 1822484465, 1275230702
790----------------------------------------
791#
792/testem/phys/setCuts 10 um
793#
794/testem/det/setAbsMat G4_POLYETHYLENE
795/testem/det/setAbsThick 4.74 mm
796/testem/det/setAbsYZ 50 mm
797/testem/det/update
798
799 Material: Galactic     density:  0.000 kg/m3   RadL: 204727512.315 pc   Nucl.Int.Length: 113804112.800 pc   Imean:  19.200 eV   temperature:   2.73 K  pressure:   0.00 atm
800   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
801
802 Material: G4_POLYETHYLENE (C_2H_4)_N-Polyethylene   density: 940.000 mg/cm3  RadL:  47.632 cm   Nucl.Int.Length:  72.090 cm   Imean:  57.400 eV
803   --->  Element: H (H)   Z =  1.0   N =   1.0   A =   1.01 g/mole
804   --->  Isotope:    H1   Z =  1   N =   1   A =   1.01 g/mole   abundance:  99.99 %
805   --->  Isotope:    H2   Z =  1   N =   2   A =   2.01 g/mole   abundance:   0.01 %  ElmMassFraction:  14.37 %  ElmAbundance  66.67 %
806   --->  Element: C (C)   Z =  6.0   N =  12.0   A =  12.01 g/mole
807   --->  Isotope:   C12   Z =  6   N =  12   A =  12.00 g/mole   abundance:  98.93 %
808   --->  Isotope:   C13   Z =  6   N =  13   A =  13.00 g/mole   abundance:   1.07 %  ElmMassFraction:  85.63 %  ElmAbundance  33.33 %
809
810 The  WORLD   is made of 7.11 mm  of Galactic. The transverse size (YZ) of the world is 6 cm
811 The ABSORBER is made of 4.74 mm  of G4_POLYETHYLENE. The transverse size (YZ) is 5 cm
812 X position of the middle of the absorber 0 fm
813Obsolete world logical volume is removed from the default region.
814World is registered to the default region.
815/testem/gun/setDefault
816/gun/particle mu+
817/gun/energy 96.2 MeV
818/run/beamOn 50000
819
820phot:   for  gamma    SubType= 12
821      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
822       PhotoElectric :     Emin=          0 eV         Emax=   10 TeV
823
824compt:   for  gamma    SubType= 13
825      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
826      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
827       Klein-Nishina :     Emin=          0 eV         Emax=   10 TeV
828
829conv:   for  gamma    SubType= 14
830      Lambda tables from 1.0219978 MeV to 10 TeV in 220 bins, spline: 1
831      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
832       Bethe-Heitler :     Emin=          0 eV         Emax=   10 TeV
833
834Rayl:   for  gamma    SubType= 11
835      Lambda tables from 100 eV  to 10 TeV in 200 bins, spline: 1
836      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
837   LivermoreRayleigh :     Emin=          0 eV         Emax=   100 GeV
838
839msc:   for e-    SubType= 10
840      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
841      RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
842      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
843          UrbanMsc93 :     Emin=          0 eV         Emax=   10 TeV
844
845eIoni:   for  e-    SubType= 2
846      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
847      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
848      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
849      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
850        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
851
852eBrem:   for  e-    SubType= 3
853      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
854      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
855      LPM flag: 1 for E > 1 GeV
856      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
857               eBrem :     Emin=          0 eV         Emax=   1 GeV
858            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
859
860eIoni:   for  e+    SubType= 2
861      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
862      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
863      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
864      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
865        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
866
867eBrem:   for  e+    SubType= 3
868      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
869      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
870      LPM flag: 1 for E > 1 GeV
871      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
872               eBrem :     Emin=          0 eV         Emax=   1 GeV
873            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
874
875annihil:   for  e+    SubType= 5
876      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
877      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
878            eplus2gg :     Emin=          0 eV         Emax=   10 TeV
879
880msc:   for proton    SubType= 10
881      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
882      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
883      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
884          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
885
886hIoni:   for  proton    SubType= 2
887      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
888      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
889      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
890      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
891               Bragg :     Emin=          0 eV         Emax=   2 MeV
892          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
893
894hBrems:   for  proton    SubType= 3
895      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
896      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
897      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
898               hBrem :     Emin=          0 eV         Emax=   10 TeV
899
900hPairProd:   for  proton    SubType= 4
901      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
902      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
903      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
904           hPairProd :     Emin=          0 eV         Emax=   10 TeV
905
906msc:   for GenericIon    SubType= 10
907      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
908      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
909          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
910
911ionIoni:   for  GenericIon    SubType= 2
912      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
913      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
914      finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
915      Stopping Power data for 17 ion/material pairs, nuclearStopping: 1
916      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
917         ParamICRU73 :     Emin=          0 eV         Emax=   10 TeV
918
919nuclearStopping:   for  GenericIon    SubType= 1
920      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
921   ICRU49NucStopping :     Emin=          0 eV         Emax=   10 TeV
922
923nuclearStopping:   for  alpha    SubType= 1
924      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
925   ICRU49NucStopping :     Emin=          0 eV         Emax=   10 TeV
926
927hIoni:   for  anti_proton    SubType= 2
928      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
929      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
930      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
931      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
932            ICRU73QO :     Emin=          0 eV         Emax=   2 MeV
933          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
934
935hIoni:   for  kaon+    SubType= 2
936      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
937      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
938      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
939      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
940               Bragg :     Emin=          0 eV         Emax=   1.052311 MeV
941          BetheBloch :     Emin=   1.052311 MeV        Emax=   10 TeV
942
943hBrems:   for  kaon+    SubType= 3
944      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
945      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
946      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
947               hBrem :     Emin=          0 eV         Emax=   10 TeV
948
949hPairProd:   for  kaon+    SubType= 4
950      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
951      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
952      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
953           hPairProd :     Emin=          0 eV         Emax=   10 TeV
954
955hIoni:   for  kaon-    SubType= 2
956      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
957      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
958      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
959      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
960            ICRU73QO :     Emin=          0 eV         Emax=   1.052311 MeV
961          BetheBloch :     Emin=   1.052311 MeV        Emax=   10 TeV
962
963hBrems:   for  kaon-    SubType= 3
964      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
965      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
966      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
967               hBrem :     Emin=          0 eV         Emax=   10 TeV
968
969hPairProd:   for  kaon-    SubType= 4
970      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
971      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
972      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
973           hPairProd :     Emin=          0 eV         Emax=   10 TeV
974
975muMsc:   for mu+    SubType= 10
976      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
977      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.459156
978      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
979        WentzelVIUni :     Emin=          0 eV         Emax=   10 TeV
980
981muIoni:   for  mu+    SubType= 2
982      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
983      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
984      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
985      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
986               Bragg :     Emin=          0 eV         Emax=   200 keV
987          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
988        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
989
990muBrems:   for  mu+    SubType= 3
991      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
992      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
993      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
994              MuBrem :     Emin=          0 eV         Emax=   10 TeV
995
996muPairProd:   for  mu+    SubType= 4
997      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
998      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
999      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1000          muPairProd :     Emin=          0 eV         Emax=   10 TeV
1001
1002CoulombScat:   for  mu+    SubType= 1
1003      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1004      11.459156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123
1005      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1006  eCoulombScattering :     Emin=          0 eV         Emax=   10 TeV
1007
1008muIoni:   for  mu-    SubType= 2
1009      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1010      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1011      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1012      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1013            ICRU73QO :     Emin=          0 eV         Emax=   200 keV
1014          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
1015        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
1016
1017muBrems:   for  mu-    SubType= 3
1018      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1019      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1020      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1021              MuBrem :     Emin=          0 eV         Emax=   10 TeV
1022
1023muPairProd:   for  mu-    SubType= 4
1024      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1025      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1026      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1027          muPairProd :     Emin=          0 eV         Emax=   10 TeV
1028
1029CoulombScat:   for  mu-    SubType= 1
1030      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1031      11.459156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123
1032      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1033  eCoulombScattering :     Emin=          0 eV         Emax=   10 TeV
1034
1035hIoni:   for  pi+    SubType= 2
1036      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1037      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1038      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1039      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1040               Bragg :     Emin=          0 eV         Emax=   297.50434 keV
1041          BetheBloch :     Emin=  297.50434 keV        Emax=   10 TeV
1042
1043hBrems:   for  pi+    SubType= 3
1044      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1045      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1046      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1047               hBrem :     Emin=          0 eV         Emax=   10 TeV
1048
1049hPairProd:   for  pi+    SubType= 4
1050      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1051      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1052      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1053           hPairProd :     Emin=          0 eV         Emax=   10 TeV
1054
1055msc:   for pi-    SubType= 10
1056      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1057      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1058      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1059          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
1060
1061hIoni:   for  pi-    SubType= 2
1062      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1063      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1064      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1065      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1066            ICRU73QO :     Emin=          0 eV         Emax=   297.50434 keV
1067          BetheBloch :     Emin=  297.50434 keV        Emax=   10 TeV
1068
1069hBrems:   for  pi-    SubType= 3
1070      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1071      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1072      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1073               hBrem :     Emin=          0 eV         Emax=   10 TeV
1074
1075hPairProd:   for  pi-    SubType= 4
1076      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1077      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1078      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1079           hPairProd :     Emin=          0 eV         Emax=   10 TeV
1080
1081Region <DefaultRegionForTheWorld> -- appears in <World> world volume
1082 Root logical volume(s) : World
1083 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
1084 Materials : Galactic G4_POLYETHYLENE
1085 Production cuts :   gamma 10 um      e- 10 um      e+ 10 um  proton 0 fm
1086
1087========= Table of registered couples ==============================
1088
1089Index : 0     used in the geometry : Yes     recalculation needed : No
1090 Material : Galactic
1091 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
1092 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
1093 Region(s) which use this couple :
1094    DefaultRegionForTheWorld
1095
1096Index : 1     used in the geometry : No      recalculation needed : No
1097 Material : Silicon
1098 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
1099 Energy thresholds :  gamma  990 eV     e-  31.951563 keV    e+  31.737567 keV proton 0 eV
1100
1101Index : 2     used in the geometry : No      recalculation needed : No
1102 Material : Gold
1103 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
1104 Energy thresholds :  gamma  8.2770687 keV    e-  80.768699 keV    e+  79.156691 keV proton 0 eV
1105
1106Index : 3     used in the geometry : Yes     recalculation needed : No
1107 Material : G4_POLYETHYLENE
1108 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
1109 Energy thresholds :  gamma  990 eV     e-  14.174755 keV    e+  13.754766 keV proton 0 eV
1110 Region(s) which use this couple :
1111    DefaultRegionForTheWorld
1112
1113====================================================================
1114
1115Start closing geometry.
1116G4GeometryManager::ReportVoxelStats -- Voxel Statistics
1117
1118    Total memory consumed for geometry optimisation:   0 kByte
1119    Total CPU time elapsed for geometry optimisation: 0 seconds
1120### Run 2 start.
1121
1122--------- Ranecu engine status ---------
1123 Initial seed (index) = 0
1124 Current couple of seeds = 1822484465, 1275230702
1125----------------------------------------
1126Start Run processing.
1127
1128---> Begin of Event: 0
1129
1130---> Begin of Event: 10000
1131
1132---> Begin of Event: 20000
1133
1134---> Begin of Event: 30000
1135
1136---> Begin of Event: 40000
1137Run terminated.
1138Run Summary
1139  Number of events processed : 50000
1140  User=6.73s Real=9.07s Sys=1.29s
1141
1142 ======================== run summary ======================
1143
1144 The run was 50000 mu+ of 96.2 MeV through 4.74 mm  of G4_POLYETHYLENE (density: 940 mg/cm3)
1145
1146 Total energy deposit in absorber per event = 1.033 MeV +- 1.01 keV
1147
1148 -----> Mean dE/dx = 2.18 MeV/cm        (2.319 MeV*cm2/g)
1149
1150 From formulas :
1151   restricted dEdx = 1.751 MeV/cm       (1.863 MeV*cm2/g)
1152   full dEdx       = 2.263 MeV/cm       (2.407 MeV*cm2/g)
1153
1154 Leakage :  primary = 95.14 MeV +- 2.335 keV   secondaries = 30.95 keV +- 3.642 keV
1155 Energy balance :  edep + eleak = 96.2 MeV
1156
1157 Total track length (charged) in absorber per event = 5.106 mm  +- 3.914 um
1158 Total track length (neutral) in absorber per event = 31.28 um  +- 2.947 um
1159
1160 Number of steps (charged) in absorber per event = 11.22 +- 0.0002798
1161 Number of steps (neutral) in absorber per event = 0.01258 +- 0.0002798
1162
1163 Number of secondaries per event : Gammas = 0.01166;   electrons = 3.988;   positrons = 2e-05
1164
1165 Number of events with the primary particle transmitted = 100 %
1166 Number of events with at least  1 particle transmitted (same charge as primary) = 100 %
1167
1168 Number of events with the primary particle reflected = 0 %
1169 Number of events with at least  1 particle reflected (same charge as primary) = 0 %
1170
1171 MultipleScattering:
1172  rms proj angle of transmit primary particle = 7.29 mrad (central part only)
1173  computed theta0 (Highland formula)          = 7.636 mrad
1174  central part defined as +- 22.91 mrad;   Tail ratio = 1.284 %
1175
1176--------- Ranecu engine status ---------
1177 Initial seed (index) = 0
1178 Current couple of seeds = 1953603820, 877017596
1179----------------------------------------
1180#
1181/testem/phys/setCuts 1 mm
1182#
1183/testem/stack/killSecondaries
1184#
1185/testem/det/setAbsMat Aluminium
1186/testem/det/setAbsThick 8.004 mm
1187/testem/det/update
1188
1189 Material: Galactic     density:  0.000 kg/m3   RadL: 204727512.315 pc   Nucl.Int.Length: 113804112.800 pc   Imean:  19.200 eV   temperature:   2.73 K  pressure:   0.00 atm
1190   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
1191
1192 Material: Aluminium     density:  2.700 g/cm3   RadL:   8.893 cm   Nucl.Int.Length:  38.879 cm   Imean: 166.000 eV
1193   --->  Element: Aluminium ( )   Z = 13.0   N =  27.0   A =  26.98 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
1194
1195 The  WORLD   is made of 1.2006 cm  of Galactic. The transverse size (YZ) of the world is 6 cm
1196 The ABSORBER is made of 8.004 mm  of Aluminium. The transverse size (YZ) is 5 cm
1197 X position of the middle of the absorber 0 fm
1198Obsolete world logical volume is removed from the default region.
1199World is registered to the default region.
1200/testem/gun/setDefault
1201/gun/particle proton
1202/gun/energy 174.664 GeV
1203/run/beamOn 50000
1204
1205phot:   for  gamma    SubType= 12
1206      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1207       PhotoElectric :     Emin=          0 eV         Emax=   10 TeV
1208
1209compt:   for  gamma    SubType= 13
1210      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1211      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1212       Klein-Nishina :     Emin=          0 eV         Emax=   10 TeV
1213
1214conv:   for  gamma    SubType= 14
1215      Lambda tables from 1.0219978 MeV to 10 TeV in 220 bins, spline: 1
1216      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1217       Bethe-Heitler :     Emin=          0 eV         Emax=   10 TeV
1218
1219Rayl:   for  gamma    SubType= 11
1220      Lambda tables from 100 eV  to 10 TeV in 200 bins, spline: 1
1221      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1222   LivermoreRayleigh :     Emin=          0 eV         Emax=   100 GeV
1223
1224msc:   for e-    SubType= 10
1225      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1226      RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
1227      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1228          UrbanMsc93 :     Emin=          0 eV         Emax=   10 TeV
1229
1230eIoni:   for  e-    SubType= 2
1231      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1232      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1233      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1234      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1235        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
1236
1237eBrem:   for  e-    SubType= 3
1238      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1239      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1240      LPM flag: 1 for E > 1 GeV
1241      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1242               eBrem :     Emin=          0 eV         Emax=   1 GeV
1243            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
1244
1245eIoni:   for  e+    SubType= 2
1246      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1247      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1248      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1249      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1250        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
1251
1252eBrem:   for  e+    SubType= 3
1253      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1254      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1255      LPM flag: 1 for E > 1 GeV
1256      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1257               eBrem :     Emin=          0 eV         Emax=   1 GeV
1258            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
1259
1260annihil:   for  e+    SubType= 5
1261      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1262      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1263            eplus2gg :     Emin=          0 eV         Emax=   10 TeV
1264
1265msc:   for proton    SubType= 10
1266      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1267      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1268      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1269          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
1270
1271hIoni:   for  proton    SubType= 2
1272      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1273      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1274      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1275      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1276               Bragg :     Emin=          0 eV         Emax=   2 MeV
1277          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
1278
1279hBrems:   for  proton    SubType= 3
1280      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1281      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1282      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1283               hBrem :     Emin=          0 eV         Emax=   10 TeV
1284
1285hPairProd:   for  proton    SubType= 4
1286      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1287      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1288      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1289           hPairProd :     Emin=          0 eV         Emax=   10 TeV
1290
1291msc:   for GenericIon    SubType= 10
1292      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
1293      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1294          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
1295
1296ionIoni:   for  GenericIon    SubType= 2
1297      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1298      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1299      finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
1300      Stopping Power data for 17 ion/material pairs, nuclearStopping: 1
1301      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1302         ParamICRU73 :     Emin=          0 eV         Emax=   10 TeV
1303
1304nuclearStopping:   for  GenericIon    SubType= 1
1305      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1306   ICRU49NucStopping :     Emin=          0 eV         Emax=   10 TeV
1307
1308nuclearStopping:   for  alpha    SubType= 1
1309      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1310   ICRU49NucStopping :     Emin=          0 eV         Emax=   10 TeV
1311
1312hIoni:   for  anti_proton    SubType= 2
1313      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1314      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1315      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1316      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1317            ICRU73QO :     Emin=          0 eV         Emax=   2 MeV
1318          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
1319
1320hIoni:   for  kaon+    SubType= 2
1321      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1322      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1323      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1324      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1325               Bragg :     Emin=          0 eV         Emax=   1.052311 MeV
1326          BetheBloch :     Emin=   1.052311 MeV        Emax=   10 TeV
1327
1328hBrems:   for  kaon+    SubType= 3
1329      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1330      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1331      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1332               hBrem :     Emin=          0 eV         Emax=   10 TeV
1333
1334hPairProd:   for  kaon+    SubType= 4
1335      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1336      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1337      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1338           hPairProd :     Emin=          0 eV         Emax=   10 TeV
1339
1340hIoni:   for  kaon-    SubType= 2
1341      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1342      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1343      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1344      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1345            ICRU73QO :     Emin=          0 eV         Emax=   1.052311 MeV
1346          BetheBloch :     Emin=   1.052311 MeV        Emax=   10 TeV
1347
1348hBrems:   for  kaon-    SubType= 3
1349      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1350      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1351      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1352               hBrem :     Emin=          0 eV         Emax=   10 TeV
1353
1354hPairProd:   for  kaon-    SubType= 4
1355      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1356      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1357      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1358           hPairProd :     Emin=          0 eV         Emax=   10 TeV
1359
1360muMsc:   for mu+    SubType= 10
1361      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1362      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.459156
1363      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1364        WentzelVIUni :     Emin=          0 eV         Emax=   10 TeV
1365
1366muIoni:   for  mu+    SubType= 2
1367      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1368      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1369      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1370      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1371               Bragg :     Emin=          0 eV         Emax=   200 keV
1372          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
1373        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
1374
1375muBrems:   for  mu+    SubType= 3
1376      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1377      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1378      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1379              MuBrem :     Emin=          0 eV         Emax=   10 TeV
1380
1381muPairProd:   for  mu+    SubType= 4
1382      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1383      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1384      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1385          muPairProd :     Emin=          0 eV         Emax=   10 TeV
1386
1387CoulombScat:   for  mu+    SubType= 1
1388      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1389      11.459156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123
1390      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1391  eCoulombScattering :     Emin=          0 eV         Emax=   10 TeV
1392
1393muIoni:   for  mu-    SubType= 2
1394      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1395      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1396      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1397      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1398            ICRU73QO :     Emin=          0 eV         Emax=   200 keV
1399          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
1400        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
1401
1402muBrems:   for  mu-    SubType= 3
1403      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1404      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1405      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1406              MuBrem :     Emin=          0 eV         Emax=   10 TeV
1407
1408muPairProd:   for  mu-    SubType= 4
1409      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1410      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1411      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1412          muPairProd :     Emin=          0 eV         Emax=   10 TeV
1413
1414CoulombScat:   for  mu-    SubType= 1
1415      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1416      11.459156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123
1417      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1418  eCoulombScattering :     Emin=          0 eV         Emax=   10 TeV
1419
1420hIoni:   for  pi+    SubType= 2
1421      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1422      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1423      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1424      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1425               Bragg :     Emin=          0 eV         Emax=   297.50434 keV
1426          BetheBloch :     Emin=  297.50434 keV        Emax=   10 TeV
1427
1428hBrems:   for  pi+    SubType= 3
1429      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1430      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1431      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1432               hBrem :     Emin=          0 eV         Emax=   10 TeV
1433
1434hPairProd:   for  pi+    SubType= 4
1435      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1436      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1437      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1438           hPairProd :     Emin=          0 eV         Emax=   10 TeV
1439
1440msc:   for pi-    SubType= 10
1441      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1442      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1443      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1444          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
1445
1446hIoni:   for  pi-    SubType= 2
1447      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1448      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1449      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1450      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1451            ICRU73QO :     Emin=          0 eV         Emax=   297.50434 keV
1452          BetheBloch :     Emin=  297.50434 keV        Emax=   10 TeV
1453
1454hBrems:   for  pi-    SubType= 3
1455      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1456      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1457      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1458               hBrem :     Emin=          0 eV         Emax=   10 TeV
1459
1460hPairProd:   for  pi-    SubType= 4
1461      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1462      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1463      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1464           hPairProd :     Emin=          0 eV         Emax=   10 TeV
1465
1466Region <DefaultRegionForTheWorld> -- appears in <World> world volume
1467 Root logical volume(s) : World
1468 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
1469 Materials : Galactic Aluminium
1470 Production cuts :   gamma 1 mm      e- 1 mm      e+ 1 mm  proton 0 fm
1471
1472========= Table of registered couples ==============================
1473
1474Index : 0     used in the geometry : Yes     recalculation needed : No
1475 Material : Galactic
1476 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1477 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
1478 Region(s) which use this couple :
1479    DefaultRegionForTheWorld
1480
1481Index : 1     used in the geometry : No      recalculation needed : No
1482 Material : Silicon
1483 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1484 Energy thresholds :  gamma  6.9501782 keV    e-  548.29137 keV    e+  526.62382 keV proton 0 eV
1485
1486Index : 2     used in the geometry : No      recalculation needed : No
1487 Material : Gold
1488 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1489 Energy thresholds :  gamma  112.64401 keV    e-  2.3097622 MeV    e+  2.1308135 MeV proton 0 eV
1490
1491Index : 3     used in the geometry : No      recalculation needed : No
1492 Material : G4_POLYETHYLENE
1493 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1494 Energy thresholds :  gamma  2.2626058 keV    e-  351.8774 keV    e+  342.54485 keV proton 0 eV
1495
1496Index : 4     used in the geometry : Yes     recalculation needed : No
1497 Material : Aluminium
1498 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1499 Energy thresholds :  gamma  6.9036292 keV    e-  598.34507 keV    e+  570.85042 keV proton 0 eV
1500 Region(s) which use this couple :
1501    DefaultRegionForTheWorld
1502
1503====================================================================
1504
1505Start closing geometry.
1506G4GeometryManager::ReportVoxelStats -- Voxel Statistics
1507
1508    Total memory consumed for geometry optimisation:   0 kByte
1509    Total CPU time elapsed for geometry optimisation: 0 seconds
1510### Run 3 start.
1511
1512--------- Ranecu engine status ---------
1513 Initial seed (index) = 0
1514 Current couple of seeds = 1953603820, 877017596
1515----------------------------------------
1516Start Run processing.
1517
1518---> Begin of Event: 0
1519
1520---> Begin of Event: 10000
1521
1522---> Begin of Event: 20000
1523
1524---> Begin of Event: 30000
1525
1526---> Begin of Event: 40000
1527Run terminated.
1528Run Summary
1529  Number of events processed : 50000
1530  User=2.22s Real=3.49s Sys=1.25s
1531
1532 ======================== run summary ======================
1533
1534 The run was 50000 proton of 175 GeV through 8 mm  of Aluminium (density: 2.7 g/cm3 )
1535
1536 Total energy deposit in absorber per event = 4.999 MeV +- 290.4 keV
1537
1538 -----> Mean dE/dx = 6.245 MeV/cm       (2.313 MeV*cm2/g)
1539
1540 From formulas :
1541   restricted dEdx = 4.024 MeV/cm       (1.491 MeV*cm2/g)
1542   full dEdx       = 6.049 MeV/cm       (2.24 MeV*cm2/g)
1543
1544 Leakage :  primary = 174.7 GeV +- 290.4 keV   secondaries = 0 eV  +- 0 eV
1545 Energy balance :  edep + eleak = 174.7 GeV
1546
1547 Total track length (charged) in absorber per event = 8.004 mm  +- 0 fm
1548 Total track length (neutral) in absorber per event = 0 fm  +- 0 fm
1549
1550 Number of steps (charged) in absorber per event = 1.268 +- 0
1551 Number of steps (neutral) in absorber per event = 0 +- 0
1552
1553 Number of secondaries per event : Gammas = 0;   electrons = 0.268;   positrons = 0.0009
1554
1555 Number of events with the primary particle transmitted = 100 %
1556 Number of events with at least  1 particle transmitted (same charge as primary) = 100 %
1557
1558 Number of events with the primary particle reflected = 0 %
1559 Number of events with at least  1 particle reflected (same charge as primary) = 0 %
1560
1561 MultipleScattering:
1562  rms proj angle of transmit primary particle = 0.02108 mrad (central part only)
1563  computed theta0 (Highland formula)          = 0.02111 mrad
1564  central part defined as +- 0.06333 mrad;   Tail ratio = 1.836 %
1565
1566--------- Ranecu engine status ---------
1567 Initial seed (index) = 0
1568 Current couple of seeds = 1162542043, 318101973
1569----------------------------------------
1570#
1571/testem/det/setAbsMat Iron
1572/testem/det/setAbsThick 1 m
1573/testem/det/update
1574
1575 Material: Galactic     density:  0.000 kg/m3   RadL: 204727512.315 pc   Nucl.Int.Length: 113804112.800 pc   Imean:  19.200 eV   temperature:   2.73 K  pressure:   0.00 atm
1576   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
1577
1578 Material:     Iron     density:  7.870 g/cm3   RadL:   1.759 cm   Nucl.Int.Length:  16.999 cm   Imean: 286.000 eV
1579   --->  Element: Iron ( )   Z = 26.0   N =  55.8   A =  55.85 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
1580
1581 The  WORLD   is made of 1.5 m   of Galactic. The transverse size (YZ) of the world is 6 cm
1582 The ABSORBER is made of 1 m   of Iron. The transverse size (YZ) is 5 cm
1583 X position of the middle of the absorber 0 fm
1584Obsolete world logical volume is removed from the default region.
1585World is registered to the default region.
1586/testem/gun/setDefault
1587/gun/particle mu+
1588/gun/energy 100 GeV
1589/run/beamOn 50000
1590
1591phot:   for  gamma    SubType= 12
1592      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1593       PhotoElectric :     Emin=          0 eV         Emax=   10 TeV
1594
1595compt:   for  gamma    SubType= 13
1596      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1597      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1598       Klein-Nishina :     Emin=          0 eV         Emax=   10 TeV
1599
1600conv:   for  gamma    SubType= 14
1601      Lambda tables from 1.0219978 MeV to 10 TeV in 220 bins, spline: 1
1602      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1603       Bethe-Heitler :     Emin=          0 eV         Emax=   10 TeV
1604
1605Rayl:   for  gamma    SubType= 11
1606      Lambda tables from 100 eV  to 10 TeV in 200 bins, spline: 1
1607      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1608   LivermoreRayleigh :     Emin=          0 eV         Emax=   100 GeV
1609
1610msc:   for e-    SubType= 10
1611      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1612      RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
1613      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1614          UrbanMsc93 :     Emin=          0 eV         Emax=   10 TeV
1615
1616eIoni:   for  e-    SubType= 2
1617      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1618      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1619      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1620      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1621        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
1622
1623eBrem:   for  e-    SubType= 3
1624      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1625      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1626      LPM flag: 1 for E > 1 GeV
1627      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1628               eBrem :     Emin=          0 eV         Emax=   1 GeV
1629            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
1630
1631eIoni:   for  e+    SubType= 2
1632      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1633      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1634      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1635      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1636        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
1637
1638eBrem:   for  e+    SubType= 3
1639      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1640      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1641      LPM flag: 1 for E > 1 GeV
1642      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1643               eBrem :     Emin=          0 eV         Emax=   1 GeV
1644            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
1645
1646annihil:   for  e+    SubType= 5
1647      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1648      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1649            eplus2gg :     Emin=          0 eV         Emax=   10 TeV
1650
1651msc:   for proton    SubType= 10
1652      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1653      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1654      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1655          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
1656
1657hIoni:   for  proton    SubType= 2
1658      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1659      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1660      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1661      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1662               Bragg :     Emin=          0 eV         Emax=   2 MeV
1663          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
1664
1665hBrems:   for  proton    SubType= 3
1666      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1667      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1668      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1669               hBrem :     Emin=          0 eV         Emax=   10 TeV
1670
1671hPairProd:   for  proton    SubType= 4
1672      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1673      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1674      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1675           hPairProd :     Emin=          0 eV         Emax=   10 TeV
1676
1677msc:   for GenericIon    SubType= 10
1678      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
1679      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1680          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
1681
1682ionIoni:   for  GenericIon    SubType= 2
1683      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1684      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1685      finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
1686      Stopping Power data for 17 ion/material pairs, nuclearStopping: 1
1687      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1688         ParamICRU73 :     Emin=          0 eV         Emax=   10 TeV
1689
1690nuclearStopping:   for  GenericIon    SubType= 1
1691      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1692   ICRU49NucStopping :     Emin=          0 eV         Emax=   10 TeV
1693
1694nuclearStopping:   for  alpha    SubType= 1
1695      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1696   ICRU49NucStopping :     Emin=          0 eV         Emax=   10 TeV
1697
1698hIoni:   for  anti_proton    SubType= 2
1699      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1700      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1701      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1702      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1703            ICRU73QO :     Emin=          0 eV         Emax=   2 MeV
1704          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
1705
1706hIoni:   for  kaon+    SubType= 2
1707      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1708      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1709      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1710      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1711               Bragg :     Emin=          0 eV         Emax=   1.052311 MeV
1712          BetheBloch :     Emin=   1.052311 MeV        Emax=   10 TeV
1713
1714hBrems:   for  kaon+    SubType= 3
1715      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1716      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1717      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1718               hBrem :     Emin=          0 eV         Emax=   10 TeV
1719
1720hPairProd:   for  kaon+    SubType= 4
1721      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1722      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1723      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1724           hPairProd :     Emin=          0 eV         Emax=   10 TeV
1725
1726hIoni:   for  kaon-    SubType= 2
1727      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1728      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1729      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1730      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1731            ICRU73QO :     Emin=          0 eV         Emax=   1.052311 MeV
1732          BetheBloch :     Emin=   1.052311 MeV        Emax=   10 TeV
1733
1734hBrems:   for  kaon-    SubType= 3
1735      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1736      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1737      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1738               hBrem :     Emin=          0 eV         Emax=   10 TeV
1739
1740hPairProd:   for  kaon-    SubType= 4
1741      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1742      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1743      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1744           hPairProd :     Emin=          0 eV         Emax=   10 TeV
1745
1746muMsc:   for mu+    SubType= 10
1747      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1748      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.459156
1749      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1750        WentzelVIUni :     Emin=          0 eV         Emax=   10 TeV
1751
1752muIoni:   for  mu+    SubType= 2
1753      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1754      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1755      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1756      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1757               Bragg :     Emin=          0 eV         Emax=   200 keV
1758          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
1759        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
1760
1761muBrems:   for  mu+    SubType= 3
1762      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1763      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1764      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1765              MuBrem :     Emin=          0 eV         Emax=   10 TeV
1766
1767muPairProd:   for  mu+    SubType= 4
1768      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1769      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1770      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1771          muPairProd :     Emin=          0 eV         Emax=   10 TeV
1772
1773CoulombScat:   for  mu+    SubType= 1
1774      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1775      11.459156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123
1776      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1777  eCoulombScattering :     Emin=          0 eV         Emax=   10 TeV
1778
1779muIoni:   for  mu-    SubType= 2
1780      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1781      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1782      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1783      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1784            ICRU73QO :     Emin=          0 eV         Emax=   200 keV
1785          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
1786        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
1787
1788muBrems:   for  mu-    SubType= 3
1789      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1790      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1791      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1792              MuBrem :     Emin=          0 eV         Emax=   10 TeV
1793
1794muPairProd:   for  mu-    SubType= 4
1795      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1796      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1797      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1798          muPairProd :     Emin=          0 eV         Emax=   10 TeV
1799
1800CoulombScat:   for  mu-    SubType= 1
1801      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1802      11.459156 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.13953123
1803      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1804  eCoulombScattering :     Emin=          0 eV         Emax=   10 TeV
1805
1806hIoni:   for  pi+    SubType= 2
1807      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1808      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1809      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1810      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1811               Bragg :     Emin=          0 eV         Emax=   297.50434 keV
1812          BetheBloch :     Emin=  297.50434 keV        Emax=   10 TeV
1813
1814hBrems:   for  pi+    SubType= 3
1815      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1816      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1817      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1818               hBrem :     Emin=          0 eV         Emax=   10 TeV
1819
1820hPairProd:   for  pi+    SubType= 4
1821      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1822      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1823      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1824           hPairProd :     Emin=          0 eV         Emax=   10 TeV
1825
1826msc:   for pi-    SubType= 10
1827      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1828      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1829      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1830          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
1831
1832hIoni:   for  pi-    SubType= 2
1833      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1834      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1835      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1836      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1837            ICRU73QO :     Emin=          0 eV         Emax=   297.50434 keV
1838          BetheBloch :     Emin=  297.50434 keV        Emax=   10 TeV
1839
1840hBrems:   for  pi-    SubType= 3
1841      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1842      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1843      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1844               hBrem :     Emin=          0 eV         Emax=   10 TeV
1845
1846hPairProd:   for  pi-    SubType= 4
1847      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1848      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1849      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1850           hPairProd :     Emin=          0 eV         Emax=   10 TeV
1851
1852Region <DefaultRegionForTheWorld> -- appears in <World> world volume
1853 Root logical volume(s) : World
1854 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
1855 Materials : Galactic Iron
1856 Production cuts :   gamma 1 mm      e- 1 mm      e+ 1 mm  proton 0 fm
1857
1858========= Table of registered couples ==============================
1859
1860Index : 0     used in the geometry : Yes     recalculation needed : No
1861 Material : Galactic
1862 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1863 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
1864 Region(s) which use this couple :
1865    DefaultRegionForTheWorld
1866
1867Index : 1     used in the geometry : No      recalculation needed : No
1868 Material : Silicon
1869 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1870 Energy thresholds :  gamma  6.9501782 keV    e-  548.29137 keV    e+  526.62382 keV proton 0 eV
1871
1872Index : 2     used in the geometry : No      recalculation needed : No
1873 Material : Gold
1874 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1875 Energy thresholds :  gamma  112.64401 keV    e-  2.3097622 MeV    e+  2.1308135 MeV proton 0 eV
1876
1877Index : 3     used in the geometry : No      recalculation needed : No
1878 Material : G4_POLYETHYLENE
1879 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1880 Energy thresholds :  gamma  2.2626058 keV    e-  351.8774 keV    e+  342.54485 keV proton 0 eV
1881
1882Index : 4     used in the geometry : No      recalculation needed : No
1883 Material : Aluminium
1884 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1885 Energy thresholds :  gamma  6.9036292 keV    e-  598.34507 keV    e+  570.85042 keV proton 0 eV
1886
1887Index : 5     used in the geometry : Yes     recalculation needed : No
1888 Material : Iron
1889 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1890 Energy thresholds :  gamma  20.643844 keV    e-  1.295918 MeV    e+  1.2116932 MeV proton 0 eV
1891 Region(s) which use this couple :
1892    DefaultRegionForTheWorld
1893
1894====================================================================
1895
1896Start closing geometry.
1897G4GeometryManager::ReportVoxelStats -- Voxel Statistics
1898
1899    Total memory consumed for geometry optimisation:   0 kByte
1900    Total CPU time elapsed for geometry optimisation: 0 seconds
1901### Run 4 start.
1902
1903--------- Ranecu engine status ---------
1904 Initial seed (index) = 0
1905 Current couple of seeds = 1162542043, 318101973
1906----------------------------------------
1907Start Run processing.
1908
1909---> Begin of Event: 0
1910
1911---> Begin of Event: 10000
1912
1913---> Begin of Event: 20000
1914
1915---> Begin of Event: 30000
1916
1917---> Begin of Event: 40000
1918Run terminated.
1919Run Summary
1920  Number of events processed : 50000
1921  User=26.38s Real=44.83s Sys=1.28s
1922
1923 ======================== run summary ======================
1924
1925 The run was 50000 mu+ of 100 GeV through 1 m   of Iron (density: 7.87 g/cm3 )
1926
1927 Total energy deposit in absorber per event = 2.126 GeV +- 14.13 MeV
1928
1929 -----> Mean dE/dx = 21.26 MeV/cm       (2.702 MeV*cm2/g)
1930
1931 From formulas :
1932   restricted dEdx = 11.22 MeV/cm       (1.425 MeV*cm2/g)
1933   full dEdx       = 21.6 MeV/cm        (2.745 MeV*cm2/g)
1934
1935 Leakage :  primary = 97.87 GeV +- 14.15 MeV   secondaries = 0 eV  +- 0 eV
1936 Energy balance :  edep + eleak = 100 GeV
1937
1938 Total track length (charged) in absorber per event = 99 cm  +- 52.13 um
1939 Total track length (neutral) in absorber per event = 0 fm  +- 0 fm
1940
1941 Number of steps (charged) in absorber per event = 47.44 +- 0.001021
1942 Number of steps (neutral) in absorber per event = 0 +- 0
1943
1944 Number of secondaries per event : Gammas = 0.04906;   electrons = 45.84;   positrons = 1.251
1945
1946 Number of events with the primary particle transmitted = 100 %
1947 Number of events with at least  1 particle transmitted (same charge as primary) = 100 %
1948
1949 Number of events with the primary particle reflected = 0 %
1950 Number of events with at least  1 particle reflected (same charge as primary) = 0 %
1951
1952 MultipleScattering:
1953  rms proj angle of transmit primary particle = 1.121 mrad (central part only)
1954  computed theta0 (Highland formula)          = 1.182 mrad
1955  central part defined as +- 3.545 mrad;   Tail ratio = 0.289 %
1956
1957--------- Ranecu engine status ---------
1958 Initial seed (index) = 0
1959 Current couple of seeds = 1305857416, 170091693
1960----------------------------------------
1961#
1962UserDetectorConstruction deleted.
1963UserPhysicsList deleted.
1964UserRunAction deleted.
1965UserPrimaryGenerator deleted.
1966G4 kernel has come to Quit state.
1967EventManager deleted.
1968UImanager deleted.
1969Units table cleared.
1970StateManager deleted.
1971RunManagerKernel is deleted.
1972RunManager is deleting.
1973
Note: See TracBrowser for help on using the repository browser.