source: trunk/examples/extended/electromagnetic/TestEm8/TestEm8.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: 19.3 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/run/verbose 2
14#/tracking/verbose 2
15#
16#  select EM PhysicsList
17#
18/testem/phys/addPhysics pai
19#/testem/phys/addPhysics  pai_photon
20#
21/testem/phys/setCuts 1 mm
22#/testem/stepMax         10 mm
23#
24#  PLOT: energy deposit distribution in absorber
25#
26/plots/setnbinEn 20
27 Nb of bins in Edep plot = 20
28/plots/setEnlow 0.0 MeV
29 Elow  in the  Edep plot = 0
30/plots/setEnhigh 0.060 MeV
31 Ehigh in the  Edep plot = 0.06
32/plots/sethistName g4.p11
33 hist file = g4.p11
34#
35/run/initialize
36userDetector->Construct() start.
37
38 The  WORLD   is made of 23.1022mm of Mylar, the transverse size (R) of the world is 100 mm.
39 The ABSORBER is made of 23mm of XeCH4C3H8, the transverse size (R) is 100 mm.
40 Z position of the (middle of the) absorber 0  mm.
41
42World is registered to the default region.
43physicsList->Construct() start.
44physicsList->Construct() start.
45physicsList->CheckParticleList() start.
46physicsList->setCut() start.
47PhysicsList::SetCuts:CutLength : 1 mm
48#
49/run/setCut 0.5 mm
50/run/setCutForRegion GasDetector 1.8 mm
51/run/particle/dumpCutValues
52#
53#/gun/particle e-
54#
55/gun/particle proton
56#
57#/gun/particle mu-
58#
59/gun/energy 200000.00 MeV
60#
61#/gun/energy 200.0 MeV
62#
63#/tracking/verbose 1
64#
65/event/printModulo 100
66/run/beamOn 1000
67G4PAIModel::Initialise for e+
68G4PAIModel::Initialise for e-
69G4PAIModel::Initialise for mu+
70G4PAIModel::Initialise for mu-
71G4PAIModel::Initialise for pi+
72G4PAIModel::Initialise for pi-
73G4PAIModel::Initialise for proton
74
75phot:   for  gamma    SubType= 12
76      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
77       PhotoElectric :     Emin=          0 eV         Emax=   10 TeV
78
79compt:   for  gamma    SubType= 13
80      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
81      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
82       Klein-Nishina :     Emin=          0 eV         Emax=   10 TeV
83
84conv:   for  gamma    SubType= 14
85      Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1
86      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
87       Bethe-Heitler :     Emin=          0 eV         Emax=   10 TeV
88
89msc:   for e-    SubType= 10
90      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
91      RangeFactor= 0.04, stepLimitType: 1, latDisplacement: 1
92      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
93          UrbanMsc93 :     Emin=          0 eV         Emax=   10 TeV
94
95eIoni:   for  e-    SubType= 2
96      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
97      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
98      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
99      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
100        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
101      ===== EM models for the G4Region  GasDetector ======
102            PAIModel :     Emin=          0 eV         Emax=   10 TeV
103
104eBrem:   for  e-    SubType= 3
105      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
106      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
107      LPM flag: 1 for E > 1 GeV
108      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
109               eBrem :     Emin=          0 eV         Emax=   1 GeV
110            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
111
112eIoni:   for  e+    SubType= 2
113      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
114      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
115      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
116      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
117        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
118      ===== EM models for the G4Region  GasDetector ======
119            PAIModel :     Emin=          0 eV         Emax=   10 TeV
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
129annihil:   for  e+    SubType= 5
130      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
131      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
132            eplus2gg :     Emin=          0 eV         Emax=   10 TeV
133
134msc:   for proton    SubType= 10
135      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
136      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
137      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
138          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
139
140hIoni:   for  proton    SubType= 2
141      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
142      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
143      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
144      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
145               Bragg :     Emin=          0 eV         Emax=   2 MeV
146          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
147      ===== EM models for the G4Region  GasDetector ======
148            PAIModel :     Emin=          0 eV         Emax=   10 TeV
149
150hBrems:   for  proton    SubType= 3
151      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
152      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
153      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
154               hBrem :     Emin=          0 eV         Emax=   10 TeV
155
156hPairProd:   for  proton    SubType= 4
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
160           hPairProd :     Emin=          0 eV         Emax=   10 TeV
161
162msc:   for GenericIon    SubType= 10
163      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
164      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
165          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
166
167ionIoni:   for  GenericIon    SubType= 2
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      finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
171      Stopping Power data for 17 ion/material pairs, nuclearStopping: 1
172      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
173            BraggIon :     Emin=          0 eV         Emax=   2 MeV
174          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
175
176hIoni:   for  anti_proton    SubType= 2
177      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
178      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
179      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
180      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
181            ICRU73QO :     Emin=          0 eV         Emax=   2 MeV
182          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
183
184hIoni:   for  kaon+    SubType= 2
185      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
186      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
187      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
188      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
189               Bragg :     Emin=          0 eV         Emax=   1.05231 MeV
190          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
191
192hBrems:   for  kaon+    SubType= 3
193      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
194      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
195      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
196               hBrem :     Emin=          0 eV         Emax=   10 TeV
197
198hPairProd:   for  kaon+    SubType= 4
199      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
200      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
201      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
202           hPairProd :     Emin=          0 eV         Emax=   10 TeV
203
204hIoni:   for  kaon-    SubType= 2
205      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
206      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
207      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
208      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
209            ICRU73QO :     Emin=          0 eV         Emax=   1.05231 MeV
210          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
211
212hBrems:   for  kaon-    SubType= 3
213      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
214      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
215      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
216               hBrem :     Emin=          0 eV         Emax=   10 TeV
217
218hPairProd:   for  kaon-    SubType= 4
219      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
220      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
221      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
222           hPairProd :     Emin=          0 eV         Emax=   10 TeV
223
224muMsc:   for mu+    SubType= 10
225      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
226      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0
227      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
228          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
229
230muIoni:   for  mu+    SubType= 2
231      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
232      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
233      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
234      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
235               Bragg :     Emin=          0 eV         Emax=   200 keV
236          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
237        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
238      ===== EM models for the G4Region  GasDetector ======
239            PAIModel :     Emin=          0 eV         Emax=   10 TeV
240
241muBrems:   for  mu+    SubType= 3
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
245              MuBrem :     Emin=          0 eV         Emax=   10 TeV
246
247muPairProd:   for  mu+    SubType= 4
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
251          muPairProd :     Emin=          0 eV         Emax=   10 TeV
252
253muIoni:   for  mu-    SubType= 2
254      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
255      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
256      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
257      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
258            ICRU73QO :     Emin=          0 eV         Emax=   200 keV
259          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
260        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
261      ===== EM models for the G4Region  GasDetector ======
262            PAIModel :     Emin=          0 eV         Emax=   10 TeV
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      ===== EM models for the G4Region  GasDetector ======
284            PAIModel :     Emin=          0 eV         Emax=   10 TeV
285
286hBrems:   for  pi+    SubType= 3
287      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
288      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
289      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
290               hBrem :     Emin=          0 eV         Emax=   10 TeV
291
292hPairProd:   for  pi+    SubType= 4
293      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
294      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
295      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
296           hPairProd :     Emin=          0 eV         Emax=   10 TeV
297
298msc:   for pi-    SubType= 10
299      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
300      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
301      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
302          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
303
304hIoni:   for  pi-    SubType= 2
305      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
306      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
307      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
308      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
309            ICRU73QO :     Emin=          0 eV         Emax=   297.504 keV
310          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
311      ===== EM models for the G4Region  GasDetector ======
312            PAIModel :     Emin=          0 eV         Emax=   10 TeV
313
314hBrems:   for  pi-    SubType= 3
315      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
316      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
317      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
318               hBrem :     Emin=          0 eV         Emax=   10 TeV
319
320hPairProd:   for  pi-    SubType= 4
321      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
322      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
323      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
324           hPairProd :     Emin=          0 eV         Emax=   10 TeV
325
326Region <DefaultRegionForTheWorld> -- appears in <World> world volume
327 Root logical volume(s) : World
328 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
329 Materials : Mylar
330 Production cuts :   gamma 500 um      e- 500 um      e+ 500 um  proton 500 um
331
332Region <GasDetector> -- appears in <World> world volume
333 Root logical volume(s) : Absorber
334 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
335 Materials : XeCH4C3H8
336 Production cuts :   gamma 1.8 mm      e- 1.8 mm      e+ 1.8 mm  proton 1.8 mm
337
338========= Table of registered couples ==============================
339
340Index : 0     used in the geometry : Yes     recalculation needed : No
341 Material : Mylar
342 Range cuts        :  gamma  500 um     e-  500 um     e+  500 um  proton 500 um
343 Energy thresholds :  gamma  2.27786 keV    e-  261.805 keV    e+  258.31 keV proton 50 keV
344 Region(s) which use this couple :
345    DefaultRegionForTheWorld
346
347Index : 1     used in the geometry : Yes     recalculation needed : No
348 Material : XeCH4C3H8
349 Range cuts        :  gamma  1.8 mm     e-  1.8 mm     e+  1.8 mm  proton 1.8 mm
350 Energy thresholds :  gamma  990 eV     e-  5.57533 keV    e+  5.40743 keV proton 180 keV
351 Region(s) which use this couple :
352    GasDetector
353
354====================================================================
355
356Start closing geometry.
357G4GeometryManager::ReportVoxelStats -- Voxel Statistics
358
359    Total memory consumed for geometry optimisation:   0 kByte
360    Total CPU time elapsed for geometry optimisation: 0 seconds
361### Run 0 start.
362
363--------- Ranecu engine status ---------
364 Initial seed (index) = 0
365 Current couple of seeds = 9876, 54321
366----------------------------------------
367Start Run processing.
368
369---> Begin of Event: 0
370
371---> Begin of Event: 100
372
373---> Begin of Event: 200
374
375---> Begin of Event: 300
376
377---> Begin of Event: 400
378
379---> Begin of Event: 500
380
381---> Begin of Event: 600
382
383---> Begin of Event: 700
384
385---> Begin of Event: 800
386
387---> Begin of Event: 900
388Run terminated.
389Run Summary
390  Number of events processed : 1000
391  User=0.11s Real=0.1s Sys=0s
392 ================== run summary =====================
393 end of Run TotNbofEvents = 1000
394    mean charged track length   in absorber=24.2704 +- 0.310817  mm 
395
396            mean energy deposit in absorber=0.0160189 +- 0.000378691  MeV
397
398 mean number of steps in absorber (charged) =0.898 +- 0.0848269     
399 mean number of steps in absorber (neutral) =0.002 +- 0.0014128     
400
401   mean number of charged secondaries = 0.257 +- 0.0163997
402
403   mean number of neutral secondaries = 0.002 +- 0.0014128
404
405   mean number of e-s =0.257  and e+s =0
406
407(number) transmission coeff=1  reflection coeff=0
408
409 energy deposit distribution
410#entries=1000    #underflows=0    #overflows=12
411 bin nb      Elow      entries     normalized
412    0         0           0           0
413    1         3           4       0.004
414    2         6          95       0.095
415    3         9         270        0.27
416    4        12         251       0.251
417    5        15         146       0.146
418    6        18          90        0.09
419    7        21          49       0.049
420    8        24          30        0.03
421    9        27          16       0.016
422   10        30          13       0.013
423   11        33           5       0.005
424   12        36           5       0.005
425   13        39           9       0.009
426   14        42           1       0.001
427   15        45           2       0.002
428   16        48           1       0.001
429   17        51           0           0
430   18        54           0           0
431   19        57           1       0.001
432
433 Emp =           0.009   width=          0.006  MeV
434
435
436--------- Ranecu engine status ---------
437 Initial seed (index) = 0
438 Current couple of seeds = 1481926969, 852281434
439----------------------------------------
440#
441UserDetectorConstruction deleted.
442UserPhysicsList deleted.
443UserRunAction deleted.
444UserPrimaryGenerator deleted.
445G4 kernel has come to Quit state.
446G4SDManager deleted.
447EventManager deleted.
448UImanager deleted.
449Units table cleared.
450StateManager deleted.
451RunManagerKernel is deleted.
452RunManager is deleting.
453
Note: See TracBrowser for help on using the repository browser.