************************************************************* Geant4 version Name: geant4-08-01-ref-00 (30-June-2006) Copyright : Geant4 Collaboration Reference : NIM A 506 (2003), 250-303 WWW : http://cern.ch/geant4 ************************************************************* <<< Geant4 Physics List engine packaging library: PACK 5.1 <<< Geant4 Physics List simulation engine: QGSP 3.1 Visualization Manager instantiating... Visualization Manager initialising... Registering graphics systems... You have successfully registered the following graphics systems. Current available graphics systems are: ASCIITree (ATree) DAWNFILE (DAWNFILE) GAGTree (GAGTree) G4HepRep (HepRepXML) G4HepRepFile (HepRepFile) RayTracer (RayTracer) VRML1FILE (VRML1FILE) VRML2FILE (VRML2FILE) FukuiRenderer (DAWN) OpenGLImmediateX (OGLIX) OpenGLStoredX (OGLSX) OpenGLImmediateXm (OGLIXm) OpenGLStoredXm (OGLSXm) OpenInventorXt (OIX) RayTracerX (RayTracerX) VRML1 (VRML1) VRML2 (VRML2) Registering model factories... You have successfully registered the following model factories. Registered model factories: generic drawByCharge drawByParticleID drawByOriginVolume Registered models: None Registered filter factories: chargeFilter particleFilter originVolumeFilter Registered filters: None current generator type: pythia 1 ****************************************************************************** ****************************************************************************** ** ** ** ** ** *......* Welcome to the Lund Monte Carlo! ** ** *:::!!:::::::::::* ** ** *::::::!!::::::::::::::* PPP Y Y TTTTT H H III A ** ** *::::::::!!::::::::::::::::* P P Y Y T H H I A A ** ** *:::::::::!!:::::::::::::::::* PPP Y T HHHHH I AAAAA ** ** *:::::::::!!:::::::::::::::::* P Y T H H I A A ** ** *::::::::!!::::::::::::::::*! P Y T H H III A A ** ** *::::::!!::::::::::::::* !! ** ** !! *:::!!:::::::::::* !! This is PYTHIA version 6.205 ** ** !! !* -><- * !! Last date of change: 1 Mar 2002 ** ** !! !! !! ** ** !! !! !! Now is 0 Jan 2000 at 0:00:00 ** ** !! !! ** ** !! lh !! Disclaimer: this program comes ** ** !! !! without any guarantees. Beware ** ** !! hh !! of errors and use common sense ** ** !! ll !! when interpreting results. ** ** !! !! ** ** !! Copyright T. Sjostrand (2001) ** ** ** ** An archive of program versions and documentation is found on the web: ** ** http://www.thep.lu.se/~torbjorn/Pythia.html ** ** ** ** When you cite this program, currently the official reference is ** ** T. Sjostrand, P. Eden, C. Friberg, L. Lonnblad, G. Miu, S. Mrenna and ** ** E. Norrbin, Computer Physics Commun. 135 (2001) 238. ** ** The large manual is ** ** T. Sjostrand, L. Lonnblad and S. Mrenna, LU TP 01-21 [hep-ph/0108264]. ** ** Also remember that the program, to a large extent, represents original ** ** physics research. Other publications of special relevance to your ** ** studies may therefore deserve separate mention. ** ** ** ** Main author: Torbjorn Sjostrand; Department of Theoretical Physics 2, ** ** Lund University, Solvegatan 14A, S-223 62 Lund, Sweden; ** ** phone: + 46 - 46 - 222 48 16; e-mail: torbjorn@thep.lu.se ** ** Author: Leif Lonnblad; Department of Theoretical Physics 2, ** ** Lund University, Solvegatan 14A, S-223 62 Lund, Sweden; ** ** phone: + 46 - 46 - 222 77 80; e-mail: leif@thep.lu.se ** ** Author: Stephen Mrenna; Computing Division, Simulations Group, ** ** Fermi National Accelerator Laboratory, MS 234, Batavia, IL 60510, USA; ** ** phone: + 1 - 630 - 840 - 2556; e-mail: mrenna@fnal.gov ** ** Author: Peter Skands; Department of Theoretical Physics 2, ** ** Lund University, Solvegatan 14A, S-223 62 Lund, Sweden; ** ** phone: + 46 - 46 - 222 31 92; e-mail: zeiler@thep.lu.se ** ** ** ** ** ****************************************************************************** ****************************************************************************** PMAS(6,1) changed from 175.00000 to 174.30000 PMAS(25,1) changed from 115.00000 to 200.00000 MSEL changed from 1 to 0 MSUB(102) changed from 0 to 1 MSUB(123) changed from 0 to 1 MSUB(124) changed from 0 to 1 MDME(210,1) changed from 1 to 0 MDME(211,1) changed from 1 to 0 MDME(212,1) changed from 1 to 0 MDME(213,1) changed from 1 to 0 MDME(214,1) changed from 1 to 0 MDME(215,1) changed from 1 to 0 MDME(218,1) changed from 1 to 0 MDME(219,1) changed from 1 to 0 MDME(220,1) changed from 1 to 0 MDME(222,1) changed from 1 to 0 MDME(223,1) changed from 1 to 0 MDME(224,1) changed from 1 to 0 MDME(225,1) changed from 1 to 1 MDME(226,1) changed from 1 to 0 MDME(216,1) changed from -1 to -1 MDME(217,1) changed from -1 to -1 MDME(218,1) changed from 0 to -1 MDME(174,1) changed from 1 to 0 MDME(175,1) changed from 1 to 0 MDME(176,1) changed from 1 to 0 MDME(177,1) changed from 1 to 0 MDME(178,1) changed from 1 to 0 MDME(179,1) changed from 1 to 0 MDME(182,1) changed from 1 to 0 MDME(183,1) changed from 1 to 0 MDME(184,1) changed from 1 to 1 MDME(185,1) changed from 1 to 0 MDME(186,1) changed from 1 to 0 MDME(187,1) changed from 1 to 0 MDME(180,1) changed from -1 to -1 MDME(181,1) changed from -1 to -1 MDME(188,1) changed from -1 to -1 MDME(189,1) changed from -1 to -1 MSTP(61) changed from 2 to 1 MSTP(71) changed from 1 to 1 MSTP(81) changed from 1 to 1 MSTP(111) changed from 1 to 1 1****************** PYINIT: initialization of PYTHIA routines ***************** ============================================================================== I I I PYTHIA will be initialized for a p on p collider I I at 14000.000 GeV center-of-mass energy I I I ============================================================================== ******** PYMAXI: summary of differential cross-section maximum search ******** ========================================================== I I I I ISUB Subprocess name I Maximum value I I I I ========================================================== I I I I 96 Semihard QCD 2 -> 2 I 1.8679E+03 I I 102 g + g -> h0 I 5.1617E-12 I I 123 f + f' -> f + f' + h0 I 4.0000E-12 I I 124 f + f' -> f" + f"' + h0 I 8.9419E-12 I I I I ========================================================== ********************** PYINIT: initialization completed ********************** phot: Total cross sections from Sandia parametrisation. compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV Sampling according Klein-Nishina model tables are built for gamma Lambda tables from 100 eV to 100 GeV in 90 bins. conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; sampling secondary e+e- according to the Bethe-Heitler model tables are built for gamma Lambda tables from 1.022 MeV to 100 GeV in 100 bins. msc: Model variant of multiple scattering for e- Lambda tables from 100 eV to 100 TeV in 120 bins. Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 eIoni: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections from Moller+Bhabha, good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 eBrem: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections from a parametrisation based on the EEDL data library. Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. eIoni: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections from Moller+Bhabha, good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 eBrem: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections from a parametrisation based on the EEDL data library. Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. annihil: Heilter model of formula of annihilation into 2 photons tables are built for e+ Lambda tables from 100 eV to 100 TeV in 120 bins. msc: Model variant of multiple scattering for proton Lambda tables from 100 eV to 100 TeV in 120 bins. Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 hIoni: tables are built for proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used to proton dE/dx and range Bether-Bloch model for Escaled > 2 MeV, ICRU49 parametrisation for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 msc: Model variant of multiple scattering for GenericIon Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 ionIoni: tables are built for GenericIon dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used to proton dE/dx and range Bether-Bloch model for Escaled > 2 MeV, ICRU49 parametrisation for alpha particles below. Stopping Power data for 8 ion/material pairs are used. Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1 hIoni: tables are built for anti_proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used to proton dE/dx and range Bether-Bloch model for Escaled > 2 MeV, ICRU49 parametrisation for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 msc: Model variant of multiple scattering for mu+ Lambda tables from 100 eV to 100 TeV in 120 bins. Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 muIoni: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 muBrems: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model muPairProd: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model muIoni: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 muBrems: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model muPairProd: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model hIoni: tables are built for pi+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used to proton dE/dx and range Bether-Bloch model for Escaled > 0.297504 MeV, ICRU49 parametrisation for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 msc: Model variant of multiple scattering for pi- Lambda tables from 100 eV to 100 TeV in 120 bins. Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 hIoni: tables are built for pi- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used to proton dE/dx and range Bether-Bloch model for Escaled > 0.297504 MeV, ICRU49 parametrisation for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 ### Run 0 start. Stage 0->1 : 1 hits found in the muon chamber. Stage 1->2 : 0 isolated muon found. ++++++++ event aborted >>> Event 0 213 hits are stored in ExN04TrackerHitsCollection. 21 hits are stored in ExN04CalorimeterHitsCollection. Total energy deposition in calorimeter : 0.49951823 (GeV) 1 hits are stored in ExN04MuonHitsCollection. Stage 0->1 : 0 hits found in the muon chamber. ++++++++ event aborted >>> Event 1 0 hits are stored in ExN04TrackerHitsCollection. 14 hits are stored in ExN04CalorimeterHitsCollection. Total energy deposition in calorimeter : 0.20789644 (GeV) 0 hits are stored in ExN04MuonHitsCollection. Stage 0->1 : 1 hits found in the muon chamber. Stage 1->2 : 1 isolated muon found. >>> Event 2 129 hits are stored in ExN04TrackerHitsCollection. 45 hits are stored in ExN04CalorimeterHitsCollection. Total energy deposition in calorimeter : 1.808844 (GeV) 1 hits are stored in ExN04MuonHitsCollection. Stage 0->1 : 1 hits found in the muon chamber. Stage 1->2 : 0 isolated muon found. ++++++++ event aborted >>> Event 3 166 hits are stored in ExN04TrackerHitsCollection. 21 hits are stored in ExN04CalorimeterHitsCollection. Total energy deposition in calorimeter : 0.45672857 (GeV) 1 hits are stored in ExN04MuonHitsCollection. Stage 0->1 : 1 hits found in the muon chamber. Stage 1->2 : 1 isolated muon found. >>> Event 4 177 hits are stored in ExN04TrackerHitsCollection. 35 hits are stored in ExN04CalorimeterHitsCollection. Total energy deposition in calorimeter : 0.78347475 (GeV) 1 hits are stored in ExN04MuonHitsCollection. 1********* PYSTAT: Statistics on Number of Events and Cross-sections ********* ============================================================================== I I I I I Subprocess I Number of points I Sigma I I I I I I----------------------------------I----------------------------I (mb) I I I I I I N:o Type I Generated Tried I I I I I I ============================================================================== I I I I I 0 All included subprocesses I 5 10 I 6.034E-12 I I 102 g + g -> h0 I 2 4 I 2.874E-12 I I 123 f + f' -> f + f' + h0 I 1 1 I 6.634E-13 I I 124 f + f' -> f" + f"' + h0 I 2 5 I 2.497E-12 I I I I I ============================================================================== ********* Fraction of events that fail fragmentation cuts = 0.00000 ********* Graphics systems deleted. Visualization Manager deleting...