Arxiv:0803.0883V3 [Hep-Ph] 2 Dec 2008 Herwig++ Physics and Manual
Total Page:16
File Type:pdf, Size:1020Kb
CERN-PH-TH/2008-038 Cavendish-HEP-08/03 KA-TP-05-2008 DCPT/08/22 IPPP/08/11 CP3-08-05 Herwig++ Physics and Manual M. B¨ahr1, S. Gieseke1, M.A. Gigg2, D. Grellscheid2, K. Hamilton3, O. Latunde- Dada4, S. Pl¨atzer1, P. Richardson2, M.H. Seymour5,6, A. Sherstnev4, J. Tully2, B.R. Webber4 1 Institut f¨ur Theoretische Physik, Universit¨at Karlsruhe. 2 Department of Physics, Durham University. 3Centre for Particle Physics and Phenomenology, Universit´e Catholique de Louvain. 4 Cavendish Laboratory, University of Cambridge. 5 School of Physics and Astronomy, University of Manchester. 6 Physics Department, CERN. Authors’ E-mail: [email protected] Abstract In this paper we describe Herwig++ version 2.3, a general-purpose Monte Carlo event generator for the simulation of hard lepton-lepton, lepton-hadron and hadron-hadron col- lisions. A number of important hard scattering processes are available, together with an interface via the Les Houches Accord to specialized matrix element generators for addi- tional processes. The simulation of Beyond the Standard Model (BSM) physics includes a range of models and allows new models to be added by encoding the Feynman rules of the model. The parton-shower approach is used to simulate initial- and final-state QCD radia- tion, including colour coherence effects, with special emphasis on the correct description of arXiv:0803.0883v3 [hep-ph] 2 Dec 2008 radiation from heavy particles. The underlying event is simulated using an eikonal multiple parton-parton scattering model. The formation of hadrons from the quarks and gluons produced in the parton shower is described using the cluster hadronization model. Hadron decays are simulated using matrix elements, where possible including spin correlations and off-shell effects. Contents 1 Introduction 4 2 Technical Details 7 3 Matrix Elements 10 3.1 Leadingordermatrixelements. ....... 10 3.2 Next-to-leading-order matrix elements . ........... 12 3.3 LesHouchesinterface. .. .. .. ..... 13 3.4 Processes with incoming photons . ....... 13 3.5 Codestructure ................................... 14 3.5.1 StandardEventHandler . 14 3.5.2 LesHouchesEventHandler . .... 15 3.5.3 Kinematiccuts ................................. 15 4 Perturbative Decays and Spin Correlations 17 4.1 Spincorrelations ................................ .... 17 4.2 StandardModeldecays. .. .. .. .... 20 4.3 QEDradiation .................................... 21 4.4 Codestructure ................................... 21 5 Physics Beyond The Standard Model 22 5.1 Hardprocess ..................................... 22 5.2 Decays.......................................... 23 5.3 Off-ShellEffects.................................. 24 5.4 Modeldescriptions ............................... .... 25 5.4.1 StandardModel................................. 25 5.4.2 Minimal Supersymmetric Standard Model . ...... 26 5.4.3 Randall-SundrumModel . 26 5.4.4 Minimal Universal Extra Dimensions Model . ....... 27 5.5 Codestructure ................................... 27 6 Parton Showers 33 6.1 Showerkinematics ................................ 33 6.2 Showerdynamics .................................. 37 6.3 Initialconditions ............................... ..... 40 6.3.1 Final-finalcolourconnection . ...... 41 6.3.2 Initial-initial colour connection . ......... 42 6.3.3 Initial-final colour connection in the hard process . ............ 42 6.3.4 Initial-final colour connection in decays . .......... 43 6.4 Final-stateradiation . ...... 44 6.4.1 Evolution .................................... 44 6.4.2 Kinematicreconstruction. ..... 46 6.5 Initial-stateradiation . ....... 47 6.5.1 Evolution .................................... 47 1 6.5.2 Kinematicreconstruction. ..... 50 6.5.3 Forcedsplitting ............................... 53 6.6 Radiationinparticledecays . ....... 54 6.6.1 Evolution .................................... 55 6.6.2 Kinematicreconstruction. ..... 56 6.7 The running coupling constant αS ........................... 57 6.7.1 The argument of αS .............................. 57 6.7.2 The Monte Carlo scheme for αS ........................ 59 6.7.3 Options for the treatment of αS inpartonshowers . 59 6.8 Matrixelementcorrections . ....... 61 6.8.1 Softmatrixelementcorrections . ...... 61 6.8.2 Hardmatrixelementcorrections. ...... 62 6.8.3 Using Herwig++ matrix element corrections . ........ 63 6.9 Showering in the POWHEG scheme . .... 63 6.10Codestructure .................................. 65 7 Hadronization 68 7.1 Gluon splitting and cluster formation . .......... 68 7.2 Clusterfission................................... 68 7.3 Clusterdecays................................... 70 7.3.1 Mixingweights ................................. 72 7.4 HadronizationinBSMmodels . ..... 73 7.4.1 Stable strongly interacting particles . .......... 74 7.4.2 Baryonnumberviolation. 74 7.5 Codestructure ................................... 74 8 Underlying Event and Beam Remnants 77 8.1 Semi-hardpartonicscatters . ....... 77 8.2 Softpartonicscatters. ...... 79 8.2.1 MonteCarloimplementation. .... 81 8.3 Connection to different simulation phases . ........... 84 8.3.1 Parton showers and hard matrix elements . ...... 84 8.3.2 Minimumbiasprocess .. .. .. 84 8.3.3 Hadronization................................. 84 8.4 UA5parametrization .............................. .... 85 8.5 Codestructure ................................... 86 9 Hadron Decays 89 9.1 Particleproperties .............................. ..... 89 9.2 Lineshapes ...................................... 93 9.3 Taudecays....................................... 94 9.3.1 Hadroniccurrents.............................. 94 9.4 Strong and electromagnetic hadron decays . .......... 97 9.4.1 Scalarmesons.................................. 97 9.4.2 Vectormesons.................................. 99 9.4.3 Tensormesons.................................. 101 2 9.4.4 BaryonDecays ................................. 102 9.4.5 Inclusive strong and electromagnetic decays . .......... 104 9.5 Weakhadronicdecays .............................. 104 9.5.1 Exclusive semi-leptonic decays . ....... 104 9.5.2 Exclusive hadronic decays . .... 108 9.5.3 Weak inclusive decays . 109 9.5.4 Leptonicdecays................................ 109 9.5.5 b sγ ...................................... 110 → 9.6 Codestructure ................................... 110 10 Summary 113 A Repository Commands 115 A.1 Example......................................... 117 B Examples 119 B.1 Switching parts of the simulation off . ........ 119 B.2 Setupforminimumbiasruns . .... 120 B.3 Simulation of several hard processes in one event . ............. 121 B.4 Changingparticleproperties . ........ 122 B.5 Changingsomesimplecuts. ..... 123 B.6 Setting up an AnalysisHandler . ....... 123 B.7 UsageofROOT..................................... 125 B.7.1 Roothistograms ................................ 126 B.7.2 rtuplewithTTree............................... 127 B.8 UsingBSMmodels .................................. 129 B.8.1 MSSM...................................... 129 B.8.2 MUED...................................... 131 B.8.3 RSModel .................................... 131 B.8.4 Disabling Selected Decay Modes . ..... 132 B.9 Intrinsic pT ....................................... 132 B.10 LesHouchesEventHandler . ....... 133 B.11UseofLHAPDF .................................... 134 B.12UseofasimplesaturationmodelforPDFs . ......... 135 C Contrib 136 D Tuning 137 3 1 Introduction Herwig++ is a general-purpose event generator for the simulation of high-energy lepton-lepton, lepton-hadron and hadron-hadron collisions with special emphasis on the accurate simulation of QCD radiation. It builds upon the heritage of the HERWIG program [1–6], while providing a much more flexible structure for further development. It already includes several features more advanced than the last FORTRAN version. Herwig++ provides a full simulation of high energy collisions with the following special features: Initial- and final-state QCD jet evolution taking account of soft gluon interference via • angular ordering; A detailed treatment of the suppression of QCD radiation from massive particles, the dead- • cone effect [7]; The simulation of BSM physics including correlations between the production and decay • of the BSM particles together with the ability to add new models by simply encoding the Feynman rules; An eikonal model for multiple partonic scatterings to describe the underlying event [8]; • A cluster model of the hadronization of jets based on non-perturbative gluon splitting; • A sophisticated model of hadron and tau decays using matrix elements to give the momenta • of the decay products for many modes and including a detailed treatment of off-shell effects and spin correlations. Some of these features were already present in the first version of Herwig++ [9]. However, there have been many improvements to both the physics and structure of the simulation following this first release, most notably the extension to hadron-hadron collisions. Given the significant differences between the current version of the program, 2.3, and that described in [9] we will describe all of the features of the program in this paper. A number of other generators are also being (re-)written for the LHC era. The PYTHIA event generator is currently being rewritten as PYTHIA8 [10]. The rewrite of ARIADNE [11] is in progress as well. Like Herwig++, this is built on the platform of ThePEG [12], which we describe below. SHERPA [13] is a completely new event generator project. It is useful to start by recalling