source: trunk/examples/extended/electromagnetic/TestEm10/TestEm10.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.8 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 hist file = g4.p11
14 Nb of bins in Edep plot = 65
15 Elow  in the  Edep plot = 0
16 Ehigh in the  Edep plot = 0.1
17Material is selected: Mylar
18Material is selected: Air
19Material is selected: Al
20Material is selected: Xe15CO2
21
22 The  WORLD   is made of 4000mm of Air, the transverse size (R) of the world is 200 mm.
23 The ABSORBER is made of 38.3mm of Xe15CO2, the transverse size (R) is 100 mm.
24 Z position of the (middle of the) absorber 293.401  mm.
25fRadZ = 129.58 mm
26fStartZ = 100 mm
27fRadThick = 59.16 mm
28fFoilNumber = 220
29fRadiatorMat = radiatorMat
30WorldMaterial = Air
31
32XTR model = transpM
33### G4VXTRenergyLoss: the number of TR radiator plates = 220
34total radiator thickness = 5.94 cm
35plate material = Mylar
36gas material = Air
37plate plasma energy = 24.5068 eV
38gas plasma energy = 0.731967 eV
39Regular transparent X-ray TR  radiator EM process is called
40Em10PhysicsList::SetCuts:CutLength for e-, e+ and gamma is: 1 mm
41Radiator gamma cut    = 1 mm
42Radiator electron cut = 1 mm
43Radiator positron cut = 1 mm
44Radiator cuts are set
45Detector gamma cut    = 1 mm
46Detector electron cut = 1 mm
47Detector positron cut = 1 mm
48Detector cuts are set
49G4PAIModel::Initialise for e+
50G4PAIModel::Initialise for e-
51G4PAIModel::Initialise for kaon+
52G4PAIModel::Initialise for kaon-
53G4PAIModel::Initialise for mu+
54G4PAIModel::Initialise for mu-
55G4PAIModel::Initialise for pi+
56G4PAIModel::Initialise for pi-
57G4PAIModel::Initialise for proton
58
59phot:   for  gamma, applyCuts: 1     SubType= 12
60      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
61       PhotoElectric :     Emin=          0 eV         Emax=   10 TeV
62
63compt:   for  gamma, applyCuts: 1     SubType= 13
64      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
65      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
66       Klein-Nishina :     Emin=          0 eV         Emax=   10 TeV
67
68conv:   for  gamma, applyCuts: 1     SubType= 14
69      Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1
70      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
71       Bethe-Heitler :     Emin=          0 eV         Emax=   10 TeV
72
73msc:   for e-    SubType= 10
74      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
75      RangeFactor= 0.04, stepLimitType: 1, latDisplacement: 1
76      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
77          UrbanMsc93 :     Emin=          0 eV         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.2, integral: 1, fluct: 1, linLossLimit= 0.01
83      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
84        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
85      ===== EM models for the G4Region  XTRdEdxDetector ======
86            PAIModel :     Emin=          0 eV         Emax=   10 TeV
87
88eBrem:   for  e-    SubType= 3
89      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
90      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
91      LPM flag: 1 for E > 1 GeV
92      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
93               eBrem :     Emin=          0 eV         Emax=   1 GeV
94            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
95
96Lorentz Factor  XTR photon number
97
98107.6   0.0003511
99123.4   0.0006267
100141.5   0.00108
101162.3   0.001512
102186.2   0.002884
103213.7   0.004693
104245.2   0.006089
105281.3   0.009623
106322.9   0.01836
107370.5   0.03094
108425.3   0.04304
109488.2   0.05196
110560.3   0.06382
111643.2   0.08555
112738.3   0.1271
113847.6   0.1977
114973     0.3108
1151117    0.4511
1161282    0.5943
1171472    0.74
1181690    1.109
1191940    1.509
1202228    1.863
1212558    2.155
1222936    2.388
1233371    2.571
1243871    2.71
1254444    2.816
1265102    2.897
1275858    2.958
1286726    3.004
1297722    3.039
1308866    3.066
1311.018e+04       3.086
1321.169e+04       3.102
1331.342e+04       3.113
1341.541e+04       3.122
1351.769e+04       3.129
1362.031e+04       3.134
1372.332e+04       3.138
1382.677e+04       3.141
1393.074e+04       3.143
1403.529e+04       3.145
1414.052e+04       3.146
1424.652e+04       3.147
1435.342e+04       3.148
1446.133e+04       3.148
1457.042e+04       3.149
1468.085e+04       3.149
1479.283e+04       3.149
148
149total time for build X-ray TR energy loss tables = 0.1 s
150
151SynRad:  Incoherent Synchrotron Radiation
152      good description for long magnets at all energies
153
154eIoni:   for  e+    SubType= 2
155      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
156      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
157      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
158      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
159        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
160      ===== EM models for the G4Region  XTRdEdxDetector ======
161            PAIModel :     Emin=          0 eV         Emax=   10 TeV
162
163eBrem:   for  e+    SubType= 3
164      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
165      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
166      LPM flag: 1 for E > 1 GeV
167      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
168               eBrem :     Emin=          0 eV         Emax=   1 GeV
169            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
170
171annihil:   for  e+, applyCuts: 1     SubType= 5
172      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
173      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
174            eplus2gg :     Emin=          0 eV         Emax=   10 TeV
175
176Lorentz Factor  XTR photon number
177
178107.6   0.0003511
179123.4   0.0006267
180141.5   0.00108
181162.3   0.001512
182186.2   0.002884
183213.7   0.004693
184245.2   0.006089
185281.3   0.009623
186322.9   0.01836
187370.5   0.03094
188425.3   0.04304
189488.2   0.05196
190560.3   0.06382
191643.2   0.08555
192738.3   0.1271
193847.6   0.1977
194973     0.3108
1951117    0.4511
1961282    0.5943
1971472    0.74
1981690    1.109
1991940    1.509
2002228    1.863
2012558    2.155
2022936    2.388
2033371    2.571
2043871    2.71
2054444    2.816
2065102    2.897
2075858    2.958
2086726    3.004
2097722    3.039
2108866    3.066
2111.018e+04       3.086
2121.169e+04       3.102
2131.342e+04       3.113
2141.541e+04       3.122
2151.769e+04       3.129
2162.031e+04       3.134
2172.332e+04       3.138
2182.677e+04       3.141
2193.074e+04       3.143
2203.529e+04       3.145
2214.052e+04       3.146
2224.652e+04       3.147
2235.342e+04       3.148
2246.133e+04       3.148
2257.042e+04       3.149
2268.085e+04       3.149
2279.283e+04       3.149
228
229total time for build X-ray TR energy loss tables = 0.1 s
230
231msc:   for proton    SubType= 10
232      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
233      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
234      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
235          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
236
237hIoni:   for  proton    SubType= 2
238      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
239      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
240      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
241      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
242               Bragg :     Emin=          0 eV         Emax=   2 MeV
243          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
244      ===== EM models for the G4Region  XTRdEdxDetector ======
245            PAIModel :     Emin=          0 eV         Emax=   10 TeV
246
247hIoni:   for  kaon+    SubType= 2
248      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
249      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
250      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
251      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
252               Bragg :     Emin=          0 eV         Emax=   1.05231 MeV
253          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
254      ===== EM models for the G4Region  XTRdEdxDetector ======
255            PAIModel :     Emin=          0 eV         Emax=   10 TeV
256
257hIoni:   for  kaon-    SubType= 2
258      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
259      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
260      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
261      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
262            ICRU73QO :     Emin=          0 eV         Emax=   1.05231 MeV
263          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
264      ===== EM models for the G4Region  XTRdEdxDetector ======
265            PAIModel :     Emin=          0 eV         Emax=   10 TeV
266
267muMsc:   for mu+    SubType= 10
268      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
269      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0
270      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
271          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
272
273muIoni:   for  mu+    SubType= 2
274      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
275      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
276      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
277      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
278               Bragg :     Emin=          0 eV         Emax=   200 keV
279          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
280        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
281      ===== EM models for the G4Region  XTRdEdxDetector ======
282            PAIModel :     Emin=          0 eV         Emax=   10 TeV
283
284muBrems:   for  mu+    SubType= 3
285      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
286      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
287      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
288              MuBrem :     Emin=          0 eV         Emax=   10 TeV
289
290muPairProd:   for  mu+    SubType= 4
291      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
292      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
293      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
294          muPairProd :     Emin=          0 eV         Emax=   10 TeV
295
296muIoni:   for  mu-    SubType= 2
297      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
298      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
299      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
300      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
301            ICRU73QO :     Emin=          0 eV         Emax=   200 keV
302          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
303        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
304      ===== EM models for the G4Region  XTRdEdxDetector ======
305            PAIModel :     Emin=          0 eV         Emax=   10 TeV
306
307muBrems:   for  mu-    SubType= 3
308      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
309      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
310      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
311              MuBrem :     Emin=          0 eV         Emax=   10 TeV
312
313muPairProd:   for  mu-    SubType= 4
314      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
315      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
316      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
317          muPairProd :     Emin=          0 eV         Emax=   10 TeV
318
319hIoni:   for  pi+    SubType= 2
320      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
321      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
322      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
323      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
324               Bragg :     Emin=          0 eV         Emax=   297.504 keV
325          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
326      ===== EM models for the G4Region  XTRdEdxDetector ======
327            PAIModel :     Emin=          0 eV         Emax=   10 TeV
328
329msc:   for pi-    SubType= 10
330      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
331      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
332      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
333          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
334
335hIoni:   for  pi-    SubType= 2
336      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
337      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
338      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
339      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
340            ICRU73QO :     Emin=          0 eV         Emax=   297.504 keV
341          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
342      ===== EM models for the G4Region  XTRdEdxDetector ======
343            PAIModel :     Emin=          0 eV         Emax=   10 TeV
344
345========= Table of registered couples ==============================
346
347Index : 0     used in the geometry : Yes     recalculation needed : No
348 Material : Air
349 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
350 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
351 Region(s) which use this couple :
352    DefaultRegionForTheWorld
353
354Index : 1     used in the geometry : Yes     recalculation needed : No
355 Material : radiatorMat
356 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
357 Energy thresholds :  gamma  1.07314 keV    e-  83.5287 keV    e+  82.4135 keV proton 0 eV
358 Region(s) which use this couple :
359    XTRradiator
360
361Index : 2     used in the geometry : Yes     recalculation needed : No
362 Material : Xe15CO2
363 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
364 Energy thresholds :  gamma  990 eV     e-  2.08861 keV    e+  2.0697 keV proton 0 eV
365 Region(s) which use this couple :
366    XTRdEdxDetector
367
368====================================================================
369
370### Run 0 start.
371
372--------- Ranecu engine status ---------
373 Initial seed (index) = 0
374 Current couple of seeds = 9876, 54321
375----------------------------------------
376Start Run processing.
377
378---> Begin of Event: 0
379Run terminated.
380Run Summary
381  Number of events processed : 1000
382  User=1.02s Real=1.46s Sys=0s
383 ================== run summary =====================
384 end of Run TotNbofEvents = 1000
385    mean charged track length   in absorber=44.0381 +- 0.536306  mm 
386
387            mean energy deposit in absorber=0.0524177 +- 0.00090763  MeV
388
389 mean number of steps in absorber (charged) =8.968 +- 0.267748     
390 mean number of steps in absorber (neutral) =2.813 +- 0.0525741     
391
392   mean number of charged secondaries = 3.505 +- 0.076681
393
394   mean number of neutral secondaries = 0.037 +- 0.00613441
395
396   mean number of e-s =3.505  and e+s =0
397
398(number) transmission coeff=0.514  reflection coeff=0.02
399
400 energy deposit distribution
401#entries=1000    #underflows=0    #overflows=61
402 bin nb      Elow      entries     normalized
403    0         0           0           0
404    1   1.53846           0           0
405    2   3.07692           0           0
406    3   4.61538           0           0
407    4   6.15385           0           0
408    5   7.69231           0           0
409    6   9.23077           0           0
410    7   10.7692           0           0
411    8   12.3077           0           0
412    9   13.8462           1       0.001
413   10   15.3846           6       0.006
414   11   16.9231          18       0.018
415   12   18.4615          14       0.014
416   13        20          19       0.019
417   14   21.5385          25       0.025
418   15   23.0769          11       0.011
419   16   24.6154          23       0.023
420   17   26.1538          34       0.034
421   18   27.6923          28       0.028
422   19   29.2308          34       0.034
423   20   30.7692          24       0.024
424   21   32.3077          36       0.036
425   22   33.8462          17       0.017
426   23   35.3846          38       0.038
427   24   36.9231          43       0.043
428   25   38.4615          29       0.029
429   26        40          34       0.034
430   27   41.5385          20        0.02
431   28   43.0769          29       0.029
432   29   44.6154          24       0.024
433   30   46.1538          29       0.029
434   31   47.6923          22       0.022
435   32   49.2308          23       0.023
436   33   50.7692          17       0.017
437   34   52.3077          30        0.03
438   35   53.8462          24       0.024
439   36   55.3846          22       0.022
440   37   56.9231          18       0.018
441   38   58.4615          19       0.019
442   39        60          13       0.013
443   40   61.5385          16       0.016
444   41   63.0769          13       0.013
445   42   64.6154          16       0.016
446   43   66.1538          10        0.01
447   44   67.6923          19       0.019
448   45   69.2308           9       0.009
449   46   70.7692          10        0.01
450   47   72.3077          11       0.011
451   48   73.8462          10        0.01
452   49   75.3846          10        0.01
453   50   76.9231           7       0.007
454   51   78.4615           9       0.009
455   52        80           8       0.008
456   53   81.5385           7       0.007
457   54   83.0769           4       0.004
458   55   84.6154           7       0.007
459   56   86.1538           8       0.008
460   57   87.6923           2       0.002
461   58   89.2308           9       0.009
462   59   90.7692           5       0.005
463   60   92.3077           6       0.006
464   61   93.8462           4       0.004
465   62   95.3846           3       0.003
466   63   96.9231           8       0.008
467   64   98.4615           4       0.004
468
469 Emp =       0.0369231   width=      0.0338462  MeV
470
471
472--------- Ranecu engine status ---------
473 Initial seed (index) = 0
474 Current couple of seeds = 368886938, 2074936418
475----------------------------------------
476
477========= Table of registered couples ==============================
478
479Index : 0     used in the geometry : Yes     recalculation needed : No
480 Material : Air
481 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
482 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
483 Region(s) which use this couple :
484    DefaultRegionForTheWorld
485
486Index : 1     used in the geometry : Yes     recalculation needed : No
487 Material : radiatorMat
488 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
489 Energy thresholds :  gamma  1.0731416 keV    e-  83.52866 keV    e+  82.413535 keV proton 0 eV
490 Region(s) which use this couple :
491    XTRradiator
492
493Index : 2     used in the geometry : Yes     recalculation needed : No
494 Material : Xe15CO2
495 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
496 Energy thresholds :  gamma  990 eV     e-  2.0886137 keV    e+  2.0697013 keV proton 0 eV
497 Region(s) which use this couple :
498    XTRdEdxDetector
499
500====================================================================
501
502G4 kernel has come to Quit state.
503
Note: See TracBrowser for help on using the repository browser.