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

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

update ti head
File size: 95.8 KB
Line 
1
2        ############################################
3        !!! WARNING - FPE detection is activated !!!
4        ############################################
5
6*************************************************************
7 Geant4 version Name: geant4-09-03-ref-09    (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: 114561548.020 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.761 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        BetheHeitler :   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  AngularGenUrban
92            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
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  AngularGenUrban
107            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
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
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
169msc:   for kaon+    SubType= 10
170      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
171      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
172      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
173          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
174
175hIoni:   for  kaon+    SubType= 2
176      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
177      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
178      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
179      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
180               Bragg :   Emin=        0 eV       Emax=    1.05231 MeV
181          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
182
183hBrems:   for  kaon+    SubType= 3
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               hBrem :   Emin=        0 eV       Emax=         10 TeV
188
189hPairProd:   for  kaon+    SubType= 4
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
193           hPairProd :   Emin=        0 eV       Emax=         10 TeV
194
195hIoni:   for  kaon-    SubType= 2
196      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
197      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
198      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
199      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
200            ICRU73QO :   Emin=        0 eV       Emax=    1.05231 MeV
201          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
202
203hBrems:   for  kaon-    SubType= 3
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               hBrem :   Emin=        0 eV       Emax=         10 TeV
208
209hPairProd:   for  kaon-    SubType= 4
210      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
211      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
212      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
213           hPairProd :   Emin=        0 eV       Emax=         10 TeV
214
215muMsc:   for mu+    SubType= 10
216      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
217      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592
218      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
219        WentzelVIUni :   Emin=        0 eV       Emax=         10 TeV
220
221muIoni:   for  mu+    SubType= 2
222      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
223      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
224      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
225      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
226               Bragg :   Emin=        0 eV       Emax=        200 keV
227          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
228        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
229
230muBrems:   for  mu+    SubType= 3
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              MuBrem :   Emin=        0 eV       Emax=         10 TeV
235
236muPairProd:   for  mu+    SubType= 4
237      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
238      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
239      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
240          muPairProd :   Emin=        0 eV       Emax=         10 TeV
241
242CoulombScat:   for  mu+    SubType= 1
243      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
244      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
245      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
246  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
247
248muIoni:   for  mu-    SubType= 2
249      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
250      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
251      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
252      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
253            ICRU73QO :   Emin=        0 eV       Emax=        200 keV
254          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
255        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
256
257muBrems:   for  mu-    SubType= 3
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              MuBrem :   Emin=        0 eV       Emax=         10 TeV
262
263muPairProd:   for  mu-    SubType= 4
264      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
265      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
266      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
267          muPairProd :   Emin=        0 eV       Emax=         10 TeV
268
269CoulombScat:   for  mu-    SubType= 1
270      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
271      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
272      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
273  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
274
275hIoni:   for  pi+    SubType= 2
276      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
277      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
278      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
279      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
280               Bragg :   Emin=        0 eV       Emax=    297.505 keV
281          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
282
283hBrems:   for  pi+    SubType= 3
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               hBrem :   Emin=        0 eV       Emax=         10 TeV
288
289hPairProd:   for  pi+    SubType= 4
290      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
291      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
292      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
293           hPairProd :   Emin=        0 eV       Emax=         10 TeV
294
295msc:   for pi-    SubType= 10
296      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
297      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
298      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
299          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
300
301hIoni:   for  pi-    SubType= 2
302      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
303      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
304      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
305      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
306            ICRU73QO :   Emin=        0 eV       Emax=    297.505 keV
307          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
308
309hBrems:   for  pi-    SubType= 3
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               hBrem :   Emin=        0 eV       Emax=         10 TeV
314
315hPairProd:   for  pi-    SubType= 4
316      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
317      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
318      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
319           hPairProd :   Emin=        0 eV       Emax=         10 TeV
320
321Region <DefaultRegionForTheWorld> -- appears in <World> world volume
322 Root logical volume(s) : World
323 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
324 Materials : Galactic Silicon
325 Production cuts :   gamma 7 um      e- 7 um      e+ 7 um  proton 0 fm
326
327========= Table of registered couples ==============================
328
329Index : 0     used in the geometry : Yes     recalculation needed : No
330 Material : Galactic
331 Range cuts        :  gamma  7 um     e-  7 um     e+  7 um  proton 0 fm
332 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
333 Region(s) which use this couple :
334    DefaultRegionForTheWorld
335
336Index : 1     used in the geometry : Yes     recalculation needed : No
337 Material : Silicon
338 Range cuts        :  gamma  7 um     e-  7 um     e+  7 um  proton 0 fm
339 Energy thresholds :  gamma  990 eV     e-  21.4886 keV    e+  20.8954 keV proton 0 eV
340 Region(s) which use this couple :
341    DefaultRegionForTheWorld
342
343====================================================================
344
345Start closing geometry.
346G4GeometryManager::ReportVoxelStats -- Voxel Statistics
347
348    Total memory consumed for geometry optimisation:   0 kByte
349    Total CPU time elapsed for geometry optimisation: 0 seconds
350### Run 0 start.
351
352--------- Ranecu engine status ---------
353 Initial seed (index) = 0
354 Current couple of seeds = 9876, 54321
355----------------------------------------
356Start Run processing.
357
358---> Begin of Event: 0
359
360---> Begin of Event: 10000
361
362---> Begin of Event: 20000
363
364---> Begin of Event: 30000
365
366---> Begin of Event: 40000
367
368---> Begin of Event: 50000
369
370---> Begin of Event: 60000
371
372---> Begin of Event: 70000
373
374---> Begin of Event: 80000
375
376---> Begin of Event: 90000
377Run terminated.
378Run Summary
379  Number of events processed : 100000
380  User=4.2s Real=6.41s Sys=2.21s
381
382 ======================== run summary ======================
383
384 The run was 100000 pi+ of 5 GeV through 20 um  of Silicon (density: 2.33 g/cm3 )
385
386 Total energy deposit in absorber per event = 6.355 keV +- 16.41 eV
387
388 -----> Mean dE/dx = 3.177 MeV/cm       (1.364 MeV*cm2/g)
389
390 From formulas :
391   restricted dEdx = 2.963 MeV/cm       (1.272 MeV*cm2/g)
392   full dEdx       = 4.726 MeV/cm       (2.028 MeV*cm2/g)
393
394 Leakage :  primary = 5 GeV +- 1.274 keV   secondaries = 3.64 keV +- 1.273 keV
395 Energy balance :  edep + eleak = 5 GeV
396
397 Total track length (charged) in absorber per event = 20.26 um  +- 7.821 nm
398 Total track length (neutral) in absorber per event = 5.528 Ang +- 3.51 Ang
399
400 Number of steps (charged) in absorber per event = 1.159 +- 8.843e-06
401 Number of steps (neutral) in absorber per event = 3e-05 +- 8.843e-06
402
403 Number of secondaries per event : Gammas = 3e-05;   electrons = 0.01661;   positrons = 0
404
405 Number of events with the primary particle transmitted = 100 %
406 Number of events with at least  1 particle transmitted (same charge as primary) = 100 %
407
408 Number of events with the primary particle reflected = 0 %
409 Number of events with at least  1 particle reflected (same charge as primary) = 0 %
410
411 MultipleScattering:
412  rms proj angle of transmit primary particle = 0.02446 mrad (central part only)
413  computed theta0 (Highland formula)          = 0.02627 mrad
414  central part defined as +- 0.0788 mrad;   Tail ratio = 1.553 %
415
416--------- Ranecu engine status ---------
417 Initial seed (index) = 0
418 Current couple of seeds = 1096967505, 715451002
419----------------------------------------
420#
421/testem/det/setAbsMat Gold
422/testem/det/setAbsThick 9.658 um
423/testem/det/setAbsYZ 1 mm
424/testem/det/update
425
426 Material: Galactic     density:  0.000 kg/m3   RadL: 204727512.315 pc   Nucl.Int.Length: 114561548.020 pc   Imean:  19.200 eV   temperature:   2.73 K  pressure:   0.00 atm
427   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
428
429 Material:     Gold     density: 19.320 g/cm3   RadL:   3.344 mm   Nucl.Int.Length:  10.539 cm   Imean: 790.000 eV
430   --->  Element: Gold ( )   Z = 79.0   N = 197.0   A = 196.97 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
431
432 The  WORLD   is made of 14.487 um  of Galactic. The transverse size (YZ) of the world is 1.2 mm
433 The ABSORBER is made of 9.658 um  of Gold. The transverse size (YZ) is 1 mm
434 X position of the middle of the absorber 0 fm
435Obsolete world logical volume is removed from the default region.
436World is registered to the default region.
437/testem/gun/setDefault
438/gun/particle e-
439/gun/energy 15.7 MeV
440/run/beamOn 50000
441
442phot:   for  gamma    SubType= 12
443      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
444       PhotoElectric :   Emin=        0 eV       Emax=         10 TeV
445
446compt:   for  gamma    SubType= 13
447      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
448      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
449       Klein-Nishina :   Emin=        0 eV       Emax=         10 TeV
450
451conv:   for  gamma    SubType= 14
452      Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1
453      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
454        BetheHeitler :   Emin=        0 eV       Emax=         10 TeV
455
456Rayl:   for  gamma    SubType= 11
457      Lambda tables from 100 eV  to 10 TeV in 200 bins, spline: 1
458      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
459   LivermoreRayleigh :   Emin=        0 eV       Emax=        100 GeV
460
461msc:   for e-    SubType= 10
462      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
463      RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
464      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
465          UrbanMsc93 :   Emin=        0 eV       Emax=         10 TeV
466
467eIoni:   for  e-    SubType= 2
468      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
469      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
470      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
471      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
472        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
473
474eBrem:   for  e-    SubType= 3
475      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
476      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
477      LPM flag: 1 for E > 1 GeV
478      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
479               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
480            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
481
482eIoni:   for  e+    SubType= 2
483      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
484      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
485      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
486      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
487        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
488
489eBrem:   for  e+    SubType= 3
490      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
491      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
492      LPM flag: 1 for E > 1 GeV
493      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
494               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
495            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
496
497annihil:   for  e+    SubType= 5
498      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
499      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
500            eplus2gg :   Emin=        0 eV       Emax=         10 TeV
501
502msc:   for proton    SubType= 10
503      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
504      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
505      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
506          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
507
508hIoni:   for  proton    SubType= 2
509      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
510      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
511      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
512      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
513               Bragg :   Emin=        0 eV       Emax=          2 MeV
514          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
515
516hBrems:   for  proton    SubType= 3
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               hBrem :   Emin=        0 eV       Emax=         10 TeV
521
522hPairProd:   for  proton    SubType= 4
523      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
524      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
525      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
526           hPairProd :   Emin=        0 eV       Emax=         10 TeV
527
528msc:   for GenericIon    SubType= 10
529      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
530      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
531          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
532
533ionIoni:   for  GenericIon    SubType= 2
534      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
535      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
536      finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
537      Stopping Power data for 17 ion/material pairs
538      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
539         ParamICRU73 :   Emin=        0 eV       Emax=         10 TeV
540
541nuclearStopping:   for  GenericIon    SubType= 1
542      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
543   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
544
545nuclearStopping:   for  alpha    SubType= 1
546      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
547   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
548
549hIoni:   for  anti_proton    SubType= 2
550      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
551      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
552      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
553      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
554            ICRU73QO :   Emin=        0 eV       Emax=          2 MeV
555          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
556
557msc:   for kaon+    SubType= 10
558      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
559      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
560      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
561          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
562
563hIoni:   for  kaon+    SubType= 2
564      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
565      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
566      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
567      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
568               Bragg :   Emin=        0 eV       Emax=    1.05231 MeV
569          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
570
571hBrems:   for  kaon+    SubType= 3
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
575               hBrem :   Emin=        0 eV       Emax=         10 TeV
576
577hPairProd:   for  kaon+    SubType= 4
578      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
579      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
580      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
581           hPairProd :   Emin=        0 eV       Emax=         10 TeV
582
583hIoni:   for  kaon-    SubType= 2
584      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
585      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
586      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
587      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
588            ICRU73QO :   Emin=        0 eV       Emax=    1.05231 MeV
589          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
590
591hBrems:   for  kaon-    SubType= 3
592      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
593      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
594      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
595               hBrem :   Emin=        0 eV       Emax=         10 TeV
596
597hPairProd:   for  kaon-    SubType= 4
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
601           hPairProd :   Emin=        0 eV       Emax=         10 TeV
602
603muMsc:   for mu+    SubType= 10
604      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
605      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592
606      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
607        WentzelVIUni :   Emin=        0 eV       Emax=         10 TeV
608
609muIoni:   for  mu+    SubType= 2
610      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
611      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
612      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
613      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
614               Bragg :   Emin=        0 eV       Emax=        200 keV
615          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
616        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
617
618muBrems:   for  mu+    SubType= 3
619      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
620      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
621      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
622              MuBrem :   Emin=        0 eV       Emax=         10 TeV
623
624muPairProd:   for  mu+    SubType= 4
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
628          muPairProd :   Emin=        0 eV       Emax=         10 TeV
629
630CoulombScat:   for  mu+    SubType= 1
631      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
632      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
633      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
634  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
635
636muIoni:   for  mu-    SubType= 2
637      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
638      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
639      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
640      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
641            ICRU73QO :   Emin=        0 eV       Emax=        200 keV
642          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
643        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
644
645muBrems:   for  mu-    SubType= 3
646      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
647      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
648      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
649              MuBrem :   Emin=        0 eV       Emax=         10 TeV
650
651muPairProd:   for  mu-    SubType= 4
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
655          muPairProd :   Emin=        0 eV       Emax=         10 TeV
656
657CoulombScat:   for  mu-    SubType= 1
658      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
659      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
660      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
661  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
662
663hIoni:   for  pi+    SubType= 2
664      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
665      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
666      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
667      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
668               Bragg :   Emin=        0 eV       Emax=    297.505 keV
669          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
670
671hBrems:   for  pi+    SubType= 3
672      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
673      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
674      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
675               hBrem :   Emin=        0 eV       Emax=         10 TeV
676
677hPairProd:   for  pi+    SubType= 4
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
681           hPairProd :   Emin=        0 eV       Emax=         10 TeV
682
683msc:   for pi-    SubType= 10
684      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
685      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
686      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
687          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
688
689hIoni:   for  pi-    SubType= 2
690      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
691      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
692      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
693      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
694            ICRU73QO :   Emin=        0 eV       Emax=    297.505 keV
695          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
696
697hBrems:   for  pi-    SubType= 3
698      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
699      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
700      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
701               hBrem :   Emin=        0 eV       Emax=         10 TeV
702
703hPairProd:   for  pi-    SubType= 4
704      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
705      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
706      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
707           hPairProd :   Emin=        0 eV       Emax=         10 TeV
708
709Region <DefaultRegionForTheWorld> -- appears in <World> world volume
710 Root logical volume(s) : World
711 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
712 Materials : Galactic Gold
713 Production cuts :   gamma 7 um      e- 7 um      e+ 7 um  proton 0 fm
714
715========= Table of registered couples ==============================
716
717Index : 0     used in the geometry : Yes     recalculation needed : No
718 Material : Galactic
719 Range cuts        :  gamma  7 um     e-  7 um     e+  7 um  proton 0 fm
720 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
721 Region(s) which use this couple :
722    DefaultRegionForTheWorld
723
724Index : 1     used in the geometry : No      recalculation needed : No
725 Material : Silicon
726 Range cuts        :  gamma  7 um     e-  7 um     e+  7 um  proton 0 fm
727 Energy thresholds :  gamma  990 eV     e-  21.4886 keV    e+  20.8954 keV proton 0 eV
728
729Index : 2     used in the geometry : Yes     recalculation needed : No
730 Material : Gold
731 Range cuts        :  gamma  7 um     e-  7 um     e+  7 um  proton 0 fm
732 Energy thresholds :  gamma  6.54229 keV    e-  65.1406 keV    e+  63.8405 keV proton 0 eV
733 Region(s) which use this couple :
734    DefaultRegionForTheWorld
735
736====================================================================
737
738Start closing geometry.
739G4GeometryManager::ReportVoxelStats -- Voxel Statistics
740
741    Total memory consumed for geometry optimisation:   0 kByte
742    Total CPU time elapsed for geometry optimisation: 0 seconds
743### Run 1 start.
744
745--------- Ranecu engine status ---------
746 Initial seed (index) = 0
747 Current couple of seeds = 1096967505, 715451002
748----------------------------------------
749Start Run processing.
750
751---> Begin of Event: 0
752
753---> Begin of Event: 10000
754
755---> Begin of Event: 20000
756
757---> Begin of Event: 30000
758
759---> Begin of Event: 40000
760Run terminated.
761Run Summary
762  Number of events processed : 50000
763  User=3.35s Real=4.56s Sys=1.21s
764
765 ======================== run summary ======================
766
767 The run was 50000 e- of 15.7 MeV through 9.66 um  of Gold (density: 19.3 g/cm3 )
768
769 Total energy deposit in absorber per event = 18.54 keV +- 65.36 eV
770
771 -----> Mean dE/dx = 19.2 MeV/cm        (0.9938 MeV*cm2/g)
772
773 From formulas :
774   restricted dEdx = 17.9 MeV/cm        (0.9266 MeV*cm2/g)
775   full dEdx       = 59.7 MeV/cm        (3.09 MeV*cm2/g)
776
777 Leakage :  primary = 15.63 MeV +- 2.365 keV   secondaries = 47.06 keV +- 2.361 keV
778 Energy balance :  edep + eleak = 15.7 MeV
779
780 Total track length (charged) in absorber per event = 9.938 um  +- 12.01 nm
781 Total track length (neutral) in absorber per event = 131.3 nm  +- 4.242 nm
782
783 Number of steps (charged) in absorber per event = 5.069 +- 1.55e-05
784 Number of steps (neutral) in absorber per event = 0.02996 +- 1.55e-05
785
786 Number of secondaries per event : Gammas = 0.02986;   electrons = 0.023;   positrons = 0
787
788 Number of events with the primary particle transmitted = 99.98 %
789 Number of events with at least  1 particle transmitted (same charge as primary) = 99.98 %
790
791 Number of events with the primary particle reflected = 0.016 %
792 Number of events with at least  1 particle reflected (same charge as primary) = 0.294 %
793
794 MultipleScattering:
795  rms proj angle of transmit primary particle = 32.82 mrad (central part only)
796  computed theta0 (Highland formula)          = 35.1 mrad
797  central part defined as +- 105.3 mrad;   Tail ratio = 2.774 %
798
799--------- Ranecu engine status ---------
800 Initial seed (index) = 0
801 Current couple of seeds = 1228101883, 1478657312
802----------------------------------------
803#
804/testem/phys/setCuts 10 um
805#
806/testem/det/setAbsMat G4_POLYETHYLENE
807/testem/det/setAbsThick 4.74 mm
808/testem/det/setAbsYZ 50 mm
809/testem/det/update
810
811 Material: Galactic     density:  0.000 kg/m3   RadL: 204727512.315 pc   Nucl.Int.Length: 114561548.020 pc   Imean:  19.200 eV   temperature:   2.73 K  pressure:   0.00 atm
812   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
813
814 Material: G4_POLYETHYLENE (C_2H_4)_N-Polyethylene   density: 940.000 mg/cm3  RadL:  47.632 cm   Nucl.Int.Length:  72.122 cm   Imean:  57.400 eV
815   --->  Element: H (H)   Z =  1.0   N =   1.0   A =   1.01 g/mole
816   --->  Isotope:    H1   Z =  1   N =   1   A =   1.01 g/mole   abundance:  99.99 %
817   --->  Isotope:    H2   Z =  1   N =   2   A =   2.01 g/mole   abundance:   0.01 %  ElmMassFraction:  14.37 %  ElmAbundance  66.67 %
818   --->  Element: C (C)   Z =  6.0   N =  12.0   A =  12.01 g/mole
819   --->  Isotope:   C12   Z =  6   N =  12   A =  12.00 g/mole   abundance:  98.93 %
820   --->  Isotope:   C13   Z =  6   N =  13   A =  13.00 g/mole   abundance:   1.07 %  ElmMassFraction:  85.63 %  ElmAbundance  33.33 %
821
822 The  WORLD   is made of 7.11 mm  of Galactic. The transverse size (YZ) of the world is 6 cm
823 The ABSORBER is made of 4.74 mm  of G4_POLYETHYLENE. The transverse size (YZ) is 5 cm
824 X position of the middle of the absorber 0 fm
825Obsolete world logical volume is removed from the default region.
826World is registered to the default region.
827/testem/gun/setDefault
828/gun/particle mu+
829/gun/energy 96.2 MeV
830/run/beamOn 50000
831
832phot:   for  gamma    SubType= 12
833      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
834       PhotoElectric :   Emin=        0 eV       Emax=         10 TeV
835
836compt:   for  gamma    SubType= 13
837      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
838      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
839       Klein-Nishina :   Emin=        0 eV       Emax=         10 TeV
840
841conv:   for  gamma    SubType= 14
842      Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1
843      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
844        BetheHeitler :   Emin=        0 eV       Emax=         10 TeV
845
846Rayl:   for  gamma    SubType= 11
847      Lambda tables from 100 eV  to 10 TeV in 200 bins, spline: 1
848      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
849   LivermoreRayleigh :   Emin=        0 eV       Emax=        100 GeV
850
851msc:   for e-    SubType= 10
852      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
853      RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
854      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
855          UrbanMsc93 :   Emin=        0 eV       Emax=         10 TeV
856
857eIoni:   for  e-    SubType= 2
858      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
859      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
860      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
861      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
862        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
863
864eBrem:   for  e-    SubType= 3
865      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
866      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
867      LPM flag: 1 for E > 1 GeV
868      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
869               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
870            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
871
872eIoni:   for  e+    SubType= 2
873      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
874      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
875      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
876      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
877        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
878
879eBrem:   for  e+    SubType= 3
880      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
881      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
882      LPM flag: 1 for E > 1 GeV
883      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
884               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
885            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
886
887annihil:   for  e+    SubType= 5
888      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
889      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
890            eplus2gg :   Emin=        0 eV       Emax=         10 TeV
891
892msc:   for proton    SubType= 10
893      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
894      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
895      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
896          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
897
898hIoni:   for  proton    SubType= 2
899      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
900      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
901      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
902      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
903               Bragg :   Emin=        0 eV       Emax=          2 MeV
904          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
905
906hBrems:   for  proton    SubType= 3
907      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
908      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
909      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
910               hBrem :   Emin=        0 eV       Emax=         10 TeV
911
912hPairProd:   for  proton    SubType= 4
913      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
914      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
915      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
916           hPairProd :   Emin=        0 eV       Emax=         10 TeV
917
918msc:   for GenericIon    SubType= 10
919      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
920      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
921          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
922
923ionIoni:   for  GenericIon    SubType= 2
924      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
925      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
926      finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
927      Stopping Power data for 17 ion/material pairs
928      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
929         ParamICRU73 :   Emin=        0 eV       Emax=         10 TeV
930
931nuclearStopping:   for  GenericIon    SubType= 1
932      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
933   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
934
935nuclearStopping:   for  alpha    SubType= 1
936      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
937   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
938
939hIoni:   for  anti_proton    SubType= 2
940      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
941      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
942      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
943      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
944            ICRU73QO :   Emin=        0 eV       Emax=          2 MeV
945          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
946
947msc:   for kaon+    SubType= 10
948      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
949      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
950      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
951          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
952
953hIoni:   for  kaon+    SubType= 2
954      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
955      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
956      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
957      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
958               Bragg :   Emin=        0 eV       Emax=    1.05231 MeV
959          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
960
961hBrems:   for  kaon+    SubType= 3
962      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
963      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
964      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
965               hBrem :   Emin=        0 eV       Emax=         10 TeV
966
967hPairProd:   for  kaon+    SubType= 4
968      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
969      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
970      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
971           hPairProd :   Emin=        0 eV       Emax=         10 TeV
972
973hIoni:   for  kaon-    SubType= 2
974      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
975      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
976      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
977      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
978            ICRU73QO :   Emin=        0 eV       Emax=    1.05231 MeV
979          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
980
981hBrems:   for  kaon-    SubType= 3
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
985               hBrem :   Emin=        0 eV       Emax=         10 TeV
986
987hPairProd:   for  kaon-    SubType= 4
988      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
989      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
990      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
991           hPairProd :   Emin=        0 eV       Emax=         10 TeV
992
993muMsc:   for mu+    SubType= 10
994      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
995      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592
996      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
997        WentzelVIUni :   Emin=        0 eV       Emax=         10 TeV
998
999muIoni:   for  mu+    SubType= 2
1000      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1001      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1002      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1003      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1004               Bragg :   Emin=        0 eV       Emax=        200 keV
1005          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
1006        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
1007
1008muBrems:   for  mu+    SubType= 3
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1012              MuBrem :   Emin=        0 eV       Emax=         10 TeV
1013
1014muPairProd:   for  mu+    SubType= 4
1015      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1016      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1017      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1018          muPairProd :   Emin=        0 eV       Emax=         10 TeV
1019
1020CoulombScat:   for  mu+    SubType= 1
1021      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1022      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
1023      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1024  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
1025
1026muIoni:   for  mu-    SubType= 2
1027      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1028      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1029      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1030      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1031            ICRU73QO :   Emin=        0 eV       Emax=        200 keV
1032          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
1033        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
1034
1035muBrems:   for  mu-    SubType= 3
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1039              MuBrem :   Emin=        0 eV       Emax=         10 TeV
1040
1041muPairProd:   for  mu-    SubType= 4
1042      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1043      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1044      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1045          muPairProd :   Emin=        0 eV       Emax=         10 TeV
1046
1047CoulombScat:   for  mu-    SubType= 1
1048      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1049      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
1050      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1051  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
1052
1053hIoni:   for  pi+    SubType= 2
1054      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1055      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1056      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1057      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1058               Bragg :   Emin=        0 eV       Emax=    297.505 keV
1059          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
1060
1061hBrems:   for  pi+    SubType= 3
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1065               hBrem :   Emin=        0 eV       Emax=         10 TeV
1066
1067hPairProd:   for  pi+    SubType= 4
1068      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1069      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1070      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1071           hPairProd :   Emin=        0 eV       Emax=         10 TeV
1072
1073msc:   for pi-    SubType= 10
1074      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1075      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1076      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1077          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
1078
1079hIoni:   for  pi-    SubType= 2
1080      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1081      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1082      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1083      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1084            ICRU73QO :   Emin=        0 eV       Emax=    297.505 keV
1085          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
1086
1087hBrems:   for  pi-    SubType= 3
1088      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1089      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1090      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1091               hBrem :   Emin=        0 eV       Emax=         10 TeV
1092
1093hPairProd:   for  pi-    SubType= 4
1094      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1095      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1096      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1097           hPairProd :   Emin=        0 eV       Emax=         10 TeV
1098
1099Region <DefaultRegionForTheWorld> -- appears in <World> world volume
1100 Root logical volume(s) : World
1101 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
1102 Materials : Galactic G4_POLYETHYLENE
1103 Production cuts :   gamma 10 um      e- 10 um      e+ 10 um  proton 0 fm
1104
1105========= Table of registered couples ==============================
1106
1107Index : 0     used in the geometry : Yes     recalculation needed : No
1108 Material : Galactic
1109 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
1110 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
1111 Region(s) which use this couple :
1112    DefaultRegionForTheWorld
1113
1114Index : 1     used in the geometry : No      recalculation needed : No
1115 Material : Silicon
1116 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
1117 Energy thresholds :  gamma  990 eV     e-  31.9516 keV    e+  31.7376 keV proton 0 eV
1118
1119Index : 2     used in the geometry : No      recalculation needed : No
1120 Material : Gold
1121 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
1122 Energy thresholds :  gamma  8.27707 keV    e-  80.7687 keV    e+  79.1567 keV proton 0 eV
1123
1124Index : 3     used in the geometry : Yes     recalculation needed : No
1125 Material : G4_POLYETHYLENE
1126 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
1127 Energy thresholds :  gamma  990 eV     e-  14.1748 keV    e+  13.7548 keV proton 0 eV
1128 Region(s) which use this couple :
1129    DefaultRegionForTheWorld
1130
1131====================================================================
1132
1133Start closing geometry.
1134G4GeometryManager::ReportVoxelStats -- Voxel Statistics
1135
1136    Total memory consumed for geometry optimisation:   0 kByte
1137    Total CPU time elapsed for geometry optimisation: 0 seconds
1138### Run 2 start.
1139
1140--------- Ranecu engine status ---------
1141 Initial seed (index) = 0
1142 Current couple of seeds = 1228101883, 1478657312
1143----------------------------------------
1144Start Run processing.
1145
1146---> Begin of Event: 0
1147
1148---> Begin of Event: 10000
1149
1150---> Begin of Event: 20000
1151
1152---> Begin of Event: 30000
1153
1154---> Begin of Event: 40000
1155Run terminated.
1156Run Summary
1157  Number of events processed : 50000
1158  User=6.36s Real=7.53s Sys=1.16s
1159
1160 ======================== run summary ======================
1161
1162 The run was 50000 mu+ of 96.2 MeV through 4.74 mm  of G4_POLYETHYLENE (density: 940 mg/cm3)
1163
1164 Total energy deposit in absorber per event = 1.031 MeV +- 1.011 keV
1165
1166 -----> Mean dE/dx = 2.175 MeV/cm       (2.314 MeV*cm2/g)
1167
1168 From formulas :
1169   restricted dEdx = 1.747 MeV/cm       (1.859 MeV*cm2/g)
1170   full dEdx       = 2.259 MeV/cm       (2.403 MeV*cm2/g)
1171
1172 Leakage :  primary = 95.14 MeV +- 1.348 keV   secondaries = 26.07 keV +- 707.1 eV
1173 Energy balance :  edep + eleak = 96.2 MeV
1174
1175 Total track length (charged) in absorber per event = 5.102 mm  +- 3.895 um
1176 Total track length (neutral) in absorber per event = 26.38 um  +- 2.481 um
1177
1178 Number of steps (charged) in absorber per event = 11.27 +- 0.0002791
1179 Number of steps (neutral) in absorber per event = 0.01202 +- 0.0002791
1180
1181 Number of secondaries per event : Gammas = 0.01128;   electrons = 3.989;   positrons = 0
1182
1183 Number of events with the primary particle transmitted = 100 %
1184 Number of events with at least  1 particle transmitted (same charge as primary) = 100 %
1185
1186 Number of events with the primary particle reflected = 0 %
1187 Number of events with at least  1 particle reflected (same charge as primary) = 0 %
1188
1189 MultipleScattering:
1190  rms proj angle of transmit primary particle = 7.275 mrad (central part only)
1191  computed theta0 (Highland formula)          = 7.636 mrad
1192  central part defined as +- 22.91 mrad;   Tail ratio = 1.259 %
1193
1194--------- Ranecu engine status ---------
1195 Initial seed (index) = 0
1196 Current couple of seeds = 1091561076, 776987475
1197----------------------------------------
1198#
1199/testem/phys/setCuts 1 mm
1200#
1201/testem/stack/killSecondaries
1202#
1203/testem/det/setAbsMat Aluminium
1204/testem/det/setAbsThick 8.004 mm
1205/testem/det/update
1206
1207 Material: Galactic     density:  0.000 kg/m3   RadL: 204727512.315 pc   Nucl.Int.Length: 114561548.020 pc   Imean:  19.200 eV   temperature:   2.73 K  pressure:   0.00 atm
1208   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
1209
1210 Material: Aluminium     density:  2.700 g/cm3   RadL:   8.893 cm   Nucl.Int.Length:  38.860 cm   Imean: 166.000 eV
1211   --->  Element: Aluminium ( )   Z = 13.0   N =  27.0   A =  26.98 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
1212
1213 The  WORLD   is made of 1.2006 cm  of Galactic. The transverse size (YZ) of the world is 6 cm
1214 The ABSORBER is made of 8.004 mm  of Aluminium. The transverse size (YZ) is 5 cm
1215 X position of the middle of the absorber 0 fm
1216Obsolete world logical volume is removed from the default region.
1217World is registered to the default region.
1218/testem/gun/setDefault
1219/gun/particle proton
1220/gun/energy 174.664 GeV
1221/run/beamOn 50000
1222
1223phot:   for  gamma    SubType= 12
1224      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1225       PhotoElectric :   Emin=        0 eV       Emax=         10 TeV
1226
1227compt:   for  gamma    SubType= 13
1228      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1229      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1230       Klein-Nishina :   Emin=        0 eV       Emax=         10 TeV
1231
1232conv:   for  gamma    SubType= 14
1233      Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1
1234      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1235        BetheHeitler :   Emin=        0 eV       Emax=         10 TeV
1236
1237Rayl:   for  gamma    SubType= 11
1238      Lambda tables from 100 eV  to 10 TeV in 200 bins, spline: 1
1239      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1240   LivermoreRayleigh :   Emin=        0 eV       Emax=        100 GeV
1241
1242msc:   for e-    SubType= 10
1243      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1244      RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
1245      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1246          UrbanMsc93 :   Emin=        0 eV       Emax=         10 TeV
1247
1248eIoni:   for  e-    SubType= 2
1249      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1250      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1251      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1252      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1253        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
1254
1255eBrem:   for  e-    SubType= 3
1256      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1257      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1258      LPM flag: 1 for E > 1 GeV
1259      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1260               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
1261            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
1262
1263eIoni:   for  e+    SubType= 2
1264      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1265      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1266      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1267      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1268        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
1269
1270eBrem:   for  e+    SubType= 3
1271      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1272      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1273      LPM flag: 1 for E > 1 GeV
1274      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1275               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
1276            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
1277
1278annihil:   for  e+    SubType= 5
1279      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1280      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1281            eplus2gg :   Emin=        0 eV       Emax=         10 TeV
1282
1283msc:   for proton    SubType= 10
1284      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1285      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1286      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1287          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
1288
1289hIoni:   for  proton    SubType= 2
1290      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1291      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1292      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1293      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1294               Bragg :   Emin=        0 eV       Emax=          2 MeV
1295          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
1296
1297hBrems:   for  proton    SubType= 3
1298      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1299      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1300      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1301               hBrem :   Emin=        0 eV       Emax=         10 TeV
1302
1303hPairProd:   for  proton    SubType= 4
1304      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1305      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1306      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1307           hPairProd :   Emin=        0 eV       Emax=         10 TeV
1308
1309msc:   for GenericIon    SubType= 10
1310      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
1311      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1312          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
1313
1314ionIoni:   for  GenericIon    SubType= 2
1315      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1316      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1317      finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
1318      Stopping Power data for 17 ion/material pairs
1319      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1320         ParamICRU73 :   Emin=        0 eV       Emax=         10 TeV
1321
1322nuclearStopping:   for  GenericIon    SubType= 1
1323      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1324   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
1325
1326nuclearStopping:   for  alpha    SubType= 1
1327      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1328   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
1329
1330hIoni:   for  anti_proton    SubType= 2
1331      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1332      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1333      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1334      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1335            ICRU73QO :   Emin=        0 eV       Emax=          2 MeV
1336          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
1337
1338msc:   for kaon+    SubType= 10
1339      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1340      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1341      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1342          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
1343
1344hIoni:   for  kaon+    SubType= 2
1345      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1346      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1347      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1348      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1349               Bragg :   Emin=        0 eV       Emax=    1.05231 MeV
1350          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
1351
1352hBrems:   for  kaon+    SubType= 3
1353      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1354      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1355      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1356               hBrem :   Emin=        0 eV       Emax=         10 TeV
1357
1358hPairProd:   for  kaon+    SubType= 4
1359      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1360      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1361      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1362           hPairProd :   Emin=        0 eV       Emax=         10 TeV
1363
1364hIoni:   for  kaon-    SubType= 2
1365      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1366      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1367      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1368      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1369            ICRU73QO :   Emin=        0 eV       Emax=    1.05231 MeV
1370          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
1371
1372hBrems:   for  kaon-    SubType= 3
1373      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1374      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1375      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1376               hBrem :   Emin=        0 eV       Emax=         10 TeV
1377
1378hPairProd:   for  kaon-    SubType= 4
1379      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1380      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1381      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1382           hPairProd :   Emin=        0 eV       Emax=         10 TeV
1383
1384muMsc:   for mu+    SubType= 10
1385      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1386      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592
1387      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1388        WentzelVIUni :   Emin=        0 eV       Emax=         10 TeV
1389
1390muIoni:   for  mu+    SubType= 2
1391      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1392      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1393      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1394      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1395               Bragg :   Emin=        0 eV       Emax=        200 keV
1396          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
1397        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
1398
1399muBrems:   for  mu+    SubType= 3
1400      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1401      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1402      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1403              MuBrem :   Emin=        0 eV       Emax=         10 TeV
1404
1405muPairProd:   for  mu+    SubType= 4
1406      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1407      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1408      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1409          muPairProd :   Emin=        0 eV       Emax=         10 TeV
1410
1411CoulombScat:   for  mu+    SubType= 1
1412      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1413      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
1414      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1415  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
1416
1417muIoni:   for  mu-    SubType= 2
1418      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1419      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1420      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1421      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1422            ICRU73QO :   Emin=        0 eV       Emax=        200 keV
1423          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
1424        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
1425
1426muBrems:   for  mu-    SubType= 3
1427      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1428      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1429      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1430              MuBrem :   Emin=        0 eV       Emax=         10 TeV
1431
1432muPairProd:   for  mu-    SubType= 4
1433      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1434      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1435      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1436          muPairProd :   Emin=        0 eV       Emax=         10 TeV
1437
1438CoulombScat:   for  mu-    SubType= 1
1439      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1440      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
1441      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1442  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
1443
1444hIoni:   for  pi+    SubType= 2
1445      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1446      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1447      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1448      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1449               Bragg :   Emin=        0 eV       Emax=    297.505 keV
1450          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
1451
1452hBrems:   for  pi+    SubType= 3
1453      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1454      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1455      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1456               hBrem :   Emin=        0 eV       Emax=         10 TeV
1457
1458hPairProd:   for  pi+    SubType= 4
1459      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1460      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1461      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1462           hPairProd :   Emin=        0 eV       Emax=         10 TeV
1463
1464msc:   for pi-    SubType= 10
1465      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1466      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1467      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1468          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
1469
1470hIoni:   for  pi-    SubType= 2
1471      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1472      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1473      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1474      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1475            ICRU73QO :   Emin=        0 eV       Emax=    297.505 keV
1476          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
1477
1478hBrems:   for  pi-    SubType= 3
1479      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1480      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1481      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1482               hBrem :   Emin=        0 eV       Emax=         10 TeV
1483
1484hPairProd:   for  pi-    SubType= 4
1485      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1486      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1487      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1488           hPairProd :   Emin=        0 eV       Emax=         10 TeV
1489
1490Region <DefaultRegionForTheWorld> -- appears in <World> world volume
1491 Root logical volume(s) : World
1492 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
1493 Materials : Galactic Aluminium
1494 Production cuts :   gamma 1 mm      e- 1 mm      e+ 1 mm  proton 0 fm
1495
1496========= Table of registered couples ==============================
1497
1498Index : 0     used in the geometry : Yes     recalculation needed : No
1499 Material : Galactic
1500 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1501 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
1502 Region(s) which use this couple :
1503    DefaultRegionForTheWorld
1504
1505Index : 1     used in the geometry : No      recalculation needed : No
1506 Material : Silicon
1507 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1508 Energy thresholds :  gamma  6.95018 keV    e-  548.291 keV    e+  526.624 keV proton 0 eV
1509
1510Index : 2     used in the geometry : No      recalculation needed : No
1511 Material : Gold
1512 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1513 Energy thresholds :  gamma  112.644 keV    e-  2.30976 MeV    e+  2.13081 MeV proton 0 eV
1514
1515Index : 3     used in the geometry : No      recalculation needed : No
1516 Material : G4_POLYETHYLENE
1517 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1518 Energy thresholds :  gamma  2.26261 keV    e-  351.877 keV    e+  342.545 keV proton 0 eV
1519
1520Index : 4     used in the geometry : Yes     recalculation needed : No
1521 Material : Aluminium
1522 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1523 Energy thresholds :  gamma  6.90363 keV    e-  598.345 keV    e+  570.85 keV proton 0 eV
1524 Region(s) which use this couple :
1525    DefaultRegionForTheWorld
1526
1527====================================================================
1528
1529Start closing geometry.
1530G4GeometryManager::ReportVoxelStats -- Voxel Statistics
1531
1532    Total memory consumed for geometry optimisation:   0 kByte
1533    Total CPU time elapsed for geometry optimisation: 0 seconds
1534### Run 3 start.
1535
1536--------- Ranecu engine status ---------
1537 Initial seed (index) = 0
1538 Current couple of seeds = 1091561076, 776987475
1539----------------------------------------
1540Start Run processing.
1541
1542---> Begin of Event: 0
1543
1544---> Begin of Event: 10000
1545
1546---> Begin of Event: 20000
1547
1548---> Begin of Event: 30000
1549
1550---> Begin of Event: 40000
1551Run terminated.
1552Run Summary
1553  Number of events processed : 50000
1554  User=2.09s Real=3.2s Sys=1.12s
1555
1556 ======================== run summary ======================
1557
1558 The run was 50000 proton of 175 GeV through 8 mm  of Aluminium (density: 2.7 g/cm3 )
1559
1560 Total energy deposit in absorber per event = 4.881 MeV +- 246.9 keV
1561
1562 -----> Mean dE/dx = 6.098 MeV/cm       (2.259 MeV*cm2/g)
1563
1564 From formulas :
1565   restricted dEdx = 4.025 MeV/cm       (1.491 MeV*cm2/g)
1566   full dEdx       = 6.049 MeV/cm       (2.24 MeV*cm2/g)
1567
1568 Leakage :  primary = 174.7 GeV +- 246.9 keV   secondaries = 0 eV  +- 0 eV
1569 Energy balance :  edep + eleak = 174.7 GeV
1570
1571 Total track length (charged) in absorber per event = 8.004 mm  +- 0 fm
1572 Total track length (neutral) in absorber per event = 0 fm  +- 0 fm
1573
1574 Number of steps (charged) in absorber per event = 1.266 +- 0
1575 Number of steps (neutral) in absorber per event = 0 +- 0
1576
1577 Number of secondaries per event : Gammas = 0;   electrons = 0.266;   positrons = 0.00074
1578
1579 Number of events with the primary particle transmitted = 100 %
1580 Number of events with at least  1 particle transmitted (same charge as primary) = 100 %
1581
1582 Number of events with the primary particle reflected = 0 %
1583 Number of events with at least  1 particle reflected (same charge as primary) = 0 %
1584
1585 MultipleScattering:
1586  rms proj angle of transmit primary particle = 0.02106 mrad (central part only)
1587  computed theta0 (Highland formula)          = 0.02111 mrad
1588  central part defined as +- 0.06333 mrad;   Tail ratio = 1.842 %
1589
1590--------- Ranecu engine status ---------
1591 Initial seed (index) = 0
1592 Current couple of seeds = 1430401471, 1601608505
1593----------------------------------------
1594#
1595/testem/det/setAbsMat Iron
1596/testem/det/setAbsThick 1 m
1597/testem/det/update
1598
1599 Material: Galactic     density:  0.000 kg/m3   RadL: 204727512.315 pc   Nucl.Int.Length: 114561548.020 pc   Imean:  19.200 eV   temperature:   2.73 K  pressure:   0.00 atm
1600   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
1601
1602 Material:     Iron     density:  7.870 g/cm3   RadL:   1.759 cm   Nucl.Int.Length:  16.969 cm   Imean: 286.000 eV
1603   --->  Element: Iron ( )   Z = 26.0   N =  55.8   A =  55.85 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
1604
1605 The  WORLD   is made of 1.5 m   of Galactic. The transverse size (YZ) of the world is 6 cm
1606 The ABSORBER is made of 1 m   of Iron. The transverse size (YZ) is 5 cm
1607 X position of the middle of the absorber 0 fm
1608Obsolete world logical volume is removed from the default region.
1609World is registered to the default region.
1610/testem/gun/setDefault
1611/gun/particle mu+
1612/gun/energy 100 GeV
1613/run/beamOn 50000
1614
1615phot:   for  gamma    SubType= 12
1616      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1617       PhotoElectric :   Emin=        0 eV       Emax=         10 TeV
1618
1619compt:   for  gamma    SubType= 13
1620      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1621      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1622       Klein-Nishina :   Emin=        0 eV       Emax=         10 TeV
1623
1624conv:   for  gamma    SubType= 14
1625      Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1
1626      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1627        BetheHeitler :   Emin=        0 eV       Emax=         10 TeV
1628
1629Rayl:   for  gamma    SubType= 11
1630      Lambda tables from 100 eV  to 10 TeV in 200 bins, spline: 1
1631      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1632   LivermoreRayleigh :   Emin=        0 eV       Emax=        100 GeV
1633
1634msc:   for e-    SubType= 10
1635      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1636      RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
1637      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1638          UrbanMsc93 :   Emin=        0 eV       Emax=         10 TeV
1639
1640eIoni:   for  e-    SubType= 2
1641      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1642      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1643      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1644      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1645        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
1646
1647eBrem:   for  e-    SubType= 3
1648      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1649      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1650      LPM flag: 1 for E > 1 GeV
1651      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1652               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
1653            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
1654
1655eIoni:   for  e+    SubType= 2
1656      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1657      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1658      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1659      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1660        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
1661
1662eBrem:   for  e+    SubType= 3
1663      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1664      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1665      LPM flag: 1 for E > 1 GeV
1666      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1667               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
1668            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
1669
1670annihil:   for  e+    SubType= 5
1671      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1672      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1673            eplus2gg :   Emin=        0 eV       Emax=         10 TeV
1674
1675msc:   for proton    SubType= 10
1676      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1677      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1678      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1679          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
1680
1681hIoni:   for  proton    SubType= 2
1682      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1683      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1684      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1685      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1686               Bragg :   Emin=        0 eV       Emax=          2 MeV
1687          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
1688
1689hBrems:   for  proton    SubType= 3
1690      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1691      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1692      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1693               hBrem :   Emin=        0 eV       Emax=         10 TeV
1694
1695hPairProd:   for  proton    SubType= 4
1696      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1697      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1698      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1699           hPairProd :   Emin=        0 eV       Emax=         10 TeV
1700
1701msc:   for GenericIon    SubType= 10
1702      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
1703      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1704          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
1705
1706ionIoni:   for  GenericIon    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.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
1710      Stopping Power data for 17 ion/material pairs
1711      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1712         ParamICRU73 :   Emin=        0 eV       Emax=         10 TeV
1713
1714nuclearStopping:   for  GenericIon    SubType= 1
1715      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1716   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
1717
1718nuclearStopping:   for  alpha    SubType= 1
1719      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1720   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
1721
1722hIoni:   for  anti_proton    SubType= 2
1723      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1724      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1725      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1726      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1727            ICRU73QO :   Emin=        0 eV       Emax=          2 MeV
1728          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
1729
1730msc:   for kaon+    SubType= 10
1731      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1732      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1733      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1734          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
1735
1736hIoni:   for  kaon+    SubType= 2
1737      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1738      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1739      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1740      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1741               Bragg :   Emin=        0 eV       Emax=    1.05231 MeV
1742          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
1743
1744hBrems:   for  kaon+    SubType= 3
1745      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1746      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1747      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1748               hBrem :   Emin=        0 eV       Emax=         10 TeV
1749
1750hPairProd:   for  kaon+    SubType= 4
1751      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1752      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1753      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1754           hPairProd :   Emin=        0 eV       Emax=         10 TeV
1755
1756hIoni:   for  kaon-    SubType= 2
1757      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1758      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1759      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1760      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1761            ICRU73QO :   Emin=        0 eV       Emax=    1.05231 MeV
1762          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
1763
1764hBrems:   for  kaon-    SubType= 3
1765      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1766      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1767      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1768               hBrem :   Emin=        0 eV       Emax=         10 TeV
1769
1770hPairProd:   for  kaon-    SubType= 4
1771      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1772      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1773      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1774           hPairProd :   Emin=        0 eV       Emax=         10 TeV
1775
1776muMsc:   for mu+    SubType= 10
1777      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1778      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592
1779      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1780        WentzelVIUni :   Emin=        0 eV       Emax=         10 TeV
1781
1782muIoni:   for  mu+    SubType= 2
1783      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1784      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1785      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1786      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1787               Bragg :   Emin=        0 eV       Emax=        200 keV
1788          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
1789        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
1790
1791muBrems:   for  mu+    SubType= 3
1792      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1793      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1794      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1795              MuBrem :   Emin=        0 eV       Emax=         10 TeV
1796
1797muPairProd:   for  mu+    SubType= 4
1798      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1799      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1800      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1801          muPairProd :   Emin=        0 eV       Emax=         10 TeV
1802
1803CoulombScat:   for  mu+    SubType= 1
1804      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1805      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
1806      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1807  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
1808
1809muIoni:   for  mu-    SubType= 2
1810      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1811      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1812      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1813      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1814            ICRU73QO :   Emin=        0 eV       Emax=        200 keV
1815          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
1816        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
1817
1818muBrems:   for  mu-    SubType= 3
1819      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1820      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1821      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1822              MuBrem :   Emin=        0 eV       Emax=         10 TeV
1823
1824muPairProd:   for  mu-    SubType= 4
1825      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1826      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1827      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1828          muPairProd :   Emin=        0 eV       Emax=         10 TeV
1829
1830CoulombScat:   for  mu-    SubType= 1
1831      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1832      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
1833      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1834  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
1835
1836hIoni:   for  pi+    SubType= 2
1837      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1838      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1839      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1840      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1841               Bragg :   Emin=        0 eV       Emax=    297.505 keV
1842          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
1843
1844hBrems:   for  pi+    SubType= 3
1845      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1846      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1847      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1848               hBrem :   Emin=        0 eV       Emax=         10 TeV
1849
1850hPairProd:   for  pi+    SubType= 4
1851      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1852      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1853      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1854           hPairProd :   Emin=        0 eV       Emax=         10 TeV
1855
1856msc:   for pi-    SubType= 10
1857      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
1858      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1859      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1860          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
1861
1862hIoni:   for  pi-    SubType= 2
1863      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1864      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1865      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1866      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1867            ICRU73QO :   Emin=        0 eV       Emax=    297.505 keV
1868          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
1869
1870hBrems:   for  pi-    SubType= 3
1871      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1872      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1873      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1874               hBrem :   Emin=        0 eV       Emax=         10 TeV
1875
1876hPairProd:   for  pi-    SubType= 4
1877      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
1878      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1879      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
1880           hPairProd :   Emin=        0 eV       Emax=         10 TeV
1881
1882Region <DefaultRegionForTheWorld> -- appears in <World> world volume
1883 Root logical volume(s) : World
1884 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
1885 Materials : Galactic Iron
1886 Production cuts :   gamma 1 mm      e- 1 mm      e+ 1 mm  proton 0 fm
1887
1888========= Table of registered couples ==============================
1889
1890Index : 0     used in the geometry : Yes     recalculation needed : No
1891 Material : Galactic
1892 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1893 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
1894 Region(s) which use this couple :
1895    DefaultRegionForTheWorld
1896
1897Index : 1     used in the geometry : No      recalculation needed : No
1898 Material : Silicon
1899 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1900 Energy thresholds :  gamma  6.95018 keV    e-  548.291 keV    e+  526.624 keV proton 0 eV
1901
1902Index : 2     used in the geometry : No      recalculation needed : No
1903 Material : Gold
1904 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1905 Energy thresholds :  gamma  112.644 keV    e-  2.30976 MeV    e+  2.13081 MeV proton 0 eV
1906
1907Index : 3     used in the geometry : No      recalculation needed : No
1908 Material : G4_POLYETHYLENE
1909 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1910 Energy thresholds :  gamma  2.26261 keV    e-  351.877 keV    e+  342.545 keV proton 0 eV
1911
1912Index : 4     used in the geometry : No      recalculation needed : No
1913 Material : Aluminium
1914 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1915 Energy thresholds :  gamma  6.90363 keV    e-  598.345 keV    e+  570.85 keV proton 0 eV
1916
1917Index : 5     used in the geometry : Yes     recalculation needed : No
1918 Material : Iron
1919 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
1920 Energy thresholds :  gamma  20.6438 keV    e-  1.29592 MeV    e+  1.21169 MeV proton 0 eV
1921 Region(s) which use this couple :
1922    DefaultRegionForTheWorld
1923
1924====================================================================
1925
1926Start closing geometry.
1927G4GeometryManager::ReportVoxelStats -- Voxel Statistics
1928
1929    Total memory consumed for geometry optimisation:   0 kByte
1930    Total CPU time elapsed for geometry optimisation: 0 seconds
1931### Run 4 start.
1932
1933--------- Ranecu engine status ---------
1934 Initial seed (index) = 0
1935 Current couple of seeds = 1430401471, 1601608505
1936----------------------------------------
1937Start Run processing.
1938
1939---> Begin of Event: 0
1940
1941---> Begin of Event: 10000
1942
1943---> Begin of Event: 20000
1944
1945---> Begin of Event: 30000
1946
1947---> Begin of Event: 40000
1948Run terminated.
1949Run Summary
1950  Number of events processed : 50000
1951  User=22.69s Real=27.71s Sys=1.28s
1952
1953 ======================== run summary ======================
1954
1955 The run was 50000 mu+ of 100 GeV through 1 m   of Iron (density: 7.87 g/cm3 )
1956
1957 Total energy deposit in absorber per event = 2.128 GeV +- 14.93 MeV
1958
1959 -----> Mean dE/dx = 21.28 MeV/cm       (2.704 MeV*cm2/g)
1960
1961 From formulas :
1962   restricted dEdx = 11.22 MeV/cm       (1.425 MeV*cm2/g)
1963   full dEdx       = 21.6 MeV/cm        (2.745 MeV*cm2/g)
1964
1965 Leakage :  primary = 97.87 GeV +- 14.93 MeV   secondaries = 0 eV  +- 0 eV
1966 Energy balance :  edep + eleak = 100 GeV
1967
1968 Total track length (charged) in absorber per event = 99.01 cm  +- 52.39 um
1969 Total track length (neutral) in absorber per event = 0 fm  +- 0 fm
1970
1971 Number of steps (charged) in absorber per event = 47.4 +- 0.001024
1972 Number of steps (neutral) in absorber per event = 0 +- 0
1973
1974 Number of secondaries per event : Gammas = 0.04828;   electrons = 45.8;   positrons = 1.249
1975
1976 Number of events with the primary particle transmitted = 100 %
1977 Number of events with at least  1 particle transmitted (same charge as primary) = 100 %
1978
1979 Number of events with the primary particle reflected = 0 %
1980 Number of events with at least  1 particle reflected (same charge as primary) = 0 %
1981
1982 MultipleScattering:
1983  rms proj angle of transmit primary particle = 1.126 mrad (central part only)
1984  computed theta0 (Highland formula)          = 1.182 mrad
1985  central part defined as +- 3.545 mrad;   Tail ratio = 0.308 %
1986
1987--------- Ranecu engine status ---------
1988 Initial seed (index) = 0
1989 Current couple of seeds = 1153020195, 500937775
1990----------------------------------------
1991#
1992UserDetectorConstruction deleted.
1993UserPhysicsList deleted.
1994UserRunAction deleted.
1995UserPrimaryGenerator deleted.
1996G4 kernel has come to Quit state.
1997EventManager deleted.
1998UImanager deleted.
1999Units table cleared.
2000StateManager deleted.
2001RunManagerKernel is deleted.
2002RunManager is deleting.
2003
Note: See TracBrowser for help on using the repository browser.