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

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

update ti head

File size: 19.3 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
14***** Table : Nb of materials = 8 *****
15
16 Material:    G4_Al     density:  2.699 g/cm3   RadL:   8.896 cm   Nucl.Int.Length:  38.877 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.753 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.959 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.334 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.306 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.261 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.261 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.517 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        BetheHeitler :   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  AngularGenUrban
127            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
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  AngularGenUrban
143            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
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
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
208msc:   for kaon+    SubType= 10
209      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
210      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
211      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
212          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
213
214hIoni:   for  kaon+    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               Bragg :   Emin=        0 eV       Emax=    1.05231 MeV
220          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
221      CSDA range table up to 1 GeV in 160 bins
222
223hIoni:   for  kaon-    SubType= 2
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      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
227      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
228            ICRU73QO :   Emin=        0 eV       Emax=    1.05231 MeV
229          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
230      CSDA range table up to 1 GeV in 160 bins
231
232muMsc:   for mu+    SubType= 10
233      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
234      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0
235      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
236          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
237
238muIoni:   for  mu+    SubType= 2
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      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
242      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
243               Bragg :   Emin=        0 eV       Emax=        200 keV
244          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
245        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
246      CSDA range table up to 1 GeV in 160 bins
247
248muBrems:   for  mu+    SubType= 3
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              MuBrem :   Emin=        0 eV       Emax=         10 TeV
253
254muPairProd:   for  mu+    SubType= 4
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
258          muPairProd :   Emin=        0 eV       Emax=         10 TeV
259
260muIoni:   for  mu-    SubType= 2
261      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
262      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
263      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
264      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
265            ICRU73QO :   Emin=        0 eV       Emax=        200 keV
266          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
267        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
268      CSDA range table up to 1 GeV in 160 bins
269
270muBrems:   for  mu-    SubType= 3
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              MuBrem :   Emin=        0 eV       Emax=         10 TeV
275
276muPairProd:   for  mu-    SubType= 4
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
280          muPairProd :   Emin=        0 eV       Emax=         10 TeV
281
282hIoni:   for  pi+    SubType= 2
283      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
284      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
285      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
286      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
287               Bragg :   Emin=        0 eV       Emax=    297.505 keV
288          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
289      CSDA range table up to 1 GeV in 160 bins
290
291hBrems:   for  pi+    SubType= 3
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               hBrem :   Emin=        0 eV       Emax=         10 TeV
296
297hPairProd:   for  pi+    SubType= 4
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
301           hPairProd :   Emin=        0 eV       Emax=         10 TeV
302
303msc:   for pi-    SubType= 10
304      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
305      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
306      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
307          UrbanMsc90 :   Emin=        0 eV       Emax=         10 TeV
308
309hIoni:   for  pi-    SubType= 2
310      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
311      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
312      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
313      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
314            ICRU73QO :   Emin=        0 eV       Emax=    297.505 keV
315          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
316      CSDA range table up to 1 GeV in 160 bins
317
318hBrems:   for  pi-    SubType= 3
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               hBrem :   Emin=        0 eV       Emax=         10 TeV
323
324hPairProd:   for  pi-    SubType= 4
325      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
326      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
327      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
328           hPairProd :   Emin=        0 eV       Emax=         10 TeV
329
330Region <DefaultRegionForTheWorld> -- appears in <Absorber> world volume
331 Root logical volume(s) : Absorber
332 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
333 Materials : G4_WATER
334 Production cuts :   gamma 1 um      e- 1 um      e+ 1 um  proton 0 fm
335
336========= Table of registered couples ==============================
337
338Index : 0     used in the geometry : Yes     recalculation needed : No
339 Material : G4_WATER
340 Range cuts        :  gamma  1 um     e-  1 um     e+  1 um  proton 0 fm
341 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 0 eV
342 Region(s) which use this couple :
343    DefaultRegionForTheWorld
344
345====================================================================
346
347Start closing geometry.
348G4GeometryManager::ReportVoxelStats -- Voxel Statistics
349
350    Total memory consumed for geometry optimisation:   0 kByte
351    Total CPU time elapsed for geometry optimisation: 0 seconds
352### Run 0 start.
353
354--------- Ranecu engine status ---------
355 Initial seed (index) = 0
356 Current couple of seeds = 9876, 54321
357----------------------------------------
358
359---> stepMax from HistoManager = 200 um
360Start Run processing.
361
362---> Begin of Event: 0
363
364---> Begin of Event: 1000
365
366---> Begin of Event: 2000
367
368---> Begin of Event: 3000
369
370---> Begin of Event: 4000
371
372---> Begin of Event: 5000
373
374---> Begin of Event: 6000
375
376---> Begin of Event: 7000
377
378---> Begin of Event: 8000
379
380---> Begin of Event: 9000
381Run terminated.
382Run Summary
383  Number of events processed : 10000
384  User=39.49s Real=39.73s Sys=0.22s
385
386 ======================== run summary ======================
387
388 The run consists of 10000 e- of 4.00 MeV through 3.00 cm  of G4_WATER (density: 1.00 g/cm3 )
389
390 ============================================================
391
392 Total Energy deposited        = 3.937 MeV +- 225.962 keV
393
394 Track length of primary track = 2.038 cm  +- 2.841 mm
395 Range from EmCalculator       = 2.048 cm  (from full dE/dx)
396
397 Projected range               = 1.306 cm  +- 4.241 mm
398
399 Nb of steps of primary track  = 240.82 +- 31.65         Step size= 84.619 um  +- 4.544 um
400
401--------- Ranecu engine status ---------
402 Initial seed (index) = 0
403 Current couple of seeds = 1150369497, 2116195687
404----------------------------------------
405UserDetectorConstruction deleted.
406UserPhysicsList deleted.
407UserRunAction deleted.
408UserPrimaryGenerator deleted.
409G4 kernel has come to Quit state.
410EventManager deleted.
411UImanager deleted.
412Units table cleared.
413StateManager deleted.
414RunManagerKernel is deleted.
415RunManager is deleting.
416
Note: See TracBrowser for help on using the repository browser.