[0.4cm] for vector-boson pair production at NNLL+NNLO
Massimiliano Grazzini***On leave of absence from INFN, Sezione di Firenze, Sesto Fiorentino, Florence, Italy., Stefan Kallweit,
[0.2cm] Dirk Rathlev and Marius Wiesemann
Physik-Institut, Universität Zürich, CH-8057 Zürich, Switzerland
PRISMA Cluster of Excellence, Institute of Physics,
[0.1cm] Johannes Gutenberg University, D-55099 Mainz, Germany
We consider the transverse-momentum () distribution of and boson pairs produced in hadron collisions. At small , the logarithmically enhanced contributions due to multiple soft-gluon emission are resummed to all orders in QCD perturbation theory. At intermediate and large values of , we consistently combine resummation with the known fixed-order results. We exploit the most advanced perturbative information that is available at present: next-to-next-to-leading logarithmic resummation combined with the next-to-next-to-leading fixed-order calculation. After integration over , we recover the known next-to-next-to-leading order result for the inclusive cross section. We present numerical results at the LHC, together with an estimate of the corresponding uncertainties. We also study the rapidity dependence of the spectrum and we consider efficiencies at different orders of resummed and fixed-order perturbation theory.
Run 1 of the Large Hadron Collider (LHC) has been a great success for the Standard Model (SM). The collected data are in good agreement with the theoretical predictions so far and led to the discovery [1, 2] of a resonance at a mass of GeV, which appears to be fully consistent with the SM Higgs boson. Among the most important reactions at hadron colliders is the production of vector-boson pairs. This class of processes gives access to the vector-boson trilinear couplings which may be modified in a large set of Beyond the Standard Model (BSM) theories. Even small deviations in both the production rate and the shape of distributions could be a signal of new physics. Anomalous couplings related to vector-boson pair production have been constrained first by LEP2, and later by the Tevatron for larger invariant masses. ATLAS and CMS will continue to tighten the bounds on anomalous couplings, especially with increasing sensitivity during Run 2 of the LHC.†††See Ref. and references therein.
On the other hand, vector-boson pair production constitutes an irreducible background to new-physics searches as well as Higgs studies. Particularly important are the off-shell effects below the and thresholds, relevant for the Higgs signal region, and the high-mass tail used to extract the width of the Higgs boson [4, 5, 6]. Furthermore, Higgs boson measurements, in particular in the channel, strongly rely on the background rejection through specific categories based on the transverse momenta of final-state particles, such as the classification into jet bins or in the Higgs transverse momentum and related variables. An accurate modelling of the respective observables for both signal and backgrounds is crucial for such analyses.
The first precise predictions for production at hadron colliders in the SM were obtained at the next-to-leading order (NLO) already more than 20 years ago for stable bosons [7, 8], and the leptonic decays were added in Ref.. The full spin correlations and off-shell effects at the NLO were first included in Refs.[10, 11] using the corresponding one-loop helicity amplitudes . An important loop-induced contribution proceeds through gluon fusion; it is enhanced by the gluon densities and was first computed for on-shell bosons in Refs.[13, 14], while leptonic decays were included later [15, 16, 17]. All the contributions to production discussed so far are implemented in the numerical program MCFM . Electroweak (EW) corrections were evaluated in Ref.[19, 20], while on-shell jet production is known through NLO QCD [21, 22]. Recently, the inclusive cross section for the production of at the next-to-next-to-leading order (NNLO) has been presented .
The production of pairs constitutes the largest cross section among all massive vector-boson pair production modes. On the other hand, its leptonic decay () embodies the most challenging experimental signature since the presence of two neutrinos prohibits a reconstruction of mass peaks. Therefore, a precise understanding of both the signal and the background is required.
The cross sections for on-shell bosons at the NLO [24, 25] as well as the gluon-fusion component [13, 26] have been known for decades. Also in this case, spin and off-shell effects were included in the NLO prediction [10, 11] after the relevant one-loop helicity amplitudes had been computed . Leptonic decays for the loop-induced gluon-fusion contribution were considered in Refs.[27, 28]. More recently, also interference effects with the Higgs production mode through gluon fusion were determined . Analogously to production, all the contributions to production discussed so far are implemented in MCFM . Furthermore, EW corrections have been evaluated [30, 20, 31]. On-shell production in association with one jet has been studied through NLO QCD in Refs.[32, 33, 34]. Detailed Monte Carlo simulations of production in association with up to one jet at NLO have been presented in Ref. . Recently, the first NNLO results for the inclusive cross section have been obtained .
Due to the very recent computation of the helicity amplitudes at two-loop order [37, 38], the inclusion of the off-shell effects and the leptonic decays in the NNLO cross section is expected in the near future. In the meanwhile, also the calculation of helicity amplitudes has been performed at two-loop order [39, 40]. This renders the evaluation of NLO QCD corrections to the gluon-fusion channel feasible.
The production of pairs in hadron collisions has been measured extensively at the Tevatron and the LHC (see Refs.[41, 42, 43, 44, 45, 46, 47] for some recent results). The cross section has also been measured already at the Tevatron, see e.g. Ref., and at the LHC both at 7 TeV [49, 50] and 8 TeV [46, 51, 52]. The ATLAS collaboration recently reported an excess  with respect to the SM prediction, which has drawn a lot of attention on the process, since the final state is a typical signature in many BSM scenarios . In the meanwhile, the excess has been alleviated to a significant degree by the recent NNLO computation . The more recent measurement by the CMS collaboration  is in good agreement with the NNLO prediction.
The transverse-momentum distributions of and pairs are among the most important differential observables for these processes. The spectrum has already been measured in the case of production  at the LHC. Transverse-momentum resummation for and production has been studied in Refs.[54, 55, 56, 57, 58]. In all these calculations the resummed computation is essentially performed up to next-to-leading logarithmic (NLL) accuracy (next-to-next-to-leading logarithmic (NNLL) effects in the Sudakov exponent are considered in Refs.[55, 57, 58]) and matched to the fixed-order result up to the first order in the QCD coupling .
In this paper we consider transverse-momentum resummation for the production of both and pairs. We use the formalism of Ref. to perform the first computation of the spectrum at full NNLL accuracy matched to the fixed-order result valid at large . Although our focus is on the inclusive spectrum of the and system, our computation is fully differential in the degrees of freedom of the vector bosons and allows us to include their EW decays, once the helicity amplitudes are implemented‡‡‡The analogous computations for Higgs, single vector-boson production, and diphoton production are presented in Refs.,  and , respectively.. The -resummation formalism of Ref. is closely related to the subtraction method of Ref., which was used to compute the NNLO cross section for these processes [23, 36]. For production we employ the four-flavour scheme (4FS) and remove all contributions with final-state bottom quarks from our computation of the transverse-momentum distribution in order to eliminate the contamination from and production. The difference to the prediction in the five-flavour scheme (5FS), where such terms have to be consistently subtracted, has been shown to be small for the NNLO inclusive cross section . Furthermore, we neglect the loop-induced gluon-fusion contribution throughout this paper, since, up to NNLO, it contributes only at .
We note that in the CMS measurement reported in Ref. an approximate NNLL prediction  of the spectrum has been used to correct the spectrum from the Monte Carlo simulation. Along these lines, the computation reported in this paper will be useful in order to validate the predictions obtained from Monte Carlo simulations for both and production, as done in the case of Higgs boson production with the calculation of Ref..
The manuscript is organized as follows. In Section 2 we review the transverse-momentum resummation formalism applied to vector-boson pair production. In Section 3 we report our numerical results, starting with remarks on the choice of the resummation scale in Section 3.1. In Section 3.2 we present our numerical predictions for the inclusive spectrum and study the ensuing uncertainties. In Section 3.3 we analyse the behaviour of the spectrum at different rapidities of the vector-boson pair. In Section 3.4 we investigate efficiencies at different orders in resummed and fixed-order perturbation theory. In Section 4 we summarize our results.
2 Transverse-momentum resummation for vector-boson
We consider the inclusive hard-scattering process
where the collision of the two hadrons and with momenta and produces the two vector bosons of momenta and . In the center-of-mass frame the momentum of the vector-boson pair is fully specified by the invariant mass , the rapidity , and the transverse-momentum vector .
The kinematics of the vector bosons is fully determined by the vector-boson pair momentum (with and ), and by the additional and independent variables that specify the angular distribution of the vector bosons with respect to . Throughout this paper we always consider quantities which are inclusive over these angular variables.
According to the QCD factorization theorem, the differential cross section can be written as
where () are the density functions of parton in hadron at the factorization scale ; is the renormalization scale§§§Throughout the paper we use parton densities in the factorization scheme and is the QCD running coupling in the renormalization scheme.; is the partonic cross section. The rapidity and the center-of-mass energy of the partonic scattering process are related to the corresponding hadronic variables and by
When the transverse momentum of the vector-boson pair is of the same order as the invariant mass , the QCD perturbative expansion is controlled by a single expansion parameter, , and fixed-order calculations can be safely applied. In this region, QCD radiative corrections are known to [21, 22, 32, 33, 34]. When the convergence of the perturbative expansion is spoiled by the presence of large logarithmic terms of the form , that need to be resummed to all orders.
The resummation is performed at the level of the partonic cross section, which is decomposed as
The first term on the right-hand side of Eq. (4) contains all the logarithmically-enhanced contributions at small and has to be evaluated by resumming them to all orders. The second term is instead free of such contributions and can thus be evaluated at fixed order in perturbation theory.
Resummation is based on the factorization of soft and collinear radiation and is viable in the impact-parameter () space, where the kinematical constraint of momentum conservation and the factorization of the phase space can be consistently taken into account [66, 67, 68]. Using the Bessel transformation between the conjugate variables and , the resummed component is expressed as [68, 59]
where is the -order Bessel function. In the case of fully inclusive resummation, the rapidity dependence is integrated out in Eq. (5). In that case it is convenient to consider Mellin moments with respect to the variable . However, in order to retain the rapidity dependence in the resummed cross section we take the ‘double’ Mellin moments with respect to the variables and at fixed ,
and organize the structure of in the following exponential form ,
where we have defined the logarithmic expansion parameter as
and ( is the Euler number). The scale appearing in Eqs. (2, 8), named resummation scale in Ref. , parametrizes the arbitrariness in the resummation procedure, and has to be chosen of the order of the hard scale . Variations of around its reference value can be exploited to estimate the size of yet uncalculated higher-order logarithmic contributions. Therefore, plays a role very similar to and for missing perturbative terms. The function does not depend on the impact parameter , and therefore includes all the perturbative terms that behave as constants as . It can thus be expanded in powers of :
where is the partonic leading-order (LO) cross section. The exponent includes the complete dependence on and, in particular, it contains all the terms that order-by-order in are logarithmically divergent as . The logarithmic expansion of reads
where the term collects the LL contributions, the function includes the NLL contributions, controls the NNLL terms and so forth.
The resummation of the large logarithmic terms carried out in Eq. (2), after transforming back to space, allows us to obtain a well behaved transverse-momentum spectrum as . However, the logarithmic expansion parameter in Eq. (8) is divergent as . This implies that the resummation produces higher-order contributions also in the high- region, which is conjugated to after Fourier transformation. In this region the fixed-order cross section is perfectly viable and any resummation effect is necessarily artificial. To reduce the impact of such contributions, the logarithmic variable is replaced by 
The variables and are equivalent when but they have a very different behaviour as . When , and . Moreover, since the behaviour of the Sudakov form factor at is related to the integral over , the replacement in Eq. (11) allows us to enforce a unitarity constraint such that the fixed-order prediction is recovered upon integration over .
A well known property of the formalism of Ref.  is that the process dependence (as well as the factorization scale and scheme dependence) is fully encoded in the hard function . In other words, the functions are universal: they depend only on the channel in which the process occurs at Born level ( annihilation in the case of vector-boson pair production). Their explicit expressions up to are given in Ref.  in terms of the universal perturbative coefficients , , , , . In particular, the LL function depends on the coefficient , the NLL function also depends on and  and the NNLL function also depends on  and [71, 72, 73].
The hard coefficients depend on the process we want to consider. The first order coefficients are known since long time [72, 73]: they can be obtained from the one-loop scattering amplitudes by using a process independent relation. By exploiting the expressions of for single Higgs  and vector-boson  production, in Ref. such a relation has been extended to : this implies that the two-loop amplitude for is the only process-dependent information needed to obtain the coefficient .
We now turn to the finite component of the transverse-momentum spectrum, i.e. the second term on the right hand side of Eq. (4). Since does not contain large logarithmic terms in the small- region, it can be evaluated by truncating the perturbative series at a given fixed order. In practice, the finite component is computed starting from the customary fixed-order (f.o.) perturbative truncation of the partonic cross section and subtracting the expansion of the resummed cross section in Eq. (5) at the same perturbative order:
At least formally, this matching procedure between resummed and finite contributions guarantees to achieve a uniform theoretical accuracy over the entire range of transverse momenta. At large values of , the resummation procedure cannot improve the fixed-order result, and the resummation (and matching) procedure is eventually superseded by the customary fixed-order calculations.
In summary, the inclusion of the functions , , in the resummed component, together with the evaluation of the finite component at NLO (i.e. ), allows us to perform the resummation at NLL+NLO accuracy. This is the theoretical accuracy of the calculations of Refs.[55, 56]. Including also the functions and , together with the finite component at NNLO (i.e. ) leads to full NNLL+NNLO accuracy. Using the recently computed two-loop amplitudes for , and the process independent relation of Ref., we are now able to present the complete result for the transverse-momentum distribution of the vector-boson pair up to NNLL+NNLO accuracy. We point out that the NNLL+NNLO (NLL+NLO) result includes the full NNLO (NLO) perturbative contribution in the entire range. In particular, the NNLO (NLO) result for the double differential cross section is exactly recovered upon integration over of the differential cross section at NNLL+NNLO (NLL+NLO) accuracy.
We conclude this Section by adding a few comments on the way in which our calculation is actually performed. The practical implementation of Eq. (4) is done in the numerical program Matrix¶¶¶Matrix is the abbreviation of “Munich Automates qT subtraction and Resummation to Integrate X-sections”, by M. Grazzini, S. Kallweit, D. Rathlev, M. Wiesemann. In preparation., which is an extension of the numerical program applied in the NNLO calculations of Refs.[76, 23, 36, 77] and based on a combination of the -subtraction formalism  with the Munich∥∥∥Munich is the abbreviation of “MUlti-chaNnel Integrator at Swiss (CH) precision”—an automated parton level NLO generator by S. Kallweit. In preparation. code. Since already performed within the -subtraction formalism, the extension of these calculations to compute the resummed cross section is conceptually quite straightforward, and is obtained by replacing the hard-collinear terms in the fixed-order computation by the proper all-order resummation formula of Eq. (5). This procedure is the same that was applied to perform the NLL+NLO calculations for  and  production, and the NNLL+NNLO calculation for Higgs boson production of Ref..
To obtain the numerical results presented here, the resummed component of Eq. (5) is evaluated with an extension of the numerical program used for the calculation of Higgs production , based on the earlier computations of Refs.[59, 64]. The hard-collinear coefficients are obtained by exploiting the implementation of the corresponding virtual amplitudes for the production of on-shell and pairs in Ref. and Ref., respectively, and the knowledge of the collinear coefficients relevant to quark-initiated processes .******We note that an independent computation of these coefficients in the framework of Soft Collinear Effective Theory has been presented in Refs.[78, 79].
The finite component of Eq. (12) is obtained from an NLO calculation of +jet, computed with the Munich code, which provides a fully automated implementation of the Catani–Seymour dipole formalism [80, 81] as well as an interface to the one-loop generator OpenLoops  to obtain all required (spin and color-correlated) tree-level and one-loop amplitudes. For the numerically stable evaluation of tensor integrals we rely on the Collier library , which is based on the Denner–Dittmaier reduction techniques [84, 85] and the scalar integrals of . To deal with problematic phase-space points, OpenLoops provides a rescue system using the quadruple-precision implementation of the OPP method in CutTools , involving scalar integrals from OneLOop .
In this Section we present our results for the resummed transverse-momentum distributions of and pairs. We compare our NNLL+NNLO predictions to the results at the NLL+NLO, and discuss the corresponding theoretical uncertainties. Additionally, we also study the rapidity dependence of the cross section as well as the -veto efficiency.
For the EW couplings we use the so-called scheme, where the input parameters are , , . In particular we set GeV, GeV, GeV. We use the NNPDF3.0 sets of parton distribution functions (PDFs)  with . At NLL+NLO and NNLL+NNLO the running of is evaluated at two- and three-loop order, respectively. For production we consider massless quarks/antiquarks. For production we make use of the 4FS, which allows us to split off all contributions related to bottom-quark final states in order to remove the and contamination from our computation. This is straightforward in the 4FS, because the process is separately finite.
We consider proton-proton collisions at TeV. The central values of the factorization and renormalization scales are set to . The choice of the central resummation scale is discussed in the next subsection.
3.1 Choice of the central resummation scale
As discussed in Section 2, the resummation scale is the scale entering the large logarithmic terms we are resumming (see Eq. (8)), and it plays the role of the scale up to which resummation is effective. In on-shell Higgs  and vector-boson  production, the scale is typically chosen equal to half the mass of the heavy boson (i.e. for Higgs production and in the case of single vector-boson production). Higher values of the scale lead to a worse matching at high . The natural extension of this choice for vector-boson pair production is a dynamical resummation scale , since is the hardness of the process. This is indeed the choice that was adopted in the calculations of Refs. [55, 56].
The following considerations apply both to and production, so we will focus on production from now on. In Fig. 1 we consider the invariant mass distribution of the pair at NLO and NNLO. We see that the distribution is strongly peaked in the threshold region, and that it quickly decreases as increases. As a consequence, for most of the events, .
We can compare the transverse-momentum distributions obtained with a dynamical resummation scale , and a fixed resummation scale . In Fig. 2 we show the ratio (blue, solid curve) of the result over the result. The bands are obtained by varying the resummation scale around the central value by a factor of two. Considering the ratio of the central curves for GeV, the differences between a fixed and a dynamical scale are extremely small and remain at the 1-2% level over the whole range. In this region of transverse momenta the uncertainty bands obtained with the two choices overlap and are similar in size. In fact, since leads to slightly larger uncertainties, it appears to be the more conservative choice. Therefore, we can conclude that either choice of the resummation scale is perfectly valid and indeed consistent with each other as expected from the discussion of the invariant mass distribution.
Looking further at the comparison of the high- tails in Fig. 2 ( GeV), we observe a very well known feature [59, 90, 60, 91, 92] of the applied matching procedure, namely the fact that for large values of the resummation scale the fixed-order cross section (black dotted curve) is not recovered in the tail of the distribution. It is important to recall that transverse-momentum resummation is supposed to improve the perturbative expansion in the low- region. At large , any large dependence on the resummation scale is necessarily artificial and an unwanted remnant of the matching procedure. This behaviour is precisely what we observe for the dynamical scale choice in Fig. 2. With this choice, in fact, the resummed result loses predictivity, as its uncertainty becomes increasingly large. By contrast, a fixed resummation scale , which is always smaller than , eventually leads to a more consistent high- behaviour of the resummed prediction.
Based on the above results, we make our default choice of the resummation scale in what follows.
3.2 Inclusive transverse-momentum distribution
We now present our resummed predictions for the inclusive transverse-momentum spectrum of the vector-boson pair and compare them with the corresponding fixed-order results. We concentrate on production since we observe no saliently different features in the case. For completeness, we provide the corresponding reference prediction with uncertainties for below.
Before presenting our resummed predictions, we recall the well known fixed-order results at and [32, 33, 34]. In Fig. 3 we show the NLO and NNLO distributions together with their perturbative uncertainties. The uncertainty bands are obtained by varying and in the range with the constraint . The lower inset shows the same results normalized to the central NLO curve. The NNLO effects range from about at GeV to about at GeV. The NLO (NNLO) uncertainty ranges from about () at GeV to about () at GeV. We note that the NLO and NNLO bands do not overlap in the region where GeV. This implies that, in this region of transverse momenta, the size of the band obtained through scale variations at NLO definitely underestimates the theoretical uncertainty.
We now move on to the resummed results. In Fig. 5 (a) the NLL+NLO spectrum is compared to the fixed-order NLO result and to the finite component of the resummed cross section (see Eq. (4)) in the region between and GeV. As expected, the NLO diverges to as , while the resummation provides a physically well behaved spectrum down to low values of , which exhibits a kinematical peak at GeV. The finite component contributes less than in the peak region, where the result is dominated by resummation, and it increases to at GeV. The lower inset shows the NLL+NLO result normalized to NLO. In Fig. 5 (b) the region between and GeV is displayed. We see that even at large values of the NLL+NLO resummed result does not match very well the fixed-order NLO result, with a difference of about %.
The analogous results at NNLL+NNLO are shown in Fig. 5. The NNLO has an unphysical (divergent) behaviour as , whereas the resummed spectrum is well behaved, with a slightly harder peak with respect to the NLL+NLO. The finite component contributes less than in the peak region, increasing to at GeV. Comparing the right panels of Fig. 5 and Fig. 5, we see that the quality of the matching at high is significantly improved when going from NLL+NLO to NNLL+NNLO, and we find that this behaviour is indeed preserved up to very high transverse momenta. The NNLL+NNLO result thus gives a prediction with uniform accuracy from small to very large transverse momenta and, in fact, provides a sufficiently large region where a hard switching to the fixed-order result is feasible. We point out that, thanks to our unitarity constraint, both at NLL+NLO and at NNLL+NNLO the integral of the resummed spectrum is in excellent agreement with the respective total cross sections; the differences are at the few-permille level.
We now turn to the scale uncertainties of our resummed results. We start our discussion by separately considering factorization and renormalization scale variations. In Fig. 6 we compare the NLL+NLO (red, dashed) and NNLL+NNLO (blue, solid) predictions with their uncertainty bands from and variations (left and right panel, respectively). In both cases, the bands are obtained by varying the factorization (renormalization) scale by a factor of two around its central value, while keeping the other scales at their default values. First of all, we notice that when going from NLL+NLO to NNLL+NNLO the spectrum becomes harder. Comparing with the results of Ref. , where the NNLL resummation was implemented without matching, we see that the increased hardness of the spectrum is a combined effect of both features, i.e. NNLL resummation and NNLO matching at high .
We note that neither in the case of the factorization scale, nor in the case of the renormalization scale, the NLL+NLO and NNLL+NNLO bands overlap. Actually, in the case of the factorization scale, there is no reduction in scale dependence when going from NLL+NLO to NNLL+NNLO, and the uncertainty slightly increases with the perturbative order, even if it is always well below , except at very low . The renormalization scale dependence instead exhibits the expected reduction when going from NLL+NLO to NNLL+NNLO.
In Fig. 7 we present our resummed predictions with uncertainty bands obtained from simultaneous variations of and (left panel) and the variation of the resummation scale (right panel). In the left panel the uncertainty bands are obtained by varying and as in Fig. 3. In the right panel the resummation scale is varied in the range . As in Fig. 6 we see that the uncertainty bands do not overlap. The uncertainty from and variations is at NLL+NLO and is reduced to at NNLL+NNLO. At NLL+NLO the resummation scale uncertainty is generally about except in the region of GeV, where it shrinks to smaller values. We find that at NNLL+NNLO the resummation scale uncertainty is reduced roughly by a factor of two in the region of transverse momenta considered in the figure.
In Figs. 9 and 9 we show our reference resummed prediction for and , respectively, with an estimate of their full perturbative uncertainty. In order to obtain a combined uncertainty from , and variations, we follow Ref.  and independently vary , and in the ranges and with the constraints and . We recall that the constraint on , which is the same as applied in Fig. 3 and Fig. 7 (left), has the purpose of avoiding large logarithmic contributions from the evolution of parton densities. Analogously, the constraint on avoids large logarithmic contributions in the expansion of the Sudakov form factor.
For production the perturbative uncertainty at NNLL+NNLO (NLL+NLO) is about () at the peak, it decreases to about () at GeV, and it increases again to () at GeV. In the high- region, the difference between the NNLL+NNLO and NLL+NLO predictions is driven by the NNLO effects, which increase the NLO result by about .
For production the uncertainties have essentially the same pattern in the small- and intermediate- region, while at high they are larger than for production, reaching about at NNLL+NNLO for GeV. We have checked that this effect is entirely driven by the resummation-scale dependence. As previously pointed out, this behaviour is not particularly worrying since, in the large- region, the resummed results should be replaced by the corresponding fixed-order prediction. Also in the case the large enhancement of the NNLL+NNLO distribution in the high- tail stems from the fixed-order cross section.
3.3 Rapidity dependence of the transverse-momentum distribution
So far, we only considered spectra for on-shell and production that are inclusive in the kinematics of the vector-boson pair. Our numerical program, however, allows us to compute arbitrary observables that are differential with respect to the phase space.
In the following we study the behaviour of the transverse-momentum spectrum in different rapidity regions of the vector-boson pair. In Fig. 10 we study the shape of the NNLL+NNLO transverse-momentum distribution, i.e. normalized such that its integral yields one, for (red, solid), (blue, dashed), (black, dotted), (magenta, dash-dotted) and (orange, dash-double dotted). The right panel shows the same results normalized to the fully inclusive distribution. We clearly see that the shapes become softer as the rapidity increases. In the central region () the distributions are still quite insensitive to the specific value of the rapidity and only slightly harder than the inclusive spectrum. In the forward rapidity region, on the other hand, the shapes become increasingly softer.
The observed pattern can be understood in the following way: rapidity and transverse momentum are two not completely independent phase-space variables. Indeed, they affect their mutual upper integration bounds. At higher rapidities the kinematically allowed range of transverse momenta is reduced: this squeezes the spectrum which consequently becomes softer. This effect has been observed also in previous studies in the case of Higgs boson production .
3.4 The cross section and -veto efficiencies
The excess in the production cross section measured by ATLAS  with respect to the SM prediction has drawn a lot of attention to the process, since the signature appears in many new physics scenarios . The inclusion of the recently computed NNLO corrections  considerably reduces the significance of the excess. However, particular attention must be payed to the modelling of the jet veto [93, 58, 52] when extrapolating from the fiducial region to obtain the inclusive cross section. Effects of jet veto resummation have been considered in Refs.[94, 95], though still matching to the fixed-order result.
In this paper we are dealing with transverse-momentum spectra, and we perform a resummation on a different variable with respect to the jet . However, the vector-boson pair and the jet are clearly related variables (actually, at they indeed coincide). We will therefore study the -veto efficiency in production at different orders in resummed and fixed-order perturbation theory. We define the -veto efficiency as
In Fig. 11 we show at the NNLL+NNLO (blue, solid), approximate NNLL+NLO (magenta, dash-double dotted), NLL+NLO (red, dashed), NNLO (black, dotted) and NLO (grey, dash-dotted). The lower inset shows the same curves normalized to our reference prediction at NNLL+NNLO. Our approximate NNLL+NLO is obtained by simply adding the function in the Sudakov exponent in Eq. (2) at NLL+NLO, and corresponds to the approximation considered in Refs.[55, 58].
For reference, the corresponding numerical values of the efficiencies are given in Tab. 1 for - GeV. The uncertainty bands are obtained by a combined variation of resummation, factorization and renormalization scales as in Fig. 9. The first thing we observe is that the NLO result appears to be well above the others and cannot be really considered a reliable prediction for the efficiency. This is because it is essentially a LO prediction at finite values of . We also note that in the small- region (say below GeV) the fixed-order NLO and NNLO predictions diverge and cannot be trusted. Comparing further the fixed-order results among each other and the resummed results among each other, we observe that higher-order corrections in fixed-order and resummed perturbation theory reduce the -veto efficiency.
Both effects can be easily understood in the light of the results presented up to now. As seen in Fig. 3, the inclusion of the NNLO corrections make the distribution harder. Furthermore, resummation effects generally harden the spectrum. A qualitatively similar result is obtained when going from NLL+NLO to NNLL+NNLO (see Fig. 9).
It is interesting to compare the approximated NNLL+NLO result with the NNLO and NNLL+NNLO predictions. For values of GeV we see that the approximated result is in between the NNLO one and our best NNLL+NNLO prediction. This means that the effect of NNLL resummation obtained by the inclusion of the function in the Sudakov exponent in Eq. (2) is quantitatively important. Nonetheless, the efficiency obtained within this approximation is still about higher than the NNLL+NNLO prediction. We also notice that in this region of , the NNLO and NLL+NLO results differ by less than .
Comparing the NNLL+NNLO and NLL+NLO results, we find that they are compatible within the corresponding uncertainties.
We add few comments on the recent measurement of the cross section carried out by the CMS collaboration . The result shows good agreement with the NNLO prediction of Ref.. The corresponding analysis, however, is based on a reweighting procedure of the spectrum of the pair. The events generated with POWHEG  plus Pythia6  were reweighted by using the calculation of Ref., which corresponds to our NNLL+NLO approximation, and includes neither the second-order hard-collinear coefficient in Eq. (2), nor the NNLO matching. The results in Fig. 11 show that the NNLL+NNLO -veto efficiency is lower than the efficiency obtained with the approximated NNLL+NLO calculation. As a consequence, a reweighting to the full NNLL+NNLO prediction for the spectrum would most likely lead to a decrease of the jet-veto efficiency.
In this paper we have studied the transverse-momentum distribution of vector-boson pairs in hadronic collisions. We presented a computation of the spectrum in which the logarithmically enhanced contributions at small are resummed up to NNLL accuracy and the ensuing result is combined with state-of-the-art (NNLO) predictions valid at large . We presented numerical results for and production at the LHC together with a study of their perturbative uncertainties.
We found that, up to relatively large transverse momenta, when scale variations are studied around the fixed resummation scale we obtain results which are fully consistent with those obtained using the dynamical choice . At large transverse momenta the fixed-order result in the tail of the distribution is nicely recovered with the fixed-scale choice. Our new NNLL+NNLO results significantly reduce the theoretical uncertainties obtained through scale variations compared to lower orders in both the peak region and the tail of the distribution.
We have also studied the rapidity dependence of the resummed transverse-momentum distribution. The rapidity dependence at NNLL+NNLO is quite flat in the central region (), but signals a substantially softer spectrum in the forward region. Due to phase-space suppression, the effect on the inclusive transverse-momentum distribution is very moderate though.
Finally, we have studied the -veto efficiency at different orders in resummed and fixed-order perturbation theory. Both NNLL resummation and the NNLO effects turned out to be important to obtain an accurate prediction for this quantity. We observed that the veto efficiency at NNLL+NNLO is % lower ( in absolute terms) with respect to the approximate NNLL+NLO calculation used in the CMS analysis of Ref.. This result suggests that our NNLL+NNLO predictions will be useful to validate the transverse-momentum spectra obtained from Monte Carlo event generators, similarly to what was done for the NNLL+NNLO calculation of Ref. in the case of Higgs boson production.
In this paper we considered the spectrum of stable vector-boson pairs. Exploiting the two-loop helicity amplitudes for leptons [37, 38] will allow us to extend the calculation to include the leptonic decay of the vector bosons and off-shell effects. The computation of the transverse-momentum spectrum with realistic experimental cuts will then become possible.
Acknowledgements. We would like to thank Andreas von Manteuffel and Lorenzo Tancredi for providing us with their private code to evaluate the helicity-averaged on-shell amplitudes in the equal-mass case. We would like to thank Giancarlo Ferrera for comments on the manuscript. This research was supported in part by the Swiss National Science Foundation (SNF) under contracts CRSII2-141847, 200021-156585 and by the Research Executive Agency (REA) of the European Union under the Grant Agreement number PITN-GA-2012-316704 (Higgstools).
-  ATLAS Collaboration, Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC, Phys.Lett. B716 (2012) 1–29, [arXiv:1207.7214].
-  CMS Collaboration, Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC, Phys.Lett. B716 (2012) 30–61, [arXiv:1207.7235].
-  ATLAS, D0, CDF, CMS Collaboration, Diboson Production at LHC and Tevatron, Int.J.Mod.Phys.Conf.Ser. 31 (2014) 1460279, [arXiv:1403.1415].
-  N. Kauer and G. Passarino, Inadequacy of zero-width approximation for a light Higgs boson signal, JHEP 1208 (2012) 116, [arXiv:1206.4803].
-  F. Caola and K. Melnikov, Constraining the Higgs boson width with production at the LHC, Phys.Rev. D88 (2013) 054024, [arXiv:1307.4935].
-  J. M. Campbell, R. K. Ellis, and C. Williams, Bounding the Higgs width at the LHC: Complementary results from , Phys.Rev. D89 (2014), no. 5 053011, [arXiv:1312.1628].
-  J. Ohnemus and J. Owens, An Order calculation of hadronic production, Phys.Rev. D43 (1991) 3626–3639.
-  B. Mele, P. Nason, and G. Ridolfi, QCD radiative corrections to Z boson pair production in hadronic collisions, Nucl.Phys. B357 (1991) 409–438.
-  J. Ohnemus, Hadronic , , and production with QCD corrections and leptonic decays, Phys.Rev. D50 (1994) 1931–1945, [hep-ph/9403331].
-  J. M. Campbell and R. K. Ellis, An Update on vector boson pair production at hadron colliders, Phys.Rev. D60 (1999) 113006, [hep-ph/9905386].
-  L. J. Dixon, Z. Kunszt, and A. Signer, Vector boson pair production in hadronic collisions at order : Lepton correlations and anomalous couplings, Phys.Rev. D60 (1999) 114037, [hep-ph/9907305].
-  L. J. Dixon, Z. Kunszt, and A. Signer, Helicity amplitudes for O() production of , , , , or pairs at hadron colliders, Nucl.Phys. B531 (1998) 3–23, [hep-ph/9803250].
-  D. A. Dicus, C. Kao, and W. Repko, Gluon Production of Gauge Bosons, Phys.Rev. D36 (1987) 1570.
-  J. van der Bij and E. N. Glover, Photon boson pair production via gluon fusion, Phys.Lett. B206 (1988) 701.
-  T. Matsuura and J. van der Bij, Characteristics of leptonic signals for Z boson pairs at hadron colliders, Z.Phys. C51 (1991) 259–266.
-  C. Zecher, T. Matsuura, and J. van der Bij, Leptonic signals from off-shell Z boson pairs at hadron colliders, Z.Phys. C64 (1994) 219–226, [hep-ph/9404295].
-  T. Binoth, N. Kauer, and P. Mertsch, Gluon-induced QCD corrections to , arXiv:0807.0024.
-  J. M. Campbell, R. K. Ellis, and C. Williams, Vector boson pair production at the LHC, JHEP 1107 (2011) 018, [arXiv:1105.0020].
-  A. Bierweiler, T. Kasprzik, and J. H. Kühn, Vector-boson pair production at the LHC to accuracy, JHEP 1312 (2013) 071, [arXiv:1305.5402].
-  J. Baglio, L. D. Ninh, and M. M. Weber, Massive gauge boson pair production at the LHC: a next-to-leading order story, Phys.Rev. D88 (2013) 113005, [arXiv:1307.4331].
-  T. Binoth, T. Gleisberg, S. Karg, N. Kauer, and G. Sanguinetti, NLO QCD corrections to +jet production at hadron colliders, Phys.Lett. B683 (2010) 154–159, [arXiv:0911.3181].
-  SM and NLO Multileg Working Group Collaboration, T. Binoth et al., The SM and NLO Multileg Working Group: Summary report, arXiv:1003.1241.
-  F. Cascioli, T. Gehrmann, M. Grazzini, S. Kallweit, P. Maierhöfer, A. von Manteuffel, S. Pozzorini, D. Rathlev, L. Tancredi, and E. Weihs, production at hadron colliders in NNLO QCD, Phys.Lett. B735 (2014) 311–313, [arXiv:1405.2219].
-  J. Ohnemus, An Order calculation of hadronic production, Phys.Rev. D44 (1991) 1403–1414.
-  S. Frixione, A Next-to-leading order calculation of the cross-section for the production of pairs in hadronic collisions, Nucl.Phys. B410 (1993) 280–324.
-  E. N. Glover and J. van der Bij, Vector boson pair production via gluon fusion, Phys.Lett. B219 (1989) 488.
-  T. Binoth, M. Ciccolini, N. Kauer, and M. Krämer, Gluon-induced background to Higgs boson searches at the LHC, JHEP 0503 (2005) 065, [hep-ph/0503094].
-  T. Binoth, M. Ciccolini, N. Kauer, and M. Krämer, Gluon-induced W-boson pair production at the LHC, JHEP 0612 (2006) 046, [hep-ph/0611170].
-  J. M. Campbell, R. K. Ellis, and C. Williams, Gluon-Gluon Contributions to Production and Higgs Interference Effects, JHEP 1110 (2011) 005, [arXiv:1107.5569].
-  A. Bierweiler, T. Kasprzik, J. H. Kühn, and S. Uccirati, Electroweak corrections to W-boson pair production at the LHC, JHEP 1211 (2012) 093, [arXiv:1208.3147].
-  M. Billoni, S. Dittmaier, B. Jäger, and C. Speckner, Next-to-leading order electroweak corrections to leptons at the LHC in double-pole approximation, JHEP 1312 (2013) 043, [arXiv:1310.1564].
-  S. Dittmaier, S. Kallweit, and P. Uwer, NLO QCD corrections to +jet production at hadron colliders, Phys.Rev.Lett. 100 (2008) 062003, [arXiv:0710.1577].
-  J. M. Campbell, R. K. Ellis, and G. Zanderighi, Next-to-leading order predictions for jet distributions at the LHC, JHEP 0712 (2007) 056, [arXiv:0710.1832].
-  S. Dittmaier, S. Kallweit, and P. Uwer, NLO QCD corrections to +jet+X including leptonic -boson decays, Nucl.Phys. B826 (2010) 18–70, [arXiv:0908.4124].
-  F. Cascioli, S. Höche, F. Krauss, P. Maierhöfer, S. Pozzorini, and F. Siegert, Precise Higgs-background predictions: merging NLO QCD and squared quark-loop corrections to four-lepton + 0,1 jet production, JHEP 1401 (2014) 046, [arXiv:1309.0500].
-  T. Gehrmann, M. Grazzini, S. Kallweit, , P. Maierhöfer, A. von Manteuffel, S. Pozzorini, D. Rathlev, and L. Tancredi, Production at Hadron Colliders in Next to Next to Leading Order QCD, Phys.Rev.Lett. 113 (2014) 212001, [arXiv:1408.5243].
-  F. Caola, J. M. Henn, K. Melnikov, A. V. Smirnov, and V. A. Smirnov, Two-loop helicity amplitudes for the production of two off-shell electroweak bosons in quark-antiquark collisions, JHEP 1411 (2014) 041, [arXiv:1408.6409].
-  T. Gehrmann, A. von Manteuffel, and L. Tancredi, The two-loop helicity amplitudes for , arXiv:1503.04812.
-  F. Caola, J. M. Henn, K. Melnikov, A. V. Smirnov, and V. A. Smirnov, Two-loop helicity amplitudes for the production of two off-shell electroweak bosons in gluon fusion, JHEP 1506 (2015) 129, [arXiv:1503.08759].
-  A. von Manteuffel and L. Tancredi, The two-loop helicity amplitudes for , JHEP 1506 (2015) 197, [arXiv:1503.08835].
-  CDF Collaboration, Measurement of production in leptonic final states at of 1.96 TeV at CDF, Phys.Rev.Lett. 108 (2012) 101801, [arXiv:1112.2978].
-  D0 Collaboration, A measurement of the and production cross sections using leptonic final states in 8.6 fb of collisions, Phys.Rev. D85 (2012) 112005, [arXiv:1201.5652].
-  ATLAS Collaboration, Measurement of production in collisions at TeV and limits on anomalous and couplings with the ATLAS detector, JHEP 1303 (2013) 128, [arXiv:1211.6096].
-  CMS Collaboration, Measurement of the production cross section and search for anomalous couplings in 2l2l’ final states in collisions at TeV, JHEP 1301 (2013) 063, [arXiv:1211.4890].
-  ATLAS Collaboration, Measurement of the total production cross section in proton-proton collisions at in 20 fb with the ATLAS detector, ATLAS-CONF-2013-020, ATLAS-COM-CONF-2013-020.
-  CMS Collaboration, Measurement of and production cross sections in pp collisions at TeV, Phys.Lett. B721 (2013) 190–211, [arXiv:1301.4698].
-  CMS Collaboration, Measurement of the production cross section and constraints on anomalous triple gauge couplings in four-lepton final states at 8 TeV, Phys.Lett. B740 (2015) 250–272, [arXiv:1406.0113].
-  D0 Collaboration, Measurements of and production in + jets final states in collisions, Phys.Rev.Lett. 108 (2012) 181803, [arXiv:1112.0536].
-  ATLAS Collaboration, Measurement of production in collisions at =7 TeV with the ATLAS detector and limits on anomalous and couplings, Phys.Rev. D87 (2013), no. 11 112001, [arXiv:1210.2979].
-  CMS Collaboration, Measurement of the Cross section in Collisions at TeV and Limits on Anomalous and couplings, Eur.Phys.J. C73 (2013), no. 10 2610, [arXiv:1306.1126].
-  ATLAS Collaboration, Measurement of the production cross section in proton-proton collisions at TeV with the ATLAS detector, ATLAS-CONF-2014-033, ATLAS-COM-CONF-2014-045.
-  CMS Collaboration, Measurement of the cross section in pp collisions at sqrt(s) = 8 TeV and limits on anomalous gauge couplings, CMS-PAS-SMP-14-016.
-  D. E. Morrissey, T. Plehn, and T. M. Tait, Physics searches at the LHC, Phys.Rept. 515 (2012) 1–113, [arXiv:0912.3259].
-  C. Balazs and C. Yuan, Higgs boson production at hadron colliders with soft gluon effects. 1. Backgrounds, Phys.Rev. D59 (1999) 114007, [hep-ph/9810319].
-  M. Grazzini, Soft-gluon effects in production at hadron colliders, JHEP 0601 (2006) 095, [hep-ph/0510337].
-  R. Frederix and M. Grazzini, Higher-order QCD effects in the search channel at the LHC, Phys.Lett. B662 (2008) 353–359, [arXiv:0801.2229].
-  Y. Wang, C. S. Li, Z. L. Liu, D. Y. Shao, and H. T. Li, Transverse-Momentum Resummation for Gauge Boson Pair Production at the Hadron Collider, Phys.Rev. D88 (2013) 114017, [arXiv:1307.7520].
-  P. Meade, H. Ramani, and M. Zeng, Transverse momentum resummation effects in measurements, Phys.Rev. D90 (2014), no. 11 114006, [arXiv:1407.4481].
-  G. Bozzi, S. Catani, D. de Florian, and M. Grazzini, Transverse-momentum resummation and the spectrum of the Higgs boson at the LHC, Nucl.Phys. B737 (2006) 73–120, [hep-ph/0508068].
-  D. de Florian, G. Ferrera, M. Grazzini, and D. Tommasini, Higgs boson production at the LHC: transverse momentum resummation effects in the , and decay modes, JHEP 1206 (2012) 132, [arXiv:1203.6321].
-  S. Catani, D. de Florian, G. Ferrera, and M. Grazzini, Vector boson production at hadron colliders: transverse-momentum resummation and leptonic decay, arXiv:1507.06937.
-  L. Cieri, F. Coradeschi, and D. de Florian, Diphoton production at hadron colliders: transverse-momentum resummation at next-to-next-to-leading logarithmic accuracy, arXiv:1505.03162.
-  S. Catani and M. Grazzini, An NNLO subtraction formalism in hadron collisions and its application to Higgs boson production at the LHC, Phys.Rev.Lett. 98 (2007) 222002, [hep-ph/0703012].
-  G. Bozzi, S. Catani, D. de Florian, and M. Grazzini, Higgs boson production at the LHC: Transverse-momentum resummation and rapidity dependence, Nucl.Phys. B791 (2008) 1–19, [arXiv:0705.3887].
-  S. Catani, L. Cieri, D. de Florian, G. Ferrera, and M. Grazzini, Universality of transverse-momentum resummation and hard factors at the NNLO, Nucl.Phys. B881 (2014) 414–443, [arXiv:1311.1654].
-  G. Parisi and R. Petronzio, Small Transverse Momentum Distributions in Hard Processes, Nucl.Phys. B154 (1979) 427.
-  G. Curci, M. Greco, and Y. Srivastava, QCD Jets From Coherent States, Nucl.Phys. B159 (1979) 451.
-  J. C. Collins, D. E. Soper, and G. F. Sterman, Transverse Momentum Distribution in Drell-Yan Pair and W and Z Boson Production, Nucl.Phys. B250 (1985) 199.
-  J. Kodaira and L. Trentadue, Summing Soft Emission in QCD, Phys.Lett. B112 (1982) 66.
-  T. Becher and M. Neubert, Drell-Yan Production at Small , Transverse Parton Distributions and the Collinear Anomaly, Eur.Phys.J. C71 (2011) 1665, [arXiv:1007.4005].
-  C. Davies and W. J. Stirling, Nonleading Corrections to the Drell-Yan Cross-Section at Small Transverse Momentum, Nucl.Phys. B244 (1984) 337.
-  D. de Florian and M. Grazzini, Next-to-next-to-leading logarithmic corrections at small transverse momentum in hadronic collisions, Phys.Rev.Lett. 85 (2000) 4678–4681, [hep-ph/0008152].
-  D. de Florian and M. Grazzini, The Structure of large logarithmic corrections at small transverse momentum in hadronic collisions, Nucl.Phys. B616 (2001) 247–285, [hep-ph/0108273].
-  S. Catani and M. Grazzini, Higgs Boson Production at Hadron Colliders: Hard-Collinear Coefficients at the NNLO, Eur.Phys.J. C72 (2012) 2013, [arXiv:1106.4652].
-  S. Catani, L. Cieri, D. de Florian, G. Ferrera, and M. Grazzini, Vector boson production at hadron colliders: hard-collinear coefficients at the NNLO, Eur.Phys.J. C72 (2012) 2195, [arXiv:1209.0158].
-  M. Grazzini, S. Kallweit, D. Rathlev, and A. Torre, production at hadron colliders in NNLO QCD, Phys.Lett. B731 (2014) 204–207, [arXiv:1309.7000].
-  M. Grazzini, S. Kallweit, and D. Rathlev, and production at the LHC in NNLO QCD, arXiv:1504.01330.
-  T. Gehrmann, T. Lübbert, and L. L. Yang, Transverse parton distribution functions at next-to-next-to-leading order: the quark-to-quark case, Phys.Rev.Lett. 109 (2012) 242003, [arXiv:1209.0682].
-  T. Gehrmann, T. Lübbert, and L. L. Yang, Calculation of the transverse parton distribution functions at next-to-next-to-leading order, JHEP 1406 (2014) 155, [arXiv:1403.6451].
-  S. Catani and M. Seymour, The Dipole formalism for the calculation of QCD jet cross-sections at next-to-leading order, Phys.Lett. B378 (1996) 287–301, [hep-ph/9602277].
-  S. Catani and M. Seymour, A General algorithm for calculating jet cross-sections in NLO QCD, Nucl.Phys. B485 (1997) 291–419, [hep-ph/9605323].
-  F. Cascioli, P. Maierhöfer, and S. Pozzorini, Scattering Amplitudes with Open Loops, Phys.Rev.Lett. 108 (2012) 111601, [arXiv:1111.5206].
-  A. Denner, S. Dittmaier, and L. Hofer, COLLIER - A fortran-library for one-loop integrals, PoS LL2014 (2014) 071, [arXiv:1407.0087].
-  A. Denner and S. Dittmaier, Reduction of one loop tensor five point integrals, Nucl.Phys. B658 (2003) 175–202, [hep-ph/0212259].
-  A. Denner and S. Dittmaier, Reduction schemes for one-loop tensor integrals, Nucl.Phys. B734 (2006) 62–115, [hep-ph/0509141].
-  A. Denner and S. Dittmaier, Scalar one-loop 4-point integrals, Nucl.Phys. B844 (2011) 199–242, [arXiv:1005.2076].
-  G. Ossola, C. G. Papadopoulos, and R. Pittau, CutTools: A Program implementing the OPP reduction method to compute one-loop amplitudes, JHEP 0803 (2008) 042, [arXiv:0711.3596].
-  A. van Hameren, OneLOop: For the evaluation of one-loop scalar functions, Comput.Phys.Commun. 182 (2011) 2427–2438, [arXiv:1007.4716].
-  NNPDF Collaboration, R. D. Ball et al., Parton distributions for the LHC Run II, JHEP 1504 (2015) 040, [arXiv:1410.8849].
-  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, Phys.Lett. B696 (2011) 207–213, [arXiv:1007.2351].
-  R. V. Harlander, A. Tripathi, and M. Wiesemann, Higgs production in bottom quark annihilation: Transverse momentum distribution at NNLONNLL, Phys.Rev. D90 (2014), no. 1 015017, [arXiv:1403.7196].
-  R. V. Harlander, H. Mantler, and M. Wiesemann, Transverse momentum resummation for Higgs production via gluon fusion in the MSSM, JHEP 1411 (2014) 116, [arXiv:1409.0531].
-  P. F. Monni and G. Zanderighi, On the excess in the inclusive cross section, JHEP 1505 (2015) 013, [arXiv:1410.4745].
-  P. Jaiswal and T. Okui, Explanation of the excess at the LHC by jet-veto resummation, Phys.Rev. D90 (2014), no. 7 073009, [arXiv:1407.4537].
-  T. Becher, R. Frederix, M. Neubert, and L. Rothen, Automated NNLL NLO resummation for jet-veto cross sections, Eur.Phys.J. C75 (2015), no. 4 154, [arXiv:1412.8408].
-  P. Nason, A New method for combining NLO QCD with shower Monte Carlo algorithms, JHEP 0411 (2004) 040, [hep-ph/0409146].
-  T. Sjostrand, S. Mrenna, and P. Z. Skands, PYTHIA 6.4 Physics and Manual, JHEP 0605 (2006) 026, [hep-ph/0603175].