source: trunk/examples/extended/electromagnetic/TestEm15/TestEm15.out @ 1309

Last change on this file since 1309 was 807, checked in by garnier, 16 years ago

update

File size: 17.3 KB
Line 
1
2*************************************************************
3 Geant4 version Name: global-V09-00-03    (9-May-2008)
4                      Copyright : Geant4 Collaboration
5                      Reference : NIM A 506 (2003), 250-303
6                            WWW : http://cern.ch/geant4
7*************************************************************
8
9
10***** Table : Nb of materials = 16 *****
11
12 Material:      Air     density:  1.290 mg/cm3  RadL: 285.161 m    Imean:  85.684 eV   temperature: 273.15 K  pressure:   1.00 atm
13   --->  Element: Nitrogen (N)   Z =  7.0   N =  14.0   A =  14.01 g/mole  ElmMassFraction:  70.00 %  ElmAbundance  72.71 %
14   --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  30.00 %  ElmAbundance  27.29 %
15
16 Material: H2liquid     density: 70.800 mg/cm3  RadL:   8.923 m    Imean:  21.800 eV
17   --->  Element: Hydrogen (H)   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
18
19 Material:    Water H_2O   density:  1.000 g/cm3   RadL:  36.092 cm   Imean:  75.000 eV
20   --->  Element: Hydrogen (H)   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction:  11.21 %  ElmAbundance  66.67 %
21   --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  88.79 %  ElmAbundance  33.33 %
22
23 Material: liquidArgon     density:  1.390 g/cm3   RadL:  14.065 cm   Imean: 187.200 eV
24   --->  Element: liquidArgon ( )   Z = 18.0   N =  40.0   A =  39.95 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
25
26 Material:   Carbon     density:  2.267 g/cm3   RadL:  18.833 cm   Imean:  78.000 eV
27   --->  Element: Carbon ( )   Z =  6.0   N =  12.0   A =  12.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
28
29 Material: Aluminium     density:  2.700 g/cm3   RadL:   8.893 cm   Imean: 166.400 eV
30   --->  Element: Aluminium ( )   Z = 13.0   N =  27.0   A =  26.98 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
31
32 Material:  Silicon     density:  2.330 g/cm3   RadL:   9.368 cm   Imean: 173.600 eV
33   --->  Element: Silicon ( )   Z = 14.0   N =  28.1   A =  28.09 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
34
35 Material: Germanium     density:  5.323 g/cm3   RadL:   2.301 cm   Imean: 332.800 eV
36   --->  Element: Germanium ( )   Z = 32.0   N =  72.6   A =  72.61 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
37
38 Material:      NaI     density:  3.670 g/cm3   RadL:   2.586 cm   Imean: 452.000 eV
39   --->  Element: Sodium (Na)   Z = 11.0   N =  23.0   A =  22.99 g/mole  ElmMassFraction:  15.34 %  ElmAbundance  50.00 %
40   --->  Element: Iodine (I)   Z = 53.0   N = 126.9   A = 126.90 g/mole  ElmMassFraction:  84.66 %  ElmAbundance  50.00 %
41
42 Material:   Iodine     density:  4.930 g/cm3   RadL:   1.720 cm   Imean: 492.900 eV
43   --->  Element: Iodine (I)   Z = 53.0   N = 126.9   A = 126.90 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
44
45 Material:      BGO     density:  7.100 g/cm3   RadL:   1.123 cm   Imean: 469.142 eV
46   --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  15.41 %  ElmAbundance  63.16 %
47   --->  Element: Germanium (Ge)   Z = 32.0   N =  72.6   A =  72.59 g/mole  ElmMassFraction:  17.48 %  ElmAbundance  15.79 %
48   --->  Element: Bismuth (Bi)   Z = 83.0   N = 209.0   A = 208.98 g/mole  ElmMassFraction:  67.10 %  ElmAbundance  21.05 %
49
50 Material:     Iron     density:  7.870 g/cm3   RadL:   1.759 cm   Imean: 286.000 eV
51   --->  Element: Iron ( )   Z = 26.0   N =  55.9   A =  55.85 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
52
53 Material: Tungsten     density: 19.300 g/cm3   RadL:   3.504 mm   Imean: 725.200 eV
54   --->  Element: Tungsten ( )   Z = 74.0   N = 183.8   A = 183.85 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
55
56 Material:     Lead     density: 11.350 g/cm3   RadL:   5.612 mm   Imean: 820.000 eV
57   --->  Element: Lead ( )   Z = 82.0   N = 207.2   A = 207.19 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
58
59 Material:  Uranium     density: 18.950 g/cm3   RadL:   3.166 mm   Imean: 892.400 eV
60   --->  Element: Uranium ( )   Z = 92.0   N = 238.0   A = 238.03 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
61
62 Material: Galactic     density:  0.000 mg/cm3  RadL: 204322175.595 pc   Imean:  21.800 eV   temperature:   2.73 K  pressure:   0.00 atm
63   --->  Element: Galactic ( )   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
64
65
66Visualization Manager instantiating...
67Visualization Manager initialising...
68Registering graphics systems...
69
70You have successfully registered the following graphics systems.
71Current available graphics systems are:
72  ASCIITree (ATree)
73  DAWNFILE (DAWNFILE)
74  G4HepRep (HepRepXML)
75  G4HepRepFile (HepRepFile)
76  RayTracer (RayTracer)
77  VRML1FILE (VRML1FILE)
78  VRML2FILE (VRML2FILE)
79  OpenGLImmediateX (OGLIX)
80  OpenGLStoredX (OGLSX)
81
82Registering model factories...
83
84You have successfully registered the following model factories.
85Registered model factories:
86  generic
87  drawByCharge
88  drawByParticleID
89  drawByOriginVolume
90  drawByAttribute
91
92Registered filter factories:
93  chargeFilter
94  particleFilter
95  originVolumeFilter
96  attributeFilter
97
98/run/verbose 2
99#
100/testem/det/setMat Water
101#
102/run/initialize
103userDetector->Construct() start.
104
105 The Box is 100 m   of Water
106World is registered to the default region.
107physicsList->Construct() start.
108physicsList->Construct() start.
109physicsList->setCut() start.
110#
111/process/inactivate eIoni
112/run/physicsModified
113/process/inactivate eBrem
114/run/physicsModified
115#
116/gun/particle e-
117#
118##/tracking/verbose 1
119##/run/beamOn 1
120#
121/gun/energy 5 MeV
122/run/beamOn 10000
123
124phot:  Total cross sections from Sandia parametrisation.
125      Sampling according PhotoElectric model
126
127compt:  Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV
128      Sampling according Klein-Nishina model
129      tables are built for  gamma
130      Lambda tables from 100 eV  to 100 GeV in 90 bins.
131
132conv:  Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z;
133      sampling secondary e+e- according Bethe-Heitler model
134      tables are built for  gamma
135      Lambda tables from 1.022 MeV to 100 GeV in 100 bins.
136
137msc:  Model variant of multiple scattering for e-
138      Lambda tables from 100 eV  to 100 TeV in 120 bins.
139      LateralDisplacementFlag=  1   Skin= 0
140      Boundary/stepping algorithm is active with RangeFactor= 0.02  Step limit type 1
141
142eIoni:   tables are built for  e-
143      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
144      Lambda tables from threshold to 100 TeV in 120 bins.
145      Delta cross sections and sampling from MollerBhabha model
146      Good description from 1 KeV to 100 GeV.
147
148eBrem:   tables are built for  e-
149      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
150      Lambda tables from threshold to 100 TeV in 120 bins.
151      Total cross sections and sampling from StandBrem model (based on the EEDL data library)
152      Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
153      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
154
155eIoni:   tables are built for  e+
156      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
157      Lambda tables from threshold to 100 TeV in 120 bins.
158      Delta cross sections and sampling from MollerBhabha model
159      Good description from 1 KeV to 100 GeV.
160
161eBrem:   tables are built for  e+
162      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
163      Lambda tables from threshold to 100 TeV in 120 bins.
164      Total cross sections and sampling from StandBrem model (based on the EEDL data library)
165      Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
166      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
167
168annihil:       Sampling according eplus2gg model
169      tables are built for  e+
170      Lambda tables from 100 eV  to 100 TeV in 120 bins.
171
172msc:  Model variant of multiple scattering for proton
173      Lambda tables from 100 eV  to 100 TeV in 120 bins.
174      LateralDisplacementFlag=  1   Skin= 0
175      Boundary/stepping algorithm is active with RangeFactor= 0.02  Step limit type 1
176
177hIoni:   tables are built for  proton
178      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
179      Lambda tables from threshold to 100 TeV in 120 bins.
180      Scaling relation is used from proton dE/dx and range.
181      Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
182      Parametrisation from Bragg for protons below. NuclearStopping= 1
183      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
184
185msc:  Model variant of multiple scattering for GenericIon
186      LateralDisplacementFlag=  0   Skin= 0
187      Boundary/stepping algorithm is active with RangeFactor= 0.2  Step limit type 1
188
189ionIoni:   tables are built for  GenericIon
190      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
191      Lambda tables from threshold to 100 TeV in 120 bins.
192      Scaling relation is used from proton dE/dx and range.
193      Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
194      Parametrisation from BraggIon for protons below. NuclearStopping= 1
195
196      Stopping Power data for 8 ion/material pairs are used.
197      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
198
199hIoni:   tables are built for  anti_proton
200      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
201      Lambda tables from threshold to 100 TeV in 120 bins.
202      Scaling relation is used from proton dE/dx and range.
203      Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
204      Parametrisation from Bragg for protons below. NuclearStopping= 1
205      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
206
207msc:  Model variant of multiple scattering for mu+
208      Lambda tables from 100 eV  to 100 TeV in 120 bins.
209      LateralDisplacementFlag=  1   Skin= 0
210      Boundary/stepping algorithm is active with RangeFactor= 0.02  Step limit type 1
211
212muIoni:   tables are built for  mu+
213      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
214      Lambda tables from threshold to 100 TeV in 120 bins.
215      Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
216      radiative corrections for E > 1 GeV
217      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
218
219muBrems:   tables are built for  mu+
220      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
221      Lambda tables from threshold to 100 TeV in 120 bins.
222      Parametrised model
223
224muPairProd:   tables are built for  mu+
225      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
226      Lambda tables from threshold to 100 TeV in 120 bins.
227      Parametrised model
228
229muIoni:   tables are built for  mu-
230      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
231      Lambda tables from threshold to 100 TeV in 120 bins.
232      Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
233      radiative corrections for E > 1 GeV
234      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
235
236muBrems:   tables are built for  mu-
237      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
238      Lambda tables from threshold to 100 TeV in 120 bins.
239      Parametrised model
240
241muPairProd:   tables are built for  mu-
242      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
243      Lambda tables from threshold to 100 TeV in 120 bins.
244      Parametrised model
245
246hIoni:   tables are built for  pi+
247      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
248      Lambda tables from threshold to 100 TeV in 120 bins.
249      Scaling relation is used from proton dE/dx and range.
250      Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
251      Parametrisation from Bragg for protons below. NuclearStopping= 1
252      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
253
254msc:  Model variant of multiple scattering for pi-
255      Lambda tables from 100 eV  to 100 TeV in 120 bins.
256      LateralDisplacementFlag=  1   Skin= 0
257      Boundary/stepping algorithm is active with RangeFactor= 0.02  Step limit type 1
258
259hIoni:   tables are built for  pi-
260      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
261      Lambda tables from threshold to 100 TeV in 120 bins.
262      Scaling relation is used from proton dE/dx and range.
263      Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
264      Parametrisation from Bragg for protons below. NuclearStopping= 1
265      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
266
267Region <DefaultRegionForTheWorld> -- appears in <World> world volume
268 Root logical volume(s) : World
269 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
270 Materials : Galactic Water
271 Production cuts :  gamma 1 mm     e- 1 mm     e+ 1 mm
272
273========= Table of registered couples ==============================
274
275Index : 0     used in the geometry : Yes     recalculation needed : No
276 Material : Galactic
277 Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
278 Energy thresholds :  gamma 990 eV     e- 990 eV     e+ 990 eV
279 Region(s) which use this couple :
280    DefaultRegionForTheWorld
281
282Index : 1     used in the geometry : Yes     recalculation needed : No
283 Material : Water
284 Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
285 Energy thresholds :  gamma 2.90186 keV    e- 347.138 keV    e+ 338.695 keV
286 Region(s) which use this couple :
287    DefaultRegionForTheWorld
288
289====================================================================
290
291Start closing geometry.
292G4GeometryManager::ReportVoxelStats -- Voxel Statistics
293
294    Total memory consumed for geometry optimisation:   0 kByte
295    Total CPU time elapsed for geometry optimisation: 0 seconds
296### Run 0 start.
297
298--------- Ranecu engine status ---------
299 Initial seed (index) = 0
300 Current couple of seeds = 9876, 54321
301----------------------------------------
302Start Run processing.
303
304---> Begin of Event: 0
305Run terminated.
306Run Summary
307  Number of events processed : 10000
308  User=0.25s Real=0.25s Sys=0s
309
310 The run consists of 10000 e- of 5 MeV through 100 m   of Water (density: 1 g/cm3 )
311
312 Process calls frequency --->   msc = 10000
313
314 truePathLength :       2.7521 cm  +- 3.269 nm
315 geomPathLength :       2.0696 cm  +- 11.392 nm
316 lateralDisplac :       1.1482 cm  +- 4.6902 nm
317 Psi            :       506.49 mrad +- 0.00019571 mrad  (29.02 deg +- 1.1213e-05 deg)
318
319 Theta_plane    :       1806.6 mrad  (103.51 deg)
320 phi correlation:       0.2624 +- 0.059273  (std::cos(phi_pos - phi_dir))
321
322 Verification from G4EmCalculator.
323
324 transport mean free path :     8.346 cm
325 range from restrict dE/dx:     2.7521 cm
326 ---> effective facRange  :     1
327
328 compute theta0 from Highland : 620.13 mrad  (35.531 deg)
329
330--------- Ranecu engine status ---------
331 Initial seed (index) = 0
332 Current couple of seeds = 397063712, 552340631
333----------------------------------------
334#
335/gun/energy 100 keV
336/run/beamOn 10000
337
338Region <DefaultRegionForTheWorld> -- appears in <World> world volume
339 Root logical volume(s) : World
340 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
341 Materials : Galactic Water
342 Production cuts :  gamma 1 mm     e- 1 mm     e+ 1 mm
343
344========= Table of registered couples ==============================
345
346Index : 0     used in the geometry : Yes     recalculation needed : No
347 Material : Galactic
348 Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
349 Energy thresholds :  gamma 9.9e+02 eV     e- 9.9e+02 eV     e+ 9.9e+02 eV
350 Region(s) which use this couple :
351    DefaultRegionForTheWorld
352
353Index : 1     used in the geometry : Yes     recalculation needed : No
354 Material : Water
355 Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
356 Energy thresholds :  gamma 2.9 keV    e- 3.5e+02 keV    e+ 3.4e+02 keV
357 Region(s) which use this couple :
358    DefaultRegionForTheWorld
359
360====================================================================
361
362### Run 1 start.
363
364--------- Ranecu engine status ---------
365 Initial seed (index) = 0
366 Current couple of seeds = 397063712, 552340631
367----------------------------------------
368Start Run processing.
369
370---> Begin of Event: 0
371Run terminated.
372Run Summary
373  Number of events processed : 10000
374  User=0.23s Real=0.23s Sys=0s
375
376 The run consists of 10000 e- of 100 keV through 100 m   of Water (density: 1 g/cm3 )
377
378 Process calls frequency --->   msc = 10000
379
380 truePathLength :       6.3051 um  +- 5404.2 fm
381 geomPathLength :       6.1807 um  +- 878.92 fm
382 lateralDisplac :       923.68 nm  +- 256.83 fm
383 Psi            :       148.35 mrad +- 8.7703e-05 mrad  (8.4998 deg +- 5.025e-06 deg)
384
385 Theta_plane    :       258.19 mrad  (14.793 deg)
386 phi correlation:       0.077745 +- 0.12942  (std::cos(phi_pos - phi_dir))
387
388 Verification from G4EmCalculator.
389
390 transport mean free path :     157.63 um
391 range from restrict dE/dx:     143.16 um
392 ---> effective facRange  :     0.04
393
394 compute theta0 from Highland : 180.68 mrad  (10.352 deg)
395
396--------- Ranecu engine status ---------
397 Initial seed (index) = 0
398 Current couple of seeds = 707823211, 530749218
399----------------------------------------
400Graphics systems deleted.
401Visualization Manager deleting...
402UserDetectorConstruction deleted.
403UserPhysicsList deleted.
404UserRunAction deleted.
405UserPrimaryGenerator deleted.
406G4 kernel has come to Quit state.
407EventManager deleted.
408Default detector region deleted.
409UImanager deleted.
410Units table cleared.
411StateManager deleted.
412RunManagerKernel is deleted.
413RunManager is deleting.
Note: See TracBrowser for help on using the repository browser.