The NNLO QCD corrections to Z boson production at large transverse momentum

Abstract

Abstract

The transverse momentum distribution of massive neutral vector bosons can be measured to high accuracy at hadron colliders. The transverse momentum is caused by a partonic recoil, and is determined by QCD dynamics. We compute the single and double-differential transverse momentum distributions for fully inclusive Z/γ production including leptonic decay to next-to-next-to-leading order (NNLO) in perturbative QCD. We also compute the same distributions normalised to the cross sections for inclusive Z/γ production, i.e. integrated over the transverse momentum of the lepton pair. We compare our predictions for the fiducial cross sections to the 8 TeV data set from the ATLAS and CMS collaborations, which both observed a tension between data and NLO theory predictions, using the experimental cuts and binning. We find that the inclusion of the NNLO QCD effects does not fully resolve the tension with the data for the unnormalised p T ^Z distribution. However, we observe that normalising the NNLO Z-boson transverse momentum distribution by the NNLO Drell-Yan cross section substantially improves the agreement between experimental data and theory, and opens the way for precision QCD studies of this observable.

Published for SISSA by Springer

Received: May 19, 2016 Accepted: July 19, 2016 Published: July 27, 2016

A. Gehrmann-De Ridder,a,b,c T. Gehrmann,b,c E.W.N. Glover,d A. Hussa and T.A. Morgand

aInstitute for Theoretical Physics, ETH,

CH-8093 Zurich, Switzerland

bDepartment of Physics, University of Zurich,

CH-8057 Zurich, Switzerland

cKavli Institute for Theoretical Physics, UC Santa Barbara, Santa Barbara, U.S.A.

dInstitute for Particle Physics Phenomenology, Department of Physics, University of Durham, Durham, DH1 3LE, U.K.

E-mail: mailto:[email protected]

Web End [email protected] , mailto:[email protected]

Web End [email protected] , mailto:[email protected]

Web End [email protected] , mailto:[email protected]

Web End [email protected] , mailto:[email protected]

Web End [email protected]

Abstract: The transverse momentum distribution of massive neutral vector bosons can be measured to high accuracy at hadron colliders. The transverse momentum is caused by a partonic recoil, and is determined by QCD dynamics. We compute the single and double-di erential transverse momentum distributions for fully inclusive Z/ production including leptonic decay to next-to-next-to-leading order (NNLO) in perturbative QCD. We also compute the same distributions normalised to the cross sections for inclusive Z/ production, i.e. integrated over the transverse momentum of the lepton pair. We compare our predictions for the ducial cross sections to the 8 TeV data set from the ATLAS and CMS collaborations, which both observed a tension between data and NLO theory predictions, using the experimental cuts and binning. We nd that the inclusion of the NNLO QCD e ects does not fully resolve the tension with the data for the unnormalised pZT distribution. However, we observe that normalising the NNLO Z-boson transverse momentum distribution by the NNLO Drell-Yan cross section substantially improves the agreement between experimental data and theory, and opens the way for precision QCD studies of this observable.

Keywords: QCD Phenomenology

ArXiv ePrint: 1605.04295

Open Access, c

[circlecopyrt] The Authors.

Article funded by SCOAP3. doi:http://dx.doi.org/10.1007/JHEP07(2016)133

Web End =10.1007/JHEP07(2016)133

The NNLO QCD corrections to Z boson production at large transverse momentum

JHEP07(2016)133

Contents

1 Introduction 1

2 Setup 3

3 The Z boson transverse momentum distribution 4

4 Double-di erential distributions 6

5 Summary and conclusions 13

JHEP07(2016)133

1 Introduction

The Drell-Yan production of lepton pairs mediated through a virtual photon or a Z boson is an important process at hadron colliders such as the Large Hadron Collider (LHC). This process has a clean signature with a large event rate, leading to small experimental errors over a wide range of energies. It is a key process in probing Standard Model physics, for example tting parton distribution functions (PDFs) and measuring the strong coupling constant s. Furthermore, it is an important background in searches for signatures of physics beyond the Standard Model, such as supersymmetry and dark matter.

There has been a signi cant amount of e ort to ensure that theoretical predictions match the high precision data available from the LHC. The inclusive cross section for Z/ production is known at NNLO accuracy in QCD [1{9]. Corrections beyond this order have been studied in the soft-virtual approximation [10{15]. The NNLO QCD corrections have been combined with a resummation of next-to-next-to-leading logarithmic e ects (NNLL) [16{19] which is necessary to predict the transverse momentum distribution of the Z boson at small pT and matched with parton showers [20, 21]. The NLO electroweak corrections [22{25] and the mixed QCD-EW corrections [26, 27] have also been computed.

As illustrated in gure 1, the transverse momentum of the Z boson is generated by the emission of QCD radiation so that the xed order calculation at O( 2s), which is NNLO

for the inclusive cross section, corresponds to only an NLO-accurate prediction of the transverse momentum distribution. The NLO EW corrections to the Z-boson transverse momentum distribution [28{31] are typically negative and at the level of a few percent at small to moderate pZT, but are enhanced by the well-known electroweak Sudakov logarithms when pZT MZ. At small pZT, the electroweak corrections approach the corrections to the

total Z production cross section, such that they largely cancel when normalising the pZT distribution to the total cross section.

Recently the O( 3s) NNLO QCD corrections to Z+jet production were computed [32{

36]. Building upon the calculation presented in ref. [32, 33], we exploit our highly exible

{ 1 {

pZT [negationslash]= 0

Figure 1. A schematic diagram demonstrating the Z boson recoiling against hard radiation.

and numerically stable code to compute the transverse momentum distribution of the Z boson at nite transverse momentum at NNLO precision. To achieve this we relax the requirement of observing a nal state jet and instead impose a low transverse momentum cut on the Z boson. This transverse momentum cut ensures the infrared niteness of the NNLO calculation, since it enforces the presence of nal-state partons to balance the transverse momentum of the Z boson.

The production of Z bosons (or, more generally, of lepton pairs with given invariant mass) at large transverse momentum has been studied extensively at the LHC by the ATLAS [37, 38], CMS [39, 40] and LHCb [41, 42] experiments. In order to reduce the systematic uncertainty on the measurement, the transverse momentum distribution is commonly normalised to the pT -inclusive Z-boson production cross section. ATLAS and CMS both observed a tension between their measurements and existing NLO QCD predictions, highlighting the potential importance of higher order corrections to this process.

Both experiments present their measurements in the form of ducial cross sections for a restricted kinematical range of the nal state leptons (in invariant mass, transverse momentum and rapidity). In view of a comparison between data and theory, this form of presenting the experimental data is preferable over a cross section that is fully inclusive in the lepton kinematics (requiring a theory-based extrapolation into phase space regions outside the detector coverage). Consequently, the theoretical calculation must take proper account of these restrictions in the nal state lepton kinematics.

The unnormalised pZT distribution represents an absolute cross section measurement based on event counting rates. As with any absolute measurement, it has the disadvantage of being sensitive to the proper modelling of acceptance corrections, and of relying on the absolute determination of the integrated luminosity of the data sample under consideration. At the LHC the luminosity uncertainty alone amounts to about 3%. In order to reduce the luminosity uncertainty, the data can be normalised to the Drell-Yan cross section for the corresponding ducial phase space. This is obtained from the cross section for Z boson production with the same transverse momentum and rapidity cuts on the individual leptons, but integrated over all possible transverse momenta of the Z boson. On the theoretical side, this amounts to normalising the distribution to the NNLO pp ! [lscript]+[lscript]+X cross section in which the ducial cuts are applied to the leptons, but which is fully inclusive on the transverse momentum of the lepton-pair.

In this paper, we compute the NNLO QCD corrections to the transverse momentum distribution for Z production at large transverse momentum, fully inclusive on the accom-

{ 2 {

Z

JHEP07(2016)133

panying hadronic nal state. We present the corrections using the nal state cuts used by the ATLAS [38] and CMS [40] collaborations. We study the NNLO corrections to both the absolute (unnormalised) di erential cross section,

d dpZT

[vextendsingle][vextendsingle][vextendsingle][vextendsingle]

pZT>20 GeV

dZJLO dpZT

+ dZJNLO dpZT

+ dZJNNLO dpZT

(1.1)

and the di erential cross section normalised to the relevant Drell-Yan cross section,

1 [notdef]

d dpZT

[vextendsingle][vextendsingle][vextendsingle][vextendsingle]

pZT>20 GeV

JHEP07(2016)133

with

[integraldisplay]

d dpZT

1

0

=

dpZT ZLO + ZNLO + ZNNLO. (1.2)

The NNLO corrections lead to a substantially better theoretical description of the experimental data. We nd that the inclusion of the NNLO QCD e ects does not fully resolve the tension with the data for the unnormalised pZT distribution. However, we observe that normalising the NNLO Z-boson transverse momentum distribution by the NNLO Drell-Yan cross section leads to an excellent agreement with the experimental distributions.

2 Setup

The NNLO corrections to the production of a Z boson at nite transverse momentum receive contributions from three types of parton-level processes: (a) the two-loop corrections to Z-boson-plus-three-parton processes [43{45], (b) the one-loop corrections to Z-boson-plus-four-parton processes [46{51] and (c) the tree-level Z-boson-plus- ve-parton processes [50{54]. The ultraviolet renormalised matrix elements for these processes are integrated over the nal state phase space appropriate to Z boson nal states, including a cut on the pZT. All three types of contributions are infrared-divergent and only their sum is nite.

In this calculation we employ the antenna subtraction method [55{63] to isolate the infrared singularities in the di erent contributions to enable their cancellation prior to the numerical implementation. The construction of the calculation and the subtraction terms is exactly as described in ref. [32, 33]. All partonic channels relevant to lepton pair production through the exchange of an on-shell or o -shell Z boson or a virtual photon are included. Our calculation is implemented in a newly developed parton-level Monte Carlo generator NNLOJET. This program provides the necessary infrastructure for the antenna subtraction of hadron collider processes at NNLO and performs the integration of all contributing subprocesses at this order. Components of it have also been used in other NNLO QCD calculations [64{69] using the antenna subtraction method. Other processes can be added to NNLOJET provided the matrix elements are available.

To describe the normalised distributions, we also implemented the NNLO QCD corrections to inclusive Z/ production including leptonic decays in NNLOJET, and validated this implementation against the publicly available codes FEWZ [8] and Vrap [9].

{ 3 {

For our numerical computations, we use the NNPDF3.0 parton distribution functions (PDFs) [70] with the value of s(MZ) = 0.118 at NNLO, and MZ = 91.1876 GeV. Note that we systematically use the same set of PDFs and the same value of s(MZ) for the

NLO and NNLO predictions. The factorisation and renormalisation scales are chosen dynamically on an event-by-event basis as,

[notdef] [notdef]R = [notdef]F =

qm2[lscript][lscript] + (pZT)2, (2.1)

where m[lscript][lscript] and pZT are the invariant mass and the transverse momentum of the nal state lepton pair respectively. The theoretical uncertainty is estimated by varying the scale choice by a factor in the range [1/2, 2]. The electroweak coupling constant is derived from the Fermi constant in the G[notdef] scheme, which absorbs large logarithms of the light fermion masses induced by the running of the coupling constant from the Thomson limit (Q2 = 0) to the electroweak scale into the tree-level coupling. We also impose a cut on the transverse momentum of the Z boson, pZT > 20 GeV. In the low transverse momentum region, large logarithmic corrections of the form logn(pT /MZ) appear at each order in the perturbative expansion in s, spoiling its convergence. A reliable theoretical prediction in this region can only be obtained by resummation [16{19] of these logarithms to all orders in perturbation theory. The cut on the transverse momentum also ensures the applicability of our approach, originally developed for Z+jet production at NNLO. In the computation of the inclusive lepton pair production cross section used to normalise the transverse momentum distributions, we choose the same scale (2.1) but varied independently over the same range of scale variation. The inclusive cross section in this case is, however, dominated by the regime where [notdef]2 m2[lscript][lscript].

3 The Z boson transverse momentum distribution

The experimental measurements of Z-boson production at nite transverse momentum are presented in the form of ducial cross sections over a restricted phase space for the nal state leptons, which is fully contained in the detectors coverage. The NNLO corrections to the transverse momentum distribution can be compared to data by considering the same cuts to the lepton kinematics as presented in the ATLAS [38] and CMS [40] analyses using data from Run 1 of the LHC with ps = 8 TeV, which are summarised in table 1. In this section we will focus on the ATLAS measurement in the ducial region de ned by a broad dilepton invariant mass window around the Z resonance, 66 GeV < m[lscript][lscript] < 116 GeV, and

compare both the absolute and normalised pZT distribution to the experimental data.

Figure 2 shows a comparison of data from the ATLAS analysis within this ducial region. At low transverse momentum, the NNLO correction increases the cross section by about 6% (compared to NLO) and signi cantly reduces the scale uncertainty. However, there is still signi cant tension between the ATLAS data and the NNLO prediction.

The data presented in gure 2 does not include the error on the integrated luminosity, which amounts to an overall normalisation uncertainty of 2.8% on all data points. To cancel the luminosity uncertainty from the measured data, it is more appropriate to normalise

{ 4 {

JHEP07(2016)133

ATLAS CMS leading lepton [notdef] [lscript]1[notdef] < 2.4 [notdef] [lscript]1[notdef] < 2.1 p[lscript]1T > 20 GeV p[lscript]1T > 25 GeV sub-leading lepton [notdef] [lscript]2[notdef] < 2.4 [notdef] [lscript]2[notdef] < 2.4 p[lscript]2T > 20 GeV p[lscript]2T,2 > 10 GeV

Table 1. Kinematical cuts used to de ne the ducial phase space for the nal state leptons in the measurements of ATLAS [38] and CMS [40].

Figure 2. The unnormalised Z-boson transverse momentum distribution for the cuts given in table 1 and 66 GeV < m[lscript][lscript] < 116 GeV. ATLAS data is taken from ref. [38]. The luminosity error is not shown. The green bands denote the NLO prediction with scale uncertainty and the blue bands show the NNLO prediction with scale uncertainty.

the data by the measured values for the inclusive lepton pair cross section in this ducial bin. The cross section for this mass window was measured to be [38],

exp(66 GeV < m[lscript][lscript] < 116 GeV) = 537.10 [notdef] 0.45% (sys.) [notdef] 2.80% (lumi.) pb.

For the NNLO prediction of the cross section in this ducial region, we obtain,

NNLO(66 GeV < m[lscript][lscript] < 116 GeV) = 507.9+2.40.7 pb.

We see that there is some tension between the measured cross section compared to the theoretical result. The normalised distribution is presented in gure 3, where excellent agreement is observed between the normalised NNLO prediction and the experimental data across a wide range of transverse momentum. The scale bands on the normalised theory predictions are obtained by independently varying the scale in the numerator and denominator, where we use the NNLO prediction for the inclusive Drell-Yan cross section in the normalisation throughout.

{ 5 {

JHEP07(2016)133

Figure 3. The normalised Z-boson transverse momentum distribution for the cuts given in table 1 and 66 GeV < m[lscript][lscript] < 116 GeV. ATLAS data is taken from ref. [38]. The green bands denote the

NLO prediction with scale uncertainty and the blue bands show the NNLO prediction with scale uncertainty.

The electroweak corrections to this distribution are known to be small at moderate transverse momenta but become more sizeable in the high-pZT tail where they can reach

15% for pZT [greaterorsimilar] 600 GeV [28{31]. It can be therefore expected that the inclusion of

electroweak corrections will further improve the prediction of the shape in the tail. However, this region of the distribution is still dominated by the statistical uncertainties of the experimental data and a fully consistent inclusion of electroweak corrections is beyond the scope of this work.

4 Double-di erential distributions

We now compare our theoretical predictions for dilepton production at large transverse momentum for di erent ranges of (a) the dilepton mass and (b) the dilepton rapidity. These double-di erential distributions have also been studied by ATLAS [38] and CMS [40] with Run 1 data using the lepton rapidity and transverse momentum cuts summarised in table 1.

We rst consider the ATLAS measurement, which covers a broader kinematical range. In gure 4 we present the unnormalised double-di erential distribution with respect to the transverse momentum of the Z boson and the invariant mass of the lepton pair, m[lscript][lscript],

normalised to the NLO prediction and compare it to the ATLAS data from ref. [38]. We observe that there is good agreement between the NNLO prediction and the data in the low m[lscript][lscript] mass windows. For values of m[lscript][lscript] close to the Z boson mass, where the statistical accuracy of the data is highest, the NNLO prediction is below the data by about 5-8%.

To obtain normalised distributions, these data sets are divided by the inclusive dilepton cross sections for each ducial bin, de ned by the lepton cuts given in table 1 and the

{ 6 {

JHEP07(2016)133

JHEP07(2016)133

Figure 4. The unnormalised double-di erential transverse momentum distribution for the Z boson in windows of invariant mass of the leptons, m[lscript][lscript], with a rapidity cut on the Z boson of [notdef]yZ[notdef] < 2.4.

The ATLAS data is taken from ref. [38]. The luminosity error is not shown. The green bands denote the NLO prediction with scale uncertainty and the blue bands show the NNLO prediction with scale uncertainty.

Figure 5. The inclusive dilepton cross section for the same m[lscript][lscript] bins as in gure 4 and with a rapidity cut on the Z boson of [notdef]yZ[notdef] < 2.4. The experimental data is taken from the ATLAS

analysis in ref. [38]. The ticks on the vertical error bands denote the systematic uncertainty from the measurement, the vertical bars without the ticks are the luminosity uncertainty only. The blue bands show the NNLO prediction with scale uncertainty.

{ 7 {

JHEP07(2016)133

Figure 6. The normalised double-di erential transverse momentum distribution for the Z boson in windows of invariant mass of the leptons, m[lscript][lscript], with a rapidity cut on the Z boson of [notdef]yZ[notdef] < 2.4.

The ATLAS data is taken from ref. [38]. The green bands denote the NLO prediction with scale uncertainty and the blue bands show the NNLO prediction w ith scale uncertainty.

appropriate dilepton invariant mass cut. We computed these ducial cross sections to NNLO accuracy (O( 2s)) and compare them to the measured values quoted by ATLAS [38]

in gure 5. We clearly observe that the central value of the NNLO prediction falls below the experimental data in all mass bins. The large scale dependence in the three lowest mass bins can be explained from the interplay between the cuts on the lepton pair invariant mass and on the single lepton transverse momentum. The latter cuts forbid events at low transverse momentum of the lepton pair down to m[lscript][lscript] = 40 GeV, such that the two lowest bins receive no leading order contribution, and the third bin only a very small one. All three bins are populated only from NLO onwards, with events containing a low-mass lepton pair recoiling against a parton at high transverse momentum. So our NNLO prediction for the inclusive cross section in these mass bins is e ectively only NLO accurate, with consequently larger scale dependence. In the three bins with larger m[lscript][lscript], the scale uncertainty on the NNLO prediction is below 0.7%, which results in tension between data and theory at the level of two standard deviations.

Combining together the unnormalised di erential distribution with the inclusive cross sections, we obtain the normalised distributions shown in gure 6. Because of the large scale uncertainty in the inclusive cross section, the theoretical errors dominate the low m[lscript][lscript]

bins. At large m[lscript][lscript], the tension between the data and NNLO theory is largely relieved. At

{ 8 {

JHEP07(2016)133

Figure 7. The unnormalised double-di erential transverse momentum distribution for the Z boson in windows of rapidity of the Z boson, yZ, with an invariant mass cut on the nal state leptons of 66 GeV < m[lscript][lscript] < 116 GeV. The ATLAS data is taken from ref. [38]. The luminosity error is not shown. The green bands denote the NLO prediction with scale uncertainty and the blue bands show the NNLO prediction.

the highest values of pZT, the tendency of the data to fall below the theory prediction may be an indication of the onset of electroweak corrections [28{31], which are negative in this region. Any remaining tension for medium values of pZT could potentially be accounted for revisiting the parton distribution functions (especially the gluon distribution) in the kinematical region relevant to this measurement.

The same tension between NNLO theory and ATLAS data for the unnormalised distribution is visible in gure 7, which shows the unnormalised double-di erential distribution with respect to the transverse momentum of the Z boson for 66 GeV < m[lscript][lscript] < 116 GeV in

di erent ranges of the rapidity of the Z boson, normalised to the NLO prediction. The NNLO corrections are uniform in rapidity and transverse momentum at the level of about 5{8%, and they decrease the residual theoretical uncertainty to 2{4%. The data, which correspond to the invariant mass bin containing the Z resonance, are consistently above the NNLO prediction.

In order to eliminate the overall luminosity uncertainty, one again normalises the data to the inclusive ducial dilepton cross section in the respective bins in yZ. Figure 8 shows the ratio of ATLAS data and the NNLO cross section compared to the NLO prediction.1

1Note that the ATLAS data is extracted from ref. [38] by summing the bins of the pZT distributions in the various yZ bins. We have checked that this procedure applied to the di erent m[lscript][lscript] bins shown in gure 5 correctly reproduces the data.

{ 9 {

JHEP07(2016)133

Figure 8. The inclusive dilepton cross section for the same [notdef]yZ[notdef] bins as in gure 7 and an invariant

mass cut on nal state leptons of 66 GeV < m[lscript][lscript] < 116 GeV. The ATLAS data is extracted from ref. [38] by summing up the transverse momentum distributions in the respective bins. The vertical error bars are given by the luminosity uncertainty. The blue bands show the NNLO prediction with scale uncertainty.

Figure 9. The normalised double-di erential transverse momentum distribution for the Z boson in windows of rapidity of the Z boson, yZ, with an invariant mass cut on nal state leptons of66 GeV < m[lscript][lscript] < 116 GeV. The ATLAS data is taken from ref. [38]. The green bands denote the

NLO prediction with scale uncertainty and the blue bands show the NNLO prediction.

{ 10 {

JHEP07(2016)133

Figure 10. The unnormalised double-di erential transverse momentum distribution for the Z boson in windows of rapidity of the Z boson, yZ, with an invariant mass cut on nal state leptons of 81 GeV < m[lscript][lscript] < 101 GeV. The CMS data is taken from ref. [40]. The experimental uncertainties include the luminosity uncertainty. The green bands denote the NLO prediction with scale uncertainty and the blue bands show the NNLO prediction.

As for the rapidity-integrated cross section in this mass bin, gure 5, we observe the NNLO prediction to be accurate to 0.9%, and to underestimate the ATLAS data by about two standard deviations. The normalised double-di erential distribution is shown in gure 9. There is excellent agreement between the normalised NNLO prediction and the ATLAS data for all rapidity bins, again o ering the possibility of further constraining the PDFs in this kinematic region.

The CMS measurement of the Z-boson transverse momentum distribution at 8 TeV [40] concentrates on the dilepton invariant mass range 81 GeV < m[lscript][lscript] < 101 GeV around the Z resonance. The ducial region of this measurement, de ned in table 1, applies asymmetric cuts on the leading and sub-leading lepton, disregarding the lepton charge. To ensure the correct implementation of these cuts in NNLOJET, and to maintain local cancellations between matrix elements and subtraction terms in the antenna subtraction method, the lepton identi cation as leading and sub-leading had to be passed alongside the lepton charge through all subtraction terms and into the event analysis routines. The double-di erential measurement is performed in pZT and yZ. Comparing the absolute double-di erential distributions of gure 10 to our NNLO predictions, we observe the same features as for the ATLAS measurement, with positive NNLO corrections at the level of 5{8% and

{ 11 {

JHEP07(2016)133

Figure 11. The normalised double-di erential transverse momentum distribution for the Z boson in windows of rapidity of the Z boson, yZ, with an invariant mass cut on nal state leptons of81 GeV < m[lscript][lscript] < 101 GeV. The CMS data is taken from ref. [40]. The green bands denote the NLO prediction with scale uncertainty and the blue bands show the NNLO prediction.

an NNLO scale uncertainty of 2{4%. Compared to NLO, inclusion of NNLO corrections brings the theoretical prediction closer to the experimental data, which are however still about 5-8% larger than expected from theory.

In the CMS study [40], the normalised cross section is obtained by dividing the distributions by the ducial cross section integrated over all bins in yZ (in contrast to ATLAS, which normalised to the ducial cross sections only integrated over the respective bin in yZ). At NNLO, we obtain for the ducial cross section with CMS cuts:

NNLO(81 GeV < m[lscript][lscript] < 101 GeV) = 450.6+2.71.6 pb.

The normalised distributions of CMS are compared to theory in gure 11, where excellent agreement is observed upon inclusion of the NNLO corrections.

Compared to NNLO theory, both ATLAS and CMS ducial cross section measurements of the Z-boson transverse momentum display a similar pattern of disagreement for the absolute distributions while having excellent agreement for the normalised distributions. The inclusion of the newly derived NNLO corrections to the transverse momentum distribution are crucial for a meaningful comparison between data and theory: (a) they reduce the theory uncertainty to a level that rmly establishes the discrepancy on the absolute cross sections, and (b) they modify the central value and the shape of the theory prediction to better agree with the data on the normalised distributions.

{ 12 {

5 Summary and conclusions

We have derived the NNLO QCD corrections to the production of Z/ bosons decaying to lepton pairs at large transverse momentum, inclusive over the hadronic nal state. This observable has typically been very easy to measure but di cult to predict theoretically because of the necessity to be fully inclusive with respect to all QCD radiation. Our calculation is performed using the parton-level Monte Carlo generator NNLOJET which

implements the antenna subtraction method for NNLO calculations of hadron collider observables. It extends our earlier calculation of Z/ +jet production at this order [32, 33]

and now also includes inclusive Z/ production. We have performed a thorough comparison of theory predictions to the 8 TeV Run 1 data of the LHC for cross sections de ned over a ducial region of lepton kinematics from the ATLAS [38] and CMS [40] collaborations. We observe that the NNLO corrections to the unnormalised distributions are moderate and positive, resulting in theory predictions that are closer to the experimental observations than at NLO. However, this improvement is much more dramatic when one considers distributions normalised to the relevant dilepton cross section. This is true for both single and double-di erential distributions. While the normalised distributions for low m[lscript][lscript] slices in comparison to ATLAS are dominated by the theoretical scale uncertainty on the inclusive dilepton cross section, for larger m[lscript][lscript] and all [notdef]yZ[notdef] bins the agreement between NNLO

theory and LHC data is excellent. With the reduction of the theoretical scale uncertainty at NNLO, a careful reinvestigation of the parametric ingredients to the NNLO theory predictions (parton distributions, strong coupling constant, electroweak parameters) appears now to be mandatory. Our calculation enables exactly these precision studies, and will allow a consistent inclusion of the precision data on the Z transverse momentum distribution into NNLO determinations of parton distributions, leading to precision constraints at this order on the gluon distribution over a large momentum range.

Acknowledgments

We would like to thank Daniel Froidevaux, Claire Gwenlan and Elzbieta Richter-Was for stimulating discussions concerning the ATLAS data, Dieter Zeppenfeld for enlightening comments that helped to shape the phenomenological study in this paper and Alexander Muck for helpful comments concerning the impact of EW corrections to these observables. We also thank Xuan Chen, Juan Cruz-Martinez, James Currie, Jan Niehues and Joo Pires for useful discussions and their many contributions to the NNLOJET code. The authors would like to acknowledge the support provided by the GridPP Collaboration. This research was supported in part by the National Science Foundation under Grant NSF PHY11-25915, in part by the Swiss National Science Foundation (SNF) under contracts 200020-162487 and CRSII2-160814, in part by the U.K. Science and Technology Facilities Council, in part by the Research Executive Agency (REA) of the European Union under the Grant Agreement PITN-GA-2012-316704 (\HiggsTools") and the ERC Advanced Grant MC@NNLO (340983).

{ 13 {

JHEP07(2016)133

Open Access. This article is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/

Web End =CC-BY 4.0 ), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.

References

[1] R. Hamberg, W.L. van Neerven and T. Matsuura, A complete calculation of the order 2s

correction to the Drell-Yan K factor, http://dx.doi.org/10.1016/S0550-3213(02)00814-3

Web End =Nucl. Phys. B 359 (1991) 343 [Erratum ibid. B 644 (2002) 403] [http://inspirehep.net/search?p=find+J+%22Nucl.Phys.,B359,343%22

Web End =INSPIRE ].

[2] W.L. van Neerven and E.B. Zijlstra, The O( 2s) corrected Drell-Yan K factor in the DIS and MS scheme, http://dx.doi.org/10.1016/j.nuclphysb.2003.12.019

Web End =Nucl. Phys. B 382 (1992) 11 [Erratum ibid. B 680 (2004) 513] [http://inspirehep.net/search?p=find+J+%22Nucl.Phys.,B382,11%22

Web End =INSPIRE ].

[3] C. Anastasiou, L.J. Dixon, K. Melnikov and F. Petriello, Dilepton rapidity distribution in the Drell-Yan process at NNLO in QCD, http://dx.doi.org/10.1103/PhysRevLett.91.182002

Web End =Phys. Rev. Lett. 91 (2003) 182002 [http://arxiv.org/abs/hep-ph/0306192

Web End =hep-ph/0306192 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0306192

Web End =INSPIRE ].

[4] K. Melnikov and F. Petriello, The W boson production cross section at the LHC through O( 2s), http://dx.doi.org/10.1103/PhysRevLett.96.231803

Web End =Phys. Rev. Lett. 96 (2006) 231803 [http://arxiv.org/abs/hep-ph/0603182

Web End =hep-ph/0603182 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0603182

Web End =INSPIRE ].

[5] K. Melnikov and F. Petriello, Electroweak gauge boson production at hadron colliders through O( 2s), http://dx.doi.org/10.1103/PhysRevD.74.114017

Web End =Phys. Rev. D 74 (2006) 114017 [http://arxiv.org/abs/hep-ph/0609070

Web End =hep-ph/0609070 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0609070

Web End =INSPIRE ].

[6] S. Catani, G. Ferrera and M. Grazzini, W boson production at hadron colliders: the lepton charge asymmetry in NNLO QCD, http://dx.doi.org/10.1007/JHEP05(2010)006

Web End =JHEP 05 (2010) 006 [http://arxiv.org/abs/1002.3115

Web End =arXiv:1002.3115 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1002.3115

Web End =INSPIRE ].

[7] S. Catani, L. Cieri, G. Ferrera, D. de Florian and M. Grazzini, Vector boson production at hadron colliders: a fully exclusive QCD calculation at NNLO, http://dx.doi.org/10.1103/PhysRevLett.103.082001

Web End =Phys. Rev. Lett. 103 (2009) http://dx.doi.org/10.1103/PhysRevLett.103.082001

Web End =082001 [http://arxiv.org/abs/0903.2120

Web End =arXiv:0903.2120 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0903.2120

Web End =INSPIRE ].

[8] R. Gavin, Y. Li, F. Petriello and S. Quackenbush, FEWZ 2.0: a code for hadronic Z production at next-to-next-to-leading order, http://dx.doi.org/10.1016/j.cpc.2011.06.008

Web End =Comput. Phys. Commun. 182 (2011) 2388 [http://arxiv.org/abs/1011.3540

Web End =arXiv:1011.3540 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1011.3540

Web End =INSPIRE ].

[9] C. Anastasiou, L.J. Dixon, K. Melnikov and F. Petriello, High precision QCD at hadron colliders: Electroweak gauge boson rapidity distributions at NNLO, http://dx.doi.org/10.1103/PhysRevD.69.094008

Web End =Phys. Rev. D 69 (2004) http://dx.doi.org/10.1103/PhysRevD.69.094008

Web End =094008 [http://arxiv.org/abs/hep-ph/0312266

Web End =hep-ph/0312266 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0312266

Web End =INSPIRE ].

[10] S. Moch and A. Vogt, Higher-order soft corrections to lepton pair and Higgs boson production, http://dx.doi.org/10.1016/j.physletb.2005.09.061

Web End =Phys. Lett. B 631 (2005) 48 [http://arxiv.org/abs/hep-ph/0508265

Web End =hep-ph/0508265 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0508265

Web End =INSPIRE ].

[11] E. Laenen and L. Magnea, Threshold resummation for electroweak annihilation from DIS data, http://dx.doi.org/10.1016/j.physletb.2005.10.038

Web End =Phys. Lett. B 632 (2006) 270 [http://arxiv.org/abs/hep-ph/0508284

Web End =hep-ph/0508284 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0508284

Web End =INSPIRE ].

[12] A. Idilbi, X.-d. Ji and F. Yuan, Resummation of threshold logarithms in e ective eld theory for DIS, Drell-Yan and Higgs production, http://dx.doi.org/10.1016/j.nuclphysb.2006.07.002

Web End =Nucl. Phys. B 753 (2006) 42 [http://arxiv.org/abs/hep-ph/0605068

Web End =hep-ph/0605068 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0605068

Web End =INSPIRE ].

[13] V. Ravindran and J. Smith, Threshold corrections to rapidity distributions of Z and W bosons beyond N2 LO at hadron colliders, http://dx.doi.org/10.1103/PhysRevD.76.114004

Web End =Phys. Rev. D 76 (2007) 114004 [http://arxiv.org/abs/0708.1689

Web End =arXiv:0708.1689 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0708.1689

Web End =INSPIRE ].

[14] T. Ahmed, M. Mahakhud, N. Rana and V. Ravindran, Drell-Yan production at threshold to third order in QCD, http://dx.doi.org/10.1103/PhysRevLett.113.112002

Web End =Phys. Rev. Lett. 113 (2014) 112002 [http://arxiv.org/abs/1404.0366

Web End =arXiv:1404.0366 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1404.0366

Web End =INSPIRE ].

{ 14 {

JHEP07(2016)133

[15] D. Bonocore, E. Laenen, L. Magnea, L. Vernazza and C.D. White, The method of regions and next-to-soft corrections in Drell-Yan production, http://dx.doi.org/10.1016/j.physletb.2015.02.008

Web End =Phys. Lett. B 742 (2015) 375 [http://arxiv.org/abs/1410.6406

Web End =arXiv:1410.6406 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1410.6406

Web End =INSPIRE ].

[16] J.C. Collins, D.E. Soper and G.F. Sterman, Transverse momentum distribution in Drell-Yan pair and W and Z boson production, http://dx.doi.org/10.1016/0550-3213(85)90479-1

Web End =Nucl. Phys. B 250 (1985) 199 [http://inspirehep.net/search?p=find+J+%22Nucl.Phys.,B250,199%22

Web End =INSPIRE ].

[17] G. Bozzi, S. Catani, G. Ferrera, D. de Florian and M. Grazzini, Transverse-momentum resummation: a perturbative study of Z production at the Tevatron, http://dx.doi.org/10.1016/j.nuclphysb.2009.02.014

Web End =Nucl. Phys. B 815 http://dx.doi.org/10.1016/j.nuclphysb.2009.02.014

Web End =(2009) 174 [http://arxiv.org/abs/0812.2862

Web End =arXiv:0812.2862 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0812.2862

Web End =INSPIRE ].

[18] G. Bozzi, S. Catani, G. Ferrera, D. de Florian and M. Grazzini, Production of Drell-Yan lepton pairs in hadron collisions: Transverse-momentum resummation at next-to-next-to-leading logarithmic accuracy, http://dx.doi.org/10.1016/j.physletb.2010.12.024

Web End =Phys. Lett. B 696 (2011) 207 [http://arxiv.org/abs/1007.2351

Web End =arXiv:1007.2351 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1007.2351

Web End =INSPIRE ].

[19] T. Becher and M. Neubert, Drell-Yan production at small qT , transverse parton distributions and the collinear anomaly, http://dx.doi.org/10.1140/epjc/s10052-011-1665-7

Web End =Eur. Phys. J. C 71 (2011) 1665 [http://arxiv.org/abs/1007.4005

Web End =arXiv:1007.4005 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1007.4005

Web End =INSPIRE ].

[20] A. Karlberg, E. Re and G. Zanderighi, NNLOPS accurate Drell-Yan production, http://dx.doi.org/10.1007/JHEP09(2014)134

Web End =JHEP 09 http://dx.doi.org/10.1007/JHEP09(2014)134

Web End =(2014) 134 [http://arxiv.org/abs/1407.2940

Web End =arXiv:1407.2940 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1407.2940

Web End =INSPIRE ].

[21] S. Hoche, Y. Li and S. Prestel, Drell-Yan lepton pair production at NNLO QCD with parton showers, http://dx.doi.org/10.1103/PhysRevD.91.074015

Web End =Phys. Rev. D 91 (2015) 074015 [http://arxiv.org/abs/1405.3607

Web End =arXiv:1405.3607 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1405.3607

Web End =INSPIRE ].

[22] U. Baur, O. Brein, W. Hollik, C. Schappacher and D. Wackeroth, Electroweak radiative corrections to neutral current Drell-Yan processes at hadron colliders, http://dx.doi.org/10.1103/PhysRevD.65.033007

Web End =Phys. Rev. D 65 http://dx.doi.org/10.1103/PhysRevD.65.033007

Web End =(2002) 033007 [http://arxiv.org/abs/hep-ph/0108274

Web End =hep-ph/0108274 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0108274

Web End =INSPIRE ].

[23] V.A. Zykunov, Weak radiative corrections to Drell-Yan process for large invariant mass of di-lepton pair, http://dx.doi.org/10.1103/PhysRevD.75.073019

Web End =Phys. Rev. D 75 (2007) 073019 [http://arxiv.org/abs/hep-ph/0509315

Web End =hep-ph/0509315 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0509315

Web End =INSPIRE ].

[24] C.M. Carloni Calame, G. Montagna, O. Nicrosini and A. Vicini, Precision electroweak calculation of the production of a high transverse-momentum lepton pair at hadron colliders, http://dx.doi.org/10.1088/1126-6708/2007/10/109

Web End =JHEP 10 (2007) 109 [http://arxiv.org/abs/0710.1722

Web End =arXiv:0710.1722 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0710.1722

Web End =INSPIRE ].

[25] L. Barze, G. Montagna, P. Nason, O. Nicrosini, F. Piccinini and A. Vicini, Neutral current Drell-Yan with combined QCD and electroweak corrections in the POWHEG BOX, http://dx.doi.org/10.1140/epjc/s10052-013-2474-y

Web End =Eur. http://dx.doi.org/10.1140/epjc/s10052-013-2474-y

Web End =Phys. J. C 73 (2013) 2474 [http://arxiv.org/abs/1302.4606

Web End =arXiv:1302.4606 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1302.4606

Web End =INSPIRE ].

[26] S. Dittmaier, A. Huss and C. Schwinn, Mixed QCD-electroweak O( s ) corrections to Drell-Yan processes in the resonance region: pole approximation and non-factorizable corrections, http://dx.doi.org/10.1016/j.nuclphysb.2014.05.027

Web End =Nucl. Phys. B 885 (2014) 318 [http://arxiv.org/abs/1403.3216

Web End =arXiv:1403.3216 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1403.3216

Web End =INSPIRE ].

[27] S. Dittmaier, A. Huss and C. Schwinn, Dominant mixed QCD-electroweak O( s )

corrections to Drell-Yan processes in the resonance region, http://dx.doi.org/10.1016/j.nuclphysb.2016.01.006

Web End =Nucl. Phys. B 904 (2016) 216 [http://arxiv.org/abs/1511.08016

Web End =arXiv:1511.08016 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1511.08016

Web End =INSPIRE ].

[28] J.H. Kuhn, A. Kulesza, S. Pozzorini and M. Schulze, One-loop weak corrections to hadronic production of Z bosons at large transverse momenta, http://dx.doi.org/10.1016/j.nuclphysb.2005.08.019

Web End =Nucl. Phys. B 727 (2005) 368 [http://arxiv.org/abs/hep-ph/0507178

Web End =hep-ph/0507178 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0507178

Web End =INSPIRE ].

[29] A. Denner, S. Dittmaier, T. Kasprzik and A. Muck, Electroweak corrections to dilepton + jet production at hadron colliders, http://dx.doi.org/10.1007/JHEP06(2011)069

Web End =JHEP 06 (2011) 069 [http://arxiv.org/abs/1103.0914

Web End =arXiv:1103.0914 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1103.0914

Web End =INSPIRE ].

[30] W. Hollik, B.A. Kniehl, E.S. Scherbakova and O.L. Veretin, Electroweak corrections to Z-boson hadroproduction at nite transverse momentum, http://dx.doi.org/10.1016/j.nuclphysb.2015.09.018

Web End =Nucl. Phys. B 900 (2015) 576 [http://arxiv.org/abs/1504.07574

Web End =arXiv:1504.07574 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1504.07574

Web End =INSPIRE ].

{ 15 {

JHEP07(2016)133

[31] S. Kallweit, J.M. Lindert, P. Maierhofer, S. Pozzorini and M. Schonherr, NLO QCD+EW predictions for V + jets including o -shell vector-boson decays and multijet merging, http://dx.doi.org/10.1007/JHEP04(2016)021

Web End =JHEP http://dx.doi.org/10.1007/JHEP04(2016)021

Web End =04 (2016) 021 [http://arxiv.org/abs/1511.08692

Web End =arXiv:1511.08692 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1511.08692

Web End =INSPIRE ].

[32] A. Gehrmann-De Ridder, T. Gehrmann, E.W.N. Glover, A. Huss and T.A. Morgan, Precise QCD predictions for the production of a Z boson in association with a hadronic jet, http://dx.doi.org/10.1103/PhysRevLett.117.022001

Web End =Phys. http://dx.doi.org/10.1103/PhysRevLett.117.022001

Web End =Rev. Lett. 117 (2016) 022001 [http://arxiv.org/abs/1507.02850

Web End =arXiv:1507.02850 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1507.02850

Web End =INSPIRE ].

[33] A. Gehrmann-De Ridder, T. Gehrmann, E.W.N. Glover, A. Huss and T.A. Morgan, NNLO QCD corrections for Z boson plus jet production, http://arxiv.org/abs/1601.04569

Web End =arXiv:1601.04569 [http://inspirehep.net/search?p=find+EPRINT+arXiv:1601.04569

Web End =INSPIRE ].

[34] R. Boughezal et al., Z-boson production in association with a jet at next-to-next-to-leading order in perturbative QCD, http://dx.doi.org/10.1103/PhysRevLett.116.152001

Web End =Phys. Rev. Lett. 116 (2016) 152001 [http://arxiv.org/abs/1512.01291

Web End =arXiv:1512.01291 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1512.01291

Web End =INSPIRE ].

[35] R. Boughezal, X. Liu and F. Petriello, A comparison of NNLO QCD predictions with 7 TeV ATLAS and CMS data for V +jet processes, http://dx.doi.org/10.1016/j.physletb.2016.06.032

Web End =Phys. Lett. B 760 (2016) 6 [http://arxiv.org/abs/1602.05612

Web End =arXiv:1602.05612 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1602.05612

Web End =INSPIRE ].

[36] R. Boughezal, X. Liu and F. Petriello, Phenomenology of the Z-boson plus jet process at NNLO, http://arxiv.org/abs/1602.08140

Web End =arXiv:1602.08140 [http://inspirehep.net/search?p=find+EPRINT+arXiv:1602.08140

Web End =INSPIRE ].

[37] ATLAS collaboration, Measurement of the Z/ boson transverse momentum distribution in pp collisions at ps = 7 TeV with the ATLAS detector, http://dx.doi.org/10.1007/JHEP09(2014)145

Web End =JHEP 09 (2014) 145 [http://arxiv.org/abs/1406.3660

Web End =arXiv:1406.3660 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1406.3660

Web End =INSPIRE ].

[38] ATLAS collaboration, Measurement of the transverse momentum and distributions of Drell-Yan lepton pairs in proton-proton collisions at ps = 8 TeV with the ATLAS detector,

http://dx.doi.org/10.1140/epjc/s10052-016-4070-4

Web End =Eur. Phys. J. C 76 (2016) 291 [http://arxiv.org/abs/1512.02192

Web End =arXiv:1512.02192 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1512.02192

Web End =INSPIRE ].

[39] CMS collaboration, Measurement of the rapidity and transverse momentum distributions of Z bosons in pp collisions at ps = 7 TeV, http://dx.doi.org/10.1103/PhysRevD.85.032002

Web End =Phys. Rev. D 85 (2012) 032002 [http://arxiv.org/abs/1110.4973

Web End =arXiv:1110.4973 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1110.4973

Web End =INSPIRE ].

[40] CMS collaboration, Measurement of the Z boson di erential cross section in transverse momentum and rapidity in proton-proton collisions at 8 TeV, http://dx.doi.org/10.1016/j.physletb.2015.07.065

Web End =Phys. Lett. B 749 (2015) 187 [http://arxiv.org/abs/1504.03511

Web End =arXiv:1504.03511 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1504.03511

Web End =INSPIRE ].

[41] LHCb collaboration, Measurement of the forward Z boson production cross-section in pp collisions at ps = 7 TeV, http://dx.doi.org/10.1007/JHEP08(2015)039

Web End =JHEP 08 (2015) 039 [http://arxiv.org/abs/1505.07024

Web End =arXiv:1505.07024 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1505.07024

Web End =INSPIRE ].

[42] LHCb collaboration, Measurement of forward W and Z boson production in pp collisions at ps = 8 TeV, http://dx.doi.org/10.1007/JHEP01(2016)155

Web End =JHEP 01 (2016) 155 [http://arxiv.org/abs/1511.08039

Web End =arXiv:1511.08039 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1511.08039

Web End =INSPIRE ].

[43] L.W. Garland, T. Gehrmann, E.W.N. Glover, A. Koukoutsakis and E. Remiddi, The two loop QCD matrix element for e+e ! 3 jets, http://dx.doi.org/10.1016/S0550-3213(02)00057-3

Web End =Nucl. Phys. B 627 (2002) 107

[http://arxiv.org/abs/hep-ph/0112081

Web End =hep-ph/0112081 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0112081

Web End =INSPIRE ].

[44] L.W. Garland, T. Gehrmann, E.W.N. Glover, A. Koukoutsakis and E. Remiddi, Two loop QCD helicity amplitudes for e+e ! three jets, http://dx.doi.org/10.1016/S0550-3213(02)00627-2

Web End =Nucl. Phys. B 642 (2002) 227

[http://arxiv.org/abs/hep-ph/0206067

Web End =hep-ph/0206067 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0206067

Web End =INSPIRE ].

[45] T. Gehrmann and L. Tancredi, Two-loop QCD helicity amplitudes for q

q ! W and

q

q ! Z0 , http://dx.doi.org/10.1007/JHEP02(2012)004

Web End =JHEP 02 (2012) 004 [http://arxiv.org/abs/1112.1531

Web End =arXiv:1112.1531 ] [

http://inspirehep.net/search?p=find+EPRINT+arXiv:1112.1531

Web End =INSPIRE ].

[46] E.W.N. Glover and D.J. Miller, The One loop QCD corrections for ! Q

http://dx.doi.org/10.1016/S0370-2693(97)00113-5

Web End =B 396 (1997) 257 [http://arxiv.org/abs/hep-ph/9609474

Web End =hep-ph/9609474 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/9609474

Web End =INSPIRE ].

{ 16 {

JHEP07(2016)133

q, http://dx.doi.org/10.1016/S0370-2693(97)00113-5

Web End =Phys. Lett.

[47] Z. Bern, L.J. Dixon, D.A. Kosower and S. Weinzierl, One loop amplitudes for e+e !

http://dx.doi.org/10.1016/S0550-3213(96)00703-1

Web End =Nucl. Phys. B 489 (1997) 3 [http://arxiv.org/abs/hep-ph/9610370

Web End =hep-ph/9610370 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/9610370

Web End =INSPIRE ].

[48] J.M. Campbell, E.W.N. Glover and D.J. Miller, The one loop QCD corrections for ! q

qgg, http://dx.doi.org/10.1016/S0370-2693(97)00909-X

Web End =Phys. Lett. B 409 (1997) 503 [http://arxiv.org/abs/hep-ph/9706297

Web End =hep-ph/9706297 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/9706297

Web End =INSPIRE ].

[49] Z. Bern, L.J. Dixon and D.A. Kosower, One loop amplitudes for e+e to four partons, http://dx.doi.org/10.1016/S0550-3213(97)00703-7

Web End =Nucl. http://dx.doi.org/10.1016/S0550-3213(97)00703-7

Web End =Phys. B 513 (1998) 3 [http://arxiv.org/abs/hep-ph/9708239

Web End =hep-ph/9708239 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/9708239

Web End =INSPIRE ].

[50] J.M. Campbell and R.K. Ellis, Next-to-leading order corrections to W + 2 jet and Z + 2 jet production at hadron colliders, http://dx.doi.org/10.1103/PhysRevD.65.113007

Web End =Phys. Rev. D 65 (2002) 113007 [http://arxiv.org/abs/hep-ph/0202176

Web End =hep-ph/0202176 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0202176

Web End =INSPIRE ].

[51] J.M. Campbell, R.K. Ellis and D.L. Rainwater, Next-to-leading order QCD predictions for W + 2 jet and Z + 2 jet production at the CERN LHC, http://dx.doi.org/10.1103/PhysRevD.68.094021

Web End =Phys. Rev. D 68 (2003) 094021 [http://arxiv.org/abs/hep-ph/0308195

Web End =hep-ph/0308195 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0308195

Web End =INSPIRE ].

[52] K. Hagiwara and D. Zeppenfeld, Amplitudes for multiparton processes involving a current at e+e, ep and hadron colliders, http://dx.doi.org/10.1016/0550-3213(89)90397-0

Web End =Nucl. Phys. B 313 (1989) 560 [http://inspirehep.net/search?p=find+J+%22Nucl.Phys.,B313,560%22

Web End =INSPIRE ].

[53] F.A. Berends, W.T. Giele and H. Kuijf, Exact expressions for processes involving a vector boson and up to ve partons, http://dx.doi.org/10.1016/0550-3213(89)90242-3

Web End =Nucl. Phys. B 321 (1989) 39 [http://inspirehep.net/search?p=find+J+%22Nucl.Phys.,B321,39%22

Web End =INSPIRE ].

[54] N.K. Falck, D. Graudenz and G. Kramer, Cross-section for ve jet production in e+e

annihilation, http://dx.doi.org/10.1016/0550-3213(89)90331-3

Web End =Nucl. Phys. B 328 (1989) 317 [http://inspirehep.net/search?p=find+J+%22Nucl.Phys.,B328,317%22

Web End =INSPIRE ].

[55] A. Gehrmann-De Ridder, T. Gehrmann and E.W.N. Glover, Antenna subtraction at NNLO, http://dx.doi.org/10.1088/1126-6708/2005/09/056

Web End =JHEP 09 (2005) 056 [http://arxiv.org/abs/hep-ph/0505111

Web End =hep-ph/0505111 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0505111

Web End =INSPIRE ].

[56] A. Gehrmann-De Ridder, T. Gehrmann and E.W.N. Glover, Gluon-gluon antenna functions from Higgs boson decay, http://dx.doi.org/10.1016/j.physletb.2005.03.003

Web End =Phys. Lett. B 612 (2005) 49 [http://arxiv.org/abs/hep-ph/0502110

Web End =hep-ph/0502110 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0502110

Web End =INSPIRE ].

[57] A. Gehrmann-De Ridder, T. Gehrmann and E.W.N. Glover, quark-gluon antenna functions from neutralino decay, http://dx.doi.org/10.1016/j.physletb.2005.02.039

Web End =Phys. Lett. B 612 (2005) 36 [http://arxiv.org/abs/hep-ph/0501291

Web End =hep-ph/0501291 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0501291

Web End =INSPIRE ].

[58] A. Daleo, T. Gehrmann and D. Ma^ tre, Antenna subtraction with hadronic initial states, http://dx.doi.org/10.1088/1126-6708/2007/04/016

Web End =JHEP 04 (2007) 016 [http://arxiv.org/abs/hep-ph/0612257

Web End =hep-ph/0612257 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0612257

Web End =INSPIRE ].

[59] A. Daleo, A. Gehrmann-De Ridder, T. Gehrmann and G. Luisoni, Antenna subtraction at NNLO with hadronic initial states: initial- nal con gurations, http://dx.doi.org/10.1007/JHEP01(2010)118

Web End =JHEP 01 (2010) 118 [http://arxiv.org/abs/0912.0374

Web End =arXiv:0912.0374 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0912.0374

Web End =INSPIRE ].

[60] T. Gehrmann and P.F. Monni, Antenna subtraction at NNLO with hadronic initial states: real-virtual initial-initial con gurations, http://dx.doi.org/10.1007/JHEP12(2011)049

Web End =JHEP 12 (2011) 049 [http://arxiv.org/abs/1107.4037

Web End =arXiv:1107.4037 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1107.4037

Web End =INSPIRE ].

[61] R. Boughezal, A. Gehrmann-De Ridder and M. Ritzmann, Antenna subtraction at NNLO with hadronic initial states: double real radiation for initial-initial con gurations with two quark avours, http://dx.doi.org/10.1007/JHEP02(2011)098

Web End =JHEP 02 (2011) 098 [http://arxiv.org/abs/1011.6631

Web End =arXiv:1011.6631 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1011.6631

Web End =INSPIRE ].

[62] A. Gehrmann-De Ridder, T. Gehrmann and M. Ritzmann, Antenna subtraction at NNLO with hadronic initial states: double real initial-initial con gurations, http://dx.doi.org/10.1007/JHEP10(2012)047

Web End =JHEP 10 (2012) 047 [http://arxiv.org/abs/1207.5779

Web End =arXiv:1207.5779 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1207.5779

Web End =INSPIRE ].

[63] J. Currie, E.W.N. Glover and S. Wells, Infrared structure at NNLO using antenna subtraction, http://dx.doi.org/10.1007/JHEP04(2013)066

Web End =JHEP 04 (2013) 066 [http://arxiv.org/abs/1301.4693

Web End =arXiv:1301.4693 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1301.4693

Web End =INSPIRE ].

[64] A. Gehrmann-De Ridder, T. Gehrmann, E.W.N. Glover and G. Heinrich, Infrared structure of e+e ! 3 jets at NNLO, http://dx.doi.org/10.1088/1126-6708/2007/11/058

Web End =JHEP 11 (2007) 058 [http://arxiv.org/abs/0710.0346

Web End =arXiv:0710.0346 ] [

{ 17 {

qq

QQ,

JHEP07(2016)133

http://inspirehep.net/search?p=find+EPRINT+arXiv:0710.0346

Web End =INSPIRE ].

[65] A. Gehrmann-De Ridder, T. Gehrmann, E.W.N. Glover and G. Heinrich, EERAD3: Event shapes and jet rates in electron-positron annihilation at order 3s, http://dx.doi.org/10.1016/j.cpc.2014.07.024

Web End =Comput. Phys. Commun.

http://dx.doi.org/10.1016/j.cpc.2014.07.024

Web End =185 (2014) 3331 [http://arxiv.org/abs/1402.4140

Web End =arXiv:1402.4140 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1402.4140

Web End =INSPIRE ].

[66] A. Gehrmann-De Ridder, T. Gehrmann, E.W.N. Glover and J. Pires, Second order QCD corrections to jet production at hadron colliders: the all-gluon contribution, http://dx.doi.org/10.1103/PhysRevLett.110.162003

Web End =Phys. Rev. Lett. http://dx.doi.org/10.1103/PhysRevLett.110.162003

Web End =110 (2013) 162003 [http://arxiv.org/abs/1301.7310

Web End =arXiv:1301.7310 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1301.7310

Web End =INSPIRE ].

[67] J. Currie, A. Gehrmann-De Ridder, E.W.N. Glover and J. Pires, NNLO QCD corrections to jet production at hadron colliders from gluon scattering, http://dx.doi.org/10.1007/JHEP01(2014)110

Web End =JHEP 01 (2014) 110 [http://arxiv.org/abs/1310.3993

Web End =arXiv:1310.3993 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1310.3993

Web End =INSPIRE ].

[68] G. Abelof, A. Gehrmann-De Ridder and I. Majer, Top quark pair production at NNLO in the quark-antiquark channel, http://dx.doi.org/10.1007/JHEP12(2015)074

Web End =JHEP 12 (2015) 074 [http://arxiv.org/abs/1506.04037

Web End =arXiv:1506.04037 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1506.04037

Web End =INSPIRE ].

[69] X. Chen, T. Gehrmann, E.W.N. Glover and M. Jaquier, Precise QCD predictions for the production of Higgs + jet nal states, http://dx.doi.org/10.1016/j.physletb.2014.11.021

Web End =Phys. Lett. B 740 (2015) 147 [http://arxiv.org/abs/1408.5325

Web End =arXiv:1408.5325 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1408.5325

Web End =INSPIRE ].

[70] NNPDF collaboration, R.D. Ball et al., Parton distributions for the LHC Run II, http://dx.doi.org/10.1007/JHEP04(2015)040

Web End =JHEP 04 http://dx.doi.org/10.1007/JHEP04(2015)040

Web End =(2015) 040 [http://arxiv.org/abs/1410.8849

Web End =arXiv:1410.8849 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1410.8849

Web End =INSPIRE ].

JHEP07(2016)133

{ 18 {