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

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

update ti head

File size: 17.6 KB
RevLine 
[807]1
[1337]2        ############################################
3        !!! WARNING - FPE detection is activated !!!
4        ############################################
5
[807]6*************************************************************
[1342]7 Geant4 version Name: geant4-09-03-ref-09    (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
14***** Table : Nb of materials = 15 *****
15
[1342]16 Material:      Air     density:  1.290 mg/cm3  RadL: 285.161 m    Nucl.Int.Length: 662.904 m    Imean:  85.703 eV   temperature: 273.15 K  pressure:   1.00 atm
[807]17   --->  Element: Nitrogen (N)   Z =  7.0   N =  14.0   A =  14.01 g/mole  ElmMassFraction:  70.00 %  ElmAbundance  72.71 %
18   --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  30.00 %  ElmAbundance  27.29 %
19
[1342]20 Material: H2liquid     density: 70.800 mg/cm3  RadL:   8.923 m    Nucl.Int.Length:   4.993 m    Imean:  19.200 eV
[807]21   --->  Element: Hydrogen (H)   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
22
[1342]23 Material:    Water H_2O   density:  1.000 g/cm3   RadL:  36.092 cm   Nucl.Int.Length:  75.537 cm   Imean:  75.000 eV
[807]24   --->  Element: Hydrogen (H)   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction:  11.21 %  ElmAbundance  66.67 %
25   --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  88.79 %  ElmAbundance  33.33 %
26
[1342]27 Material: liquidArgon     density:  1.390 g/cm3   RadL:  14.065 cm   Nucl.Int.Length:  86.006 cm   Imean: 188.000 eV
[807]28   --->  Element: liquidArgon ( )   Z = 18.0   N =  40.0   A =  39.95 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
29
[1342]30 Material:   Carbon     density:  2.267 g/cm3   RadL:  18.833 cm   Nucl.Int.Length:  35.376 cm   Imean:  81.000 eV
[807]31   --->  Element: Carbon ( )   Z =  6.0   N =  12.0   A =  12.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
32
[1342]33 Material: Aluminium     density:  2.700 g/cm3   RadL:   8.893 cm   Nucl.Int.Length:  38.860 cm   Imean: 166.000 eV
[807]34   --->  Element: Aluminium ( )   Z = 13.0   N =  27.0   A =  26.98 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
35
[1342]36 Material:  Silicon     density:  2.330 g/cm3   RadL:   9.368 cm   Nucl.Int.Length:  45.761 cm   Imean: 173.000 eV
[807]37   --->  Element: Silicon ( )   Z = 14.0   N =  28.1   A =  28.09 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
38
[1342]39 Material: Germanium     density:  5.323 g/cm3   RadL:   2.301 cm   Nucl.Int.Length:  27.333 cm   Imean: 350.000 eV
[807]40   --->  Element: Germanium ( )   Z = 32.0   N =  72.6   A =  72.61 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
41
[1342]42 Material:      NaI     density:  3.670 g/cm3   RadL:   2.586 cm   Nucl.Int.Length:  42.858 cm   Imean: 452.000 eV
[807]43   --->  Element: Sodium (Na)   Z = 11.0   N =  23.0   A =  22.99 g/mole  ElmMassFraction:  15.34 %  ElmAbundance  50.00 %
44   --->  Element: Iodine (I)   Z = 53.0   N = 126.9   A = 126.90 g/mole  ElmMassFraction:  84.66 %  ElmAbundance  50.00 %
45
[1342]46 Material:   Iodine     density:  4.930 g/cm3   RadL:   1.720 cm   Nucl.Int.Length:  35.657 cm   Imean: 491.000 eV
[807]47   --->  Element: Iodine (I)   Z = 53.0   N = 126.9   A = 126.90 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
48
[1342]49 Material:      BGO     density:  7.100 g/cm3   RadL:   1.123 cm   Nucl.Int.Length:  22.788 cm   Imean: 473.785 eV
[807]50   --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  15.41 %  ElmAbundance  63.16 %
51   --->  Element: Germanium (Ge)   Z = 32.0   N =  72.6   A =  72.59 g/mole  ElmMassFraction:  17.48 %  ElmAbundance  15.79 %
52   --->  Element: Bismuth (Bi)   Z = 83.0   N = 209.0   A = 208.98 g/mole  ElmMassFraction:  67.10 %  ElmAbundance  21.05 %
53
[1342]54 Material:     Iron     density:  7.870 g/cm3   RadL:   1.759 cm   Nucl.Int.Length:  16.969 cm   Imean: 286.000 eV
[1337]55   --->  Element: Iron ( )   Z = 26.0   N =  55.8   A =  55.85 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
[807]56
[1342]57 Material: Tungsten     density: 19.300 g/cm3   RadL:   3.504 mm   Nucl.Int.Length:  10.306 cm   Imean: 727.000 eV
[807]58   --->  Element: Tungsten ( )   Z = 74.0   N = 183.8   A = 183.85 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
59
[1342]60 Material:     Lead     density: 11.350 g/cm3   RadL:   5.612 mm   Nucl.Int.Length:  18.258 cm   Imean: 823.000 eV
[807]61   --->  Element: Lead ( )   Z = 82.0   N = 207.2   A = 207.19 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
62
[1342]63 Material:  Uranium     density: 18.950 g/cm3   RadL:   3.166 mm   Nucl.Int.Length:  11.447 cm   Imean: 890.000 eV
[807]64   --->  Element: Uranium ( )   Z = 92.0   N = 238.0   A = 238.03 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
65
66
67/run/verbose 2
68#
69/testem/det/setMat Water
70/testem/det/setSize 1 cm
71#
[1337]72/testem/phys/addPhysics standard
73PhysicsList::AddPhysicsList: <standard>
[807]74#
75/run/initialize
76userDetector->Construct() start.
77
78 The Box is 1 cm  of Water
79Water is registered to the default region.
80physicsList->Construct() start.
81physicsList->Construct() start.
[1337]82physicsList->CheckParticleList() start.
[807]83physicsList->setCut() start.
84#
85/gun/particle gamma
86/gun/energy 100 keV
[1337]87/run/beamOn 1000000
[807]88
[1337]89phot:   for  gamma    SubType= 12
90      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]91       PhotoElectric :   Emin=        0 eV       Emax=         10 TeV
[807]92
[1337]93compt:   for  gamma    SubType= 13
94      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
95      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]96       Klein-Nishina :   Emin=        0 eV       Emax=         10 TeV
[807]97
[1337]98conv:   for  gamma    SubType= 14
99      Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1
100      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]101        BetheHeitler :   Emin=        0 eV       Emax=         10 TeV
[807]102
[1337]103eIoni:   for  e-    SubType= 2
104      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
105      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
106      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
107      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]108        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
[807]109
[1337]110eBrem:   for  e-    SubType= 3
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      LPM flag: 1 for E > 1 GeV
114      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]115               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
116            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
[807]117
[1337]118eIoni:   for  e+    SubType= 2
119      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
120      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
121      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
122      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]123        MollerBhabha :   Emin=        0 eV       Emax=         10 TeV
[807]124
[1337]125eBrem:   for  e+    SubType= 3
126      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
127      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
128      LPM flag: 1 for E > 1 GeV
129      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]130               eBrem :   Emin=        0 eV       Emax=          1 GeV  AngularGenUrban
131            eBremLPM :   Emin=        1 GeV      Emax=         10 TeV  AngularGenUrban
[807]132
[1337]133annihil:   for  e+    SubType= 5
134      Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
135      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]136            eplus2gg :   Emin=        0 eV       Emax=         10 TeV
[807]137
[1337]138hIoni:   for  proton    SubType= 2
139      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
140      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
141      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
142      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]143               Bragg :   Emin=        0 eV       Emax=          2 MeV
144          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
[807]145
[1337]146ionIoni:   for  GenericIon    SubType= 2
147      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
148      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
149      finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 0, fluct: 1, linLossLimit= 0.02
[1342]150      Stopping Power data for 17 ion/material pairs
[1337]151      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]152            BraggIon :   Emin=        0 eV       Emax=          2 MeV
153          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
[807]154
[1337]155hIoni:   for  anti_proton    SubType= 2
156      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
157      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
158      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
159      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]160            ICRU73QO :   Emin=        0 eV       Emax=          2 MeV
161          BetheBloch :   Emin=        2 MeV      Emax=         10 TeV
[807]162
[1337]163hIoni:   for  kaon+    SubType= 2
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      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
167      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]168               Bragg :   Emin=        0 eV       Emax=    1.05231 MeV
169          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
[807]170
[1337]171hIoni:   for  kaon-    SubType= 2
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      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
175      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]176            ICRU73QO :   Emin=        0 eV       Emax=    1.05231 MeV
177          BetheBloch :   Emin=  1.05231 MeV      Emax=         10 TeV
[807]178
[1337]179muIoni:   for  mu+    SubType= 2
180      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
181      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
182      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
183      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]184               Bragg :   Emin=        0 eV       Emax=        200 keV
185          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
186        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
[807]187
[1337]188muBrems:   for  mu+    SubType= 3
189      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
190      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
191      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]192              MuBrem :   Emin=        0 eV       Emax=         10 TeV
[807]193
[1337]194muPairProd:   for  mu+    SubType= 4
195      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
196      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
197      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]198          muPairProd :   Emin=        0 eV       Emax=         10 TeV
[807]199
[1337]200muIoni:   for  mu-    SubType= 2
201      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
202      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
203      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
204      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]205            ICRU73QO :   Emin=        0 eV       Emax=        200 keV
206          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
207        MuBetheBloch :   Emin=        1 GeV      Emax=         10 TeV
[807]208
[1337]209muBrems:   for  mu-    SubType= 3
210      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
211      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
212      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]213              MuBrem :   Emin=        0 eV       Emax=         10 TeV
[807]214
[1337]215muPairProd:   for  mu-    SubType= 4
216      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
217      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
218      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]219          muPairProd :   Emin=        0 eV       Emax=         10 TeV
[807]220
[1337]221hIoni:   for  pi+    SubType= 2
222      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
223      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
224      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
225      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]226               Bragg :   Emin=        0 eV       Emax=    297.505 keV
227          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
[807]228
[1337]229hIoni:   for  pi-    SubType= 2
230      dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
231      Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
232      finalRange(mm)= 1, dRoverRange= 0.2, integral: 0, fluct: 1, linLossLimit= 0.01
233      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
[1342]234            ICRU73QO :   Emin=        0 eV       Emax=    297.505 keV
235          BetheBloch :   Emin=  297.505 keV      Emax=         10 TeV
[1337]236
[807]237Region <DefaultRegionForTheWorld> -- appears in <Water> world volume
238 Root logical volume(s) : Water
239 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
240 Materials : Water
[1337]241 Production cuts :   gamma 1 mm      e- 1 mm      e+ 1 mm  proton 0 fm
[807]242
243========= Table of registered couples ==============================
244
245Index : 0     used in the geometry : Yes     recalculation needed : No
246 Material : Water
[1337]247 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
248 Energy thresholds :  gamma  2.93991 keV    e-  349.408 keV    e+  341.455 keV proton 0 eV
[807]249 Region(s) which use this couple :
250    DefaultRegionForTheWorld
251
252====================================================================
253
254Start closing geometry.
255G4GeometryManager::ReportVoxelStats -- Voxel Statistics
256
257    Total memory consumed for geometry optimisation:   0 kByte
258    Total CPU time elapsed for geometry optimisation: 0 seconds
259### Run 0 start.
260
261--------- Ranecu engine status ---------
262 Initial seed (index) = 0
263 Current couple of seeds = 9876, 54321
264----------------------------------------
265Start Run processing.
266
267---> Begin of Event: 0
268
[1337]269---> Begin of Event: 100000
[807]270
[1337]271---> Begin of Event: 200000
[807]272
[1337]273---> Begin of Event: 300000
[807]274
[1337]275---> Begin of Event: 400000
[807]276
[1337]277---> Begin of Event: 500000
[807]278
[1337]279---> Begin of Event: 600000
[807]280
[1337]281---> Begin of Event: 700000
[807]282
[1337]283---> Begin of Event: 800000
[807]284
[1337]285---> Begin of Event: 900000
[807]286Run terminated.
287Run Summary
[1337]288  Number of events processed : 1000000
[1342]289  User=7.56s Real=7.58s Sys=0.02s
[807]290
[1337]291 The run consists of 1000000 gamma of 100 keV through 1 cm  of Water (density: 1 g/cm3 )
[807]292
[1337]293 Process calls frequency --->   Transportation = 848441 compt = 149125  phot = 2434
[807]294
[1337]295 Nb of incident particles unaltered after 1 cm  of Water : 848441 over 1000000 incident particles.  Ratio = 84.844 %
296 ---> CrossSection per volume:  0.16435 cm^-1   CrossSection per mass: 0.16435 cm2/g
[807]297
298 Verification from G4EmCalculator:
[1337]299        compt= 0.1623 cm2/g     phot= 0.0026853 cm2/g   total= 0.16499 cm2/g
[807]300        Expected ratio of transmitted particles= 84.791 %
301
302--------- Ranecu engine status ---------
303 Initial seed (index) = 0
[1337]304 Current couple of seeds = 1922923057, 1535294815
[807]305----------------------------------------
306#
307/gun/particle e-
308/gun/energy 100 MeV
[1337]309/run/beamOn 1000000
[807]310
311Region <DefaultRegionForTheWorld> -- appears in <Water> world volume
312 Root logical volume(s) : Water
313 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
314 Materials : Water
[1337]315 Production cuts :   gamma 1 mm      e- 1 mm      e+ 1 mm  proton 0 fm
[807]316
317========= Table of registered couples ==============================
318
319Index : 0     used in the geometry : Yes     recalculation needed : No
320 Material : Water
[1337]321 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm
[1342]322 Energy thresholds :  gamma  2.93991 keV    e-  349.408 keV    e+  341.455 keV proton 0 eV
[807]323 Region(s) which use this couple :
324    DefaultRegionForTheWorld
325
326====================================================================
327
328### Run 1 start.
329
330--------- Ranecu engine status ---------
331 Initial seed (index) = 0
[1337]332 Current couple of seeds = 1922923057, 1535294815
[807]333----------------------------------------
334Start Run processing.
335
336---> Begin of Event: 0
337
[1337]338---> Begin of Event: 100000
[807]339
[1337]340---> Begin of Event: 200000
[807]341
[1337]342---> Begin of Event: 300000
[807]343
[1337]344---> Begin of Event: 400000
[807]345
[1337]346---> Begin of Event: 500000
[807]347
[1337]348---> Begin of Event: 600000
[807]349
[1337]350---> Begin of Event: 700000
[807]351
[1337]352---> Begin of Event: 800000
[807]353
[1337]354---> Begin of Event: 900000
[807]355Run terminated.
356Run Summary
[1337]357  Number of events processed : 1000000
[1342]358  User=8.44s Real=8.45s Sys=0.01s
[807]359
[1337]360 The run consists of 1000000 e- of 100 MeV through 1 cm  of Water (density: 1 g/cm3 )
[807]361
[1337]362 Process calls frequency --->   Transportation = 540804 eBrem = 277603  eIoni = 181593
[807]363
[1337]364 Nb of incident particles unaltered after 1 cm  of Water : 540804 over 1000000 incident particles.  Ratio = 54.08 %
365 ---> CrossSection per volume:  0.6147 cm^-1    CrossSection per mass: 0.6147 cm2/g
[807]366
367 Verification from G4EmCalculator:
[1337]368        eBrem= 0.37212 cm2/g    eIoni= 0.24337 cm2/g    total= 0.61549 cm2/g
369        Expected ratio of transmitted particles= 54.037 %
[807]370
371--------- Ranecu engine status ---------
372 Initial seed (index) = 0
[1337]373 Current couple of seeds = 490323513, 1321056075
[807]374----------------------------------------
375UserDetectorConstruction deleted.
376UserPhysicsList deleted.
377UserRunAction deleted.
378UserPrimaryGenerator deleted.
379G4 kernel has come to Quit state.
380EventManager deleted.
381UImanager deleted.
382Units table cleared.
383StateManager deleted.
384RunManagerKernel is deleted.
385RunManager is deleting.
[1337]386
Note: See TracBrowser for help on using the repository browser.