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

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