source: trunk/examples/extended/electromagnetic/TestEm9/TestEm9.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: 18.2 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 1
14#/tracking/verbose 1
15#
16#  select EM PhysicsList
17#
18/testem/phys/addPhysics emstandard_opt1
19PhysicsList::Set emstandard_opt1 EM physics
20#
21/testem/phys/setCuts 0.3 mm
22/testem/phys/VertexCuts 0.01 mm
23/testem/phys/MuonCuts 10.0 mm
24/testem/stepMax 100 mm
25#
26/testem/det/acceptance1 0.8271 0.02974 2.5
27/testem/det/acceptance9 0.9570 0.01286 2.5
28/testem/det/acceptance25 0.9793 0.00884 2.5
29#
30#/process/eLoss/preciseRange  true
31#/process/eLoss/integral      true
32/run/initialize
33Ecal is 36 cm  of G4_CESIUM_IODIDE
34Absorber is 2 mm  of G4_Al
35Vertex is 3 cm  of 3 layers of Si of 100 um  npads= 601
36### New geometry is constructed
37#
38#/testem/histo/fileName  e_pr.root
39#/testem/histo/fileType  root
40/gun/particle e-
41/gun/energy 1 GeV
42/run/beamOn 100
43
44phot:   for  gamma, applyCuts: 1     SubType= 12
45      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
46       PhotoElectric :     Emin=          0 eV         Emax=   10 TeV
47
48compt:   for  gamma, applyCuts: 1     SubType= 13
49      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
50      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
51       Klein-Nishina :     Emin=          0 eV         Emax=   10 TeV
52
53conv:   for  gamma, applyCuts: 1     SubType= 14
54      Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1
55      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
56       Bethe-Heitler :     Emin=          0 eV         Emax=   10 TeV
57
58msc:   for e-    SubType= 10
59      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
60      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
61      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
62          UrbanMsc93 :     Emin=          0 eV         Emax=   10 TeV
63
64eIoni:   for  e-    SubType= 2
65      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
66      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
67      finalRange(mm)= 1, dRoverRange= 0.8, integral: 1, fluct: 1, linLossLimit= 0.01
68      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
69        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
70
71eBrem:   for  e-    SubType= 3
72      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
73      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
74      LPM flag: 1 for E > 1 GeV
75      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
76               eBrem :     Emin=          0 eV         Emax=   1 GeV
77            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
78
79eIoni:   for  e+    SubType= 2
80      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
81      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
82      finalRange(mm)= 1, dRoverRange= 0.8, integral: 1, fluct: 1, linLossLimit= 0.01
83      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
84        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
85
86eBrem:   for  e+    SubType= 3
87      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
88      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
89      LPM flag: 1 for E > 1 GeV
90      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
91               eBrem :     Emin=          0 eV         Emax=   1 GeV
92            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
93
94annihil:   for  e+, applyCuts: 1     SubType= 5
95      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
96      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
97            eplus2gg :     Emin=          0 eV         Emax=   10 TeV
98
99msc:   for proton    SubType= 10
100      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
101      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
102      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
103          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
104
105hIoni:   for  proton    SubType= 2
106      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
107      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
108      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
109      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
110               Bragg :     Emin=          0 eV         Emax=   2 MeV
111          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
112
113hBrems:   for  proton    SubType= 3
114      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
115      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
116      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
117               hBrem :     Emin=          0 eV         Emax=   10 TeV
118
119hPairProd:   for  proton    SubType= 4
120      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
121      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
122      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
123           hPairProd :     Emin=          0 eV         Emax=   10 TeV
124
125msc:   for GenericIon    SubType= 10
126      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
127      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
128          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
129
130ionIoni:   for  GenericIon    SubType= 2
131      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
132      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
133      finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
134      Stopping Power data for 17 ion/material pairs, nuclearStopping: 1
135      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
136            BraggIon :     Emin=          0 eV         Emax=   2 MeV
137          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
138
139hIoni:   for  anti_proton    SubType= 2
140      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
141      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
142      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
143      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
144            ICRU73QO :     Emin=          0 eV         Emax=   2 MeV
145          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
146
147hIoni:   for  kaon+    SubType= 2
148      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
149      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
150      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
151      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
152               Bragg :     Emin=          0 eV         Emax=   1.05231 MeV
153          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
154
155hBrems:   for  kaon+    SubType= 3
156      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
157      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
158      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
159               hBrem :     Emin=          0 eV         Emax=   10 TeV
160
161hPairProd:   for  kaon+    SubType= 4
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
165           hPairProd :     Emin=          0 eV         Emax=   10 TeV
166
167hIoni:   for  kaon-    SubType= 2
168      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
169      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
170      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
171      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
172            ICRU73QO :     Emin=          0 eV         Emax=   1.05231 MeV
173          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
174
175hBrems:   for  kaon-    SubType= 3
176      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
177      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
178      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
179               hBrem :     Emin=          0 eV         Emax=   10 TeV
180
181hPairProd:   for  kaon-    SubType= 4
182      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
183      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
184      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
185           hPairProd :     Emin=          0 eV         Emax=   10 TeV
186
187muMsc:   for mu+    SubType= 10
188      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
189      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0
190      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
191          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
192
193muIoni:   for  mu+    SubType= 2
194      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
195      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
196      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
197      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
198               Bragg :     Emin=          0 eV         Emax=   200 keV
199          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
200        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
201
202muBrems:   for  mu+    SubType= 3
203      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
204      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
205      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
206              MuBrem :     Emin=          0 eV         Emax=   10 TeV
207
208muPairProd:   for  mu+    SubType= 4
209      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
210      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
211      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
212          muPairProd :     Emin=          0 eV         Emax=   10 TeV
213
214muIoni:   for  mu-    SubType= 2
215      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
216      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
217      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
218      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
219            ICRU73QO :     Emin=          0 eV         Emax=   200 keV
220          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
221        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
222
223muBrems:   for  mu-    SubType= 3
224      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
225      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
226      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
227              MuBrem :     Emin=          0 eV         Emax=   10 TeV
228
229muPairProd:   for  mu-    SubType= 4
230      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
231      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
232      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
233          muPairProd :     Emin=          0 eV         Emax=   10 TeV
234
235hIoni:   for  pi+    SubType= 2
236      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
237      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
238      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
239      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
240               Bragg :     Emin=          0 eV         Emax=   297.504 keV
241          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
242
243hBrems:   for  pi+    SubType= 3
244      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
245      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
246      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
247               hBrem :     Emin=          0 eV         Emax=   10 TeV
248
249hPairProd:   for  pi+    SubType= 4
250      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
251      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
252      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
253           hPairProd :     Emin=          0 eV         Emax=   10 TeV
254
255msc:   for pi-    SubType= 10
256      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
257      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
258      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
259          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
260
261hIoni:   for  pi-    SubType= 2
262      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
263      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
264      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
265      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
266            ICRU73QO :     Emin=          0 eV         Emax=   297.504 keV
267          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
268
269hBrems:   for  pi-    SubType= 3
270      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
271      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
272      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
273               hBrem :     Emin=          0 eV         Emax=   10 TeV
274
275hPairProd:   for  pi-    SubType= 4
276      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
277      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
278      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
279           hPairProd :     Emin=          0 eV         Emax=   10 TeV
280
281========= Table of registered couples ==============================
282
283Index : 0     used in the geometry : Yes     recalculation needed : No
284 Material : G4_AIR
285 Range cuts        :  gamma  300 um     e-  300 um     e+  300 um  proton 0 fm
286 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
287 Region(s) which use this couple :
288    DefaultRegionForTheWorld
289
290Index : 1     used in the geometry : Yes     recalculation needed : No
291 Material : G4_CESIUM_IODIDE
292 Range cuts        :  gamma  300 um     e-  300 um     e+  300 um  proton 0 fm
293 Energy thresholds :  gamma  18.2919 keV    e-  291.524 keV    e+  281.891 keV proton 0 eV
294 Region(s) which use this couple :
295    DefaultRegionForTheWorld
296
297Index : 2     used in the geometry : Yes     recalculation needed : No
298 Material : G4_Al
299 Range cuts        :  gamma  300 um     e-  300 um     e+  300 um  proton 0 fm
300 Energy thresholds :  gamma  3.92577 keV    e-  258.31 keV    e+  253.154 keV proton 0 eV
301 Region(s) which use this couple :
302    DefaultRegionForTheWorld
303
304Index : 3     used in the geometry : Yes     recalculation needed : No
305 Material : G4_AIR
306 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
307 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
308 Region(s) which use this couple :
309    VertexDetector
310
311Index : 4     used in the geometry : Yes     recalculation needed : No
312 Material : G4_Al
313 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
314 Energy thresholds :  gamma  990 eV     e-  34.1725 keV    e+  33.9436 keV proton 0 eV
315 Region(s) which use this couple :
316    VertexDetector
317
318Index : 5     used in the geometry : Yes     recalculation needed : No
319 Material : G4_Si
320 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
321 Energy thresholds :  gamma  990 eV     e-  31.9516 keV    e+  31.7376 keV proton 0 eV
322 Region(s) which use this couple :
323    VertexDetector
324
325Index : 6     used in the geometry : Yes     recalculation needed : No
326 Material : G4_Fe
327 Range cuts        :  gamma  1 cm     e-  1 cm     e+  1 cm  proton 0 fm
328 Energy thresholds :  gamma  60.499 keV    e-  13.0779 MeV    e+  12.2279 MeV proton 0 eV
329 Region(s) which use this couple :
330    MuonDetector
331
332Index : 7     used in the geometry : Yes     recalculation needed : No
333 Material : G4_Al
334 Range cuts        :  gamma  1 cm     e-  1 cm     e+  1 cm  proton 0 fm
335 Energy thresholds :  gamma  20.3682 keV    e-  4.52289 MeV    e+  4.28616 MeV proton 0 eV
336 Region(s) which use this couple :
337    MuonDetector
338
339====================================================================
340
341### Run 0 start
342HistoManager: Histograms are booked and run has been started
343Start Run processing.
344Run terminated.
345Run Summary
346  Number of events processed : 100
347  User=1.85s Real=2.86s Sys=0s
348RunAction: End of run actions are started
349HistoManager: End of run actions are started
350=================================================================
351Number of events             100
352Average number of e-         404.1
353Average number of gamma      680.3
354Average number of e+         48.44
355Average number of steps      3032
356Edep 1x1 =                   0.8268 +- 0.003033  res=  3.576 %
357Edep 3x3 =                   0.9564 +- 0.00132  res=  1.36 %
358Edep 5x5 =                   0.98 +- 0.0007972  res=  0.8052 %
359===========  Mean values without trancating =====================
360Edep 1x1 =                   0.8215 +- 0.003784  res=  4.607 %
361Edep 3x3 =                   0.9552 +- 0.001452  res=  1.52 %
362Edep 5x5 =                   0.9793 +- 0.0009604  res=  0.9807 %
363===========  Ratios without trancating ===========================
364  E1/E9  =                   0.8599 +- 0.003448
365  E1/E25 =                   0.8388 +- 0.003691
366  E9/E25 =                   0.9754 +- 0.001015
367Beam Energy                  1 GeV
368==================================================================
369
370
371<<<<<ACCEPTANCE>>>>> 100 events for Crystal Calorimeter
372Edep1x1: 0.8268  delEdep1x1= -0.0003204 nrms= -0.1077
373Erms1x1: 0.02957  delErms1x1= -0.0001739 nrms= -0.05847
374Edep3x3: 0.9564  delEdep3x3= -0.0006442 nrms= -0.5009
375Erms3x3: 0.01301  delErms3x3= 0.0001451 nrms= 0.1128
376Edep5x5: 0.98  delEdep5x5= 0.0007303 nrms= 0.8261
377Erms5x5: 0.007892  delErms5x5= -0.0009484 nrms= -1.073
378<<<<<END>>>>>   IS ACCEPTED
379
380   Z  bremsstrahlung photoeffect  compton    conversion
381  13           0           0         158           0
382  53         284          41           0          45
383  55         298          48           0          49
384#
385G4 kernel has come to Quit state.
386
Note: See TracBrowser for help on using the repository browser.