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
RevLine 
[807]1
[1337]2        ############################################
3        !!! WARNING - FPE detection is activated !!!
4        ############################################
5
[807]6*************************************************************
[1342]7 Geant4 version Name: geant4-09-03-ref-09    (25-June-2010)
[807]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#
[1337]19/testem/phys/addPhysics emstandard_opt3
20PhysicsList::AddPhysicsList: <emstandard_opt3>
[807]21#
[1337]22/testem/phys/setCuts 7 um
[807]23#
24/run/initialize
25userDetector->Construct() start.
26
[1342]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
[807]28   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
29
[1342]30 Material:  Silicon     density:  2.330 g/cm3   RadL:   9.368 cm   Nucl.Int.Length:  45.761 cm   Imean: 173.000 eV
[807]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.
[1337]39physicsList->CheckParticleList() start.
[807]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
[1337]48/testem/histo/setFileType root
[807]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
[1337]54phot:   for  gamma    SubType= 12
55      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]56       PhotoElectric :   Emin=        0 eV       Emax=         10 TeV
[807]57
[1337]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 ======
[1342]61       Klein-Nishina :   Emin=        0 eV       Emax=         10 TeV
[807]62
[1337]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 ======
[1342]66        BetheHeitler :   Emin=        0 eV       Emax=         10 TeV
[807]67
[1337]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 ======
[1342]71   LivermoreRayleigh :   Emin=        0 eV       Emax=        100 GeV
[807]72
[1337]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 ======
[1342]77          UrbanMsc93 :   Emin=        0 eV       Emax=         10 TeV
[807]78
[1337]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 ======
[1342]84        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
[807]85
[1337]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 ======
[1342]91               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
92            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
[807]93
[1337]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 ======
[1342]99        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
[807]100
[1337]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 ======
[1342]106               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
107            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
[807]108
[1337]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 ======
[1342]112            eplus2gg :   Emin=        0 eV       Emax=         10 TeV
[807]113
[1337]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 ======
[1342]118          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
[807]119
[1337]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 ======
[1342]125               Bragg :   Emin=        0 eV       Emax=          2 MeV
126          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
[807]127
[1337]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 ======
[1342]132               hBrem :   Emin=        0 eV       Emax=         10 TeV
[807]133
[1337]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 ======
[1342]138           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[807]139
[1337]140msc:   for GenericIon    SubType= 10
141      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
142      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]143          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
[807]144
[1337]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
[1342]149      Stopping Power data for 17 ion/material pairs
[1337]150      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]151         ParamICRU73 :   Emin=        0 eV       Emax=         10 TeV
[807]152
[1337]153nuclearStopping:   for  GenericIon    SubType= 1
154      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]155   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
[807]156
[1337]157nuclearStopping:   for  alpha    SubType= 1
158      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]159   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
[807]160
[1337]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 ======
[1342]166            ICRU73QO :   Emin=        0 eV       Emax=          2 MeV
167          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
[807]168
[1342]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
[1337]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 ======
[1342]180               Bragg :   Emin=        0 eV       Emax=    1.05231 MeV
181          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
[807]182
[1337]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 ======
[1342]187               hBrem :   Emin=        0 eV       Emax=         10 TeV
[807]188
[1337]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 ======
[1342]193           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[807]194
[1337]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 ======
[1342]200            ICRU73QO :   Emin=        0 eV       Emax=    1.05231 MeV
201          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
[807]202
[1337]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 ======
[1342]207               hBrem :   Emin=        0 eV       Emax=         10 TeV
[807]208
[1337]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 ======
[1342]213           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[807]214
[1337]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 ======
[1342]219        WentzelVIUni :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]226               Bragg :   Emin=        0 eV       Emax=        200 keV
227          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
228        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
[1337]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 ======
[1342]234              MuBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]240          muPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]246  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]253            ICRU73QO :   Emin=        0 eV       Emax=        200 keV
254          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
255        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
[1337]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 ======
[1342]261              MuBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]267          muPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]273  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]280               Bragg :   Emin=        0 eV       Emax=    297.505 keV
281          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
[1337]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 ======
[1342]287               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]293           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]299          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]306            ICRU73QO :   Emin=        0 eV       Emax=    297.505 keV
307          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
[1337]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 ======
[1342]313               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]319           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]320
[807]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
[1337]325 Production cuts :   gamma 7 um      e- 7 um      e+ 7 um  proton 0 fm
[807]326
327========= Table of registered couples ==============================
328
329Index : 0     used in the geometry : Yes     recalculation needed : No
330 Material : Galactic
[1337]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
[807]333 Region(s) which use this couple :
334    DefaultRegionForTheWorld
335
336Index : 1     used in the geometry : Yes     recalculation needed : No
337 Material : Silicon
[1337]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
[807]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
[1342]380  User=4.2s Real=6.41s Sys=2.21s
[807]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
[1342]386 Total energy deposit in absorber per event = 6.355 keV +- 16.41 eV
[807]387
[1342]388 -----> Mean dE/dx = 3.177 MeV/cm       (1.364 MeV*cm2/g)
[807]389
390 From formulas :
[1337]391   restricted dEdx = 2.963 MeV/cm       (1.272 MeV*cm2/g)
392   full dEdx       = 4.726 MeV/cm       (2.028 MeV*cm2/g)
[807]393
[1342]394 Leakage :  primary = 5 GeV +- 1.274 keV   secondaries = 3.64 keV +- 1.273 keV
[807]395 Energy balance :  edep + eleak = 5 GeV
396
[1342]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
[807]399
[1342]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
[807]402
[1342]403 Number of secondaries per event : Gammas = 3e-05;   electrons = 0.01661;   positrons = 0
[807]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:
[1342]412  rms proj angle of transmit primary particle = 0.02446 mrad (central part only)
[807]413  computed theta0 (Highland formula)          = 0.02627 mrad
[1342]414  central part defined as +- 0.0788 mrad;   Tail ratio = 1.553 %
[807]415
416--------- Ranecu engine status ---------
417 Initial seed (index) = 0
[1342]418 Current couple of seeds = 1096967505, 715451002
[807]419----------------------------------------
[1337]420#
421/testem/det/setAbsMat Gold
422/testem/det/setAbsThick 9.658 um
423/testem/det/setAbsYZ 1 mm
424/testem/det/update
425
[1342]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
[1337]427   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
428
[1342]429 Material:     Gold     density: 19.320 g/cm3   RadL:   3.344 mm   Nucl.Int.Length:  10.539 cm   Imean: 790.000 eV
[1337]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 ======
[1342]444       PhotoElectric :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]449       Klein-Nishina :   Emin=        0 eV       Emax=         10 TeV
[1337]450
451conv:   for  gamma    SubType= 14
[1342]452      Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1
[1337]453      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]454        BetheHeitler :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]459   LivermoreRayleigh :   Emin=        0 eV       Emax=        100 GeV
[1337]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 ======
[1342]465          UrbanMsc93 :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]472        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]479               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
480            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
[1337]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 ======
[1342]487        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]494               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
495            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
[1337]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 ======
[1342]500            eplus2gg :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]506          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]513               Bragg :   Emin=        0 eV       Emax=          2 MeV
514          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
[1337]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 ======
[1342]520               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]526           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]527
528msc:   for GenericIon    SubType= 10
529      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
530      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]531          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
[1337]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
[1342]537      Stopping Power data for 17 ion/material pairs
[1337]538      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]539         ParamICRU73 :   Emin=        0 eV       Emax=         10 TeV
[1337]540
541nuclearStopping:   for  GenericIon    SubType= 1
542      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]543   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
[1337]544
545nuclearStopping:   for  alpha    SubType= 1
546      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]547   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]554            ICRU73QO :   Emin=        0 eV       Emax=          2 MeV
555          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
[1337]556
[1342]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
[1337]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 ======
[1342]568               Bragg :   Emin=        0 eV       Emax=    1.05231 MeV
569          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
[1337]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 ======
[1342]575               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]581           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]588            ICRU73QO :   Emin=        0 eV       Emax=    1.05231 MeV
589          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
[1337]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 ======
[1342]595               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]601           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]602
603muMsc:   for mu+    SubType= 10
604      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
[1342]605      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592
[1337]606      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]607        WentzelVIUni :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]614               Bragg :   Emin=        0 eV       Emax=        200 keV
615          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
616        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
[1337]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 ======
[1342]622              MuBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]628          muPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]629
630CoulombScat:   for  mu+    SubType= 1
631      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
[1342]632      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
[1337]633      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]634  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]641            ICRU73QO :   Emin=        0 eV       Emax=        200 keV
642          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
643        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
[1337]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 ======
[1342]649              MuBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]655          muPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]656
657CoulombScat:   for  mu-    SubType= 1
658      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
[1342]659      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
[1337]660      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]661  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]668               Bragg :   Emin=        0 eV       Emax=    297.505 keV
669          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
[1337]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 ======
[1342]675               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]681           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]687          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]694            ICRU73QO :   Emin=        0 eV       Emax=    297.505 keV
695          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
[1337]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 ======
[1342]701               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]707           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]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
[1342]727 Energy thresholds :  gamma  990 eV     e-  21.4886 keV    e+  20.8954 keV proton 0 eV
[1337]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
[1342]732 Energy thresholds :  gamma  6.54229 keV    e-  65.1406 keV    e+  63.8405 keV proton 0 eV
[1337]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
[1342]747 Current couple of seeds = 1096967505, 715451002
[1337]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
[1342]763  User=3.35s Real=4.56s Sys=1.21s
[1337]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
[1342]769 Total energy deposit in absorber per event = 18.54 keV +- 65.36 eV
[1337]770
[1342]771 -----> Mean dE/dx = 19.2 MeV/cm        (0.9938 MeV*cm2/g)
[1337]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
[1342]777 Leakage :  primary = 15.63 MeV +- 2.365 keV   secondaries = 47.06 keV +- 2.361 keV
[1337]778 Energy balance :  edep + eleak = 15.7 MeV
779
[1342]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
[1337]782
[1342]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
[1337]785
[1342]786 Number of secondaries per event : Gammas = 0.02986;   electrons = 0.023;   positrons = 0
[1337]787
[1342]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 %
[1337]790
[1342]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 %
[1337]793
794 MultipleScattering:
[1342]795  rms proj angle of transmit primary particle = 32.82 mrad (central part only)
[1337]796  computed theta0 (Highland formula)          = 35.1 mrad
[1342]797  central part defined as +- 105.3 mrad;   Tail ratio = 2.774 %
[1337]798
799--------- Ranecu engine status ---------
800 Initial seed (index) = 0
[1342]801 Current couple of seeds = 1228101883, 1478657312
[1337]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
[1342]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
[1337]812   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
813
[1342]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
[1337]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 ======
[1342]834       PhotoElectric :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]839       Klein-Nishina :   Emin=        0 eV       Emax=         10 TeV
[1337]840
841conv:   for  gamma    SubType= 14
[1342]842      Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1
[1337]843      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]844        BetheHeitler :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]849   LivermoreRayleigh :   Emin=        0 eV       Emax=        100 GeV
[1337]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 ======
[1342]855          UrbanMsc93 :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]862        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]869               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
870            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
[1337]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 ======
[1342]877        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]884               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
885            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
[1337]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 ======
[1342]890            eplus2gg :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]896          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]903               Bragg :   Emin=        0 eV       Emax=          2 MeV
904          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
[1337]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 ======
[1342]910               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]916           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]917
918msc:   for GenericIon    SubType= 10
919      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
920      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]921          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
[1337]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
[1342]927      Stopping Power data for 17 ion/material pairs
[1337]928      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]929         ParamICRU73 :   Emin=        0 eV       Emax=         10 TeV
[1337]930
931nuclearStopping:   for  GenericIon    SubType= 1
932      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]933   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
[1337]934
935nuclearStopping:   for  alpha    SubType= 1
936      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]937   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]944            ICRU73QO :   Emin=        0 eV       Emax=          2 MeV
945          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
[1337]946
[1342]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
[1337]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 ======
[1342]958               Bragg :   Emin=        0 eV       Emax=    1.05231 MeV
959          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
[1337]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 ======
[1342]965               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]971           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]978            ICRU73QO :   Emin=        0 eV       Emax=    1.05231 MeV
979          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
[1337]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 ======
[1342]985               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]991           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]992
993muMsc:   for mu+    SubType= 10
994      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
[1342]995      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592
[1337]996      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]997        WentzelVIUni :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1004               Bragg :   Emin=        0 eV       Emax=        200 keV
1005          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
1006        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
[1337]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 ======
[1342]1012              MuBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1018          muPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]1019
1020CoulombScat:   for  mu+    SubType= 1
1021      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
[1342]1022      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
[1337]1023      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1024  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1031            ICRU73QO :   Emin=        0 eV       Emax=        200 keV
1032          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
1033        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
[1337]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 ======
[1342]1039              MuBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1045          muPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]1046
1047CoulombScat:   for  mu-    SubType= 1
1048      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
[1342]1049      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
[1337]1050      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1051  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1058               Bragg :   Emin=        0 eV       Emax=    297.505 keV
1059          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
[1337]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 ======
[1342]1065               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1071           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1077          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1084            ICRU73QO :   Emin=        0 eV       Emax=    297.505 keV
1085          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
[1337]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 ======
[1342]1091               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1097           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]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
[1342]1117 Energy thresholds :  gamma  990 eV     e-  31.9516 keV    e+  31.7376 keV proton 0 eV
[1337]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
[1342]1122 Energy thresholds :  gamma  8.27707 keV    e-  80.7687 keV    e+  79.1567 keV proton 0 eV
[1337]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
[1342]1127 Energy thresholds :  gamma  990 eV     e-  14.1748 keV    e+  13.7548 keV proton 0 eV
[1337]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
[1342]1142 Current couple of seeds = 1228101883, 1478657312
[1337]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
[1342]1158  User=6.36s Real=7.53s Sys=1.16s
[1337]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
[1342]1164 Total energy deposit in absorber per event = 1.031 MeV +- 1.011 keV
[1337]1165
[1342]1166 -----> Mean dE/dx = 2.175 MeV/cm       (2.314 MeV*cm2/g)
[1337]1167
1168 From formulas :
[1342]1169   restricted dEdx = 1.747 MeV/cm       (1.859 MeV*cm2/g)
1170   full dEdx       = 2.259 MeV/cm       (2.403 MeV*cm2/g)
[1337]1171
[1342]1172 Leakage :  primary = 95.14 MeV +- 1.348 keV   secondaries = 26.07 keV +- 707.1 eV
[1337]1173 Energy balance :  edep + eleak = 96.2 MeV
1174
[1342]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
[1337]1177
[1342]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
[1337]1180
[1342]1181 Number of secondaries per event : Gammas = 0.01128;   electrons = 3.989;   positrons = 0
[1337]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:
[1342]1190  rms proj angle of transmit primary particle = 7.275 mrad (central part only)
[1337]1191  computed theta0 (Highland formula)          = 7.636 mrad
[1342]1192  central part defined as +- 22.91 mrad;   Tail ratio = 1.259 %
[1337]1193
1194--------- Ranecu engine status ---------
1195 Initial seed (index) = 0
[1342]1196 Current couple of seeds = 1091561076, 776987475
[1337]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
[1342]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
[1337]1208   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
1209
[1342]1210 Material: Aluminium     density:  2.700 g/cm3   RadL:   8.893 cm   Nucl.Int.Length:  38.860 cm   Imean: 166.000 eV
[1337]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 ======
[1342]1225       PhotoElectric :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1230       Klein-Nishina :   Emin=        0 eV       Emax=         10 TeV
[1337]1231
1232conv:   for  gamma    SubType= 14
[1342]1233      Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1
[1337]1234      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1235        BetheHeitler :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1240   LivermoreRayleigh :   Emin=        0 eV       Emax=        100 GeV
[1337]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 ======
[1342]1246          UrbanMsc93 :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1253        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1260               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
1261            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
[1337]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 ======
[1342]1268        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1275               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
1276            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
[1337]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 ======
[1342]1281            eplus2gg :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1287          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1294               Bragg :   Emin=        0 eV       Emax=          2 MeV
1295          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
[1337]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 ======
[1342]1301               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1307           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]1308
1309msc:   for GenericIon    SubType= 10
1310      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
1311      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1312          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
[1337]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
[1342]1318      Stopping Power data for 17 ion/material pairs
[1337]1319      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1320         ParamICRU73 :   Emin=        0 eV       Emax=         10 TeV
[1337]1321
1322nuclearStopping:   for  GenericIon    SubType= 1
1323      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1324   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
[1337]1325
1326nuclearStopping:   for  alpha    SubType= 1
1327      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1328   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1335            ICRU73QO :   Emin=        0 eV       Emax=          2 MeV
1336          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
[1337]1337
[1342]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
[1337]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 ======
[1342]1349               Bragg :   Emin=        0 eV       Emax=    1.05231 MeV
1350          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
[1337]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 ======
[1342]1356               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1362           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1369            ICRU73QO :   Emin=        0 eV       Emax=    1.05231 MeV
1370          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
[1337]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 ======
[1342]1376               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1382           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]1383
1384muMsc:   for mu+    SubType= 10
1385      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
[1342]1386      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592
[1337]1387      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1388        WentzelVIUni :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1395               Bragg :   Emin=        0 eV       Emax=        200 keV
1396          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
1397        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
[1337]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 ======
[1342]1403              MuBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1409          muPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]1410
1411CoulombScat:   for  mu+    SubType= 1
1412      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
[1342]1413      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
[1337]1414      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1415  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1422            ICRU73QO :   Emin=        0 eV       Emax=        200 keV
1423          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
1424        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
[1337]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 ======
[1342]1430              MuBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1436          muPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]1437
1438CoulombScat:   for  mu-    SubType= 1
1439      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
[1342]1440      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
[1337]1441      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1442  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1449               Bragg :   Emin=        0 eV       Emax=    297.505 keV
1450          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
[1337]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 ======
[1342]1456               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1462           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1468          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1475            ICRU73QO :   Emin=        0 eV       Emax=    297.505 keV
1476          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
[1337]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 ======
[1342]1482               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1488           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]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
[1342]1508 Energy thresholds :  gamma  6.95018 keV    e-  548.291 keV    e+  526.624 keV proton 0 eV
[1337]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
[1342]1513 Energy thresholds :  gamma  112.644 keV    e-  2.30976 MeV    e+  2.13081 MeV proton 0 eV
[1337]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
[1342]1518 Energy thresholds :  gamma  2.26261 keV    e-  351.877 keV    e+  342.545 keV proton 0 eV
[1337]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
[1342]1523 Energy thresholds :  gamma  6.90363 keV    e-  598.345 keV    e+  570.85 keV proton 0 eV
[1337]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
[1342]1538 Current couple of seeds = 1091561076, 776987475
[1337]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
[1342]1554  User=2.09s Real=3.2s Sys=1.12s
[1337]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
[1342]1560 Total energy deposit in absorber per event = 4.881 MeV +- 246.9 keV
[1337]1561
[1342]1562 -----> Mean dE/dx = 6.098 MeV/cm       (2.259 MeV*cm2/g)
[1337]1563
1564 From formulas :
[1342]1565   restricted dEdx = 4.025 MeV/cm       (1.491 MeV*cm2/g)
[1337]1566   full dEdx       = 6.049 MeV/cm       (2.24 MeV*cm2/g)
1567
[1342]1568 Leakage :  primary = 174.7 GeV +- 246.9 keV   secondaries = 0 eV  +- 0 eV
[1337]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
[1342]1574 Number of steps (charged) in absorber per event = 1.266 +- 0
[1337]1575 Number of steps (neutral) in absorber per event = 0 +- 0
1576
[1342]1577 Number of secondaries per event : Gammas = 0;   electrons = 0.266;   positrons = 0.00074
[1337]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:
[1342]1586  rms proj angle of transmit primary particle = 0.02106 mrad (central part only)
[1337]1587  computed theta0 (Highland formula)          = 0.02111 mrad
[1342]1588  central part defined as +- 0.06333 mrad;   Tail ratio = 1.842 %
[1337]1589
1590--------- Ranecu engine status ---------
1591 Initial seed (index) = 0
[1342]1592 Current couple of seeds = 1430401471, 1601608505
[1337]1593----------------------------------------
1594#
1595/testem/det/setAbsMat Iron
1596/testem/det/setAbsThick 1 m
1597/testem/det/update
1598
[1342]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
[1337]1600   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
1601
[1342]1602 Material:     Iron     density:  7.870 g/cm3   RadL:   1.759 cm   Nucl.Int.Length:  16.969 cm   Imean: 286.000 eV
[1337]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 ======
[1342]1617       PhotoElectric :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1622       Klein-Nishina :   Emin=        0 eV       Emax=         10 TeV
[1337]1623
1624conv:   for  gamma    SubType= 14
[1342]1625      Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1
[1337]1626      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1627        BetheHeitler :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1632   LivermoreRayleigh :   Emin=        0 eV       Emax=        100 GeV
[1337]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 ======
[1342]1638          UrbanMsc93 :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1645        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1652               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
1653            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
[1337]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 ======
[1342]1660        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1667               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
1668            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
[1337]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 ======
[1342]1673            eplus2gg :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1679          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1686               Bragg :   Emin=        0 eV       Emax=          2 MeV
1687          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
[1337]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 ======
[1342]1693               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1699           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]1700
1701msc:   for GenericIon    SubType= 10
1702      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
1703      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1704          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
[1337]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
[1342]1710      Stopping Power data for 17 ion/material pairs
[1337]1711      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1712         ParamICRU73 :   Emin=        0 eV       Emax=         10 TeV
[1337]1713
1714nuclearStopping:   for  GenericIon    SubType= 1
1715      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1716   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
[1337]1717
1718nuclearStopping:   for  alpha    SubType= 1
1719      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1720   ICRU49NucStopping :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1727            ICRU73QO :   Emin=        0 eV       Emax=          2 MeV
1728          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
[1337]1729
[1342]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
[1337]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 ======
[1342]1741               Bragg :   Emin=        0 eV       Emax=    1.05231 MeV
1742          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
[1337]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 ======
[1342]1748               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1754           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1761            ICRU73QO :   Emin=        0 eV       Emax=    1.05231 MeV
1762          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
[1337]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 ======
[1342]1768               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1774           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]1775
1776muMsc:   for mu+    SubType= 10
1777      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
[1342]1778      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592
[1337]1779      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1780        WentzelVIUni :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1787               Bragg :   Emin=        0 eV       Emax=        200 keV
1788          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
1789        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
[1337]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 ======
[1342]1795              MuBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1801          muPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]1802
1803CoulombScat:   for  mu+    SubType= 1
1804      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
[1342]1805      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
[1337]1806      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1807  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1814            ICRU73QO :   Emin=        0 eV       Emax=        200 keV
1815          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
1816        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
[1337]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 ======
[1342]1822              MuBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1828          muPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]1829
1830CoulombScat:   for  mu-    SubType= 1
1831      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
[1342]1832      11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
[1337]1833      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]1834  eCoulombScattering :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1841               Bragg :   Emin=        0 eV       Emax=    297.505 keV
1842          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
[1337]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 ======
[1342]1848               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1854           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1860          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1867            ICRU73QO :   Emin=        0 eV       Emax=    297.505 keV
1868          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
[1337]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 ======
[1342]1874               hBrem :   Emin=        0 eV       Emax=         10 TeV
[1337]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 ======
[1342]1880           hPairProd :   Emin=        0 eV       Emax=         10 TeV
[1337]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
[1342]1900 Energy thresholds :  gamma  6.95018 keV    e-  548.291 keV    e+  526.624 keV proton 0 eV
[1337]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
[1342]1905 Energy thresholds :  gamma  112.644 keV    e-  2.30976 MeV    e+  2.13081 MeV proton 0 eV
[1337]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
[1342]1910 Energy thresholds :  gamma  2.26261 keV    e-  351.877 keV    e+  342.545 keV proton 0 eV
[1337]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
[1342]1915 Energy thresholds :  gamma  6.90363 keV    e-  598.345 keV    e+  570.85 keV proton 0 eV
[1337]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
[1342]1920 Energy thresholds :  gamma  20.6438 keV    e-  1.29592 MeV    e+  1.21169 MeV proton 0 eV
[1337]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
[1342]1935 Current couple of seeds = 1430401471, 1601608505
[1337]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
[1342]1951  User=22.69s Real=27.71s Sys=1.28s
[1337]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
[1342]1957 Total energy deposit in absorber per event = 2.128 GeV +- 14.93 MeV
[1337]1958
[1342]1959 -----> Mean dE/dx = 21.28 MeV/cm       (2.704 MeV*cm2/g)
[1337]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
[1342]1965 Leakage :  primary = 97.87 GeV +- 14.93 MeV   secondaries = 0 eV  +- 0 eV
[1337]1966 Energy balance :  edep + eleak = 100 GeV
1967
[1342]1968 Total track length (charged) in absorber per event = 99.01 cm  +- 52.39 um
[1337]1969 Total track length (neutral) in absorber per event = 0 fm  +- 0 fm
1970
[1342]1971 Number of steps (charged) in absorber per event = 47.4 +- 0.001024
[1337]1972 Number of steps (neutral) in absorber per event = 0 +- 0
1973
[1342]1974 Number of secondaries per event : Gammas = 0.04828;   electrons = 45.8;   positrons = 1.249
[1337]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:
[1342]1983  rms proj angle of transmit primary particle = 1.126 mrad (central part only)
[1337]1984  computed theta0 (Highland formula)          = 1.182 mrad
[1342]1985  central part defined as +- 3.545 mrad;   Tail ratio = 0.308 %
[1337]1986
1987--------- Ranecu engine status ---------
1988 Initial seed (index) = 0
[1342]1989 Current couple of seeds = 1153020195, 500937775
[1337]1990----------------------------------------
1991#
[807]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.
[1337]2003
Note: See TracBrowser for help on using the repository browser.