SHERPA is a general-purpose Monte Carlo event generator for the simulation of particle collisions in high-energy collider experiments. We summarise essential features and improvements of the SHERPA 2.2 release series, which is heavily used for event generation in the analysis and interpretation of LHC Run 1 and Run 2 data. We highlight a decade of developments towards ever higher precision in the simulation of particle-collision events. Figure 1: Overview of the SHERPA 2.2 event generator framework. Interfaces/OutputsAMEGIC [6] and COMIX [7,8]. They are used for the simulation of parton-level events within the Standard Model and beyond, and for the decay of heavy resonances such as W , Z, or Higgs bosons or top quarks. Both include automated methods for efficient phase-space integration and algorithms for the subtraction of infrared divergences in calculations at next-to-leading order (NLO) in QCD [9, 10, 11] and the electroweak theory [12]. For the evaluation of virtual corrections at NLO accuracy SHERPA relies on interfaces to dedicated one-loop providers, e.g. BLACKHAT [13], OPENLOOPS [14] and RECOLA [15,16]. The default parton-showering algorithm of the SHERPA 2.2 series is the CSSHOWER [17], based on Catani-Seymour dipole factorisation [9,10,18]. As of version 2.2.0 SHERPA also features an independent second shower implementation, DIRE [19,20,21]. For the matching of NLO QCD matrix elements with parton showers SHERPA implements the MC@NLO method [22,23]. For NNLO QCD calculations the UN 2 LOPS method [24, 25] is used. The merging of multi-jet production processes at leading order [26,27,28] and next-to-leading order [29,30] is based on truncated parton showers. Multiple parton interactions are implemented via the Sjöstrand-van-Zijl model [31]. The hadronisation of partons into hadrons is modelled by a cluster fragmentation model [32]. Alternatively, in particular for uncertainty estimations, an interface to the Lund fragmentation model [33] of PYTHIA [34] is available. SHERPA provides a large library for the simulation of τ -lepton and hadron decays, including many form-factor models. Furthermore, a module for the simulation of QED final-state radiation in particle decays [35], which is accurate to first order in α for many channels is built-in. To account for spin correlations in production and subsequent decay processes the algorithm described in [36] is implemented. Events generated with SHERPA can be cast into various output formats for further processing, with the HEPMC [37] format being the most commonly used. In the specific case of parton-level events, at the leading and next-to-leading order in QCD, additional output formats are supported. They include Les Houches Event Files [38], NTUPLE files for NLO QCD events [39] and cross-section interpolation grids produced via MCGRID [40,41] in the APPLGRID [42] and FASTNLO [43,44] formats. To analyse events on-the-fly a runtime interface to the RIVET package [45] can be used conveniently.

We compute the total cross-section for Higgs boson production in bottom-quark fusion using the so-called FONLL method for the matching of a scheme in which the $b$-quark is treated as a massless parton to that in which it is treated as a massive final-state particle. We discuss the general framework for the application of the FONLL method to this process, and then we present explicit expressions for the case in which the next-to-next-to-leading-log five-flavor scheme result is combined with the leading-order $\cal O(\alpha_s^2)$ four-flavor scheme computation. We compare our results in this case to the four-and five-flavor scheme computations, and to the so-called Santander matching.Comment: 15 pages, 4 figures; three typos corrected. Final version, to be published in Phys. Lett.

We compute the total cross-section for Higgs boson production in bottom-quark fusion using the so-called FONLL method for the matching of a scheme in which the b-quark is treated as a massless parton to that in which it is treated as a massive final-state particle, and extend our previous results to the case in which the next-to-next-to-leading-log five-flavor scheme result is combined with the next-to-leading-order O(α 3 s ) four-flavor scheme computation.

Jet substructure tools have proven useful in a number of high-energy particle-physics studies. A particular case is the discrimination, or tagging, between a boosted jet originated from an electroweak boson (signal), and a standard QCD parton (background). A common way to achieve this is to cut on a measure of the radiation inside the jet, i.e. a jet shape. Over the last few years, analytic calculations of jet substructure have allowed for a deeper understanding of these tools and for the development of more efficient ones. However, analytic calculations are often limited to the region where the jet shape is small. In this paper we introduce a new approach in perturbative QCD to compute jet shapes for a generic boosted jets, waiving the above limitation. We focus on an example common in the substructure literature: the jet mass distribution after a cut on the N -subjettiness τ 21 ratio, extending previous works to the region relevant for phenomenology. We compare our analytic predictions to Monte Carlo simulations for both plain and SoftDrop-groomed jets. We use our results to construct analytically a decorrelated tagger.

We present a new calculation for the production of isolated photon pairs at the LHC with $$ {\mathrm{NNLL}}_{{\mathcal{T}}_0}^{\prime } $$ NNLL T 0 ′ +NNLO accuracy. This is the first implementation within the Geneva Monte Carlo framework of a process with a nontrivial Born-level definition which suffers from QED singularities. Throughout the computation we use a smooth-cone isolation algorithm to remove such divergences. The higher-order resummation of the 0-jettiness resolution variable $$ {\mathcal{T}}_0 $$ T 0 is based on a factorisation formula derived within Soft-Collinear Effective Theory which predicts all of the singular, virtual and real NNLO corrections. Starting from this precise parton-level prediction and by employing the Geneva method, we provide fully showered and hadronised events using Pythia8, while retaining the NNLO QCD accuracy for observables which are inclusive over the additional radiation. We compare our final predictions to LHC data at 7 TeV and find good agreement.

In this work we explore first necessary steps to contruct a fully massive version of a variable flavour number scheme. In particular we focus, as an example, on an extension of the five-flavour scheme, where instead of neglecting explicit initial state quark mass effects, we retain all massive dependence, while keeping the resummation properties of the massless five-flavour scheme. We name this scheme five-flavour-massive (5FMS) scheme. Apart from consistently modified parton distribution functions, we provide all the ingredients that are needed to implement this scheme at MC@NLO accuracy, in a Monte Carlo event generator. As proof of concept we implement this scheme in SHERPA, and perfom a comparison of the new scheme with traditional ones for the simple process of scalar particle production in bottom quark fusion.where each phase space integral by itself is infrared finite.

Reprinted with permission from the American Physical Society: Krauss, Frank, Napoletano, Davide and Schumann, Steen (2017) 'Simulating b-associated production of Z and Higgs bosons with the SHERPA event generator.', Physical review D., 95 (3). 036012 c 2017 American Physical Society by the American Physical Society. Readers may view, browse, and/or download material for temporary copying purposes only, provided these uses are for noncommercial personal purposes. Except as provided by law, this material may not be further reproduced, distributed, transmitted, modied, adapted, performed, displayed, published, or sold in whole or part, without prior written permission from the American Physical Society.Additional information: Use policyThe full-text may be used and/or reproduced, and given to third parties in any format or medium, without prior permission or charge, for personal research or study, educational, or not-for-prot purposes provided that:• a full bibliographic reference is made to the original source • a link is made to the metadata record in DRO • the full-text is not changed in any way The full-text must not be sold in any format or medium without the formal permission of the copyright holders.Please consult the full DRO policy for further details. We compare four-and five-flavor scheme predictions for b-associated production of Z and Higgs bosons. The results are obtained with SHERPA 's MC@NLO implementation for the four-flavor scheme, treating the b's as massive, and with multijet merging at leading and next-to-leading order for the fiveflavor schemes. Comparison with data for Z þ bðbÞ production at the 7 TeV LHC exhibits strengths and weaknesses of the different approaches and is used to validate predictions for b-associated Higgs-boson production at the 13 TeV Run II.

No abstract