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

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

update

File size: 18.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 = 14 *****
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 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: CO2 density: 0.001 mg/cm3 RadL: 361.873 km Imean: 90.164 eV temperature: 273.15 K pressure: 1.00 atm
24 ---> Element: Hydrogen (C) Z = 6.0 N = 12.0 A = 12.00 g/mole ElmMassFraction: 27.27 % ElmAbundance 33.33 %
25 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 72.73 % ElmAbundance 66.67 %
26
27 Material: D2_gas density: 0.036 mg/cm3 RadL: 13.184 km Imean: 41.800 eV temperature: 273.15 K pressure: 1.00 atm
28 ---> Element: D2_gas ( ) Z = 2.0 N = 2.0 A = 2.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
29
30 Material: liquidArgon density: 1.390 g/cm3 RadL: 14.065 cm Imean: 187.200 eV
31 ---> Element: liquidArgon ( ) Z = 18.0 N = 40.0 A = 39.95 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
32
33 Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Imean: 166.400 eV
34 ---> Element: Aluminium ( ) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
35
36 Material: Silicon density: 2.330 g/cm3 RadL: 9.368 cm Imean: 173.600 eV
37 ---> Element: Silicon ( ) Z = 14.0 N = 28.1 A = 28.09 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
38
39 Material: Germanium density: 5.323 g/cm3 RadL: 2.301 cm Imean: 332.800 eV
40 ---> Element: Germanium ( ) Z = 32.0 N = 72.6 A = 72.61 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
41
42 Material: BGO density: 7.100 g/cm3 RadL: 1.123 cm Imean: 469.142 eV
43 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 15.41 % ElmAbundance 63.16 %
44 ---> Element: Germanium (Ge) Z = 32.0 N = 72.6 A = 72.59 g/mole ElmMassFraction: 17.48 % ElmAbundance 15.79 %
45 ---> Element: Bismuth (Bi) Z = 83.0 N = 209.0 A = 208.98 g/mole ElmMassFraction: 67.10 % ElmAbundance 21.05 %
46
47 Material: Iron density: 7.870 g/cm3 RadL: 1.759 cm Imean: 286.000 eV
48 ---> Element: Iron ( ) Z = 26.0 N = 55.9 A = 55.85 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
49
50 Material: Tungsten density: 19.300 g/cm3 RadL: 3.504 mm Imean: 725.200 eV
51 ---> Element: Tungsten ( ) Z = 74.0 N = 183.8 A = 183.85 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
52
53 Material: Lead density: 11.350 g/cm3 RadL: 5.612 mm Imean: 820.000 eV
54 ---> Element: Lead ( ) Z = 82.0 N = 207.2 A = 207.19 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
55
56 Material: Uranium density: 18.950 g/cm3 RadL: 3.166 mm Imean: 892.400 eV
57 ---> Element: Uranium ( ) Z = 92.0 N = 238.0 A = 238.03 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
58
59
60/run/verbose 2
61#
62/testem/det/setMat Aluminium
63/testem/det/setSize 10 m
64#
65/run/initialize
66userDetector->Construct() start.
67
68 The Box is 10 m of Aluminium
69Aluminium is registered to the default region.
70physicsList->Construct() start.
71physicsList->Construct() start.
72physicsList->setCut() start.
73PhysicsList::SetCuts:CutLength : 1 mm
74#
75/testem/gun/setDefault
76/gun/particle e-
77/gun/energy 100 MeV
78#
79/testem/histo/setHisto 1 100 0 50 cm
80----> SetHisto 1: total track length of primary particle (cm); 100 bins from 0 cm to 50 cm
81/testem/histo/setHisto 2 100 0 300 none
82----> SetHisto 2: nb steps of primary particle; 100 bins from 0 none to 300 none
83/testem/histo/setHisto 3 200 0 20 mm
84----> SetHisto 3: step size of primary particle (mm); 200 bins from 0 mm to 20 mm
85#
86/testem/histo/setFileType hbook
87###/testem/histo/setFileType root
88###/testem/histo/setFileType XML
89#
90/testem/phys/setCuts 1 mm
91/testem/histo/setFileName run0
92/run/beamOn 2000
93
94phot: Total cross sections from Sandia parametrisation.
95 Sampling according PhotoElectric model
96
97compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV
98 Sampling according Klein-Nishina model
99 tables are built for gamma
100 Lambda tables from 100 eV to 100 GeV in 90 bins.
101
102conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z;
103 sampling secondary e+e- according Bethe-Heitler model
104 tables are built for gamma
105 Lambda tables from 1.022 MeV to 100 GeV in 100 bins.
106
107msc: Model variant of multiple scattering for e-
108 Lambda tables from 100 eV to 100 TeV in 120 bins.
109 LateralDisplacementFlag= 1 Skin= 2
110 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 2
111### The process eBrem is added to the list of collaborative processes of eIoni
112### The process eBrem is added to the list of collaborative processes of eIoni
113### The process muBrems is added to the list of collaborative processes of muIoni
114### The process muPairProd is added to the list of collaborative processes of muIoni
115### The process muBrems is added to the list of collaborative processes of muIoni
116### The process muPairProd is added to the list of collaborative processes of muIoni
117
118eIoni: tables are built for e-
119 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
120 Lambda tables from threshold to 100 TeV in 120 bins.
121 Delta cross sections and sampling from MollerBhabha model
122 Good description from 1 KeV to 100 GeV.
123 Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
124 CSDA range table up to 1 GeV in 70 bins.
125 Subcutoff sampling in 1 regions
126
127eBrem: tables are built for e-
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 Total cross sections and sampling from StandBrem model (based on the EEDL data library)
131 Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
132 Subcutoff sampling in 1 regions
133
134eIoni: tables are built for e+
135 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
136 Lambda tables from threshold to 100 TeV in 120 bins.
137 Delta cross sections and sampling from MollerBhabha model
138 Good description from 1 KeV to 100 GeV.
139 Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
140 CSDA range table up to 1 GeV in 70 bins.
141 Subcutoff sampling in 1 regions
142
143eBrem: tables are built for e+
144 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
145 Lambda tables from threshold to 100 TeV in 120 bins.
146 Total cross sections and sampling from StandBrem model (based on the EEDL data library)
147 Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
148 Subcutoff sampling in 1 regions
149
150annihil: Sampling according eplus2gg model
151 tables are built for e+
152 Lambda tables from 100 eV to 100 TeV in 120 bins.
153
154msc: Model variant of multiple scattering for proton
155 Lambda tables from 100 eV to 100 TeV in 120 bins.
156 LateralDisplacementFlag= 1 Skin= 2
157 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 2
158
159hIoni: tables are built for proton
160 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
161 Lambda tables from threshold to 100 TeV in 120 bins.
162 Scaling relation is used from proton dE/dx and range.
163 Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
164 Parametrisation from Bragg for protons below. NuclearStopping= 1
165 Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
166 CSDA range table up to 1 GeV in 70 bins.
167 Subcutoff sampling in 1 regions
168
169msc: Model variant of multiple scattering for GenericIon
170 LateralDisplacementFlag= 0 Skin= 0
171 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 1
172
173ionIoni: tables are built for GenericIon
174 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
175 Lambda tables from threshold to 100 TeV in 120 bins.
176 Scaling relation is used from proton dE/dx and range.
177 Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
178 Parametrisation from BraggIon for protons below. NuclearStopping= 1
179
180 Stopping Power data for 8 ion/material pairs are used.
181 Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
182 CSDA range table up to 1 GeV in 70 bins.
183 Subcutoff sampling in 1 regions
184
185hIoni: tables are built for anti_proton
186 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
187 Lambda tables from threshold to 100 TeV in 120 bins.
188 Scaling relation is used from proton dE/dx and range.
189 Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
190 Parametrisation from Bragg for protons below. NuclearStopping= 1
191 Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
192 CSDA range table up to 1 GeV in 70 bins.
193 Subcutoff sampling in 1 regions
194
195msc: Model variant of multiple scattering for mu+
196 Lambda tables from 100 eV to 100 TeV in 120 bins.
197 LateralDisplacementFlag= 1 Skin= 2
198 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 2
199
200muIoni: tables are built for mu+
201 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
202 Lambda tables from threshold to 100 TeV in 120 bins.
203 Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
204 radiative corrections for E > 1 GeV
205 Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
206 CSDA range table up to 1 GeV in 70 bins.
207 Subcutoff sampling in 1 regions
208
209muBrems: tables are built for mu+
210 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
211 Lambda tables from threshold to 100 TeV in 120 bins.
212 Parametrised model
213 Subcutoff sampling in 1 regions
214
215muPairProd: tables are built for mu+
216 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
217 Lambda tables from threshold to 100 TeV in 120 bins.
218 Parametrised model
219 Subcutoff sampling in 1 regions
220
221muIoni: tables are built for mu-
222 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
223 Lambda tables from threshold to 100 TeV in 120 bins.
224 Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
225 radiative corrections for E > 1 GeV
226 Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
227 CSDA range table up to 1 GeV in 70 bins.
228 Subcutoff sampling in 1 regions
229
230muBrems: tables are built for mu-
231 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
232 Lambda tables from threshold to 100 TeV in 120 bins.
233 Parametrised model
234 Subcutoff sampling in 1 regions
235
236muPairProd: 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 Subcutoff sampling in 1 regions
241
242hIoni: tables are built for pi+
243 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
244 Lambda tables from threshold to 100 TeV in 120 bins.
245 Scaling relation is used from proton dE/dx and range.
246 Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
247 Parametrisation from Bragg for protons below. NuclearStopping= 1
248 Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
249 CSDA range table up to 1 GeV in 70 bins.
250 Subcutoff sampling in 1 regions
251
252msc: Model variant of multiple scattering for pi-
253 Lambda tables from 100 eV to 100 TeV in 120 bins.
254 LateralDisplacementFlag= 1 Skin= 2
255 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 2
256
257hIoni: tables are built for pi-
258 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
259 Lambda tables from threshold to 100 TeV in 120 bins.
260 Scaling relation is used from proton dE/dx and range.
261 Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
262 Parametrisation from Bragg for protons below. NuclearStopping= 1
263 Step function: finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1
264 CSDA range table up to 1 GeV in 70 bins.
265 Subcutoff sampling in 1 regions
266
267Region <DefaultRegionForTheWorld> -- appears in <Aluminium> world volume
268 Root logical volume(s) : Aluminium
269 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
270 Materials : Aluminium
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 : Aluminium
277 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm
278 Energy thresholds : gamma 6.88731 keV e- 596.68 keV e+ 568.011 keV
279 Region(s) which use this couple :
280 DefaultRegionForTheWorld
281
282====================================================================
283
284Start closing geometry.
285G4GeometryManager::ReportVoxelStats -- Voxel Statistics
286
287 Total memory consumed for geometry optimisation: 0 kByte
288 Total CPU time elapsed for geometry optimisation: 0 seconds
289### Run 0 start.
290
291--------- Ranecu engine status ---------
292 Initial seed (index) = 0
293 Current couple of seeds = 9876, 54321
294----------------------------------------
295Start Run processing.
296
297---> Begin of Event: 0
298Run terminated.
299Run Summary
300 Number of events processed : 2000
301 User=13.41s Real=18.92s Sys=0.58s
302
303 ======================== run summary ======================
304
305 The run was: 2000 e- of 100 MeV through 10 m of Aluminium (density: 2.7 g/cm3 )
306
307 ============================================================
308
309 total energy deposit: 99.792 MeV
310
311 nb tracks/event neutral: 27.749 charged: 147.2
312 nb steps/event neutral: 140 charged: 470.24
313
314 nb of process calls per event:
315 eBrem eIoni compt phot msc conv annihilTransportation
316 25.228 300.25 112.25 25.379 143.34 1.41 1.41 0.9675
317
318---------------------------------------------------------
319 Primary particle :
320 true Range = 11.136 cm rms = 3.8075 cm
321 proj Range = 9.6876 cm rms = 3.5531 cm
322 proj/true = 0.86994
323 transverse dispersion at end = 1.6484 cm
324 mass true Range from simulation = 30.067 g/cm2
325 from PhysicsTable (csda range) = 31.791 g/cm2
326---------------------------------------------------------
327
328
329--------- Ranecu engine status ---------
330 Initial seed (index) = 0
331 Current couple of seeds = 1933278320, 292162522
332----------------------------------------
333#
334/process/eLoss/verbose 0
335#
336/testem/phys/setCuts 10 um
337/testem/histo/setFileName run1
338/run/beamOn 2000
339
340Region <DefaultRegionForTheWorld> -- appears in <Aluminium> world volume
341 Root logical volume(s) : Aluminium
342 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
343 Materials : Aluminium
344 Production cuts : gamma 10 um e- 10 um e+ 10 um
345
346========= Table of registered couples ==============================
347
348Index : 0 used in the geometry : Yes recalculation needed : No
349 Material : Aluminium
350 Range cuts : gamma 10 um e- 10 um e+ 10 um
351 Energy thresholds : gamma 990 eV e- 33.887552 keV e+ 33.472938 keV
352 Region(s) which use this couple :
353 DefaultRegionForTheWorld
354
355====================================================================
356
357### Run 1 start.
358
359--------- Ranecu engine status ---------
360 Initial seed (index) = 0
361 Current couple of seeds = 1933278320, 292162522
362----------------------------------------
363Start Run processing.
364
365---> Begin of Event: 0
366Run terminated.
367Run Summary
368 Number of events processed : 2000
369 User=24.74s Real=30.49s Sys=0.6s
370
371 ======================== run summary ======================
372
373 The run was: 2000 e- of 100 MeV through 10 m of Aluminium (density: 2.7 g/cm3 )
374
375 ============================================================
376
377 total energy deposit: 99.813 MeV
378
379 nb tracks/event neutral: 36.679 charged: 291.47
380 nb steps/event neutral: 146.65 charged: 905.42
381
382 nb of process calls per event:
383 eIoni eBrem msc phot compt conv annihilTransportation
384 583.88 34.748 285.38 34.35 109.97 1.412 1.412 0.9225
385
386---------------------------------------------------------
387 Primary particle :
388 true Range = 11.151 cm rms = 3.8821 cm
389 proj Range = 9.7241 cm rms = 3.6334 cm
390 proj/true = 0.87201
391 transverse dispersion at end = 1.707 cm
392 mass true Range from simulation = 30.109 g/cm2
393 from PhysicsTable (csda range) = 31.791 g/cm2
394---------------------------------------------------------
395
396
397--------- Ranecu engine status ---------
398 Initial seed (index) = 0
399 Current couple of seeds = 268076810, 1628262656
400----------------------------------------
401UserDetectorConstruction deleted.
402UserPhysicsList deleted.
403UserRunAction deleted.
404UserPrimaryGenerator deleted.
405G4 kernel has come to Quit state.
406EventManager deleted.
407Default detector region deleted.
408UImanager deleted.
409Units table cleared.
410StateManager deleted.
411RunManagerKernel is deleted.
412RunManager is deleting.
Note: See TracBrowser for help on using the repository browser.