NLO QCD corrections to W+W gamma and Z Z gamma production with leptonic decays
KA–TP–11–2009 SFB/CPP09105 FTUV–09–1101 IFUM938FT
NLO QCD corrections to and production with leptonic decays
G. Bozzi, F. Campanario, V. Hankele and D. Zeppenfeld
Dipartimento di Fisica, Università di Milano and INFN,
Sezione di Milano
Via Celoria 16, I20133 Milano, Italy
Institut für Theoretische Physik, Karlsruhe Institute of Technology,
Universität Karlsruhe, 76128 Karlsruhe, Germany
Departament de Física Teòrica and IFIC, Universitat de València  CSIC,
E46100, Burjassot, València, Spain
The computation of the QCD corrections to the cross sections for and production in hadronic collisions is presented. We consider the case of a real photon in the final state, but include full leptonic decays of the and bosons. Numerical results for the LHC and the Tevatron are obtained through a fully flexible parton level Monte Carlo based on the structure of the VBFNLO program, allowing an easy implementation of arbitrary cuts and distributions. We show the dependence on scale variations of the integrated cross sections and provide evidence that NLO QCD corrections strongly modify the LO predictions for observables at the LHC both in magnitude and in shape.
1 Introduction
The experimental precision that will be reached in cross section measurements at the CERN Large Hadron Collider (LHC) demands an effort of the theoretical community in providing accurate phenomenological predictions. NLO QCD corrections for cross sections and distributions have thus become mandatory, and many relevant processes at hadron colliders are now known to this accuracy.
Events with multiple gauge bosons in the final state provide an irreducible background to many new physics searches (see, for instance, [1] for a discussion of the relevant backgrounds in the search of New Physics at the LHC). In addition, the triple gauge couplings involved in the contributing diagrams allow for restrictive tests of the gauge sector of the Standard Model. The process + X is particularly important since it is also sensitive to four gauge boson couplings, namely the and vertices [2].
In this paper, we present the corrections for the processes
(1.1) 
Similar to previous work on triple weak boson production [3, 4, 5, 6], we find that the QCD corrections are sizeable and also modify the shape of the differential distributions for many observables: this proves that a simple rescaling of the LO results is not adequate and a full NLO Monte Carlo is needed for a quantitative determination of quartic couplings at the LHC. We have implemented our calculation within the VBFNLO framework [7], a fullyflexible parton level Monte Carlo program which allows the definition of arbitrary acceptance cuts and distributions.
The paper is organized as follows: in Section 2 we provide an example of the relevant Feynman diagrams at tree level, a short account of the strategies used to compute the real and virtual corrections and the various checks performed both internally and against other available codes. In Section 3 we show numerical results, including the scale variations of the LO and NLO integrated cross sections and some selected differential distributions. Conclusions are given in Section 4.
2 The calculation
We consider the processes (1) up to order in the limit where all fermions are massless. Among the LO diagrams (110 in the “” case, 336 in the “” case) we can distinguish three different topologies: they correspond to the cases when 1, 2, or 3 vector bosons are attached to the quark line (respectively case a), b), or c) of Fig.1, where examples for “” production are shown. The single and double vertex topologies of Fig.1 a) and b) also exist for “” production since we include offshell effects in our calculation, i.e. the photon or a boson can be radiated off a final state lepton.)
In all cases, the leptonic decays of the and bosons and combinations of subgraphs corresponding to decays of the same virtual , , or photon are factorized in the form of leptonic tensors and computed independently from the rest of the cross section, in analogy with the procedure in [8]: this reduces computational time since the same decay current may appear in many Feynman diagrams and in different subprocesses and can, thus, be computed once for each phasespace point and saved for later use. For the computation of the matrix elements, we use the helicity method introduced in [9].
The cross section at order contains contributions from real emission and virtual corrections from the interference of oneloop diagrams with the Born amplitude. We use the CataniSeymour dipole subtraction method [10] to handle the cancellation of infrared divergences between real emission and virtual corrections. Factorization of initial state collinear singularities into the parton distribution functions leads to additional “finite collinear terms”.
The NLO real corrections are generated by the emission of a gluon off the quark line or by the emission of a quark through a splitting in the initial state. The strategy of computing leptonic tensors separately proves to be particularly useful in this case, since we have a total of 272 (896) different Feynman diagrams contributing to the “” (“”) process when attaching an additional gluon to the Born diagrams.
In the case of real emissions, the definition of an isolated onshell photon requires further discussion. Photon emission collinear to a massless quark results in an additional singularity, which needs to be regularized. The simple rejection of events containing partons inside a cone of fixed width drawn around the photon is not a permissible choice since it will suppress the phase space for soft gluon emissions and spoil the cancellation of infrared divergences. In this paper we avoid the complications of quark fragmentation to photons by using an effective method which has been proposed by Frixione [11]. The precise definition of this photon isolation cut is discussed in Section 3.
The NLO virtual corrections come from the interference of leading order diagrams with oneloop diagrams obtained by attaching a gluon to the quarkantiquark line. In the case where only one vector boson is attached to the quark line, i.e. topology a), one only has vertex corrections, while diagrams belonging to topologies b) and c) give rise to boxes and pentagons. The tensor coefficients of the loop integrals have been computed à la PassarinoVeltman [12] up to box level, while the DennerDittmaier reduction procedure [13] has been used for the pentagons, which avoids the problem of small Gram determinants occurring in pentagons with planar configurations of the external momenta. DennerDittmaier reduction is now fully implemented in the public version of the VBFNLO code. The overall contribution from virtual corrections can be written as
(2.1) 
where is the Born amplitude and is the partonic centerofmass energy, i.e. the invariant mass of the final state system. The term is the finite part of the virtual corrections to 2 and 3 weak boson amplitudes as in Fig.1b) and c) which we call boxline and pentline contributions in the following, so named since boxes and pentagons constitute the most complex loop diagrams, respectively.
We have verified analytically and numerically that in order to obtain the above factorization formula of the infrared divergences against the born amplitude, the transversality property of the photon () must be used, otherwise, additional IR singularities which are not proportional to the born amplitude do appear. After checking that there are no additional terms due to the massless onshell photon, the analytical structure of can be obtained from the massive weak boson case considered in [5]. Nevertheless, we have recomputed equivalent pentline and boxline contributions for massless particles and found that they agree at the double precision accuracy level with those of [5], once the proper scalar integrals are used.
A further reduction of computing time is obtained for “” production by implementing a trick already used in [8, 4, 5]. We write the polarization vector of the vector bosons as
(2.2) 
where is the momentum of the massive boson (). A pentagon contracted with one of the external momenta can always be reduced to a difference of box diagrams. Thus we effectively isolate a true pentagon contribution and shift the rest of the 5point contribution to the less timeconsuming boxes. The remainder is chosen in such a way that
(2.3) 
i.e., the timecomponent of the shifted polarization vector is zero in the centerofmass system of the Vpair. This choice yields particularly small pentline contributions which, hence, can be determined with lower Monte Carlo statistics. For “” production, the pentline contributions group into gauge invariant subsets which are invariant under the replacement of Eq.(2.2). Thus the trick does not speed up the code but provides a consistency check for the calculation.
We have performed a number of checks on our final results. First, we compared all tree level matrix elements (including the ones for real emission corrections) against amplitudes generated by MadGraph[14]: the agreement is at the level of machine precision (). Second, we have compared the LO cross sections both for and against HELAC[15] and found agreement at the per mill level. Third, we have compared integrated tree level cross sections for and against MadEvent and Sherpa [16] finding agreetment at the per mill level. Finally, we also performed gauge tests based on Ward identities for the NLO virtual corrections at different levels of complexity.
3 Results
Numerical results have been obtained through the implementation of our calculation into a NLO Monte Carlo program based on the structure of the VBFNLO code [7]. This and other processes with a real photon in the final state will be included in the public version of VBFNLO in a future release.
The value of the total cross section for “” production, including leptonic Zdecays, turns out to be less than a femtobarn at LHC energies: since it will probably be too hard to detect, at least in an early luminosity run, in this section we focus on the process.
For the electroweak parameters we use the and boson masses and the Fermi constant as input. From these we derive the electromagnetic coupling and the weak mixing angle via tree level relations, i.e. we use
(3.1) 
We do not consider bottom and top quark effects. The remaining quarks are assumed to be massless and we work in the approximation where the CKM matrix is the identity matrix. We choose the invariant “” mass as the central value for the factorization and renormalization scales:
(3.2) 
We use the CTEQ6L1 parton distribution function at LO and the CTEQ6M set with at NLO [17].
Since the photon can be emitted either from the initial quark line or from each of the two final state charged leptons, we divide the phase space into three separate regions to consider all the possibilities and then sum the contributions to get the total result. The regions are generated as (approximately) onshell and threebody decay and triple electroweak boson production, respectively. We impose a set of minimal cuts on the rapidity, , and the transverse momenta , , of the photon and the charged leptons which are designed to represent typical experimental requirements. Furthermore, leptons, photons and jets must be well separated in the rapidityazimuthal angle plane. Specifically, the cuts we impose are
(3.3) 
where, in our simulations, a jet is defined as a colored parton of transverse momentum GeV. For photon isolation, we implement the procedure introduced in [11]: if is a parton with transverse energy and a separation with a photon of transverse momentum , then the event is accepted only if
(3.4) 
where is a fixed separation that we set equal to 0.7. A quick look at Eq.(3.4) reveals that a sufficiently soft gluon can be arbitrarily close to the photon axis, while the energy of an exactly collinear parton must be vanishing in order to pass the isolation cut. Collinearonly events (leading to fragmentation contributions) are thus rejected, while soft emissions are retained, as desired. Additionally, for “” production we impose that the invariant mass of any combination of two charged leptons, , be larger than 15 GeV in order to avoid virtualphoton singularities in at a low . In Table 1, we give results for the integrated cross sections for “” and “” production at the LHC for the given cuts as well as for a harder cut on the photon transverse momentum of GeV. In Table 2, the numbers for the Tevatron are presented for the cuts as in Eqs.(3.3,3.4), but for less restrictive transverse momentum cuts, GeV and GeV. The centerofmass energy is set to TeV for the LHC and TeV for Tevatron collisions, respectively.
LHC  LO [fb]  NLO [fb] 

1.695  2.881  
GeV  

0.9580  1.738 
GeV  

0.07786  0.1062 
GeV  

0.03969  0.05577 
GeV 
Tevatron 
LO [fb]  NLO [fb] 

0.9015  1.271  
GeV  

0.3755  0.5342 
GeV  

0.05469  0.07681 
GeV  

0.02440  0.03441 
GeV 
In the following we show results for the process , i.e we do not consider interference effects due to identical leptons, which generally will be small. The cross sections have been multiplied by a combinatorial factor of 4, however, and correspond to the production of any combination of electrons or muons.
In order to estimate the scale uncertainty of the total cross section, we study the variation of the numerical results in the interval
(3.5) 
where is the “” or “” invariant mass. The variation of the total “” cross section at LO and NLO is shown on the left panel of Figure 2. It is evident that the LO scale variation greatly underestimates the size of NLO corrections, which give a Kfactor of 1.70 for . The NLO scale uncertainty is about 6 when varying the factorization and the renormalization scale up and down by a factor 2 around the reference scale , slightly smaller than for the LO uncertainty, and is dominated by the dependence on , which gives a negative slope with increasing energy, while it has a quite flat dependence on .
On the right panel of Figure 2, the scale dependence and relative size of the different NLO contributions are shown. One clearly sees that almost the entire scale variation of the integrated NLO cross section is accounted for by the real emission contributions, defined here as the real emission cross section minus the CataniSeymour subtraction terms plus the finite collinear terms. The virtual contributions proportional to the Born matrix element (due to the term in Eq. (2.1)) constitute the bulk of the NLO corrections, being roughly twice as large as the real emissions, while the finite virtual remainders due to boxes and pentagons only represent 3% and less than 1% of the total result, respectively, and their scale dependence is basically flat.
Similar considerations apply for “” production where a combinatorial factor of 2 has been included for . In Figure 3, one observes that the overall scale variation decreases when varying the factorization and the renormalization scale up and down by a factor 2 around the reference scale from roughly 7% at LO to 1% at NLO, and we have a Kfactor of 1.36 for . The main difference with respect to the “” case is the numerical contribution of the real corrections, which amounts to 20 % of the total contribution or less.
The NLO contributions cannot be described by a simple rescaling of the LO results since the corrections show a strong dependence on the phase space region under investigation. As practical examples, we plot two differential distributions at LO and NLO together with the associated Kfactor, defined as
(3.6) 
The observables, , considered are the transversemomentum of the photon and the separation between the softest (lowest) lepton and the photon in production at the LHC. In Figure 4 we show the transversemomentum distribution of the onshell photon: the corresponding Kfactor increases from 1.5 at GeV to 2.1 at GeV. In Figure 5, we consider the separation between the photon and the softest lepton and a similar variation of the Kfactor in the range is observed.
4 Conclusions
We have calculated the NLO QCD corrections to the processes and with full leptonic decays of the and bosons. These processes can be relevant both for New Physics searches (as a background) and for the measurement of quartic gauge couplings (as a signal) at the LHC.
Our numerical results show sizeable cross section increases with respect to the LO calculation, ranging from 40% to as much as 100% in certain regions of phase space. and production at the LHC provide additional examples of cross sections whose theoretical errors at LO are substantially underestimated by considering scale variations only: the LO factorization scale variation is much smaller than the NLO correction. Remaining NLO scale variations are at the 10% level (e.g. for the integrated production cross section at the LHC when varying by a factor of 2 around the reference scale ).
Given the size of the higherorder corrections and, in particular, their strong dependence on the observable and on different phase space regions under investigation, a fullyexclusive NLO parton Monte Carlo for and production is required to match the expected precision of the LHC measurements. We plan to incorporate this and other processes with a final state photon into the VBFNLO package in the near future.
Acknowledgments
We would like to thank Christoph Englert and Michael Rauch for helpful discussions, and Malgorzata Worek, for the comparison with HELAC. This research was supported in part by the Deutsche Forschungsgemeinschaft via the Sonderforschungsbereich/Transregio SFB/TR9 “Computational Particle Physics”. F.C. acknowledges a postdoctoral fellowship of the Generalitat Valenciana (Beca Postdoctoral d’Excel.lència) and partial support by European FEDER and Spanish MICINN under grant FPA200802878. The Feynman diagrams in this paper were drawn using Axodraw [18].
References
 J. M. Campbell, J. W. Huston and W. J. Stirling, Rept. Prog. Phys. 70, 89 (2007) [arXiv:hepph/0611148].
 S. Godfrey, arXiv:hepph/9505252; P. J. Dervan, A. Signer, W. J. Stirling and A. Werthenbach, J. Phys. G 26 (2000) 607 [arXiv:hepph/0002175]; O. J. P. Eboli, M. C. GonzalezGarcia, S. M. Lietti and S. F. Novaes, Phys. Rev. D 63 (2001) 075008 [arXiv:hepph/0009262]; P. J. Bell, arXiv:0907.5299 [hepph].
 A. Lazopoulos, K. Melnikov and F. Petriello, Phys. Rev. D 76 (2007) 014001 [arXiv:hepph/0703273].
 V. Hankele and D. Zeppenfeld, Phys. Lett. B 661 (2008) 103 [arXiv:0712.3544 [hepph]].
 F. Campanario, V. Hankele, C. Oleari, S. Prestel and D. Zeppenfeld, Phys. Rev. D 78 (2008) 094012 [arXiv:0809.0790 [hepph]].
 T. Binoth, G. Ossola, C. G. Papadopoulos and R. Pittau, JHEP 0806 (2008) 082 [arXiv:0804.0350 [hepph]].
 K. Arnold et al., Comput. Phys. Commun. 180 (2009) 1661 [arXiv:0811.4559 [hepph]].
 B. Jager, C. Oleari and D. Zeppenfeld, JHEP 0607 (2006) 015 [arXiv:hepph/0603177]; B. Jager, C. Oleari and D. Zeppenfeld, Phys. Rev. D 73, 113006 (2006) [arXiv:hepph/0604200]; G. Bozzi, B. Jager, C. Oleari and D. Zeppenfeld, Phys. Rev. D 75 (2007) 073004 [arXiv:hepph/0701105].
 K. Hagiwara and D. Zeppenfeld, Nucl. Phys. B 274 (1986) 1; K. Hagiwara and D. Zeppenfeld, Nucl. Phys. B 313 (1989) 560.
 S. Catani and M. H. Seymour, Nucl. Phys. B 485 (1997) 291 [Erratumibid. B 510 (1998) 503] [arXiv:hepph/9605323].
 S. Frixione, Phys. Lett. B 429, 369 (1998) [arXiv:hepph/9801442].
 G. Passarino and M. J. G. Veltman, Nucl. Phys. B 160 (1979) 151.
 A. Denner and S. Dittmaier, Nucl. Phys. B 658, 175 (2003) [arXiv:hepph/0212259]; A. Denner and S. Dittmaier, Nucl. Phys. B 734, 62 (2006) [arXiv:hepph/0509141].
 T. Stelzer and W. F. Long, Comput. Phys. Commun. 81 (1994) 357 [arXiv:hepph/9401258]; F. Maltoni and T. Stelzer, JHEP 0302 (2003) 027 [arXiv:hepph/0208156].
 A. Cafarella, C. G. Papadopoulos and M. Worek, Comput. Phys. Commun. 180 (2009) 1941 [arXiv:0710.2427 [hepph]]; C. G. Papadopoulos and M. Worek, Eur. Phys. J. C 50 (2007) 843 [arXiv:hepph/0512150]. A. Kanaki and C. G. Papadopoulos, Comput. Phys. Commun. 132 (2000) 306 [arXiv:hepph/0002082].
 T. Gleisberg, S. Hoche, F. Krauss, M. Schonherr, S. Schumann, F. Siegert and J. Winter, JHEP 0902 (2009) 007 [arXiv:0811.4622].
 J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadolsky and W. K. Tung, JHEP 0207 (2002) 012 [arXiv:hepph/0201195].
 J. A. M. Vermaseren, Comput. Phys. Commun. 83 (1994) 45.