1 | |
---|
2 | ************************************************************* |
---|
3 | Geant4 version Name: geant4-08-01-ref-00 (30-June-2006) |
---|
4 | Copyright : Geant4 Collaboration |
---|
5 | Reference : NIM A 506 (2003), 250-303 |
---|
6 | WWW : http://cern.ch/geant4 |
---|
7 | ************************************************************* |
---|
8 | |
---|
9 | |
---|
10 | <<< Geant4 Physics List engine packaging library: PACK 5.1 |
---|
11 | <<< Geant4 Physics List simulation engine: QGSP 3.1 |
---|
12 | |
---|
13 | Visualization Manager instantiating... |
---|
14 | Visualization Manager initialising... |
---|
15 | Registering graphics systems... |
---|
16 | |
---|
17 | You have successfully registered the following graphics systems. |
---|
18 | Current available graphics systems are: |
---|
19 | ASCIITree (ATree) |
---|
20 | DAWNFILE (DAWNFILE) |
---|
21 | GAGTree (GAGTree) |
---|
22 | G4HepRep (HepRepXML) |
---|
23 | G4HepRepFile (HepRepFile) |
---|
24 | RayTracer (RayTracer) |
---|
25 | VRML1FILE (VRML1FILE) |
---|
26 | VRML2FILE (VRML2FILE) |
---|
27 | FukuiRenderer (DAWN) |
---|
28 | OpenGLImmediateX (OGLIX) |
---|
29 | OpenGLStoredX (OGLSX) |
---|
30 | OpenGLImmediateXm (OGLIXm) |
---|
31 | OpenGLStoredXm (OGLSXm) |
---|
32 | OpenInventorXt (OIX) |
---|
33 | RayTracerX (RayTracerX) |
---|
34 | VRML1 (VRML1) |
---|
35 | VRML2 (VRML2) |
---|
36 | |
---|
37 | Registering model factories... |
---|
38 | |
---|
39 | You have successfully registered the following model factories. |
---|
40 | Registered model factories: |
---|
41 | generic |
---|
42 | drawByCharge |
---|
43 | drawByParticleID |
---|
44 | drawByOriginVolume |
---|
45 | |
---|
46 | Registered models: |
---|
47 | None |
---|
48 | |
---|
49 | Registered filter factories: |
---|
50 | chargeFilter |
---|
51 | particleFilter |
---|
52 | originVolumeFilter |
---|
53 | |
---|
54 | Registered filters: |
---|
55 | None |
---|
56 | |
---|
57 | current generator type: pythia |
---|
58 | 1 |
---|
59 | ****************************************************************************** |
---|
60 | ****************************************************************************** |
---|
61 | ** ** |
---|
62 | ** ** |
---|
63 | ** *......* Welcome to the Lund Monte Carlo! ** |
---|
64 | ** *:::!!:::::::::::* ** |
---|
65 | ** *::::::!!::::::::::::::* PPP Y Y TTTTT H H III A ** |
---|
66 | ** *::::::::!!::::::::::::::::* P P Y Y T H H I A A ** |
---|
67 | ** *:::::::::!!:::::::::::::::::* PPP Y T HHHHH I AAAAA ** |
---|
68 | ** *:::::::::!!:::::::::::::::::* P Y T H H I A A ** |
---|
69 | ** *::::::::!!::::::::::::::::*! P Y T H H III A A ** |
---|
70 | ** *::::::!!::::::::::::::* !! ** |
---|
71 | ** !! *:::!!:::::::::::* !! This is PYTHIA version 6.205 ** |
---|
72 | ** !! !* -><- * !! Last date of change: 1 Mar 2002 ** |
---|
73 | ** !! !! !! ** |
---|
74 | ** !! !! !! Now is 0 Jan 2000 at 0:00:00 ** |
---|
75 | ** !! !! ** |
---|
76 | ** !! lh !! Disclaimer: this program comes ** |
---|
77 | ** !! !! without any guarantees. Beware ** |
---|
78 | ** !! hh !! of errors and use common sense ** |
---|
79 | ** !! ll !! when interpreting results. ** |
---|
80 | ** !! !! ** |
---|
81 | ** !! Copyright T. Sjostrand (2001) ** |
---|
82 | ** ** |
---|
83 | ** An archive of program versions and documentation is found on the web: ** |
---|
84 | ** http://www.thep.lu.se/~torbjorn/Pythia.html ** |
---|
85 | ** ** |
---|
86 | ** When you cite this program, currently the official reference is ** |
---|
87 | ** T. Sjostrand, P. Eden, C. Friberg, L. Lonnblad, G. Miu, S. Mrenna and ** |
---|
88 | ** E. Norrbin, Computer Physics Commun. 135 (2001) 238. ** |
---|
89 | ** The large manual is ** |
---|
90 | ** T. Sjostrand, L. Lonnblad and S. Mrenna, LU TP 01-21 [hep-ph/0108264]. ** |
---|
91 | ** Also remember that the program, to a large extent, represents original ** |
---|
92 | ** physics research. Other publications of special relevance to your ** |
---|
93 | ** studies may therefore deserve separate mention. ** |
---|
94 | ** ** |
---|
95 | ** Main author: Torbjorn Sjostrand; Department of Theoretical Physics 2, ** |
---|
96 | ** Lund University, Solvegatan 14A, S-223 62 Lund, Sweden; ** |
---|
97 | ** phone: + 46 - 46 - 222 48 16; e-mail: torbjorn@thep.lu.se ** |
---|
98 | ** Author: Leif Lonnblad; Department of Theoretical Physics 2, ** |
---|
99 | ** Lund University, Solvegatan 14A, S-223 62 Lund, Sweden; ** |
---|
100 | ** phone: + 46 - 46 - 222 77 80; e-mail: leif@thep.lu.se ** |
---|
101 | ** Author: Stephen Mrenna; Computing Division, Simulations Group, ** |
---|
102 | ** Fermi National Accelerator Laboratory, MS 234, Batavia, IL 60510, USA; ** |
---|
103 | ** phone: + 1 - 630 - 840 - 2556; e-mail: mrenna@fnal.gov ** |
---|
104 | ** Author: Peter Skands; Department of Theoretical Physics 2, ** |
---|
105 | ** Lund University, Solvegatan 14A, S-223 62 Lund, Sweden; ** |
---|
106 | ** phone: + 46 - 46 - 222 31 92; e-mail: zeiler@thep.lu.se ** |
---|
107 | ** ** |
---|
108 | ** ** |
---|
109 | ****************************************************************************** |
---|
110 | ****************************************************************************** |
---|
111 | PMAS(6,1) changed from 175.00000 to 174.30000 |
---|
112 | PMAS(25,1) changed from 115.00000 to 200.00000 |
---|
113 | MSEL changed from 1 to 0 |
---|
114 | MSUB(102) changed from 0 to 1 |
---|
115 | MSUB(123) changed from 0 to 1 |
---|
116 | MSUB(124) changed from 0 to 1 |
---|
117 | MDME(210,1) changed from 1 to 0 |
---|
118 | MDME(211,1) changed from 1 to 0 |
---|
119 | MDME(212,1) changed from 1 to 0 |
---|
120 | MDME(213,1) changed from 1 to 0 |
---|
121 | MDME(214,1) changed from 1 to 0 |
---|
122 | MDME(215,1) changed from 1 to 0 |
---|
123 | MDME(218,1) changed from 1 to 0 |
---|
124 | MDME(219,1) changed from 1 to 0 |
---|
125 | MDME(220,1) changed from 1 to 0 |
---|
126 | MDME(222,1) changed from 1 to 0 |
---|
127 | MDME(223,1) changed from 1 to 0 |
---|
128 | MDME(224,1) changed from 1 to 0 |
---|
129 | MDME(225,1) changed from 1 to 1 |
---|
130 | MDME(226,1) changed from 1 to 0 |
---|
131 | MDME(216,1) changed from -1 to -1 |
---|
132 | MDME(217,1) changed from -1 to -1 |
---|
133 | MDME(218,1) changed from 0 to -1 |
---|
134 | MDME(174,1) changed from 1 to 0 |
---|
135 | MDME(175,1) changed from 1 to 0 |
---|
136 | MDME(176,1) changed from 1 to 0 |
---|
137 | MDME(177,1) changed from 1 to 0 |
---|
138 | MDME(178,1) changed from 1 to 0 |
---|
139 | MDME(179,1) changed from 1 to 0 |
---|
140 | MDME(182,1) changed from 1 to 0 |
---|
141 | MDME(183,1) changed from 1 to 0 |
---|
142 | MDME(184,1) changed from 1 to 1 |
---|
143 | MDME(185,1) changed from 1 to 0 |
---|
144 | MDME(186,1) changed from 1 to 0 |
---|
145 | MDME(187,1) changed from 1 to 0 |
---|
146 | MDME(180,1) changed from -1 to -1 |
---|
147 | MDME(181,1) changed from -1 to -1 |
---|
148 | MDME(188,1) changed from -1 to -1 |
---|
149 | MDME(189,1) changed from -1 to -1 |
---|
150 | MSTP(61) changed from 2 to 1 |
---|
151 | MSTP(71) changed from 1 to 1 |
---|
152 | MSTP(81) changed from 1 to 1 |
---|
153 | MSTP(111) changed from 1 to 1 |
---|
154 | 1****************** PYINIT: initialization of PYTHIA routines ***************** |
---|
155 | |
---|
156 | ============================================================================== |
---|
157 | I I |
---|
158 | I PYTHIA will be initialized for a p on p collider I |
---|
159 | I at 14000.000 GeV center-of-mass energy I |
---|
160 | I I |
---|
161 | ============================================================================== |
---|
162 | |
---|
163 | ******** PYMAXI: summary of differential cross-section maximum search ******** |
---|
164 | |
---|
165 | ========================================================== |
---|
166 | I I I |
---|
167 | I ISUB Subprocess name I Maximum value I |
---|
168 | I I I |
---|
169 | ========================================================== |
---|
170 | I I I |
---|
171 | I 96 Semihard QCD 2 -> 2 I 1.8679E+03 I |
---|
172 | I 102 g + g -> h0 I 5.1617E-12 I |
---|
173 | I 123 f + f' -> f + f' + h0 I 4.0000E-12 I |
---|
174 | I 124 f + f' -> f" + f"' + h0 I 8.9419E-12 I |
---|
175 | I I I |
---|
176 | ========================================================== |
---|
177 | |
---|
178 | ********************** PYINIT: initialization completed ********************** |
---|
179 | |
---|
180 | phot: Total cross sections from Sandia parametrisation. |
---|
181 | |
---|
182 | compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV |
---|
183 | Sampling according Klein-Nishina model |
---|
184 | tables are built for gamma |
---|
185 | Lambda tables from 100 eV to 100 GeV in 90 bins. |
---|
186 | |
---|
187 | conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; |
---|
188 | sampling secondary e+e- according to the Bethe-Heitler model |
---|
189 | tables are built for gamma |
---|
190 | Lambda tables from 1.022 MeV to 100 GeV in 100 bins. |
---|
191 | |
---|
192 | msc: Model variant of multiple scattering for e- |
---|
193 | Lambda tables from 100 eV to 100 TeV in 120 bins. |
---|
194 | Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 |
---|
195 | |
---|
196 | eIoni: tables are built for e- |
---|
197 | dE/dx and range tables from 100 eV to 100 TeV in 120 bins. |
---|
198 | Lambda tables from threshold to 100 TeV in 120 bins. |
---|
199 | Delta cross sections from Moller+Bhabha, good description from 1 KeV to 100 GeV. |
---|
200 | Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 |
---|
201 | |
---|
202 | eBrem: tables are built for e- |
---|
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 | Total cross sections from a parametrisation based on the EEDL data library. |
---|
206 | Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. |
---|
207 | |
---|
208 | eIoni: tables are built for e+ |
---|
209 | dE/dx and range tables from 100 eV to 100 TeV in 120 bins. |
---|
210 | Lambda tables from threshold to 100 TeV in 120 bins. |
---|
211 | Delta cross sections from Moller+Bhabha, good description from 1 KeV to 100 GeV. |
---|
212 | Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 |
---|
213 | |
---|
214 | eBrem: tables are built for e+ |
---|
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 | Total cross sections from a parametrisation based on the EEDL data library. |
---|
218 | Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. |
---|
219 | |
---|
220 | annihil: Heilter model of formula of annihilation into 2 photons |
---|
221 | tables are built for e+ |
---|
222 | Lambda tables from 100 eV to 100 TeV in 120 bins. |
---|
223 | |
---|
224 | msc: Model variant of multiple scattering for proton |
---|
225 | Lambda tables from 100 eV to 100 TeV in 120 bins. |
---|
226 | Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 |
---|
227 | |
---|
228 | hIoni: tables are built for proton |
---|
229 | dE/dx and range tables from 100 eV to 100 TeV in 120 bins. |
---|
230 | Lambda tables from threshold to 100 TeV in 120 bins. |
---|
231 | Scaling relation is used to proton dE/dx and range |
---|
232 | Bether-Bloch model for Escaled > 2 MeV, ICRU49 parametrisation for protons below. |
---|
233 | Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 |
---|
234 | |
---|
235 | msc: Model variant of multiple scattering for GenericIon |
---|
236 | Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 |
---|
237 | |
---|
238 | ionIoni: tables are built for GenericIon |
---|
239 | dE/dx and range tables from 100 eV to 100 TeV in 120 bins. |
---|
240 | Lambda tables from threshold to 100 TeV in 120 bins. |
---|
241 | Scaling relation is used to proton dE/dx and range |
---|
242 | Bether-Bloch model for Escaled > 2 MeV, ICRU49 parametrisation for alpha particles below. |
---|
243 | Stopping Power data for 8 ion/material pairs are used. |
---|
244 | Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1 |
---|
245 | |
---|
246 | hIoni: tables are built for anti_proton |
---|
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 to proton dE/dx and range |
---|
250 | Bether-Bloch model for Escaled > 2 MeV, ICRU49 parametrisation for protons below. |
---|
251 | Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 |
---|
252 | |
---|
253 | msc: Model variant of multiple scattering for mu+ |
---|
254 | Lambda tables from 100 eV to 100 TeV in 120 bins. |
---|
255 | Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 |
---|
256 | |
---|
257 | muIoni: tables are built for mu+ |
---|
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 | Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, |
---|
261 | radiative corrections for E > 1 GeV |
---|
262 | Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 |
---|
263 | |
---|
264 | muBrems: tables are built for mu+ |
---|
265 | dE/dx and range tables from 100 eV to 100 TeV in 120 bins. |
---|
266 | Lambda tables from threshold to 100 TeV in 120 bins. |
---|
267 | Parametrised model |
---|
268 | |
---|
269 | muPairProd: tables are built for mu+ |
---|
270 | dE/dx and range tables from 100 eV to 100 TeV in 120 bins. |
---|
271 | Lambda tables from threshold to 100 TeV in 120 bins. |
---|
272 | Parametrised model |
---|
273 | |
---|
274 | muIoni: tables are built for mu- |
---|
275 | dE/dx and range tables from 100 eV to 100 TeV in 120 bins. |
---|
276 | Lambda tables from threshold to 100 TeV in 120 bins. |
---|
277 | Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, |
---|
278 | radiative corrections for E > 1 GeV |
---|
279 | Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 |
---|
280 | |
---|
281 | muBrems: tables are built for mu- |
---|
282 | dE/dx and range tables from 100 eV to 100 TeV in 120 bins. |
---|
283 | Lambda tables from threshold to 100 TeV in 120 bins. |
---|
284 | Parametrised model |
---|
285 | |
---|
286 | muPairProd: tables are built for mu- |
---|
287 | dE/dx and range tables from 100 eV to 100 TeV in 120 bins. |
---|
288 | Lambda tables from threshold to 100 TeV in 120 bins. |
---|
289 | Parametrised model |
---|
290 | |
---|
291 | hIoni: tables are built for pi+ |
---|
292 | dE/dx and range tables from 100 eV to 100 TeV in 120 bins. |
---|
293 | Lambda tables from threshold to 100 TeV in 120 bins. |
---|
294 | Scaling relation is used to proton dE/dx and range |
---|
295 | Bether-Bloch model for Escaled > 0.297504 MeV, ICRU49 parametrisation for protons below. |
---|
296 | Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 |
---|
297 | |
---|
298 | msc: Model variant of multiple scattering for pi- |
---|
299 | Lambda tables from 100 eV to 100 TeV in 120 bins. |
---|
300 | Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 |
---|
301 | |
---|
302 | hIoni: tables are built for pi- |
---|
303 | dE/dx and range tables from 100 eV to 100 TeV in 120 bins. |
---|
304 | Lambda tables from threshold to 100 TeV in 120 bins. |
---|
305 | Scaling relation is used to proton dE/dx and range |
---|
306 | Bether-Bloch model for Escaled > 0.297504 MeV, ICRU49 parametrisation for protons below. |
---|
307 | Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 |
---|
308 | ### Run 0 start. |
---|
309 | Stage 0->1 : 1 hits found in the muon chamber. |
---|
310 | Stage 1->2 : 0 isolated muon found. |
---|
311 | ++++++++ event aborted |
---|
312 | >>> Event 0 |
---|
313 | 213 hits are stored in ExN04TrackerHitsCollection. |
---|
314 | 21 hits are stored in ExN04CalorimeterHitsCollection. |
---|
315 | Total energy deposition in calorimeter : 0.49951823 (GeV) |
---|
316 | 1 hits are stored in ExN04MuonHitsCollection. |
---|
317 | Stage 0->1 : 0 hits found in the muon chamber. |
---|
318 | ++++++++ event aborted |
---|
319 | >>> Event 1 |
---|
320 | 0 hits are stored in ExN04TrackerHitsCollection. |
---|
321 | 14 hits are stored in ExN04CalorimeterHitsCollection. |
---|
322 | Total energy deposition in calorimeter : 0.20789644 (GeV) |
---|
323 | 0 hits are stored in ExN04MuonHitsCollection. |
---|
324 | Stage 0->1 : 1 hits found in the muon chamber. |
---|
325 | Stage 1->2 : 1 isolated muon found. |
---|
326 | >>> Event 2 |
---|
327 | 129 hits are stored in ExN04TrackerHitsCollection. |
---|
328 | 45 hits are stored in ExN04CalorimeterHitsCollection. |
---|
329 | Total energy deposition in calorimeter : 1.808844 (GeV) |
---|
330 | 1 hits are stored in ExN04MuonHitsCollection. |
---|
331 | Stage 0->1 : 1 hits found in the muon chamber. |
---|
332 | Stage 1->2 : 0 isolated muon found. |
---|
333 | ++++++++ event aborted |
---|
334 | >>> Event 3 |
---|
335 | 166 hits are stored in ExN04TrackerHitsCollection. |
---|
336 | 21 hits are stored in ExN04CalorimeterHitsCollection. |
---|
337 | Total energy deposition in calorimeter : 0.45672857 (GeV) |
---|
338 | 1 hits are stored in ExN04MuonHitsCollection. |
---|
339 | Stage 0->1 : 1 hits found in the muon chamber. |
---|
340 | Stage 1->2 : 1 isolated muon found. |
---|
341 | >>> Event 4 |
---|
342 | 177 hits are stored in ExN04TrackerHitsCollection. |
---|
343 | 35 hits are stored in ExN04CalorimeterHitsCollection. |
---|
344 | Total energy deposition in calorimeter : 0.78347475 (GeV) |
---|
345 | 1 hits are stored in ExN04MuonHitsCollection. |
---|
346 | 1********* PYSTAT: Statistics on Number of Events and Cross-sections ********* |
---|
347 | |
---|
348 | ============================================================================== |
---|
349 | I I I I |
---|
350 | I Subprocess I Number of points I Sigma I |
---|
351 | I I I I |
---|
352 | I----------------------------------I----------------------------I (mb) I |
---|
353 | I I I I |
---|
354 | I N:o Type I Generated Tried I I |
---|
355 | I I I I |
---|
356 | ============================================================================== |
---|
357 | I I I I |
---|
358 | I 0 All included subprocesses I 5 10 I 6.034E-12 I |
---|
359 | I 102 g + g -> h0 I 2 4 I 2.874E-12 I |
---|
360 | I 123 f + f' -> f + f' + h0 I 1 1 I 6.634E-13 I |
---|
361 | I 124 f + f' -> f" + f"' + h0 I 2 5 I 2.497E-12 I |
---|
362 | I I I I |
---|
363 | ============================================================================== |
---|
364 | |
---|
365 | ********* Fraction of events that fail fragmentation cuts = 0.00000 ********* |
---|
366 | |
---|
367 | Graphics systems deleted. |
---|
368 | Visualization Manager deleting... |
---|