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

Last change on this file since 1309 was 807, checked in by garnier, 16 years ago

update

File size: 15.6 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### NIST DataBase for Materials is used
10G4Material: Mean excitation energy is changed for G4_AIR Iold= 85.627eV; Inew= 85.7 eV;
11G4Material: Mean excitation energy is changed for G4_Al Iold= 166.4eV; Inew= 166 eV;
12G4Material: Mean excitation energy is changed for G4_Si Iold= 173.6eV; Inew= 173 eV;
13G4Material: Mean excitation energy is changed for G4_Fe Iold= 286eV; Inew= 286 eV;
14G4Material: Mean excitation energy is changed for G4_CESIUM_IODIDE Iold= 491.166eV; Inew= 553.1 eV;
15/run/verbose 1
16#/tracking/verbose 1
17#
18# select EM PhysicsList
19#
20/testem/phys/addPhysics emstandard_opt1
21PhysicsList::Set emstandard_opt1 EM physics
22#
23/testem/phys/setCuts 0.3 mm
24/testem/phys/VertexCuts 0.01 mm
25/testem/phys/MuonCuts 10.0 mm
26/testem/stepMax 100 mm
27#
28/testem/det/acceptance1 0.8271 0.02974 2.5
29/testem/det/acceptance9 0.9570 0.01286 2.5
30/testem/det/acceptance25 0.9793 0.00884 2.5
31#
32#/process/eLoss/preciseRange true
33#/process/eLoss/integral true
34/run/initialize
35
36***** Table : Nb of materials = 5 *****
37
38 Material:   G4_AIR     density:  1.205 mg/cm3  RadL: 303.921 m    Imean:  85.700 eV   temperature: 273.15 K  pressure:   1.00 atm
39   --->  Element: C (C)   Z =  6.0   N =  12.0   A =  12.01 g/mole
40   --->  Isotope:     C   Z =  6   N =  12   A =  12.00 g/mole   abundance:  98.93 %
41   --->  Isotope:     C   Z =  6   N =  13   A =  13.00 g/mole   abundance:   1.07 %  ElmMassFraction:   0.01 %  ElmAbundance   0.02 %
42   --->  Element: N (N)   Z =  7.0   N =  14.0   A =  14.01 g/mole
43   --->  Isotope:     N   Z =  7   N =  14   A =  14.00 g/mole   abundance:  99.63 %
44   --->  Isotope:     N   Z =  7   N =  15   A =  15.00 g/mole   abundance:   0.37 %  ElmMassFraction:  75.53 %  ElmAbundance  78.44 %
45   --->  Element: O (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole
46   --->  Isotope:     O   Z =  8   N =  16   A =  15.99 g/mole   abundance:  99.76 %
47   --->  Isotope:     O   Z =  8   N =  17   A =  17.00 g/mole   abundance:   0.04 %
48   --->  Isotope:     O   Z =  8   N =  18   A =  18.00 g/mole   abundance:   0.20 %  ElmMassFraction:  23.18 %  ElmAbundance  21.07 %
49   --->  Element: Ar (Ar)   Z = 18.0   N =  40.0   A =  39.95 g/mole
50   --->  Isotope:    Ar   Z = 18   N =  36   A =  35.97 g/mole   abundance:   0.34 %
51   --->  Isotope:    Ar   Z = 18   N =  38   A =  37.96 g/mole   abundance:   0.06 %
52   --->  Isotope:    Ar   Z = 18   N =  40   A =  39.96 g/mole   abundance:  99.60 %  ElmMassFraction:   1.28 %  ElmAbundance   0.47 %
53
54 Material:    G4_Al     density:  2.699 g/cm3   RadL:   8.896 cm   Imean: 166.000 eV
55   --->  Element: Al (Al)   Z = 13.0   N =  27.0   A =  26.98 g/mole
56   --->  Isotope:    Al   Z = 13   N =  27   A =  26.98 g/mole   abundance: 100.00 %  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
57
58 Material:    G4_Si     density:  2.330 g/cm3   RadL:   9.366 cm   Imean: 173.000 eV
59   --->  Element: Si (Si)   Z = 14.0   N =  28.1   A =  28.09 g/mole
60   --->  Isotope:    Si   Z = 14   N =  28   A =  27.98 g/mole   abundance:  92.23 %
61   --->  Isotope:    Si   Z = 14   N =  29   A =  28.98 g/mole   abundance:   4.68 %
62   --->  Isotope:    Si   Z = 14   N =  30   A =  29.97 g/mole   abundance:   3.09 %  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
63
64 Material:    G4_Fe     density:  7.874 g/cm3   RadL:   1.757 cm   Imean: 286.000 eV
65   --->  Element: Fe (Fe)   Z = 26.0   N =  55.9   A =  55.85 g/mole
66   --->  Isotope:    Fe   Z = 26   N =  54   A =  53.94 g/mole   abundance:   5.84 %
67   --->  Isotope:    Fe   Z = 26   N =  56   A =  55.93 g/mole   abundance:  91.75 %
68   --->  Isotope:    Fe   Z = 26   N =  57   A =  56.94 g/mole   abundance:   2.12 %
69   --->  Isotope:    Fe   Z = 26   N =  58   A =  57.93 g/mole   abundance:   0.28 %  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
70
71 Material: G4_CESIUM_IODIDE     density:  4.510 g/cm3   RadL:   1.860 cm   Imean: 553.100 eV
72   --->  Element: I (I)   Z = 53.0   N = 127.0   A = 126.90 g/mole
73   --->  Isotope:     I   Z = 53   N = 127   A = 126.90 g/mole   abundance: 100.00 %  ElmMassFraction:  48.85 %  ElmAbundance  50.00 %
74   --->  Element: Cs (Cs)   Z = 55.0   N = 133.0   A = 132.91 g/mole
75   --->  Isotope:    Cs   Z = 55   N = 133   A = 132.91 g/mole   abundance: 100.00 %  ElmMassFraction:  51.15 %  ElmAbundance  50.00 %
76
77
78Ecal is 36 cm  of G4_CESIUM_IODIDE
79Absorber is 2 mm  of G4_Al
80Vertex is 3 cm  of 3 layers of Si of 100 um  npads= 601
81### New geometry is constructed
82PhysicsList: world cuts are set cutG= 0.3 mm    cutE= 0.3 mm
83Vertex cuts are set
84Muon cuts are set
85#
86#/testem/histo/fileName e_pr.root
87#/testem/histo/fileType root
88/gun/particle e-
89/gun/energy 1 GeV
90/run/beamOn 100
91
92phot:  Total cross sections from Sandia parametrisation.
93      Sampling according PhotoElectric model
94
95compt:  Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV
96      Sampling according Klein-Nishina model
97      tables are built for  gamma
98      Lambda tables from 100 eV  to 100 GeV in 90 bins.
99
100conv:  Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z;
101      sampling secondary e+e- according Bethe-Heitler model
102      tables are built for  gamma
103      Lambda tables from 1.022 MeV to 100 GeV in 100 bins.
104
105msc:  Model variant of multiple scattering for e-
106      Lambda tables from 100 eV  to 100 TeV in 120 bins.
107      LateralDisplacementFlag=  1   Skin= 0
108      Boundary/stepping algorithm is active with RangeFactor= 0.2  Step limit type 0
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= 0.8, integral: 1, fluct: 1
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
123eIoni:   tables are built for  e+
124      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
125      Lambda tables from threshold to 100 TeV in 120 bins.
126      Delta cross sections and sampling from MollerBhabha model
127      Good description from 1 KeV to 100 GeV.
128      Step function: finalRange(mm)= 1, dRoverRange= 0.8, integral: 1, fluct: 1
129
130eBrem:   tables are built for  e+
131      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
132      Lambda tables from threshold to 100 TeV in 120 bins.
133      Total cross sections and sampling from StandBrem model (based on the EEDL data library)
134      Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
135
136annihil:       Sampling according eplus2gg model
137      tables are built for  e+
138      Lambda tables from 100 eV  to 100 TeV in 120 bins.
139
140msc:  Model variant of multiple scattering for proton
141      Lambda tables from 100 eV  to 100 TeV in 120 bins.
142      LateralDisplacementFlag=  1   Skin= 0
143      Boundary/stepping algorithm is active with RangeFactor= 0.2  Step limit type 0
144
145hIoni:   tables are built for  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= 0.2, integral: 1, fluct: 1
152
153msc:  Model variant of multiple scattering for GenericIon
154      LateralDisplacementFlag=  0   Skin= 0
155      Boundary/stepping algorithm is active with RangeFactor= 0.2  Step limit type 1
156
157ionIoni:   tables are built for  GenericIon
158      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
159      Lambda tables from threshold to 100 TeV in 120 bins.
160      Scaling relation is used from proton dE/dx and range.
161      Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
162      Parametrisation from BraggIon for protons below. NuclearStopping= 1
163
164      Stopping Power data for 8 ion/material pairs are used.
165      Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1
166
167hIoni:   tables are built for  anti_proton
168      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
169      Lambda tables from threshold to 100 TeV in 120 bins.
170      Scaling relation is used from proton dE/dx and range.
171      Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
172      Parametrisation from Bragg for protons below. NuclearStopping= 1
173      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
174
175msc:  Model variant of multiple scattering for mu+
176      Lambda tables from 100 eV  to 100 TeV in 120 bins.
177      LateralDisplacementFlag=  1   Skin= 0
178      Boundary/stepping algorithm is active with RangeFactor= 0.2  Step limit type 0
179
180muIoni:   tables are built for  mu+
181      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
182      Lambda tables from threshold to 100 TeV in 120 bins.
183      Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
184      radiative corrections for E > 1 GeV
185      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
186
187muBrems:   tables are built for  mu+
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      Parametrised model
191
192muPairProd:   tables are built for  mu+
193      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
194      Lambda tables from threshold to 100 TeV in 120 bins.
195      Parametrised model
196
197muIoni:   tables are built for  mu-
198      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
199      Lambda tables from threshold to 100 TeV in 120 bins.
200      Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
201      radiative corrections for E > 1 GeV
202      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
203
204muBrems:   tables are built for  mu-
205      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
206      Lambda tables from threshold to 100 TeV in 120 bins.
207      Parametrised model
208
209muPairProd:   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
214hIoni:   tables are built for  pi+
215      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
216      Lambda tables from threshold to 100 TeV in 120 bins.
217      Scaling relation is used from proton dE/dx and range.
218      Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
219      Parametrisation from Bragg for protons below. NuclearStopping= 1
220      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
221
222msc:  Model variant of multiple scattering for pi-
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.2  Step limit type 0
226
227hIoni:   tables are built for  pi-
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      Scaling relation is used from proton dE/dx and range.
231      Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
232      Parametrisation from Bragg for protons below. NuclearStopping= 1
233      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
234
235========= Table of registered couples ==============================
236
237Index : 0     used in the geometry : Yes     recalculation needed : No
238 Material : G4_AIR
239 Range cuts        :  gamma 300 um     e- 300 um     e+ 300 um
240 Energy thresholds :  gamma 990 eV     e- 990 eV     e+ 990 eV
241 Region(s) which use this couple :
242    DefaultRegionForTheWorld
243
244Index : 1     used in the geometry : Yes     recalculation needed : No
245 Material : G4_CESIUM_IODIDE
246 Range cuts        :  gamma 300 um     e- 300 um     e+ 300 um
247 Energy thresholds :  gamma 18.1187 keV    e- 288.607 keV    e+ 278.143 keV
248 Region(s) which use this couple :
249    DefaultRegionForTheWorld
250
251Index : 2     used in the geometry : Yes     recalculation needed : No
252 Material : G4_Al
253 Range cuts        :  gamma 300 um     e- 300 um     e+ 300 um
254 Energy thresholds :  gamma 3.89739 keV    e- 258.339 keV    e+ 252.056 keV
255 Region(s) which use this couple :
256    DefaultRegionForTheWorld
257
258Index : 3     used in the geometry : Yes     recalculation needed : No
259 Material : G4_AIR
260 Range cuts        :  gamma 10 um     e- 10 um     e+ 10 um
261 Energy thresholds :  gamma 990 eV     e- 990 eV     e+ 990 eV
262 Region(s) which use this couple :
263    VertexDetector
264
265Index : 4     used in the geometry : Yes     recalculation needed : No
266 Material : G4_Al
267 Range cuts        :  gamma 10 um     e- 10 um     e+ 10 um
268 Energy thresholds :  gamma 990 eV     e- 33.8876 keV    e+ 33.4729 keV
269 Region(s) which use this couple :
270    VertexDetector
271
272Index : 5     used in the geometry : Yes     recalculation needed : No
273 Material : G4_Si
274 Range cuts        :  gamma 10 um     e- 10 um     e+ 10 um
275 Energy thresholds :  gamma 990 eV     e- 31.4747 keV    e+ 31.0896 keV
276 Region(s) which use this couple :
277    VertexDetector
278
279Index : 6     used in the geometry : Yes     recalculation needed : No
280 Material : G4_Fe
281 Range cuts        :  gamma 1 cm     e- 1 cm     e+ 1 cm
282 Energy thresholds :  gamma 60.4019 keV    e- 12.9518 MeV    e+ 12.0297 MeV
283 Region(s) which use this couple :
284    MuonDetector
285
286Index : 7     used in the geometry : Yes     recalculation needed : No
287 Material : G4_Al
288 Range cuts        :  gamma 1 cm     e- 1 cm     e+ 1 cm
289 Energy thresholds :  gamma 20.4949 keV    e- 4.54874 MeV    e+ 4.2772 MeV
290 Region(s) which use this couple :
291    MuonDetector
292
293====================================================================
294
295### Run 0 start
296HistoManager: Histograms are booked and run has been started
297Start Run processing.
298Run terminated.
299Run Summary
300  Number of events processed : 100
301  User=4.19s Real=4.21s Sys=0s
302RunAction: End of run actions are started
303HistoManager: End of run actions are started
304=================================================================
305Number of events             100
306Average number of e-         1320
307Average number of gamma      701.3
308Average number of e+         48.12
309Average number of steps      3980
310Edep 1x1 =                   0.8336 +- 0.002518  res=  3.021 %
311Edep 3x3 =                   0.9601 +- 0.001163  res=  1.212 %
312Edep 5x5 =                   0.9807 +- 0.0007853  res=  0.8007 %
313===========  Mean values without trancating =====================
314Edep 1x1 =                   0.8327 +- 0.002681  res=  3.219 %
315Edep 3x3 =                   0.9596 +- 0.001241  res=  1.293 %
316Edep 5x5 =                   0.9807 +- 0.0007853  res=  0.8007 %
317==================================================================
318
319
320<<<<<ACCEPTANCE>>>>> 100 events for Crystal Calorimeter
321Edep1x1: 0.8336  delEdep1x1= 0.006515 nrms= 2.191
322Erms1x1: 0.02518  delErms1x1= -0.00456 nrms= -1.533
323Edep3x3: 0.9601  delEdep3x3= 0.003087 nrms= 2.401
324Erms3x3: 0.01163  delErms3x3= -0.001226 nrms= -0.9531
325Edep5x5: 0.9807  delEdep5x5= 0.001428 nrms= 1.615
326Erms5x5: 0.007853  delErms5x5= -0.0009872 nrms= -1.117
327<<<<<END>>>>>   IS ACCEPTED
328
329#
330G4 kernel has come to Quit state.
Note: See TracBrowser for help on using the repository browser.