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

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

update
File size: 14.8 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 = 15 *****
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
63/run/verbose 2
64#
65/testem/det/setMat Water
66/testem/det/setSize 1 cm
67#
68#### /testem/phys/addPhysics livermore
69#### /testem/phys/addPhysics penelope
70#
71/run/initialize
72userDetector->Construct() start.
73
74 The Box is 1 cm  of Water
75Water is registered to the default region.
76physicsList->Construct() start.
77physicsList->Construct() start.
78physicsList->setCut() start.
79#
80/gun/particle gamma
81/gun/energy 100 keV
82/run/beamOn 100000
83
84phot:  Total cross sections from Sandia parametrisation.
85      Sampling according PhotoElectric model
86
87compt:  Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV
88      Sampling according Klein-Nishina model
89      tables are built for  gamma
90      Lambda tables from 100 eV  to 100 GeV in 90 bins.
91
92conv:  Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z;
93      sampling secondary e+e- according Bethe-Heitler model
94      tables are built for  gamma
95      Lambda tables from 1.022 MeV to 100 GeV in 100 bins.
96
97eIoni:   tables are built for  e-
98      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
99      Lambda tables from threshold to 100 TeV in 120 bins.
100      Delta cross sections and sampling from MollerBhabha model
101      Good description from 1 KeV to 100 GeV.
102      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
103
104eBrem:   tables are built for  e-
105      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
106      Lambda tables from threshold to 100 TeV in 120 bins.
107      Total cross sections and sampling from StandBrem model (based on the EEDL data library)
108      Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
109
110eIoni:   tables are built for  e+
111      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
112      Lambda tables from threshold to 100 TeV in 120 bins.
113      Delta cross sections and sampling from MollerBhabha model
114      Good description from 1 KeV to 100 GeV.
115      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
116
117eBrem:   tables are built for  e+
118      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
119      Lambda tables from threshold to 100 TeV in 120 bins.
120      Total cross sections and sampling from StandBrem model (based on the EEDL data library)
121      Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
122
123annihil:       Sampling according eplus2gg model
124      tables are built for  e+
125      Lambda tables from 100 eV  to 100 TeV in 120 bins.
126
127hIoni:   tables are built for  proton
128      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
129      Lambda tables from threshold to 100 TeV in 120 bins.
130      Scaling relation is used from proton dE/dx and range.
131      Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
132      Parametrisation from Bragg for protons below. NuclearStopping= 1
133      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
134
135ionIoni:   tables are built for  GenericIon
136      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
137      Lambda tables from threshold to 100 TeV in 120 bins.
138      Scaling relation is used from proton dE/dx and range.
139      Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
140      Parametrisation from BraggIon for protons below. NuclearStopping= 1
141
142      Stopping Power data for 8 ion/material pairs are used.
143      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
144
145hIoni:   tables are built for  anti_proton
146      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
147      Lambda tables from threshold to 100 TeV in 120 bins.
148      Scaling relation is used from proton dE/dx and range.
149      Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
150      Parametrisation from Bragg for protons below. NuclearStopping= 1
151      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
152
153muIoni:   tables are built for  mu+
154      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
155      Lambda tables from threshold to 100 TeV in 120 bins.
156      Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
157      radiative corrections for E > 1 GeV
158      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
159
160muBrems:   tables are built for  mu+
161      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
162      Lambda tables from threshold to 100 TeV in 120 bins.
163      Parametrised model
164
165muPairProd:   tables are built for  mu+
166      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
167      Lambda tables from threshold to 100 TeV in 120 bins.
168      Parametrised model
169
170muIoni:   tables are built for  mu-
171      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
172      Lambda tables from threshold to 100 TeV in 120 bins.
173      Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
174      radiative corrections for E > 1 GeV
175      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
176
177muBrems:   tables are built for  mu-
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      Parametrised model
181
182muPairProd:   tables are built for  mu-
183      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
184      Lambda tables from threshold to 100 TeV in 120 bins.
185      Parametrised model
186
187hIoni:   tables are built for  pi+
188      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
189      Lambda tables from threshold to 100 TeV in 120 bins.
190      Scaling relation is used from proton dE/dx and range.
191      Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
192      Parametrisation from Bragg for protons below. NuclearStopping= 1
193      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
194
195hIoni:   tables are built for  pi-
196      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
197      Lambda tables from threshold to 100 TeV in 120 bins.
198      Scaling relation is used from proton dE/dx and range.
199      Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
200      Parametrisation from Bragg for protons below. NuclearStopping= 1
201      Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 0, fluct: 0
202
203Region <DefaultRegionForTheWorld> -- appears in <Water> world volume
204 Root logical volume(s) : Water
205 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
206 Materials : Water
207 Production cuts :  gamma 1 mm     e- 1 mm     e+ 1 mm
208
209========= Table of registered couples ==============================
210
211Index : 0     used in the geometry : Yes     recalculation needed : No
212 Material : Water
213 Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
214 Energy thresholds :  gamma 2.90186 keV    e- 347.138 keV    e+ 338.695 keV
215 Region(s) which use this couple :
216    DefaultRegionForTheWorld
217
218====================================================================
219
220Start closing geometry.
221G4GeometryManager::ReportVoxelStats -- Voxel Statistics
222
223    Total memory consumed for geometry optimisation:   0 kByte
224    Total CPU time elapsed for geometry optimisation: 0 seconds
225### Run 0 start.
226
227--------- Ranecu engine status ---------
228 Initial seed (index) = 0
229 Current couple of seeds = 9876, 54321
230----------------------------------------
231Start Run processing.
232
233---> Begin of Event: 0
234
235---> Begin of Event: 10000
236
237---> Begin of Event: 20000
238
239---> Begin of Event: 30000
240
241---> Begin of Event: 40000
242
243---> Begin of Event: 50000
244
245---> Begin of Event: 60000
246
247---> Begin of Event: 70000
248
249---> Begin of Event: 80000
250
251---> Begin of Event: 90000
252Run terminated.
253Run Summary
254  Number of events processed : 100000
255  User=1.17s Real=1.16s Sys=0s
256
257 The run consists of 100000 gamma of 100 keV through 1 cm  of Water (density: 1 g/cm3 )
258
259 Process calls frequency --->   Transportation = 84934  compt = 14842   phot = 224
260
261 Nb of incident particles unaltered after 1 cm  of Water : 84934 over 100000 incident particles.  Ratio = 84.934 %
262 ---> CrossSection per volume:  0.1633 cm^-1    CrossSection per mass: 0.1633 cm2/g
263
264 Verification from G4EmCalculator:
265        compt= 0.1623 cm2/g     phot= 0.0026853 cm2/g   total= 0.16498 cm2/g
266        Expected ratio of transmitted particles= 84.791 %
267
268--------- Ranecu engine status ---------
269 Initial seed (index) = 0
270 Current couple of seeds = 507682210, 1650158306
271----------------------------------------
272#
273/gun/particle e-
274/gun/energy 100 MeV
275/run/beamOn 100000
276
277Region <DefaultRegionForTheWorld> -- appears in <Water> world volume
278 Root logical volume(s) : Water
279 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
280 Materials : Water
281 Production cuts :  gamma 1 mm     e- 1 mm     e+ 1 mm
282
283========= Table of registered couples ==============================
284
285Index : 0     used in the geometry : Yes     recalculation needed : No
286 Material : Water
287 Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
288 Energy thresholds :  gamma 2.9018564 keV    e- 347.13789 keV    e+ 338.6954 keV
289 Region(s) which use this couple :
290    DefaultRegionForTheWorld
291
292====================================================================
293
294### Run 1 start.
295
296--------- Ranecu engine status ---------
297 Initial seed (index) = 0
298 Current couple of seeds = 507682210, 1650158306
299----------------------------------------
300Start Run processing.
301
302---> Begin of Event: 0
303
304---> Begin of Event: 10000
305
306---> Begin of Event: 20000
307
308---> Begin of Event: 30000
309
310---> Begin of Event: 40000
311
312---> Begin of Event: 50000
313
314---> Begin of Event: 60000
315
316---> Begin of Event: 70000
317
318---> Begin of Event: 80000
319
320---> Begin of Event: 90000
321Run terminated.
322Run Summary
323  Number of events processed : 100000
324  User=1.31s Real=1.32s Sys=0s
325
326 The run consists of 100000 e- of 100 MeV through 1 cm  of Water (density: 1 g/cm3 )
327
328 Process calls frequency --->   eBrem = 29069   Transportation = 52915  eIoni = 18016
329
330 Nb of incident particles unaltered after 1 cm  of Water : 52915 over 100000 incident particles.  Ratio = 52.915 %
331 ---> CrossSection per volume:  0.63648 cm^-1   CrossSection per mass: 0.63648 cm2/g
332
333 Verification from G4EmCalculator:
334        eBrem= 0.39397 cm2/g    eIoni= 0.24497 cm2/g    total= 0.63894 cm2/g
335        Expected ratio of transmitted particles= 52.785 %
336
337--------- Ranecu engine status ---------
338 Initial seed (index) = 0
339 Current couple of seeds = 1186383716, 309205105
340----------------------------------------
341UserDetectorConstruction deleted.
342UserPhysicsList deleted.
343UserRunAction deleted.
344UserPrimaryGenerator deleted.
345G4 kernel has come to Quit state.
346EventManager deleted.
347Default detector region deleted.
348UImanager deleted.
349Units table cleared.
350StateManager deleted.
351RunManagerKernel is deleted.
352RunManager is deleting.
Note: See TracBrowser for help on using the repository browser.