source: trunk/examples/extended/electromagnetic/TestEm2/TestEm2.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.7 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
14***** Table : Nb of materials = 9 *****
15
16 Material:    Water     density:  1.000 g/cm3   RadL:  36.092 cm   Nucl.Int.Length:  75.416 cm   Imean:  75.000 eV
17   --->  Element: Hydrogen (H)   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction:  11.21 %  ElmAbundance  66.67 %
18   --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  88.79 %  ElmAbundance  33.33 %
19
20 Material: liquidArgon     density:  1.390 g/cm3   RadL:  14.065 cm   Nucl.Int.Length:  86.078 cm   Imean: 188.000 eV
21   --->  Element: liquidArgon ( )   Z = 18.0   N =  40.0   A =  39.95 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
22
23 Material: Aluminium     density:  2.700 g/cm3   RadL:   8.893 cm   Nucl.Int.Length:  38.879 cm   Imean: 166.000 eV
24   --->  Element: Aluminium ( )   Z = 13.0   N =  27.0   A =  26.98 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
25
26 Material:     Iron     density:  7.870 g/cm3   RadL:   1.759 cm   Nucl.Int.Length:  16.999 cm   Imean: 286.000 eV
27   --->  Element: Iron ( )   Z = 26.0   N =  55.8   A =  55.85 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
28
29 Material:   Copper     density:  8.960 g/cm3   RadL:   1.436 cm   Nucl.Int.Length:  15.588 cm   Imean: 322.000 eV
30   --->  Element: Copper ( )   Z = 29.0   N =  63.5   A =  63.55 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
31
32 Material: Tungsten     density: 19.350 g/cm3   RadL:   3.495 mm   Nucl.Int.Length:  10.285 cm   Imean: 727.000 eV
33   --->  Element: Tungsten ( )   Z = 74.0   N = 183.8   A = 183.84 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
34
35 Material:     Lead     density: 11.350 g/cm3   RadL:   5.612 mm   Nucl.Int.Length:  18.247 cm   Imean: 823.000 eV
36   --->  Element: Lead ( )   Z = 82.0   N = 207.2   A = 207.19 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
37
38 Material:  Uranium     density: 18.950 g/cm3   RadL:   3.166 mm   Nucl.Int.Length:  11.446 cm   Imean: 890.000 eV
39   --->  Element: Uranium ( )   Z = 92.0   N = 238.0   A = 238.03 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
40
41 Material:      BGO     density:  7.100 g/cm3   RadL:   1.123 cm   Nucl.Int.Length:  22.806 cm   Imean: 473.785 eV
42   --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  15.41 %  ElmAbundance  63.16 %
43   --->  Element: Germanium (Ge)   Z = 32.0   N =  72.6   A =  72.59 g/mole  ElmMassFraction:  17.48 %  ElmAbundance  15.79 %
44   --->  Element: Bismuth (Bi)   Z = 83.0   N = 209.0   A = 208.98 g/mole  ElmMassFraction:  67.10 %  ElmAbundance  21.05 %
45
46
47/run/verbose 1
48/testem/det/setMat G4_PbWO4
49/testem/det/setLbin 40 0.5
50/testem/det/setRbin 50 0.1
51/testem/run/acceptance 0.949 0.0082 3.0
52/testem/phys/addPhysics emstandard_opt0
53PhysicsList::AddPhysicsList: <emstandard_opt0>
54/testem/phys/setCuts 1.0 mm
55/run/initialize
56Absorber is 17.8491 cm  of G4_PbWO4
57 Material: G4_PbWO4     density:  8.280 g/cm3   RadL:   8.925 mm   Nucl.Int.Length:  20.737 cm   Imean: 542.741 eV
58   --->  Element: O (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole
59   --->  Isotope:   O16   Z =  8   N =  16   A =  15.99 g/mole   abundance:  99.76 %
60   --->  Isotope:   O17   Z =  8   N =  17   A =  17.00 g/mole   abundance:   0.04 %
61   --->  Isotope:   O18   Z =  8   N =  18   A =  18.00 g/mole   abundance:   0.20 %  ElmMassFraction:  14.06 %  ElmAbundance  66.67 %
62   --->  Element: Pb (Pb)   Z = 82.0   N = 207.2   A = 207.22 g/mole
63   --->  Isotope: Pb204   Z = 82   N = 204   A = 203.97 g/mole   abundance:   1.40 %
64   --->  Isotope: Pb206   Z = 82   N = 206   A = 205.97 g/mole   abundance:  24.10 %
65   --->  Isotope: Pb207   Z = 82   N = 207   A = 206.98 g/mole   abundance:  22.10 %
66   --->  Isotope: Pb208   Z = 82   N = 208   A = 207.98 g/mole   abundance:  52.40 %  ElmMassFraction:  45.54 %  ElmAbundance  16.67 %
67   --->  Element: W (W)   Z = 74.0   N = 183.9   A = 183.84 g/mole
68   --->  Isotope:  W180   Z = 74   N = 180   A = 179.95 g/mole   abundance:   0.12 %
69   --->  Isotope:  W182   Z = 74   N = 182   A = 181.95 g/mole   abundance:  26.50 %
70   --->  Isotope:  W183   Z = 74   N = 183   A = 182.95 g/mole   abundance:  14.31 %
71   --->  Isotope:  W184   Z = 74   N = 184   A = 183.95 g/mole   abundance:  30.64 %
72   --->  Isotope:  W186   Z = 74   N = 186   A = 185.95 g/mole   abundance:  28.43 %  ElmMassFraction:  40.40 %  ElmAbundance  16.67 %
73PhysicsList::SetCuts:CutLength : 1 mm
74/process/eLoss/verbose 1
75/gun/particle e-
76/gun/energy 5 GeV
77/testem/event/printModulo 20
78/run/beamOn 100
79
80phot:   for  gamma    SubType= 12
81      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
82       PhotoElectric :     Emin=          0 eV         Emax=   10 TeV
83
84compt:   for  gamma    SubType= 13
85      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
86      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
87       Klein-Nishina :     Emin=          0 eV         Emax=   10 TeV
88
89conv:   for  gamma    SubType= 14
90      Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1
91      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
92       Bethe-Heitler :     Emin=          0 eV         Emax=   10 TeV
93
94msc:   for e-    SubType= 10
95      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
96      RangeFactor= 0.04, stepLimitType: 1, latDisplacement: 1
97      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
98          UrbanMsc93 :     Emin=          0 eV         Emax=   10 TeV
99
100eIoni:   for  e-    SubType= 2
101      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
102      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
103      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
104      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
105        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
106
107eBrem:   for  e-    SubType= 3
108      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
109      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
110      LPM flag: 1 for E > 1 GeV
111      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
112               eBrem :     Emin=          0 eV         Emax=   1 GeV
113            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
114
115eIoni:   for  e+    SubType= 2
116      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
117      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
118      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
119      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
120        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
121
122eBrem:   for  e+    SubType= 3
123      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
124      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
125      LPM flag: 1 for E > 1 GeV
126      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
127               eBrem :     Emin=          0 eV         Emax=   1 GeV
128            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
129
130annihil:   for  e+    SubType= 5
131      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
132      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
133            eplus2gg :     Emin=          0 eV         Emax=   10 TeV
134
135msc:   for proton    SubType= 10
136      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
137      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
138      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
139          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
140
141hIoni:   for  proton    SubType= 2
142      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
143      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
144      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
145      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
146               Bragg :     Emin=          0 eV         Emax=   2 MeV
147          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
148
149hBrems:   for  proton    SubType= 3
150      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
151      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
152      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
153               hBrem :     Emin=          0 eV         Emax=   10 TeV
154
155hPairProd:   for  proton    SubType= 4
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           hPairProd :     Emin=          0 eV         Emax=   10 TeV
160
161msc:   for GenericIon    SubType= 10
162      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
163      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
164          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
165
166ionIoni:   for  GenericIon    SubType= 2
167      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
168      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
169      finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
170      Stopping Power data for 17 ion/material pairs, nuclearStopping: 1
171      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
172            BraggIon :     Emin=          0 eV         Emax=   2 MeV
173          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
174
175hIoni:   for  anti_proton    SubType= 2
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      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
179      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
180            ICRU73QO :     Emin=          0 eV         Emax=   2 MeV
181          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
182
183hIoni:   for  kaon+    SubType= 2
184      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
185      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
186      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
187      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
188               Bragg :     Emin=          0 eV         Emax=   1.05231 MeV
189          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
190
191hBrems:   for  kaon+    SubType= 3
192      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
193      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
194      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
195               hBrem :     Emin=          0 eV         Emax=   10 TeV
196
197hPairProd:   for  kaon+    SubType= 4
198      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
199      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
200      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
201           hPairProd :     Emin=          0 eV         Emax=   10 TeV
202
203hIoni:   for  kaon-    SubType= 2
204      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
205      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
206      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
207      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
208            ICRU73QO :     Emin=          0 eV         Emax=   1.05231 MeV
209          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
210
211hBrems:   for  kaon-    SubType= 3
212      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
213      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
214      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
215               hBrem :     Emin=          0 eV         Emax=   10 TeV
216
217hPairProd:   for  kaon-    SubType= 4
218      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
219      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
220      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
221           hPairProd :     Emin=          0 eV         Emax=   10 TeV
222
223muMsc:   for mu+    SubType= 10
224      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
225      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0
226      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
227          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
228
229muIoni:   for  mu+    SubType= 2
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      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
233      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
234               Bragg :     Emin=          0 eV         Emax=   200 keV
235          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
236        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
237
238muBrems:   for  mu+    SubType= 3
239      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
240      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
241      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
242              MuBrem :     Emin=          0 eV         Emax=   10 TeV
243
244muPairProd:   for  mu+    SubType= 4
245      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
246      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
247      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
248          muPairProd :     Emin=          0 eV         Emax=   10 TeV
249
250muIoni:   for  mu-    SubType= 2
251      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
252      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
253      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
254      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
255            ICRU73QO :     Emin=          0 eV         Emax=   200 keV
256          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
257        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
258
259muBrems:   for  mu-    SubType= 3
260      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
261      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
262      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
263              MuBrem :     Emin=          0 eV         Emax=   10 TeV
264
265muPairProd:   for  mu-    SubType= 4
266      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
267      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
268      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
269          muPairProd :     Emin=          0 eV         Emax=   10 TeV
270
271hIoni:   for  pi+    SubType= 2
272      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
273      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
274      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
275      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
276               Bragg :     Emin=          0 eV         Emax=   297.504 keV
277          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
278
279hBrems:   for  pi+    SubType= 3
280      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
281      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
282      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
283               hBrem :     Emin=          0 eV         Emax=   10 TeV
284
285hPairProd:   for  pi+    SubType= 4
286      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
287      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
288      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
289           hPairProd :     Emin=          0 eV         Emax=   10 TeV
290
291msc:   for pi-    SubType= 10
292      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
293      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
294      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
295          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
296
297hIoni:   for  pi-    SubType= 2
298      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
299      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
300      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
301      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
302            ICRU73QO :     Emin=          0 eV         Emax=   297.504 keV
303          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
304
305hBrems:   for  pi-    SubType= 3
306      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
307      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
308      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
309               hBrem :     Emin=          0 eV         Emax=   10 TeV
310
311hPairProd:   for  pi-    SubType= 4
312      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
313      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
314      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
315           hPairProd :     Emin=          0 eV         Emax=   10 TeV
316
317========= Table of registered couples ==============================
318
319Index : 0     used in the geometry : Yes     recalculation needed : No
320 Material : G4_PbWO4
321 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
322 Energy thresholds :  gamma  89.3347 keV    e-  1.12535 MeV    e+  1.05931 MeV proton 0 eV
323 Region(s) which use this couple :
324    DefaultRegionForTheWorld
325
326====================================================================
327
328### Run 0 start.
329
330--------- Ranecu engine status ---------
331 Initial seed (index) = 0
332 Current couple of seeds = 9876, 54321
333----------------------------------------
334Start Run processing.
335
336---> Begin of Event: 0
337
338---> Begin of Event: 20
339
340---> Begin of Event: 40
341
342---> Begin of Event: 60
343
344---> Begin of Event: 80
345Run terminated.
346Run Summary
347  Number of events processed : 100
348  User=7.07s Real=9.12s Sys=0s
349
350 SUMMARY
351
352 energy deposit :   94.69 % E0 +-    1.11 % E0
353 charged traklen:  479.10 radl +-    5.93 radl
354 neutral traklen: 4445.58 radl +-  130.36 radl
355
356 90.00 % confinement: radius = 3.10 radl  (2.76 cm )
357
358--------- Ranecu engine status ---------
359 Initial seed (index) = 0
360 Current couple of seeds = 859713455, 2085840081
361----------------------------------------
362
363<<<<ACCEPTANCE>>>> 100 events for Total Energy in Absorber
364Edep: 0.94689776  delEdep= -0.0021022368 nrms= -2.5637034
365Erms: 0.01113771  delErms= 0.0029377098 nrms= 3.5825729
366<<<<END>>>>   IS ACCEPTED
367
368G4 kernel has come to Quit state.
369
Note: See TracBrowser for help on using the repository browser.