source: trunk/examples/extended/electromagnetic/TestEm3/TestEm3.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.5 KB
RevLine 
[807]1
[1337]2        ############################################
3        !!! WARNING - FPE detection is activated !!!
4        ############################################
5
[807]6*************************************************************
[1337]7 Geant4 version Name: geant4-09-03-ref-06    (25-June-2010)
[807]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#
[1337]15/testem/phys/addPhysics emstandard_opt0
[807]16#
17/run/initialize
18userDetector->Construct() start.
19
20-------------------------------------------------------------
21 ---> The calorimeter is 50 layers of:
22                Lead:    2.3 mm
23         liquidArgon:    5.7 mm
24-------------------------------------------------------------
25
[1337]26 Material: Galactic     density:  0.000 kg/m3   RadL: 204322111.300 pc   Nucl.Int.Length: 113728944.943 pc   Imean:  19.200 eV   temperature:   2.73 K  pressure:   0.00 atm
[807]27   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
28
[1337]29 Material:     Lead     density: 11.350 g/cm3   RadL:   5.613 mm   Nucl.Int.Length:  18.247 cm   Imean: 823.000 eV
[807]30   --->  Element: Pb (Pb)   Z = 82.0   N = 207.2   A = 207.22 g/mole
[1337]31   --->  Isotope: Pb204   Z = 82   N = 204   A = 203.97 g/mole   abundance:   1.40 %
32   --->  Isotope: Pb206   Z = 82   N = 206   A = 205.97 g/mole   abundance:  24.10 %
33   --->  Isotope: Pb207   Z = 82   N = 207   A = 206.98 g/mole   abundance:  22.10 %
34   --->  Isotope: Pb208   Z = 82   N = 208   A = 207.98 g/mole   abundance:  52.40 %  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
[807]35
[1337]36 Material: liquidArgon     density:  1.390 g/cm3   RadL:  14.064 cm   Nucl.Int.Length:  86.022 cm   Imean: 188.000 eV
[807]37   --->  Element: Ar (Ar)   Z = 18.0   N =  40.0   A =  39.95 g/mole
[1337]38   --->  Isotope:  Ar36   Z = 18   N =  36   A =  35.97 g/mole   abundance:   0.34 %
39   --->  Isotope:  Ar38   Z = 18   N =  38   A =  37.96 g/mole   abundance:   0.06 %
40   --->  Isotope:  Ar40   Z = 18   N =  40   A =  39.96 g/mole   abundance:  99.60 %  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
[807]41
42-------------------------------------------------------------
43World is registered to the default region.
44physicsList->Construct() start.
45physicsList->Construct() start.
[1337]46physicsList->CheckParticleList() start.
[807]47physicsList->setCut() start.
48PhysicsList::SetCuts:CutLength : 1 mm
49/process/eLoss/verbose 1
50/process/msc/Skin 1
51/run/physicsModified
52#
53/gun/particle e-
54/gun/energy 1 GeV
55#
56/testem/histo/setFileName testem3
[1337]57/testem/histo/setFileType root
58/testem/histo/setHisto 1 100 700. 1000. MeV
[807]59----> SetHisto 1: Edep in absorber 1 (MeV);  100 bins from 700 MeV to 1000 MeV
[1337]60/testem/histo/setHisto 2 100 0. 300. MeV
[807]61----> SetHisto 2: Edep in absorber 2 (MeV);  100 bins from 0 MeV to 300 MeV
[1337]62/testem/histo/setHisto 11 52 0. 52. none
[807]63----> SetHisto 11: longit. profile of Edep (MeV/event) in absorber 1;  52 bins from 0 none to 52 none
[1337]64/testem/histo/setHisto 12 52 0. 52. none
[807]65----> SetHisto 12: longit. profile of Edep (MeV/event) in absorber 2;  52 bins from 0 none to 52 none
[1337]66/testem/histo/setHisto 21 102 0. 102. none
[807]67----> SetHisto 21: energy flow (MeV/event);  102 bins from 0 none to 102 none
[1337]68/testem/histo/setHisto 22 102 0. 102. none
[807]69----> SetHisto 22: lateral energy leak (MeV/event);  102 bins from 0 none to 102 none
70#
71/testem/event/printModulo 50
72#
73/run/beamOn 100
74
[1337]75phot:   for  gamma    SubType= 12
76      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
77       PhotoElectric :     Emin=          0 eV         Emax=   10 TeV
[807]78
[1337]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
[807]83
[1337]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
[807]88
[1337]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
[807]94
[1337]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
[807]101
[1337]102eBrem:   for  e-    SubType= 3
103      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
104      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
105      LPM flag: 1 for E > 1 GeV
106      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
107               eBrem :     Emin=          0 eV         Emax=   1 GeV
108            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
[807]109
[1337]110eIoni:   for  e+    SubType= 2
111      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
112      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
113      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
114      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
115        MollerBhabha :     Emin=          0 eV         Emax=   10 TeV
[807]116
[1337]117eBrem:   for  e+    SubType= 3
118      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
119      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
120      LPM flag: 1 for E > 1 GeV
121      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
122               eBrem :     Emin=          0 eV         Emax=   1 GeV
123            eBremRel :     Emin=          1 GeV        Emax=   10 TeV
[807]124
[1337]125annihil:   for  e+    SubType= 5
126      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
127      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
128            eplus2gg :     Emin=          0 eV         Emax=   10 TeV
[807]129
[1337]130msc:   for proton    SubType= 10
131      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
132      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
133      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
134          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
[807]135
[1337]136hIoni:   for  proton    SubType= 2
137      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
138      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
139      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
140      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
141               Bragg :     Emin=          0 eV         Emax=   2 MeV
142          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
[807]143
[1337]144hBrems:   for  proton    SubType= 3
145      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
146      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
147      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
148               hBrem :     Emin=          0 eV         Emax=   10 TeV
[807]149
[1337]150hPairProd:   for  proton    SubType= 4
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           hPairProd :     Emin=          0 eV         Emax=   10 TeV
[807]155
[1337]156msc:   for GenericIon    SubType= 10
157      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
158      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
159          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
[807]160
[1337]161ionIoni:   for  GenericIon    SubType= 2
162      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
163      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
164      finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
165      Stopping Power data for 17 ion/material pairs, nuclearStopping: 1
166      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
167            BraggIon :     Emin=          0 eV         Emax=   2 MeV
168          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
[807]169
[1337]170hIoni:   for  anti_proton    SubType= 2
171      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
172      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
173      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
174      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
175            ICRU73QO :     Emin=          0 eV         Emax=   2 MeV
176          BetheBloch :     Emin=          2 MeV        Emax=   10 TeV
[807]177
[1337]178hIoni:   for  kaon+    SubType= 2
179      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
180      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
181      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
182      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
183               Bragg :     Emin=          0 eV         Emax=   1.05231 MeV
184          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
[807]185
[1337]186hBrems:   for  kaon+    SubType= 3
187      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
188      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
189      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
190               hBrem :     Emin=          0 eV         Emax=   10 TeV
[807]191
[1337]192hPairProd:   for  kaon+    SubType= 4
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           hPairProd :     Emin=          0 eV         Emax=   10 TeV
[807]197
[1337]198hIoni:   for  kaon-    SubType= 2
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      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
202      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
203            ICRU73QO :     Emin=          0 eV         Emax=   1.05231 MeV
204          BetheBloch :     Emin=    1.05231 MeV        Emax=   10 TeV
[807]205
[1337]206hBrems:   for  kaon-    SubType= 3
207      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
208      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
209      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
210               hBrem :     Emin=          0 eV         Emax=   10 TeV
[807]211
[1337]212hPairProd:   for  kaon-    SubType= 4
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           hPairProd :     Emin=          0 eV         Emax=   10 TeV
[807]217
[1337]218muMsc:   for mu+    SubType= 10
219      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
220      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0
221      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
222          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
[807]223
[1337]224muIoni:   for  mu+    SubType= 2
225      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
226      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
227      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
228      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
229               Bragg :     Emin=          0 eV         Emax=   200 keV
230          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
231        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
[807]232
[1337]233muBrems:   for  mu+    SubType= 3
234      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
235      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
236      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
237              MuBrem :     Emin=          0 eV         Emax=   10 TeV
[807]238
[1337]239muPairProd:   for  mu+    SubType= 4
240      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
241      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
242      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
243          muPairProd :     Emin=          0 eV         Emax=   10 TeV
244
245muIoni:   for  mu-    SubType= 2
246      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
247      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
248      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
249      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
250            ICRU73QO :     Emin=          0 eV         Emax=   200 keV
251          BetheBloch :     Emin=        200 keV        Emax=   1 GeV
252        MuBetheBloch :     Emin=          1 GeV        Emax=   10 TeV
253
254muBrems:   for  mu-    SubType= 3
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              MuBrem :     Emin=          0 eV         Emax=   10 TeV
259
260muPairProd:   for  mu-    SubType= 4
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      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
264          muPairProd :     Emin=          0 eV         Emax=   10 TeV
265
266hIoni:   for  pi+    SubType= 2
267      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
268      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
269      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
270      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
271               Bragg :     Emin=          0 eV         Emax=   297.504 keV
272          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
273
274hBrems:   for  pi+    SubType= 3
275      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
276      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
277      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
278               hBrem :     Emin=          0 eV         Emax=   10 TeV
279
280hPairProd:   for  pi+    SubType= 4
281      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
282      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
283      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
284           hPairProd :     Emin=          0 eV         Emax=   10 TeV
285
286msc:   for pi-    SubType= 10
287      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
288      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
289      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
290          UrbanMsc90 :     Emin=          0 eV         Emax=   10 TeV
291
292hIoni:   for  pi-    SubType= 2
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      finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
296      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
297            ICRU73QO :     Emin=          0 eV         Emax=   297.504 keV
298          BetheBloch :     Emin=    297.504 keV        Emax=   10 TeV
299
300hBrems:   for  pi-    SubType= 3
301      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
302      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
303      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
304               hBrem :     Emin=          0 eV         Emax=   10 TeV
305
306hPairProd:   for  pi-    SubType= 4
307      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
308      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
309      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
310           hPairProd :     Emin=          0 eV         Emax=   10 TeV
311
[807]312Region <DefaultRegionForTheWorld> -- appears in <World> world volume
313 Root logical volume(s) : World
314 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
315 Materials : Galactic Lead liquidArgon
[1337]316 Production cuts :   gamma 1 mm      e- 1 mm      e+ 1 mm  proton 1 mm
[807]317
318========= Table of registered couples ==============================
319
320Index : 0     used in the geometry : Yes     recalculation needed : No
321 Material : Galactic
[1337]322 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
323 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 100 keV
[807]324 Region(s) which use this couple :
325    DefaultRegionForTheWorld
326
327Index : 1     used in the geometry : Yes     recalculation needed : No
328 Material : Lead
[1337]329 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
330 Energy thresholds :  gamma  101.843 keV    e-  1.36749 MeV    e+  1.27862 MeV proton 100 keV
[807]331 Region(s) which use this couple :
332    DefaultRegionForTheWorld
333
334Index : 2     used in the geometry : Yes     recalculation needed : No
335 Material : liquidArgon
[1337]336 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
337 Energy thresholds :  gamma  6.19986 keV    e-  344.855 keV    e+  333.46 keV proton 100 keV
[807]338 Region(s) which use this couple :
339    DefaultRegionForTheWorld
340
341====================================================================
342
343Start closing geometry.
344G4GeometryManager::ReportVoxelStats -- Voxel Statistics
345
346    Total memory consumed for geometry optimisation:   3 kByte
347    Total CPU time elapsed for geometry optimisation: 0 seconds
348
349    Voxelisation: top CPU users:
350    Percent   Total CPU    System CPU       Memory  Volume
351    -------   ----------   ----------     --------  ----------
352       0.00         0.00         0.00            3k Calorimeter
353       0.00         0.00         0.00            0k Layer
354
355    Voxelisation: top memory users:
356    Percent     Memory      Heads    Nodes   Pointers    Total CPU    Volume
357    -------   --------     ------   ------   --------   ----------    ----------
358      92.05          2k         1       50         50         0.00    Calorimeter
359       7.95          0k         1        3          4         0.00    Layer
360### Run 0 start.
361
362--------- Ranecu engine status ---------
363 Initial seed (index) = 0
364 Current couple of seeds = 9876, 54321
365----------------------------------------
366Start Run processing.
367
368---> Begin Of Event: 0
369
370---> Begin Of Event: 50
371Run terminated.
372Run Summary
373  Number of events processed : 100
[1337]374  User=5.56s Real=9.12s Sys=0s
[807]375
376------------------------------------------------------------
[1337]377      material   Edep       RMS         sqrt(E0(GeV))*rmsE/Emean     total tracklen
[807]378 
[1337]379          Lead: 781.73 MeV :  14.91 MeV     1.907 +- 0.1907 %      54.5 cm  +- 1.17 cm
380   liquidArgon: 207.32 MeV :  15.05 MeV     7.262 +- 0.7262 %      1.02 m   +- 7.76 cm
[807]381
382------------------------------------------------------------
[1337]383 Beam particle e-  E = 1 GeV
384 Mean number of gamma  570
385 Mean number of e-     960
386 Mean number of e+     54.5
387------------------------------------------------------------
[807]388
389 Energy deposition from Energy flow balance :
390  material       Total Edep
391 
[1337]392      Lead:      781.726 MeV
393liquidArgon:     207.311 MeV
[807]394
395------------------------------------------------------------
396
397
398--------- Ranecu engine status ---------
399 Initial seed (index) = 0
[1337]400 Current couple of seeds = 1438848957, 1887175185
[807]401----------------------------------------
402UserDetectorConstruction deleted.
403UserPhysicsList deleted.
404UserRunAction deleted.
405UserPrimaryGenerator deleted.
406G4 kernel has come to Quit state.
407EventManager deleted.
408UImanager deleted.
409Units table cleared.
410StateManager deleted.
411RunManagerKernel is deleted.
412RunManager is deleting.
[1337]413
Note: See TracBrowser for help on using the repository browser.