source: trunk/examples/extended/electromagnetic/TestEm9/TestEm9.out @ 1342

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

update ti head

File size: 18.5 KB
Line 
1
2        ############################################
3        !!! WARNING - FPE detection is activated !!!
4        ############################################
5
6*************************************************************
7 Geant4 version Name: geant4-09-03-ref-09    (25-June-2010)
8                      Copyright : Geant4 Collaboration
9                      Reference : NIM A 506 (2003), 250-303
10                            WWW : http://cern.ch/geant4
11*************************************************************
12
13/run/verbose 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        BetheHeitler :   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  AngularGenUrban
77            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
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  AngularGenUrban
92            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
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
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
147msc:   for kaon+    SubType= 10
148      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
149      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
150      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
151          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
152
153hIoni:   for  kaon+    SubType= 2
154      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
155      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
156      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
157      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
158               Bragg :   Emin=        0 eV       Emax=    1.05231 MeV
159          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
160
161hBrems:   for  kaon+    SubType= 3
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               hBrem :   Emin=        0 eV       Emax=         10 TeV
166
167hPairProd:   for  kaon+    SubType= 4
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
171           hPairProd :   Emin=        0 eV       Emax=         10 TeV
172
173hIoni:   for  kaon-    SubType= 2
174      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
175      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
176      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
177      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
178            ICRU73QO :   Emin=        0 eV       Emax=    1.05231 MeV
179          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
180
181hBrems:   for  kaon-    SubType= 3
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               hBrem :   Emin=        0 eV       Emax=         10 TeV
186
187hPairProd:   for  kaon-    SubType= 4
188      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
189      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
190      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
191           hPairProd :   Emin=        0 eV       Emax=         10 TeV
192
193msc:   for mu+    SubType= 10
194      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
195      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
196      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
197          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
198
199muIoni:   for  mu+    SubType= 2
200      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
201      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
202      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
203      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
204               Bragg :   Emin=        0 eV       Emax=        200 keV
205          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
206        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
207
208muBrems:   for  mu+    SubType= 3
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              MuBrem :   Emin=        0 eV       Emax=         10 TeV
213
214muPairProd:   for  mu+    SubType= 4
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
218          muPairProd :   Emin=        0 eV       Emax=         10 TeV
219
220muIoni:   for  mu-    SubType= 2
221      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
222      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
223      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
224      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
225            ICRU73QO :   Emin=        0 eV       Emax=        200 keV
226          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
227        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
228
229muBrems:   for  mu-    SubType= 3
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              MuBrem :   Emin=        0 eV       Emax=         10 TeV
234
235muPairProd:   for  mu-    SubType= 4
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
239          muPairProd :   Emin=        0 eV       Emax=         10 TeV
240
241hIoni:   for  pi+    SubType= 2
242      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
243      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
244      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
245      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
246               Bragg :   Emin=        0 eV       Emax=    297.505 keV
247          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
248
249hBrems:   for  pi+    SubType= 3
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               hBrem :   Emin=        0 eV       Emax=         10 TeV
254
255hPairProd:   for  pi+    SubType= 4
256      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
257      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
258      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
259           hPairProd :   Emin=        0 eV       Emax=         10 TeV
260
261msc:   for pi-    SubType= 10
262      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
263      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
264      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
265          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
266
267hIoni:   for  pi-    SubType= 2
268      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
269      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
270      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
271      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
272            ICRU73QO :   Emin=        0 eV       Emax=    297.505 keV
273          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
274
275hBrems:   for  pi-    SubType= 3
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               hBrem :   Emin=        0 eV       Emax=         10 TeV
280
281hPairProd:   for  pi-    SubType= 4
282      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
283      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
284      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
285           hPairProd :   Emin=        0 eV       Emax=         10 TeV
286
287========= Table of registered couples ==============================
288
289Index : 0     used in the geometry : Yes     recalculation needed : No
290 Material : G4_AIR
291 Range cuts        :  gamma  300 um     e-  300 um     e+  300 um  proton 0 fm
292 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
293 Region(s) which use this couple :
294    DefaultRegionForTheWorld
295
296Index : 1     used in the geometry : Yes     recalculation needed : No
297 Material : G4_CESIUM_IODIDE
298 Range cuts        :  gamma  300 um     e-  300 um     e+  300 um  proton 0 fm
299 Energy thresholds :  gamma  18.2919 keV    e-  291.524 keV    e+  281.891 keV proton 0 eV
300 Region(s) which use this couple :
301    DefaultRegionForTheWorld
302
303Index : 2     used in the geometry : Yes     recalculation needed : No
304 Material : G4_Al
305 Range cuts        :  gamma  300 um     e-  300 um     e+  300 um  proton 0 fm
306 Energy thresholds :  gamma  3.92577 keV    e-  258.31 keV    e+  253.154 keV proton 0 eV
307 Region(s) which use this couple :
308    DefaultRegionForTheWorld
309
310Index : 3     used in the geometry : Yes     recalculation needed : No
311 Material : G4_AIR
312 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
313 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
314 Region(s) which use this couple :
315    VertexDetector
316
317Index : 4     used in the geometry : Yes     recalculation needed : No
318 Material : G4_Al
319 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
320 Energy thresholds :  gamma  990 eV     e-  34.1725 keV    e+  33.9436 keV proton 0 eV
321 Region(s) which use this couple :
322    VertexDetector
323
324Index : 5     used in the geometry : Yes     recalculation needed : No
325 Material : G4_Si
326 Range cuts        :  gamma  10 um     e-  10 um     e+  10 um  proton 0 fm
327 Energy thresholds :  gamma  990 eV     e-  31.9516 keV    e+  31.7376 keV proton 0 eV
328 Region(s) which use this couple :
329    VertexDetector
330
331Index : 6     used in the geometry : Yes     recalculation needed : No
332 Material : G4_Fe
333 Range cuts        :  gamma  1 cm     e-  1 cm     e+  1 cm  proton 0 fm
334 Energy thresholds :  gamma  60.499 keV    e-  13.0779 MeV    e+  12.2279 MeV proton 0 eV
335 Region(s) which use this couple :
336    MuonDetector
337
338Index : 7     used in the geometry : Yes     recalculation needed : No
339 Material : G4_Al
340 Range cuts        :  gamma  1 cm     e-  1 cm     e+  1 cm  proton 0 fm
341 Energy thresholds :  gamma  20.3682 keV    e-  4.52289 MeV    e+  4.28616 MeV proton 0 eV
342 Region(s) which use this couple :
343    MuonDetector
344
345====================================================================
346
347### Run 0 start
348HistoManager: Histograms are booked and run has been started
349Start Run processing.
350Run terminated.
351Run Summary
352  Number of events processed : 100
353  User=1.88s Real=1.89s Sys=0s
354RunAction: End of run actions are started
355HistoManager: End of run actions are started
356=================================================================
357Number of events             100
358Average number of e-         405.9
359Average number of gamma      680.2
360Average number of e+         47.97
361Average number of steps      3041
362Edep 1x1 =                   0.8243 +- 0.002658  res=  3.192 %
363Edep 3x3 =                   0.957 +- 0.001221  res=  1.269 %
364Edep 5x5 =                   0.9793 +- 0.0008076  res=  0.8247 %
365===========  Mean values without trancating =====================
366Edep 1x1 =                   0.8226 +- 0.002854  res=  3.469 %
367Edep 3x3 =                   0.9566 +- 0.001255  res=  1.312 %
368Edep 5x5 =                   0.9793 +- 0.0008076  res=  0.8247 %
369===========  Ratios without trancating ===========================
370  E1/E9  =                   0.8599 +- 0.002718
371  E1/E25 =                   0.84 +- 0.002731
372  E9/E25 =                   0.9768 +- 0.0008467
373Beam Energy                  1 GeV
374==================================================================
375
376
377<<<<<ACCEPTANCE>>>>> 100 events for Crystal Calorimeter
378Edep1x1: 0.8243  delEdep1x1= -0.002786 nrms= -0.9369
379Erms1x1: 0.02631  delErms1x1= -0.003426 nrms= -1.152
380Edep3x3: 0.957  delEdep3x3= -1.314e-05 nrms= -0.01022
381Erms3x3: 0.01214  delErms3x3= -0.0007156 nrms= -0.5564
382Edep5x5: 0.9793  delEdep5x5= 3.176e-05 nrms= 0.03593
383Erms5x5: 0.008076  delErms5x5= -0.0007638 nrms= -0.864
384<<<<<END>>>>>   IS ACCEPTED
385
386   Z  bremsstrahlung photoeffect  compton    conversion
387  13           0           0         160           0
388  53         282          41           0          44
389  55         301          48           0          49
390#
391G4 kernel has come to Quit state.
392
Note: See TracBrowser for help on using the repository browser.