[807] | 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... |
---|