# NLO QCD Production of Higgs boson plus jets with GoSam

## Abstract

After reviewing the main features of the GoSam framework for automated one-loop calculations, we present a selection of recent phenomenological results obtained with it. In particular, we focus on the recent calculation of NLO QCD corrections to the production of a Higgs boson in conjunction with jets at the LHC.

NLO QCD Production of Higgs boson plus jets with GoSam
\FullConferenceThe European Physical Society Conference on High Energy Physics

18-24 July, 2013

Stockholm, Sweden

## 1 Introduction

The large amount of data accumulated by the experimental collaborations at the Large Hadron Collider (LHC) allowed for a very detailed investigation of the Standard Model (SM) of particle physics. Moreover, the discovery of a Higgs boson with mass of about 126 GeV [1] finally confirmed the validity of the electroweak symmetry breaking mechanism [2].

In all these analyses, for example to further study the properties of the recently discovered Higgs boson, theory predictions play a fundamental role. They are not only needed for the signal, but also for the modeling of the relevant background processes, which share similar experimental signatures. Further, precise theory predictions are important in order to constrain model parameters in the event that a signal of New Physics is detected. Since leading-order (LO) results are affected by large uncertainties, theory predictions are not reliable without accounting for higher orders. Therefore, it is of primary interest to provide theoretical tools which are able to perform the comparison of LHC data to theory at NLO accuracy.

In the past few years, the progress in the automation of NLO calculations for multi-particle final states has been tremendous and led to the so-called “NLO revolution” [3]. Several automated frameworks for one-loop calculations [4] have been presented, which are based on various new theoretical developments [6]. It is fascinating to witness the number and quality of advanced automated NLO calculations that have been performed with different techniques.

In this presentation, we review the main features of the GoSam framework [5] for the automated computation of one-loop amplitudes. GoSam has been recently employed in several calculations at NLO QCD accuracy [7, 8, 9, 10, 11, 12, 13, 14] related to signal and backgrounds for Higgs boson production, as well as in the context of Beyond Standard Model (BSM) scenarios [15, 16] and electroweak studies [17], and has been successfully interfaced with Monte Carlo programs to merge multiple NLO matrix elements with parton showers [18, 19].

We also briefly describe a selection of recent phenomenological results obtained with GoSam, with particular attention to the recent calculations of NLO QCD corrections to the production of Higgs boson in conjunction with jets at the LHC.

## 2 The GoSam framework

GoSam combines automated diagram generation and algebraic manipulation [20, 21, 22, 23] with integrand-level reduction techniques [24]. Amplitudes are generated via Feynman diagrams, using QGRAF [20], FORM [21], spinney [23] and haggies [22]. The individual program tasks are managed by python scripts, so that the only task required from the user is the preparation of an input file in order to launch the generation of the source code and its compiling, without having to worry about the internal details. The input file contains specific information about: i) the process, such as a list of initial and final state particles, the order in the coupling constants, and the model; ii) the scheme employed, such as the regularization and renormalization schemes; iii) the system , such as paths to libraries or compiler options; iv) optional information to control the code generation.

After the generation of all contributing diagrams, the virtual corrections are evaluated using the -dimensional integrand-level reduction method, as implemented in Samurai [26] library, which allows for the combined determination of both cut-constructible and rational terms at once. Alternatively, the tensorial decomposition provided by Golem95C [27, 28, 29] is also available. Such reduction, which is numerically stable but more time consuming, is employed as a rescue system. After the reduction, all relevant master integrals can be computed by means of Golem95C [29], QCDLoop [30], or OneLoop [31].

As a novel approach to the integrand reduction, the method proposed in [25] allows for the extraction of all the coefficients in the integrand decomposition by performing a Laurent expansion, whenever the analytic form of the numerator function is known. This method has been implemented, within the GoSam framework, in the c++ library Ninja, showing an improvement in the computational performance, both in terms of speed and precision, with respect to the standard algorithms. More details are provided in the talk of T. Peraro at this conference [32]. The new library has been recently employed in the evaluation of NLO QCD corrections to [13].

GoSam can be used to generate and evaluate one-loop corrections in both QCD and electro-weak theory. Model files for BSM theories can be generated from a Universal FeynRules Output (UFO) [33] or LanHEP [34] file.

#### Code development

New features have been recently implemented within GoSam, with respect to the current public version. In order to deal with the complexity level of calculations such as [12], the GoSam code has been enhanced. On the one side, the generation algorithm has been improved by a more efficient diagrammatic layout: Feynman diagrams are grouped according to their topologies, namely global numerators are constructed by combining diagrams that have a common set, or subset, of denominators, irrespectively of the specific particle content. On the other side, additional improvements in the performances of GoSam have been achieved by exploiting the optimized manipulation of polynomial expressions available in Form 4.0 [35]. The possibility of employing numerical polarization vectors and the option to sum diagrams sharing the same propagators algebraically during the generation of the code led to an enormous gain in code generation time and reduction of code size.

Concerning the amplitude reduction, aside from the already mentioned new integrand-level reduction via Laurent expansion [25], GoSam has been enhanced to reduce integrands that may exhibit numerators with rank larger than the number of the denominators. This is indeed the case in the presence of effective couplings [9, 12], which appear in the large top-mass approximation, or when dealing with spin-2 particles [16]. For these cases, within the context of integrand-reduction techniques, the parametrization of the residues at the multiple-cut has to be extended and the decomposition of any one-loop amplitude acquires new master integrals [25]. The extended integrand decomposition has been implemented in the samurai library [36].

The new developments regarding the improved generation and reduction algorithms will be publicly available in the GoSam 2.0 release, which is currently in preparation.

#### Interfacing with MC and BLHA

The computation of physical observables at NLO accuracy, such as cross sections and differential distributions, requires to combine the one-loop results for the virtual amplitudes obtained with GoSam, with other tools that can take care of the computation of the real emission contributions and of the subtraction terms, needed to control the cancellation of IR singularities.

This can be obtained by embedding the calculation of virtual corrections within a Monte Carlo framework (MC), that can take care of the phase-space integration, and of the combination of the different pieces of the calculation. A table with a comprehensive list of GoSam interfaces with MC programs has been recently presented in [37].

In order to facilitate the communication between the programs computing virtual one-loop amplitudes and the MC frameworks, a standard interface called the *Binoth Les Houches Accord* (BLHA) [38] has been designed. Within the BLHA, the interaction between the One-loop Program (OLP) and the Monte Carlo framework (MC) proceeds in two phases. During the first phase, called pre-runtime phase, the MC creates an order file, which contains information about the setup and the subprocesses needed
from the OLP in order to perform the computation. The OLP reads the order file, checks availability for each item, and returns a contract file telling the MC what it can provide. In the second stage, the MC requires from the OLP the values of the virtual one-loop amplitudes at specific phase-space points.

## 3 Higgs boson production in Gluon Fusion

At the LHC, the dominant Higgs production mechanism proceeds via gluon fusion (GF), where the coupling of the Higgs to the gluons is mediated by a heavy quark loop. For this reason, the calculation of higher order corrections for the GF production of a Higgs boson in association with jets has received a lot of attention in the theory community over the past decade [39].

The developments in GoSam, described in Section 2, allowed to compute the NLO QCD corrections to the production of jets [9] and, for the first time, also jets [12] in GF (in the large top-mass limit). While a fully automated BLHA interface between GoSam and Sherpa [40] has been used for , the complexity of the integration for the process forced us to employ a hybrid setup which combines GoSam, Sherpa and the MadDipole/Madgraph4/MadEvent framework [41]. This calculation is indeed challenging both on the side of real-emission contributions and of the virtual corrections, which alone involve more than one-loop Feynman diagrams with up to rank-seven hexagons.

In the calculation of jets the cteq6L1 and cteq6mE parton-distribution functions were used for LO and NLO respectively, and a minimal set of cuts based on the anti- jet algorithm with , GeV and was applied. Figure 1 shows the distributions of the three jets and of the Higgs boson, respectively. The NLO corrections enhance all distributions for values lower than GeV, whereas their contribution is negative at higher .

This study also shows that the virtual contributions for generated by GoSam is ready to be paired with available Monte Carlo programs to aim at further phenomenological studies.

## 4 Other Phenomenological results

#### Diphotons+jets

GoSam in combination with MadDipole/MadGraph4/MadEvent has been used to calculate the NLO QCD corrections to jet [10] and jets [11], where the former also includes the fragmentation component. This calculation allowed for a first reliable prediction of the absolute normalization of this process, and demonstrated that the shape of important kinematical distributions is modified by higher-order effects.

#### Beyond the Standard Model

GoSam has been used to calculate the NLO Susy-QCD corrections to the production of a pair of the lightest neutralinos plus one jet at the LHC at TeV, which appears as a monojet signature in combination with missing energy. All non-resonant diagrams have been fully included, namely without using the assumption that production and decay factorize. We observe that the NLO corrections to the missing transverse energy are large, mainly due to additional channels opening up at NLO. The detailed setup can be found in [15].

Another recent BSM result obtained with GoSam+MadDipole/MadGraph4/MadEvent is the calculation of NLO QCD corrections to the production of a graviton in association with one jet [16], where the graviton decays into a photon pair, within ADD models of large extra dimensions [42]. The calculation is quite complicated due to the tensor structure introduced by spin-2 particles, and the non-standard propagator of the graviton, coming from the summation over Kaluza-Klein modes. It is interesting to notice that the -factors of the invariant mass distribution of the photon pair emitted in the decay of the graviton are not uniform. Since the latter is used to derive exclusion limits, it is advisable to take into account the effect of NLO corrections. For details we refer to [16].

#### Higgs+Vector Boson+jet

GoSam was interfaced with the POWHEG BOX to compute the associate production of a Higgs boson, a vector boson, and one jet [18]. In this calculation, the improved MiNLO procedure [43] was used to obtain NLO accurate predictions also when the jet is not resolved. Using this interface the generation of any process is fully automated, except for the construction of the Born phase space.

#### NLO QCD corrections to

The production rate for a Higgs boson associated with a top-antitop pair () is particularly interesting to study the properties of the newly discovered Higgs boson, since it is directly proportional to the SM Yukawa coupling of the Higgs boson to the top quark. We recently presented the complete NLO QCD corrections to the process jet () at the LHC [13]. The goal of the calculation was twofold. On the one hand, it is important for the phenomenological analyses at the LHC, in particular for the high- region, where the presence of the additional jet can be relevant. On the other hand, from the technical point of view, for the presence of two mass scales (Higgs boson and top quark) and internal massive particles, together with a high number of diagrams, constitutes a challenge for many reduction algorithms. This calculation represents the first application of the novel reduction algorithm, implemented in the library Ninja, based on integrand-level reduction via Laurent expansion [25].

## 5 Conclusions

GoSam is a flexible and widely applicable tool for the automated calculation of the virtual part of multi-particle scattering amplitudes at NLO accuracy. After interfacing it with MC programs, that can perform integration over phase space and combine the contributions coming from real emission and subtraction terms as well, total cross-sections and differential distributions can be easily obtained for a variety of processes of interest at the LHC.

Boosted by state-of-the-art techniques for the reduction of the scattering amplitudes, GoSam provides a reliable answer for multi-leg amplitudes in the presence of massive internal and external legs and propagators, such as the production of a Higgs boson in conjunction with a top-quark pair, as well as in configurations with relatively high multiplicity, such as Higgs boson plus jets production. While the GoSam code will be further improved, it will be interesting to observe whether the attempt of extending integrand-level techniques to higher orders [44] will succeed and provide a comparable level of automation, at least for the calculation of the virtual parts.

Other challenges for the near future involve interfacing GoSam with MC programs for an automated generation of the full cross section including parton showering, and ultimately the production of codes and results to be used within experimental analyses. We believe that the amount of recent calculations that were produced with the GoSam framework shows, both in terms of stability and precision, that it is an ideal multi-purpose tool for studying the physics of the LHC.

### Acknowledgments

The work of G.C. was supported by DFG Sonderforschungsbereich Transregio 9, Computergestützte Theoretische Teilchenphysik and the support of the Research Executive Agency (REA) of the European Union under the Grant Agreement number PITN-GA-2010-264564 (LHCPhenoNet). H.v.D., G.L., P.M., and T.P. are supported by the Alexander von Humboldt Foundation, in the framework of the Sofja Kovaleskaja Award Project “Advanced Mathematical Methods for Particle Physics”, endowed by the German Federal Ministry of Education and Research. The work of G.O. was supported in part by the National Science Foundation under Grant PHY-1068550 and PSC-CUNY Award No. 65188-00 43. This research work benefited of computing resources from the Rechenzentrum Garching and the the CTP cluster of the New York City College of Technology.

### References

- Aad G et al. (ATLAS Collaboration) 2012 Phys.Lett. B716 1–29; Chatrchyan S et al. (CMS Collaboration) 2012 Phys.Lett. B716 30–61
- Englert F and Brout R 1964 Phys.Rev.Lett. 13 321–323; Higgs P W 1964 Phys.Lett. 12 132–133
- Salam G P 2010 PoS ICHEP2010 556; Perret-Gallix D 2013 J.Phys.Conf.Ser. 454 012051
- Berger C, Bern Z, Dixon L, Febres Cordero F, Forde D et al. 2008 Phys.Rev. D78 036003; Bevilacqua G, Czakon M, Garzelli M, van Hameren A, Kardos A et al. 2013 Comput.Phys.Commun. 184 986–997; Hirschi V, Frederix R, Frixione S, Garzelli M V, Maltoni F et al. 2011 JHEP 1105 044; Cascioli F, Maierhofer P and Pozzorini S 2012 Phys.Rev.Lett. 108 111601; Agrawal S, Hahn T and Mirabella E 2012 PoS LL2012 046; Badger S, Biedermann B, Uwer P and Yundin V 2013 Comput.Phys.Commun. 184 1981–1998; Actis S, Denner A, Hofer L, Scharf A and Uccirati S 2013 JHEP 1304 037
- Cullen G, Greiner N, Heinrich G, Luisoni G, Mastrolia P et al. 2012 Eur.Phys.J. C72 1889
- Ellis R K, Kunszt Z, Melnikov K and Zanderighi G 2012 Phys.Rept. 518 141–250; Bern Z, Dixon L and Kosower D 2007 Annals Phys. 322 1587; Ossola G 2013 (Preprint 1310.3214)
- Greiner N, Guffanti A, Reiter T and Reuter J 2011 Phys.Rev.Lett. 107 102002
- Greiner N, Heinrich G, Mastrolia P, Ossola G, Reiter T et al. 2012 Phys.Lett. B713 277–283
- van Deurzen H, Greiner N, Luisoni G, Mastrolia P, Mirabella E et al. 2013 Phys.Lett. B721 74–81
- Gehrmann T, Greiner N and Heinrich G 2013 JHEP 1306 058
- Gehrmann T, Greiner N and Heinrich G 2013 (Preprint 1308.3660)
- Cullen G, van Deurzen H, Greiner N, Luisoni G, Mastrolia P et al. 2013 Phys.Rev.Lett. 111 131801
- van Deurzen H, Luisoni G, Mastrolia P, Mirabella E, Ossola G et al. 2013 (Preprint 1307.8437)
- Dolan M J, Englert C, Greiner N and Spannowsky M 2013 (Preprint 1310.1084)
- Cullen G, Greiner N and Heinrich G 2013 Eur.Phys.J. C73 2388
- Greiner N, Heinrich G, Reichel J and von Soden-Fraunhofen J F 2013 (Preprint 1308.2194)
- Chiesa M, Montagna G, Barze‘ L, Moretti M, Nicrosini O et al. 2013 Phys.Rev.Lett. 111 121801; Mishra K, Becher T, Barze L, Chiesa M, Dittmaier S et al. 2013 (Preprint 1308.1430)
- Luisoni G, Nason P, Oleari C and Tramontano F 2013 (Preprint 1306.2542)
- Hoeche S, Huang J, Luisoni G, Schoenherr M and Winter J 2013 Phys.Rev. D88 014040
- Nogueira P 1993 J.Comput.Phys. 105 279–289
- Vermaseren J A M 2000 (Preprint math-ph/0010025)
- Reiter T 2010 Comput.Phys.Commun. 181 1301–1331
- Cullen G, Koch-Janusz M and Reiter T 2011 Comput.Phys.Commun. 182 2368–2387
- Ossola G, Papadopoulos C G and Pittau R 2007 Nucl.Phys. B763 147–169; Ossola G, Papadopoulos C G and Pittau R 2007 JHEP 0707 085; Ellis R K, Giele W T and Kunszt Z 2008 JHEP 03 003; Ossola G, Papadopoulos C G and Pittau R 2008 JHEP 0805 004; Mastrolia P, Ossola G, Papadopoulos C and Pittau R 2008 JHEP 0806 030
- Mastrolia P, Mirabella E and Peraro T 2012 JHEP 1206 095
- Mastrolia P, Ossola G, Reiter T and Tramontano F 2010 JHEP 1008 080
- Binoth T, Guillet J P, Heinrich G, Pilon E and Reiter T 2009 Comput.Phys.Commun. 180 2317–2330
- Heinrich G, Ossola G, Reiter T and Tramontano F 2010 JHEP 1010 105
- Cullen G, Guillet J, Heinrich G, Kleinschmidt T, Pilon E et al. 2011 Comput.Phys.Commun. 182 2276–2284
- van Oldenborgh G 1991 Comput.Phys.Commun. 66 1–15; Ellis R K and Zanderighi G 2008 JHEP 02 002
- van Hameren A 2011 Comput.Phys.Commun. 182 2427–2438
- Peraro T, presentation at EPS-HEP 2013, in these proceedings
- Christensen N D and Duhr C 2009 Comput. Phys. Commun. 180 1614–1641; Degrande C, Duhr C, Fuks B, Grellscheid D, Mattelaer O et al. 2011 (Preprint 1108.2040); Alloul A, Christensen N D, Degrande C, Duhr C and Fuks B 2013 (Preprint 1310.1921)
- Semenov A 2010 (Preprint 1005.1909)
- Kuipers J, Ueda T, Vermaseren J and Vollinga J 2013 Comput.Phys.Commun. 184 1453–1467
- Mastrolia P, Mirabella E, Ossola G, Peraro T and van Deurzen H 2012 PoS LL2012 028 (Preprint 1209.5678)
- Cullen G, van Deurzen H, Greiner N, Heinrich G, Luisoni G et al. 2013 (Preprint 1309.3741)
- Binoth T, Boudjema F, Dissertori G, Lazopoulos A, Denner A et al. 2010 Comput.Phys.Commun. 181 1612–1622; Alioli S, Badger S, Bellm J, Biedermann B, Boudjema F et al. 2013 (Preprint 1308.3462)
- Dittmaier S et al. (LHC Higgs Cross Section Working Group) 2011 (Preprint 1101.0593); Dittmaier S, Dittmaier S, Mariotti C, Passarino G, Tanaka R et al. 2012 (Preprint 1201.3084); Heinemeyer S et al. (The LHC Higgs Cross Section Working Group) 2013 (Preprint 1307.1347)
- Gleisberg T, Hoeche S, Krauss F, Schonherr M, Schumann S et al. 2009 JHEP 0902 007
- Frederix R, Gehrmann T and Greiner N 2008 JHEP 0809 122; Frederix R, Gehrmann T and Greiner N 2010 JHEP 1006 086; Stelzer T and Long W 1994 Comput.Phys.Commun. 81 357–371; Maltoni F and Stelzer T 2003 JHEP 0302 027; Alwall J, Demin P, de Visscher S, Frederix R, Herquet M et al. 2007 JHEP 0709 028
- Arkani-Hamed N, Dimopoulos S and Dvali G 1998 Phys.Lett. B429 263–272; Antoniadis I, Arkani-Hamed N, Dimopoulos S and Dvali G 1998 Phys.Lett. B436 257–263
- Hamilton K, Nason P, Oleari C and Zanderighi G 2013 JHEP 1305 082
- Mastrolia P and Ossola G 2011 JHEP 1111 014 (Preprint 1107.6041); Badger S, Frellesvig H and Zhang Y 2012 JHEP 1204 055; Zhang Y 2012 JHEP 1209 042; Mastrolia P, Mirabella E, Ossola G and Peraro T 2012 Phys.Lett. B718 173–177; Feng B and Huang R 2013 JHEP 1302 117; Mastrolia P, Mirabella E, Ossola G and Peraro T 2013 Phys.Rev. D87 085026; Mastrolia P, Mirabella E, Ossola G and Peraro T 2013 (Preprint 1307.5832); Badger S, Frellesvig H and Zhang Y 2013 (Preprint 1310.1051)