source: trunk/examples/extended/electromagnetic/TestEm3/TestEm3.out @ 1339

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

tag geant4.9.4 beta 1 + modifs locales

File size: 19.5 KB
Line 
1
2        ############################################
3        !!! WARNING - FPE detection is activated !!!
4        ############################################
5
6*************************************************************
7 Geant4 version Name: geant4-09-03-ref-06    (25-June-2010)
8                      Copyright : Geant4 Collaboration
9                      Reference : NIM A 506 (2003), 250-303
10                            WWW : http://cern.ch/geant4
11*************************************************************
12
13/run/verbose 2
14#
15/testem/phys/addPhysics emstandard_opt0
16#
17/run/initialize
18userDetector->Construct() start.
19
20-------------------------------------------------------------
21 ---> The calorimeter is 50 layers of:
22                Lead:    2.3 mm
23         liquidArgon:    5.7 mm
24-------------------------------------------------------------
25
26 Material: Galactic     density:  0.000 kg/m3   RadL: 204322111.300 pc   Nucl.Int.Length: 113728944.943 pc   Imean:  19.200 eV   temperature:   2.73 K  pressure:   0.00 atm
27   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
28
29 Material:     Lead     density: 11.350 g/cm3   RadL:   5.613 mm   Nucl.Int.Length:  18.247 cm   Imean: 823.000 eV
30   --->  Element: Pb (Pb)   Z = 82.0   N = 207.2   A = 207.22 g/mole
31   --->  Isotope: Pb204   Z = 82   N = 204   A = 203.97 g/mole   abundance:   1.40 %
32   --->  Isotope: Pb206   Z = 82   N = 206   A = 205.97 g/mole   abundance:  24.10 %
33   --->  Isotope: Pb207   Z = 82   N = 207   A = 206.98 g/mole   abundance:  22.10 %
34   --->  Isotope: Pb208   Z = 82   N = 208   A = 207.98 g/mole   abundance:  52.40 %  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
35
36 Material: liquidArgon     density:  1.390 g/cm3   RadL:  14.064 cm   Nucl.Int.Length:  86.022 cm   Imean: 188.000 eV
37   --->  Element: Ar (Ar)   Z = 18.0   N =  40.0   A =  39.95 g/mole
38   --->  Isotope:  Ar36   Z = 18   N =  36   A =  35.97 g/mole   abundance:   0.34 %
39   --->  Isotope:  Ar38   Z = 18   N =  38   A =  37.96 g/mole   abundance:   0.06 %
40   --->  Isotope:  Ar40   Z = 18   N =  40   A =  39.96 g/mole   abundance:  99.60 %  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
41
42-------------------------------------------------------------
43World is registered to the default region.
44physicsList->Construct() start.
45physicsList->Construct() start.
46physicsList->CheckParticleList() start.
47physicsList->setCut() start.
48PhysicsList::SetCuts:CutLength : 1 mm
49/process/eLoss/verbose 1
50/process/msc/Skin 1
51/run/physicsModified
52#
53/gun/particle e-
54/gun/energy 1 GeV
55#
56/testem/histo/setFileName testem3
57/testem/histo/setFileType root
58/testem/histo/setHisto 1 100 700. 1000. MeV
59----> SetHisto 1: Edep in absorber 1 (MeV);  100 bins from 700 MeV to 1000 MeV
60/testem/histo/setHisto 2 100 0. 300. MeV
61----> SetHisto 2: Edep in absorber 2 (MeV);  100 bins from 0 MeV to 300 MeV
62/testem/histo/setHisto 11 52 0. 52. none
63----> SetHisto 11: longit. profile of Edep (MeV/event) in absorber 1;  52 bins from 0 none to 52 none
64/testem/histo/setHisto 12 52 0. 52. none
65----> SetHisto 12: longit. profile of Edep (MeV/event) in absorber 2;  52 bins from 0 none to 52 none
66/testem/histo/setHisto 21 102 0. 102. none
67----> SetHisto 21: energy flow (MeV/event);  102 bins from 0 none to 102 none
68/testem/histo/setHisto 22 102 0. 102. none
69----> SetHisto 22: lateral energy leak (MeV/event);  102 bins from 0 none to 102 none
70#
71/testem/event/printModulo 50
72#
73/run/beamOn 100
74
75phot:   for  gamma    SubType= 12
76      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
77       PhotoElectric :     Emin=          0 eV         Emax=   10 TeV
78
79compt:   for  gamma    SubType= 13
80      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
81      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
82       Klein-Nishina :     Emin=          0 eV         Emax=   10 TeV
83
84conv:   for  gamma    SubType= 14
85      Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1
86      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
87       Bethe-Heitler :     Emin=          0 eV         Emax=   10 TeV
88
89msc:   for e-    SubType= 10
90      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
91      RangeFactor= 0.04, stepLimitType: 1, latDisplacement: 1
92      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
93          UrbanMsc93 :     Emin=          0 eV         Emax=   10 TeV
94
95eIoni:   for  e-    SubType= 2
96      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
97      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
98      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
99      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
100        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
101
102eBrem:   for  e-    SubType= 3
103      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
104      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
105      LPM flag: 1 for E > 1 GeV
106      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
107               eBrem :     Emin=          0 eV         Emax=   1 GeV
108            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
109
110eIoni:   for  e+    SubType= 2
111      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
112      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
113      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
114      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
115        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
116
117eBrem:   for  e+    SubType= 3
118      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
119      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
120      LPM flag: 1 for E > 1 GeV
121      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
122               eBrem :     Emin=          0 eV         Emax=   1 GeV
123            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
124
125annihil:   for  e+    SubType= 5
126      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
127      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
128            eplus2gg :     Emin=          0 eV         Emax=   10 TeV
129
130msc:   for proton    SubType= 10
131      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
132      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
133      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
134          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
135
136hIoni:   for  proton    SubType= 2
137      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
138      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
139      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
140      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
141               Bragg :     Emin=          0 eV         Emax=   2 MeV
142          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
143
144hBrems:   for  proton    SubType= 3
145      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
146      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
147      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
148               hBrem :     Emin=          0 eV         Emax=   10 TeV
149
150hPairProd:   for  proton    SubType= 4
151      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
152      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
153      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
154           hPairProd :     Emin=          0 eV         Emax=   10 TeV
155
156msc:   for GenericIon    SubType= 10
157      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
158      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
159          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
160
161ionIoni:   for  GenericIon    SubType= 2
162      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
163      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
164      finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
165      Stopping Power data for 17 ion/material pairs, nuclearStopping: 1
166      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
167            BraggIon :     Emin=          0 eV         Emax=   2 MeV
168          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
169
170hIoni:   for  anti_proton    SubType= 2
171      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
172      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
173      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
174      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
175            ICRU73QO :     Emin=          0 eV         Emax=   2 MeV
176          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
177
178hIoni:   for  kaon+    SubType= 2
179      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
180      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
181      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
182      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
183               Bragg :     Emin=          0 eV         Emax=   1.05231 MeV
184          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
185
186hBrems:   for  kaon+    SubType= 3
187      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
188      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
189      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
190               hBrem :     Emin=          0 eV         Emax=   10 TeV
191
192hPairProd:   for  kaon+    SubType= 4
193      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
194      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
195      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
196           hPairProd :     Emin=          0 eV         Emax=   10 TeV
197
198hIoni:   for  kaon-    SubType= 2
199      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
200      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
201      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
202      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
203            ICRU73QO :     Emin=          0 eV         Emax=   1.05231 MeV
204          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
205
206hBrems:   for  kaon-    SubType= 3
207      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
208      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
209      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
210               hBrem :     Emin=          0 eV         Emax=   10 TeV
211
212hPairProd:   for  kaon-    SubType= 4
213      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
214      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
215      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
216           hPairProd :     Emin=          0 eV         Emax=   10 TeV
217
218muMsc:   for mu+    SubType= 10
219      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
220      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0
221      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
222          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
223
224muIoni:   for  mu+    SubType= 2
225      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
226      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
227      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
228      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
229               Bragg :     Emin=          0 eV         Emax=   200 keV
230          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
231        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
232
233muBrems:   for  mu+    SubType= 3
234      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
235      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
236      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
237              MuBrem :     Emin=          0 eV         Emax=   10 TeV
238
239muPairProd:   for  mu+    SubType= 4
240      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
241      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
242      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
243          muPairProd :     Emin=          0 eV         Emax=   10 TeV
244
245muIoni:   for  mu-    SubType= 2
246      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
247      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
248      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
249      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
250            ICRU73QO :     Emin=          0 eV         Emax=   200 keV
251          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
252        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
253
254muBrems:   for  mu-    SubType= 3
255      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
256      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
257      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
258              MuBrem :     Emin=          0 eV         Emax=   10 TeV
259
260muPairProd:   for  mu-    SubType= 4
261      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
262      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
263      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
264          muPairProd :     Emin=          0 eV         Emax=   10 TeV
265
266hIoni:   for  pi+    SubType= 2
267      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
268      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
269      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
270      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
271               Bragg :     Emin=          0 eV         Emax=   297.504 keV
272          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
273
274hBrems:   for  pi+    SubType= 3
275      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
276      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
277      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
278               hBrem :     Emin=          0 eV         Emax=   10 TeV
279
280hPairProd:   for  pi+    SubType= 4
281      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
282      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
283      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
284           hPairProd :     Emin=          0 eV         Emax=   10 TeV
285
286msc:   for pi-    SubType= 10
287      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
288      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
289      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
290          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
291
292hIoni:   for  pi-    SubType= 2
293      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
294      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
295      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
296      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
297            ICRU73QO :     Emin=          0 eV         Emax=   297.504 keV
298          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
299
300hBrems:   for  pi-    SubType= 3
301      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
302      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
303      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
304               hBrem :     Emin=          0 eV         Emax=   10 TeV
305
306hPairProd:   for  pi-    SubType= 4
307      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
308      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
309      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
310           hPairProd :     Emin=          0 eV         Emax=   10 TeV
311
312Region <DefaultRegionForTheWorld> -- appears in <World> world volume
313 Root logical volume(s) : World
314 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
315 Materials : Galactic Lead liquidArgon
316 Production cuts :   gamma 1 mm      e- 1 mm      e+ 1 mm  proton 1 mm
317
318========= Table of registered couples ==============================
319
320Index : 0     used in the geometry : Yes     recalculation needed : No
321 Material : Galactic
322 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
323 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 100 keV
324 Region(s) which use this couple :
325    DefaultRegionForTheWorld
326
327Index : 1     used in the geometry : Yes     recalculation needed : No
328 Material : Lead
329 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
330 Energy thresholds :  gamma  101.843 keV    e-  1.36749 MeV    e+  1.27862 MeV proton 100 keV
331 Region(s) which use this couple :
332    DefaultRegionForTheWorld
333
334Index : 2     used in the geometry : Yes     recalculation needed : No
335 Material : liquidArgon
336 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
337 Energy thresholds :  gamma  6.19986 keV    e-  344.855 keV    e+  333.46 keV proton 100 keV
338 Region(s) which use this couple :
339    DefaultRegionForTheWorld
340
341====================================================================
342
343Start closing geometry.
344G4GeometryManager::ReportVoxelStats -- Voxel Statistics
345
346    Total memory consumed for geometry optimisation:   3 kByte
347    Total CPU time elapsed for geometry optimisation: 0 seconds
348
349    Voxelisation: top CPU users:
350    Percent   Total CPU    System CPU       Memory  Volume
351    -------   ----------   ----------     --------  ----------
352       0.00         0.00         0.00            3k Calorimeter
353       0.00         0.00         0.00            0k Layer
354
355    Voxelisation: top memory users:
356    Percent     Memory      Heads    Nodes   Pointers    Total CPU    Volume
357    -------   --------     ------   ------   --------   ----------    ----------
358      92.05          2k         1       50         50         0.00    Calorimeter
359       7.95          0k         1        3          4         0.00    Layer
360### Run 0 start.
361
362--------- Ranecu engine status ---------
363 Initial seed (index) = 0
364 Current couple of seeds = 9876, 54321
365----------------------------------------
366Start Run processing.
367
368---> Begin Of Event: 0
369
370---> Begin Of Event: 50
371Run terminated.
372Run Summary
373  Number of events processed : 100
374  User=5.56s Real=9.12s Sys=0s
375
376------------------------------------------------------------
377      material   Edep       RMS         sqrt(E0(GeV))*rmsE/Emean     total tracklen
378 
379          Lead: 781.73 MeV :  14.91 MeV     1.907 +- 0.1907 %      54.5 cm  +- 1.17 cm
380   liquidArgon: 207.32 MeV :  15.05 MeV     7.262 +- 0.7262 %      1.02 m   +- 7.76 cm
381
382------------------------------------------------------------
383 Beam particle e-  E = 1 GeV
384 Mean number of gamma  570
385 Mean number of e-     960
386 Mean number of e+     54.5
387------------------------------------------------------------
388
389 Energy deposition from Energy flow balance :
390  material       Total Edep
391 
392      Lead:      781.726 MeV
393liquidArgon:     207.311 MeV
394
395------------------------------------------------------------
396
397
398--------- Ranecu engine status ---------
399 Initial seed (index) = 0
400 Current couple of seeds = 1438848957, 1887175185
401----------------------------------------
402UserDetectorConstruction deleted.
403UserPhysicsList deleted.
404UserRunAction deleted.
405UserPrimaryGenerator deleted.
406G4 kernel has come to Quit state.
407EventManager deleted.
408UImanager deleted.
409Units table cleared.
410StateManager deleted.
411RunManagerKernel is deleted.
412RunManager is deleting.
413
Note: See TracBrowser for help on using the repository browser.