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

Last change on this file since 1036 was 807, checked in by garnier, 17 years ago

update

File size: 15.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 = 13 *****
11
12 Material: G4_Al density: 2.699 g/cm3 RadL: 8.896 cm Imean: 166.000 eV
13 ---> Element: Al (Al) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
14
15 Material: G4_Si density: 2.330 g/cm3 RadL: 9.366 cm Imean: 173.000 eV
16 ---> Element: Si (Si) Z = 14.0 N = 28.1 A = 28.09 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
17
18 Material: G4_Fe density: 7.874 g/cm3 RadL: 1.757 cm Imean: 286.000 eV
19 ---> Element: Fe (Fe) Z = 26.0 N = 55.8 A = 55.85 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
20
21 Material: G4_Ge density: 5.323 g/cm3 RadL: 2.301 cm Imean: 350.000 eV
22 ---> Element: Ge (Ge) Z = 32.0 N = 72.6 A = 72.61 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
23
24 Material: G4_Mo density: 10.220 g/cm3 RadL: 9.591 mm Imean: 424.000 eV
25 ---> Element: Mo (Mo) Z = 42.0 N = 95.9 A = 95.93 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
26
27 Material: G4_Ta density: 16.654 g/cm3 RadL: 4.094 mm Imean: 718.000 eV
28 ---> Element: Ta (Ta) Z = 73.0 N = 180.9 A = 180.95 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
29
30 Material: G4_W density: 19.300 g/cm3 RadL: 3.504 mm Imean: 727.000 eV
31 ---> Element: W (W) Z = 74.0 N = 183.8 A = 183.84 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
32
33 Material: G4_Au density: 19.320 g/cm3 RadL: 3.344 mm Imean: 790.000 eV
34 ---> Element: Au (Au) Z = 79.0 N = 197.0 A = 196.97 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
35
36 Material: G4_Pb density: 11.350 g/cm3 RadL: 5.613 mm Imean: 823.000 eV
37 ---> Element: Pb (Pb) Z = 82.0 N = 207.2 A = 207.22 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
38
39 Material: G4_PbWO4 density: 8.280 g/cm3 RadL: 8.925 mm Imean: 541.543 eV
40 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 14.06 % ElmAbundance 66.67 %
41 ---> Element: Pb (Pb) Z = 82.0 N = 207.2 A = 207.22 g/mole ElmMassFraction: 45.54 % ElmAbundance 16.67 %
42 ---> Element: W (W) Z = 74.0 N = 183.8 A = 183.84 g/mole ElmMassFraction: 40.40 % ElmAbundance 16.67 %
43
44 Material: G4_SODIUM_IODIDE density: 3.667 g/cm3 RadL: 2.588 cm Imean: 452.000 eV
45 ---> Element: Na (Na) Z = 11.0 N = 23.0 A = 22.99 g/mole ElmMassFraction: 15.34 % ElmAbundance 50.00 %
46 ---> Element: I (I) Z = 53.0 N = 126.9 A = 126.90 g/mole ElmMassFraction: 84.66 % ElmAbundance 50.00 %
47
48 Material: G4_AIR density: 1.205 mg/cm3 RadL: 303.921 m Imean: 85.700 eV temperature: 273.15 K pressure: 1.00 atm
49 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 0.01 % ElmAbundance 0.02 %
50 ---> Element: N (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 75.53 % ElmAbundance 78.44 %
51 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 23.18 % ElmAbundance 21.07 %
52 ---> Element: Ar (Ar) Z = 18.0 N = 39.9 A = 39.95 g/mole ElmMassFraction: 1.28 % ElmAbundance 0.47 %
53
54 Material: Water density: 1.000 g/cm3 RadL: 36.083 cm Imean: 75.000 eV
55 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.19 % ElmAbundance 66.67 %
56 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.81 % ElmAbundance 33.33 %
57
58
59Visualization Manager instantiating...
60Visualization Manager initialising...
61Registering graphics systems...
62
63You have successfully registered the following graphics systems.
64Current available graphics systems are:
65 ASCIITree (ATree)
66 DAWNFILE (DAWNFILE)
67 G4HepRep (HepRepXML)
68 G4HepRepFile (HepRepFile)
69 RayTracer (RayTracer)
70 VRML1FILE (VRML1FILE)
71 VRML2FILE (VRML2FILE)
72 OpenGLImmediateX (OGLIX)
73 OpenGLStoredX (OGLSX)
74
75Registering model factories...
76
77You have successfully registered the following model factories.
78Registered model factories:
79 generic
80 drawByCharge
81 drawByParticleID
82 drawByOriginVolume
83 drawByAttribute
84
85Registered filter factories:
86 chargeFilter
87 particleFilter
88 originVolumeFilter
89 attributeFilter
90
91/run/verbose 2
92#
93/testem/det/setMat G4_Si
94/testem/det/setSizeX 1 mm
95/testem/det/setSizeYZ 1 mm
96#
97/testem/phys/addPhysics standard
98PhysicsList::AddPhysicsList: <standard>
99###/testem/phys/addPhysics penelope # em physics
100#
101/run/initialize
102userDetector->Construct() start.
103
104---------------------------------------------------------
105---> The Absorber is 1 mm of G4_Si divided in 1 slices of 1 mm
106
107---------------------------------------------------------
108Absorber is registered to the default region.
109physicsList->Construct() start.
110physicsList->Construct() start.
111physicsList->setCut() start.
112PhysicsList::SetCuts:CutLength : 1 mm
113#
114/gun/particle e-
115/gun/energy 500 keV
116#
117/testem/histo/setFileType hbook
118/testem/histo/setFileName testem11
119/testem/histo/setHisto 1 100 0. 1. mm
120----> SetHisto 1: Edep (MeV/mm) along absorber (mm); 100 bins from 0 mm to 1 mm
121/testem/histo/setHisto 2 102 0. 510. keV
122----> SetHisto 2: total Energy deposited in absorber (keV); 102 bins from 0 keV to 510 keV
123/testem/histo/setHisto 3 100 0. 2. mm
124----> SetHisto 3: true track length of the primary particle (mm); 100 bins from 0 mm to 2 mm
125/testem/histo/setHisto 4 100 0. 200. um
126----> SetHisto 4: true step size of the primary particle (um); 100 bins from 0 um to 200 um
127/testem/histo/setHisto 5 100 0. 1. mm
128----> SetHisto 5: projected range of the primary particle (mm); 100 bins from 0 mm to 1 mm
129#
130/testem/event/printModulo 10000
131#
132/run/beamOn 20000
133
134phot: Total cross sections from Sandia parametrisation.
135 Sampling according PhotoElectric model
136
137compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV
138 Sampling according Klein-Nishina model
139 tables are built for gamma
140 Lambda tables from 100 eV to 100 GeV in 90 bins.
141
142conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z;
143 sampling secondary e+e- according Bethe-Heitler model
144 tables are built for gamma
145 Lambda tables from 1.022 MeV to 100 GeV in 100 bins.
146
147msc: Model variant of multiple scattering for e-
148 Lambda tables from 100 eV to 100 TeV in 120 bins.
149 LateralDisplacementFlag= 1 Skin= 0
150 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 1
151
152eIoni: tables are built for e-
153 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
154 Lambda tables from threshold to 100 TeV in 120 bins.
155 Delta cross sections and sampling from MollerBhabha model
156 Good description from 1 KeV to 100 GeV.
157 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
158 CSDA range table up to 1 GeV in 70 bins.
159
160eBrem: tables are built for e-
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 Total cross sections and sampling from StandBrem model (based on the EEDL data library)
164 Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
165
166eIoni: tables are built for e+
167 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
168 Lambda tables from threshold to 100 TeV in 120 bins.
169 Delta cross sections and sampling from MollerBhabha model
170 Good description from 1 KeV to 100 GeV.
171 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
172 CSDA range table up to 1 GeV in 70 bins.
173
174eBrem: tables are built for e+
175 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
176 Lambda tables from threshold to 100 TeV in 120 bins.
177 Total cross sections and sampling from StandBrem model (based on the EEDL data library)
178 Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
179
180annihil: Sampling according eplus2gg model
181 tables are built for e+
182 Lambda tables from 100 eV to 100 TeV in 120 bins.
183
184msc: Model variant of multiple scattering for proton
185 Lambda tables from 100 eV to 100 TeV in 120 bins.
186 LateralDisplacementFlag= 1 Skin= 0
187 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 1
188
189hIoni: tables are built for proton
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 Bragg for protons below. NuclearStopping= 1
195 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
196 CSDA range table up to 1 GeV in 70 bins.
197
198msc: Model variant of multiple scattering for GenericIon
199 LateralDisplacementFlag= 0 Skin= 0
200 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 1
201
202ionIoni: tables are built for GenericIon
203 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
204 Lambda tables from threshold to 100 TeV in 120 bins.
205 Scaling relation is used from proton dE/dx and range.
206 Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
207 Parametrisation from BraggIon for protons below. NuclearStopping= 1
208
209 Stopping Power data for 8 ion/material pairs are used.
210 Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1
211 CSDA range table up to 1 GeV in 70 bins.
212
213hIoni: tables are built for anti_proton
214 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
215 Lambda tables from threshold to 100 TeV in 120 bins.
216 Scaling relation is used from proton dE/dx and range.
217 Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
218 Parametrisation from Bragg for protons below. NuclearStopping= 1
219 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
220 CSDA range table up to 1 GeV in 70 bins.
221
222msc: Model variant of multiple scattering for mu+
223 Lambda tables from 100 eV to 100 TeV in 120 bins.
224 LateralDisplacementFlag= 1 Skin= 0
225 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 1
226
227muIoni: tables are built for mu+
228 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
229 Lambda tables from threshold to 100 TeV in 120 bins.
230 Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
231 radiative corrections for E > 1 GeV
232 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
233 CSDA range table up to 1 GeV in 70 bins.
234
235muBrems: tables are built for mu+
236 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
237 Lambda tables from threshold to 100 TeV in 120 bins.
238 Parametrised model
239
240muPairProd: tables are built for mu+
241 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
242 Lambda tables from threshold to 100 TeV in 120 bins.
243 Parametrised model
244
245muIoni: tables are built for mu-
246 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
247 Lambda tables from threshold to 100 TeV in 120 bins.
248 Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
249 radiative corrections for E > 1 GeV
250 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
251 CSDA range table up to 1 GeV in 70 bins.
252
253muBrems: tables are built for mu-
254 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
255 Lambda tables from threshold to 100 TeV in 120 bins.
256 Parametrised model
257
258muPairProd: tables are built for mu-
259 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
260 Lambda tables from threshold to 100 TeV in 120 bins.
261 Parametrised model
262
263hIoni: tables are built for pi+
264 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
265 Lambda tables from threshold to 100 TeV in 120 bins.
266 Scaling relation is used from proton dE/dx and range.
267 Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
268 Parametrisation from Bragg for protons below. NuclearStopping= 1
269 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
270 CSDA range table up to 1 GeV in 70 bins.
271
272msc: Model variant of multiple scattering for pi-
273 Lambda tables from 100 eV to 100 TeV in 120 bins.
274 LateralDisplacementFlag= 1 Skin= 0
275 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 1
276
277hIoni: tables are built for pi-
278 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
279 Lambda tables from threshold to 100 TeV in 120 bins.
280 Scaling relation is used from proton dE/dx and range.
281 Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
282 Parametrisation from Bragg for protons below. NuclearStopping= 1
283 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
284 CSDA range table up to 1 GeV in 70 bins.
285
286Region <DefaultRegionForTheWorld> -- appears in <Absorber> world volume
287 Root logical volume(s) : Absorber
288 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
289 Materials : G4_Si
290 Production cuts : gamma 1 mm e- 1 mm e+ 1 mm
291
292========= Table of registered couples ==============================
293
294Index : 0 used in the geometry : Yes recalculation needed : No
295 Material : G4_Si
296 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm
297 Energy thresholds : gamma 6.88731 keV e- 540.718 keV e+ 521.113 keV
298 Region(s) which use this couple :
299 DefaultRegionForTheWorld
300
301====================================================================
302
303Start closing geometry.
304G4GeometryManager::ReportVoxelStats -- Voxel Statistics
305
306 Total memory consumed for geometry optimisation: 0 kByte
307 Total CPU time elapsed for geometry optimisation: 0 seconds
308### Run 0 start.
309
310--------- Ranecu engine status ---------
311 Initial seed (index) = 0
312 Current couple of seeds = 9876, 54321
313----------------------------------------
314
315---> stepMax from HistoManager = 10 um
316Start Run processing.
317
318---> Begin of Event: 0
319
320---> Begin of Event: 10000
321Run terminated.
322Run Summary
323 Number of events processed : 20000
324 User=22.26s Real=71.63s Sys=5.33s
325
326 ======================== run summary ======================
327
328 The run consists of 20000 e- of 500.00 keV through 1.00 mm of G4_Si (density: 2.33 g/cm3 )
329
330 ============================================================
331
332 Total Energy deposited = 449.403 keV +- 103.444 keV
333
334 Track length of primary track = 875.191 um +- 189.287 um
335 Range from EmCalculator = 937.365 um (from full dE/dx)
336
337 Projected range = 336.299 um +- 207.386 um
338
339 Nb of steps of primary track = 88.04 +- 18.91 Step size= 9.933 um +- 86.256 nm
340
341 absorbed = 78.50 % transmit = 5.42 % reflected = 16.08 %
342
343--------- Ranecu engine status ---------
344 Initial seed (index) = 0
345 Current couple of seeds = 1622151268, 626459451
346----------------------------------------
347Graphics systems deleted.
348Visualization Manager deleting...
349UserDetectorConstruction deleted.
350UserPhysicsList deleted.
351UserRunAction deleted.
352UserPrimaryGenerator deleted.
353G4 kernel has come to Quit state.
354EventManager deleted.
355Default detector region deleted.
356UImanager deleted.
357Units table cleared.
358StateManager deleted.
359RunManagerKernel is deleted.
360RunManager is deleting.
Note: See TracBrowser for help on using the repository browser.