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

Last change on this file since 1036 was 807, checked in by garnier, 17 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.