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

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