source: trunk/examples/extended/electromagnetic/TestEm12/TestEm12.out @ 1337

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

tag geant4.9.4 beta 1 + modifs locales
File size: 19.0 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 = 8 *****
15
16 Material:    G4_Al     density:  2.699 g/cm3   RadL:   8.896 cm   Nucl.Int.Length:  38.894 cm   Imean: 166.000 eV
17   --->  Element: Al (Al)   Z = 13.0   N =  27.0   A =  26.98 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
18
19 Material:    G4_Si     density:  2.330 g/cm3   RadL:   9.366 cm   Nucl.Int.Length:  45.660 cm   Imean: 173.000 eV
20   --->  Element: Si (Si)   Z = 14.0   N =  28.1   A =  28.09 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
21
22 Material:    G4_Fe     density:  7.874 g/cm3   RadL:   1.757 cm   Nucl.Int.Length:  16.990 cm   Imean: 286.000 eV
23   --->  Element: Fe (Fe)   Z = 26.0   N =  55.8   A =  55.85 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
24
25 Material:    G4_Ge     density:  5.323 g/cm3   RadL:   2.301 cm   Nucl.Int.Length:  27.431 cm   Imean: 350.000 eV
26   --->  Element: Ge (Ge)   Z = 32.0   N =  72.6   A =  72.61 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
27
28 Material:     G4_W     density: 19.300 g/cm3   RadL:   3.504 mm   Nucl.Int.Length:  10.312 cm   Imean: 727.000 eV
29   --->  Element: W (W)   Z = 74.0   N = 183.8   A = 183.84 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
30
31 Material:    G4_Pb     density: 11.350 g/cm3   RadL:   5.613 mm   Nucl.Int.Length:  18.248 cm   Imean: 823.000 eV
32   --->  Element: Pb (Pb)   Z = 82.0   N = 207.2   A = 207.22 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
33
34 Material:   G4_AIR     density:  1.205 mg/cm3  RadL: 303.921 m    Nucl.Int.Length: 710.095 m    Imean:  85.700 eV   temperature: 273.15 K  pressure:   1.00 atm
35   --->  Element: C (C)   Z =  6.0   N =  12.0   A =  12.01 g/mole  ElmMassFraction:   0.01 %  ElmAbundance   0.02 %
36   --->  Element: N (N)   Z =  7.0   N =  14.0   A =  14.01 g/mole  ElmMassFraction:  75.53 %  ElmAbundance  78.44 %
37   --->  Element: O (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  23.18 %  ElmAbundance  21.07 %
38   --->  Element: Ar (Ar)   Z = 18.0   N =  39.9   A =  39.95 g/mole  ElmMassFraction:   1.28 %  ElmAbundance   0.47 %
39
40 Material: G4_WATER H_2O   density:  1.000 g/cm3   RadL:  36.083 cm   Nucl.Int.Length:  75.423 cm   Imean:  78.000 eV
41   --->  Element: H (H)   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction:  11.19 %  ElmAbundance  66.67 %
42   --->  Element: O (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  88.81 %  ElmAbundance  33.33 %
43
44
45/run/verbose 2
46#
47/testem/det/setMat G4_WATER
48/testem/det/setRadius 3 cm
49/testem/det/setNbOfLayers 1
50#
51/testem/phys/addPhysics local
52PhysicsList::AddPhysicsList: <local>
53#
54/run/initialize
55userDetector->Construct() start.
56
57---------------------------------------------------------
58---> The Absorber is a sphere of 3 cm  radius of G4_WATER divided in 1 slices of 3 cm
59
60---------------------------------------------------------
61Absorber is registered to the default region.
62physicsList->Construct() start.
63physicsList->Construct() start.
64physicsList->CheckParticleList() start.
65physicsList->setCut() start.
66PhysicsList::SetCuts:CutLength : 1 mm
67#
68/testem/phys/setCuts 1 um
69#
70/gun/particle e-
71/gun/energy 4 MeV
72#
73/testem/histo/setFileType root
74/testem/histo/setFileName testem12
75/testem/histo/setHisto 1 120 0. 2.4 cm
76----> SetHisto 1: dE/dr (MeV/mm) along radius (cm);  120 bins from 0 cm to 2.4 cm
77/testem/histo/setHisto 2 180 0. 4.5 MeV
78----> SetHisto 2: total Energy deposited in absorber (MeV);  180 bins from 0 MeV to 4.5 MeV
79/testem/histo/setHisto 3 100 0. 3. cm
80----> SetHisto 3: true track length of the primary particle (cm);  100 bins from 0 cm to 3 cm
81/testem/histo/setHisto 4 100 0. 500. um
82----> SetHisto 4: true step size of the primary particle (um);  100 bins from 0 um to 500 um
83/testem/histo/setHisto 5 200 -3. 3. cm
84----> SetHisto 5: projected range of the primary particle (cm);  200 bins from -3 cm to 3 cm
85/testem/histo/setcsdaRange 2.037 cm
86/testem/histo/setHisto 8 120 0. 1.2 none
87----> SetHisto 8: d(E/E0)/d(r/r0) along r/r0;  120 bins from 0 none to 1.2 none
88#
89/testem/event/printModulo 1000
90#
91/run/beamOn 10000
92
93phot:   for  gamma    SubType= 12
94      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
95       PhotoElectric :     Emin=          0 eV         Emax=   10 TeV
96
97compt:   for  gamma    SubType= 13
98      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
99      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
100       Klein-Nishina :     Emin=          0 eV         Emax=   10 TeV
101
102conv:   for  gamma    SubType= 14
103      Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1
104      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
105       Bethe-Heitler :     Emin=          0 eV         Emax=   10 TeV
106
107msc:   for e-    SubType= 10
108      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
109      RangeFactor= 0.04, stepLimitType: 1, latDisplacement: 1
110      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
111          UrbanMsc93 :     Emin=          0 eV         Emax=   10 TeV
112
113eIoni:   for  e-    SubType= 2
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      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
117      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
118        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
119      CSDA range table up to 1 GeV in 160 bins
120
121eBrem:   for  e-    SubType= 3
122      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
123      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
124      LPM flag: 1 for E > 1 GeV
125      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
126               eBrem :     Emin=          0 eV         Emax=   1 GeV
127            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
128
129eIoni:   for  e+    SubType= 2
130      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
131      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
132      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
133      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
134        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
135      CSDA range table up to 1 GeV in 160 bins
136
137eBrem:   for  e+    SubType= 3
138      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
139      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
140      LPM flag: 1 for E > 1 GeV
141      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
142               eBrem :     Emin=          0 eV         Emax=   1 GeV
143            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
144
145annihil:   for  e+    SubType= 5
146      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
147      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
148            eplus2gg :     Emin=          0 eV         Emax=   10 TeV
149
150msc:   for proton    SubType= 10
151      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
152      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
153      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
154          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
155
156hIoni:   for  proton    SubType= 2
157      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
158      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
159      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
160      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
161               Bragg :     Emin=          0 eV         Emax=   2 MeV
162          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
163      CSDA range table up to 1 GeV in 160 bins
164
165hBrems:   for  proton    SubType= 3
166      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
167      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
168      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
169               hBrem :     Emin=          0 eV         Emax=   10 TeV
170
171hPairProd:   for  proton    SubType= 4
172      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
173      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
174      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
175           hPairProd :     Emin=          0 eV         Emax=   10 TeV
176
177msc:   for GenericIon    SubType= 10
178      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
179      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
180          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
181
182ionIoni:   for  GenericIon    SubType= 2
183      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
184      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
185      finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
186      Stopping Power data for 17 ion/material pairs, nuclearStopping: 1
187      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
188         ParamICRU73 :     Emin=          0 eV         Emax=   10 TeV
189      CSDA range table up to 1 GeV in 160 bins
190
191nuclearStopping:   for  GenericIon    SubType= 1
192      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
193   ICRU49NucStopping :     Emin=          0 eV         Emax=   10 TeV
194
195nuclearStopping:   for  alpha    SubType= 1
196      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
197   ICRU49NucStopping :     Emin=          0 eV         Emax=   10 TeV
198
199hIoni:   for  anti_proton    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            ICRU73QO :     Emin=          0 eV         Emax=   2 MeV
205          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
206      CSDA range table up to 1 GeV in 160 bins
207
208hIoni:   for  kaon+    SubType= 2
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      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
212      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
213               Bragg :     Emin=          0 eV         Emax=   1.05231 MeV
214          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
215      CSDA range table up to 1 GeV in 160 bins
216
217hIoni:   for  kaon-    SubType= 2
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      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
221      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
222            ICRU73QO :     Emin=          0 eV         Emax=   1.05231 MeV
223          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
224      CSDA range table up to 1 GeV in 160 bins
225
226muMsc:   for mu+    SubType= 10
227      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
228      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0
229      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
230          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
231
232muIoni:   for  mu+    SubType= 2
233      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
234      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
235      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
236      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
237               Bragg :     Emin=          0 eV         Emax=   200 keV
238          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
239        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
240      CSDA range table up to 1 GeV in 160 bins
241
242muBrems:   for  mu+    SubType= 3
243      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
244      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
245      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
246              MuBrem :     Emin=          0 eV         Emax=   10 TeV
247
248muPairProd:   for  mu+    SubType= 4
249      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
250      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
251      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
252          muPairProd :     Emin=          0 eV         Emax=   10 TeV
253
254muIoni:   for  mu-    SubType= 2
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      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
258      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
259            ICRU73QO :     Emin=          0 eV         Emax=   200 keV
260          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
261        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
262      CSDA range table up to 1 GeV in 160 bins
263
264muBrems:   for  mu-    SubType= 3
265      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
266      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
267      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
268              MuBrem :     Emin=          0 eV         Emax=   10 TeV
269
270muPairProd:   for  mu-    SubType= 4
271      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
272      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
273      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
274          muPairProd :     Emin=          0 eV         Emax=   10 TeV
275
276hIoni:   for  pi+    SubType= 2
277      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
278      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
279      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
280      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
281               Bragg :     Emin=          0 eV         Emax=   297.504 keV
282          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
283      CSDA range table up to 1 GeV in 160 bins
284
285hBrems:   for  pi+    SubType= 3
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               hBrem :     Emin=          0 eV         Emax=   10 TeV
290
291hPairProd:   for  pi+    SubType= 4
292      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
293      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
294      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
295           hPairProd :     Emin=          0 eV         Emax=   10 TeV
296
297msc:   for pi-    SubType= 10
298      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
299      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
300      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
301          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
302
303hIoni:   for  pi-    SubType= 2
304      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
305      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
306      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
307      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
308            ICRU73QO :     Emin=          0 eV         Emax=   297.504 keV
309          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
310      CSDA range table up to 1 GeV in 160 bins
311
312hBrems:   for  pi-    SubType= 3
313      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
314      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
315      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
316               hBrem :     Emin=          0 eV         Emax=   10 TeV
317
318hPairProd:   for  pi-    SubType= 4
319      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
320      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
321      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
322           hPairProd :     Emin=          0 eV         Emax=   10 TeV
323
324Region <DefaultRegionForTheWorld> -- appears in <Absorber> world volume
325 Root logical volume(s) : Absorber
326 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
327 Materials : G4_WATER
328 Production cuts :   gamma 1 um      e- 1 um      e+ 1 um  proton 0 fm
329
330========= Table of registered couples ==============================
331
332Index : 0     used in the geometry : Yes     recalculation needed : No
333 Material : G4_WATER
334 Range cuts        :  gamma  1 um     e-  1 um     e+  1 um  proton 0 fm
335 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
336 Region(s) which use this couple :
337    DefaultRegionForTheWorld
338
339====================================================================
340
341Start closing geometry.
342G4GeometryManager::ReportVoxelStats -- Voxel Statistics
343
344    Total memory consumed for geometry optimisation:   0 kByte
345    Total CPU time elapsed for geometry optimisation: 0 seconds
346### Run 0 start.
347
348--------- Ranecu engine status ---------
349 Initial seed (index) = 0
350 Current couple of seeds = 9876, 54321
351----------------------------------------
352
353---> stepMax from HistoManager = 2e+02 um
354Start Run processing.
355
356---> Begin of Event: 0
357
358---> Begin of Event: 1000
359
360---> Begin of Event: 2000
361
362---> Begin of Event: 3000
363
364---> Begin of Event: 4000
365
366---> Begin of Event: 5000
367
368---> Begin of Event: 6000
369
370---> Begin of Event: 7000
371
372---> Begin of Event: 8000
373
374---> Begin of Event: 9000
375Run terminated.
376Run Summary
377  Number of events processed : 10000
378  User=42.93s Real=71.08s Sys=0.28s
379
380 ======================== run summary ======================
381
382 The run consists of 10000 e- of 4.00 MeV through 3.00 cm  of G4_WATER (density: 1.00 g/cm3 )
383
384 ============================================================
385
386 Total Energy deposited        = 3.932 MeV +- 234.119 keV
387
388 Track length of primary track = 2.045 cm  +- 2.840 mm
389 Range from EmCalculator       = 2.048 cm  (from full dE/dx)
390
391 Projected range               = 1.316 cm  +- 4.268 mm
392
393 Nb of steps of primary track  = 241.47 +- 31.53         Step size= 84.668 um  +- 4.553 um
394
395--------- Ranecu engine status ---------
396 Initial seed (index) = 0
397 Current couple of seeds = 343308114, 1826359441
398----------------------------------------
399UserDetectorConstruction deleted.
400UserPhysicsList deleted.
401UserRunAction deleted.
402UserPrimaryGenerator deleted.
403G4 kernel has come to Quit state.
404EventManager deleted.
405UImanager deleted.
406Units table cleared.
407StateManager deleted.
408RunManagerKernel is deleted.
409RunManager is deleting.
410
Note: See TracBrowser for help on using the repository browser.