Top-pair production at the LHC through NNLO QCD and NLO EW

A bstractIn this work we present for the first time predictions for top-quark pair differential distributions at the LHC at NNLO QCD accuracy and including EW corrections. For the latter we include not only contributions of Oαs2α$$ \mathcal{O}\left({\alpha}_s^2\alpha \right) $$, but also those of order Oαsα2$$ \mathcal{O}\left({\alpha}_s{\alpha}^2\right) $$ and Oα3$$ \mathcal{O}\left({\alpha}^3\right) $$. Besides providing phenomenological predictions for all main differential distributions with stable top quarks, we also study the following issues. 1) The effect of the photon PDF on top-pair spectra: we find it to be strongly dependent on the PDF set used — especially for the top pT distribution. 2) The difference between the additive and multiplicative approaches for combining QCD and EW corrections: with our scale choice, we find relatively small differences between the central predictions, but reduced scale dependence within the multiplicative approach. 3) The potential effect from the radiation of heavy bosons on inclusive top-pair spectra: we find it to be, typically, negligible.


Introduction
The availability of NNLO QCD predictions for stable top-pair production at the LHC, both for the total cross-section [1][2][3][4] with NNLL soft-gluon resummation [5,6] and for all the main differential distributions [7,8], has made it possible to compare Standard Model (SM) theory with LHC data at the few-percent level accuracy. Such a high precision has led, among others, to further scrutiny of the differences between LHC measurements [9] and the ability of Monte-Carlo event generators to describe hadronic tt production. As a result of these ongoing studies, new MC developments are taking place, such as the incorporation of non-resonant and interference effects [10,11], which builds upon previous works that included NLO top decay corrections through-fixed order [12][13][14][15][16] and/or showered [17][18][19] calculations.
One of the remaining ways for further improving SM theory predictions is by consistently including the so called Electro-Weak (EW) corrections on top of the NNLO QCD ones. Weak [20][21][22][23][24][25][26][27][28], QED [29] and EW (weak+QED) [30][31][32][33][34] corrections to top-quark pair production have been known for quite some time, and also EW corrections to the fully off-shell dilepton signature are nowadays available [35]. As it has been documented in the literature, although EW effects are rather small at the level of total cross-section, they can have a sizeable impact on differential distributions and also on the top-quark charge asymmetry.
The goal of this work is to consistently merge existing NNLO QCD predictions with EW corrections into a single coherent prediction and to study its phenomenological impact. This is achieved by combining the NNLO QCD predictions from ref. [8] with the complete LO and NLO contributions derived within the framework of ref. [34]. Specifically, we include the NLO EW effects of O(α 2 s α), all subleading NLO (O(α s α 2 ) and O(α 3 )) terms as well as the LO (O(α s α) and O(α 2 )) contributions.
Another motivation for this study stems from recent developments in understanding the photon content of the proton [36,37]. As shown in ref. [34], depending on the PDF set, photon-initiated contributions can be numerically significant in some regions of phase space. 1 If the photon density from the NNPDF3.0QED set [40,41] is employed, the photon-initiated contribution is large in size and of opposite sign with respect to the Sudakov EW corrections, leading to the almost complete cancellation of the two effects. Nevertheless, large PDF uncertainties from the photon PDF are still present after this cancellation. On the other hand, theoretical consensus about the correctness of the novel approach introduced in ref. [37] appears to have emerged by now. 2 The PDF set provided with ref. [37], named LUXQED, includes a photon PDF whose central value and relative uncertainty are both much smaller than in the case of NNPDF3.0QED. Thus, at variance with the NNPDF3.0QED set, neither large cancellation between Sudakov effects and photon-induced contributions nor large photon PDF uncertainty is present in LUXQED-based predictions.
In order to document the ambiguity arising from the differences between the photon densities in the available PDF sets, with the exception of sec. 2, in this work we always give predictions for top-pair differential distributions at the LHC based on the LUXQED [37] and NNPDF3.0QED [40,41] PDF sets. We believe that our findings will provide a valuable input to future PDF determinations including EW effects. This paper is organised as follows: sec. 2 is devoted to the phenomenological study of our combined QCD and EW predictions for the LHC at 13 TeV. The reasons behind some of the choices made in sec. 2 -like the choice of PDF set and combination approach -are revealed in sec. 3, where we compare in great detail two approaches for combining NNLO QCD and EW corrections in top-pair differential distributions. The so-called additive approach is discussed in sec. 3.1, while the multiplicative one in sec. 3.2. Section 4 is dedicated to studying the impact of the photon PDF on top-pair spectra. In sec. 5 we provide an estimate of the impact of inclusive Heavy Boson Radiation (HBR), namely the contribution from ttV final states with V = H, W ± , Z. While most of the notation is introduced in the main text some technical details are delegated to Appendix A.

Phenomenological predictions for the LHC at 13 TeV
In this section we present predictions for tt distributions for the LHC at 13 TeV at NNLO QCD accuracy including also EW corrections. We focus on the following distributions: the top-pair invariant mass m(tt), the top/antitop average transverse momentum (p T,avt ) and rapidity (y avt ) and the rapidity y(tt) of the tt system. The p T,avt (y avt ) distributions are calculated not on an event-by-event basis but by averaging the results of the histograms for the transverse momentum (rapidity) of the top and the antitop.
Our calculation is performed using the following input parameters m t = 173.3 GeV , m H = 125.09 GeV , m W = 80.385 GeV , m Z = 91.1876 GeV , (2.1) while all other fermion masses are set to zero. All masses are renormalised on-shell and all decay widths are set to zero. The renormalisation of α s is performed in the 5-flavour scheme while EW input parameters and the associated α renormalisation condition are in the G µ -scheme, with G µ = 1.1663787 · 10 −5 GeV −2 . (2. 2) The EW corrections have been calculated in a completely automated way via an extension of the MadGraph5 aMC@NLO code [43] that has been already validated in refs. [44,45], and in ref. [46] for the calculation of the complete NLO corrections.
We work with dynamical renormalisation (µ r ) and factorisation (µ f ) scales. Their common central value is defined as where m T,t ≡ m 2 t + p 2 T,t and m T,t ≡ m 2 t + p 2 T,t are the transverse masses of the top and antitop quarks. As already mentioned, p T,avt and y avt distributions are obtained by averaging the top and antitop distributions for the transverse momentum and rapidity, respectively.
These scale choices have been motivated and studied at length in ref. [8]. In all cases theoretical uncertainties due to missing higher orders are estimated via the 7-point variation of µ r and µ f in the interval {µ/2, 2µ} with 1/2 ≤ µ r /µ f ≤ 2. We remark that the combination of QCD and EW corrections is independently performed for each value of µ f,r .
For theoretical consistency, a set of PDFs including QED effects in the DGLAP evolution should always be preferred whenever NLO EW corrections are computed. At the moment, the only two NNLO QCD accurate PDF sets that include them are NNPDF3.0QED and LUXQED. 3 Both sets have a photon density, which induces additional contributions to tt production [29,34].
As motivated and discussed at length in sec. 3, the phenomenological predictions in this section are based on the LUXQED PDF set and on the multiplicative approach for combining QCD and EW corrections, which we will denote as QCD × EW. We invite the dσ/dp T   interested reader to consult secs. 3 and 4 where detailed comparisons between the two PDF sets as well as between the two approaches for combining QCD and EW corrections can be found.
From the plots shown in fig. 1 we conclude that the impact of the EW corrections relative to NNLO QCD depends strongly on the kinematic distribution. The smallest impact is observed in the two rapidity distributions: the relative effect for y avt is around 2 permil and is much smaller than the scale uncertainty. The y(tt) distribution is slightly more sensitive, with a relative impact of slightly above 1% for large values of y(tt). This correction is also well within the scale-variation uncertainty band. The impact of EW corrections on the m(tt) distribution is larger. Relative to NNLO QCD it varies between +2% at the absolute threshold and -6% at high energies. Still, this correction is well within the scale variation uncertainty. The small sensitivity of the y avt and y(tt) distributions to EW corrections supports the findings of ref. [9] where these two distributions were used for constraining the PDF of the gluon.
The largest correction due to EW effects is observed in the p T,avt distribution. Relative to NNLO QCD, the correction ranges from +2% at low p T,avt to -25% at p T,avt ∼ 3 TeV. The correction is significant and is comparable to the scale variation band already for p T,avt ∼ 500 GeV. Overall, the EW contribution to the p T,avt distribution is as large as the total theory uncertainty band in the full kinematic range p T,avt ≤ 3 TeV considered in this work.
The fraction of the theory uncertainty induced by PDFs is strongly dependent on kinematics. For the y avt and y(tt) distributions, the PDF error is slightly smaller than the scale uncertainty for central rapidities, but is larger in the peripheral region, especially for the y(tt) distribution. The PDF uncertainty becomes the dominant source of theory error in the p T,avt distribution for p T,avt as large as 500 GeV, while for the m(tt) distribution it begins to dominate over the scale uncertainty for m(tt) ∼ 2.5 TeV.
There are many applications for the results derived in this work. Examples are: inclusion of EW effects in PDF determinations from LHC tt distributions, high-mass LHC searches, precision SM LHC measurements and benchmarking of LHC event generators. A practical and sufficiently accurate procedure for the utilisation of our results could be as follows. One starts by deriving an analytic fit for the QCD × EW/QCD K-factor (all Kfactors in fig. 1 are available in electronic form 4 ); as evident from fig. 1 it is a very smooth function for all four differential distributions. Under the assumption that this K-factor is PDF independent, such an analytic fit could then be used to rescale the NNLO-QCDaccurate differential distributions derived with any PDF set from existing NNLO QCD [50] fastNLO [51,52] tables. Regarding the PDF error of NNLO QCD differential distributions, it can be calculated very fast with any PDF set with the help of the fastNLO tables of ref. [50]. As we show in the following the PDF error of the QCD and combined QCD and EW predictions is almost the same, especially for the LUXQED PDF set used for our phenomenological predictions.

Comparison of two approaches for combining NNLO QCD predictions and EW corrections
In this work we compare two approaches for combining QCD and EW corrections. For brevity, we will refer to them as additive and multiplicative approaches. As already mentioned, the results presented in sec. 2 have been calculated using the multiplicative approach.
In the additive approach the NNLO QCD predictions (defined as the complete set of O(α n s ) terms up to n = 4) are combined with all possible remaining LO and NLO terms arising from QCD and electroweak interactions in the Standard Model. In other words, at LO we include not only the purely QCD O(α 2 s ) contribution, but also all O(α s α) and O(α 2 ) terms. Similarly, at NLO we take into account not only the NLO QCD O(α 3 s ) contribution but also the O(α 2 s α) one, the so-called NLO EW, as well as the subleading contributions of O(α s α 2 ) and O(α 3 ). For brevity, we will denote as "EW corrections" the sum of all LO and NLO terms of the form O(α m s α n ) with n > 0. Moreover, when we will refer to "QCD" results, we will understand predictions at NNLO QCD accuracy. For a generic observable Σ in the additive approach we denote the prediction at this level of accuracy as Σ QCD+EW .
In the multiplicative approach one assumes complete factorisation of NLO QCD and NLO EW effects. This approach is presented in sec. 3.2 and is denoted as Σ QCD×EW . The precise definition of the various quantities mentioned in the text is given in appendix A where an appropriate notation for the classification of the different contributions is introduced. Here, we just state the most relevant definitions for the following discussion, where Σ denotes a generic observable in tt production and K NLO QCD is the standard NLO/LO K-factor in QCD. A variation of the multiplicative approach denoted as Σ QCD 2 ×EW will also be considered; it is defined similarly to Σ QCD×EW in eq. (3.1) but with NNLO/LO QCD K-factor.
As it has been discussed in ref. [34], the usage of different PDF sets leads to a very different impact of the photon-induced contribution on tt distributions. While in the case of NNPDF3.0QED the impact of photon-induced contributions is relatively large and with very large uncertainties, in the case of LUXQED it is expected to be negligible. For this reason in the rest of this work we always show predictions with both PDF sets.

Additive combination
Distributions for p T,avt and m(tt) are shown in fig. 2, while the y avt and y(tt) distributions are shown in fig. 3. The format of the plots for all distributions is as follows: for each observable, we show two plots side-by-side, with the same layout. The plot on the lefthand side shows predictions obtained using the LUXQED set, while for the one on the right the NNPDF3.0QED set is employed. Results at NNLO QCD accuracy are labelled as "QCD" while the combination of NNLO QCD predictions and EW corrections in the additive approach are labelled as "QCD+EW".
In each plot the three insets display ratios of different quantities 5 over the centralscale QCD result (i.e., in the case of LUXQED, the black line in the main panel of fig. 1).  In the first inset we show the scale uncertainty due to EW corrections alone (red band), without QCD contributions (Σ EW using the notation of Appendix A). This quantity can be compared to the scale uncertainty of the QCD prediction at NNLO accuracy (grey band). In the second inset we present the scale-uncertainty band (red) for the combined QCD+EW prediction. The grey band corresponds to the NNLO QCD scale-uncertainty band already shown in the first inset. The third inset is equivalent to the second one, but it shows the PDF uncertainties. We combine, for each one of the PDF members, the QCD prediction and the EW corrections into the QCD+EW result. The PDF uncertainty band of the QCD+EW prediction is shown in red while the grey band corresponds to the PDF uncertainty of the QCD prediction. For all insets, when the grey band is covered by the red one, its borders are displayed as black dashed lines. As can be seen in figs. 2 and 3, the effect of EW corrections is, in general, within the NNLO QCD scale uncertainty. A notable exception is the case of the p T,avt distribution with LUXQED. In the tail of this distribution the effect of Sudakov logarithms is large and negative, of the order of -(10-20%), and is not compensated by the photon-induced contribution. On the contrary, in the case of NNPDF3.0QED, photon-induced contributions mostly compensate the negative corrections due to Sudakov logarithms. As it has already been noted in ref. [34], with this PDF set, the effect of photon-induced contributions is not negligible also for large values of m(tt), y avt and y(tt).
As it can be seen in the first inset, in the large p T,avt regime the scale dependence of  the EW corrections alone is of the same size as, or even larger than, the scale variation at NNLO QCD accuracy. For this reason, as evident from the second inset, the scale uncertainty of the combined QCD+EW prediction is much larger than in the purely QCD case, both with the LUXQED and NNPDF3.0QED PDF sets. This feature is present only in the tail of the p T,avt distribution. The PDF uncertainties (third inset) for all distributions do not exhibit large differences between QCD and QCD+EW predictions, despite the fact that the photon-induced contribution in NNPDF3.0QED is large and has very large PDF uncertainty (relative to LUXQED).

Multiplicative combination and comparison with the additive one
The additive approach Σ QCD+EW for combining QCD and EW corrections discussed in sec. 3.1 is exact to the order at which the perturbative expansion of the production crosssection is truncated. An alternative possibility for combining QCD and EW corrections is what we already called the multiplicative approach, Σ QCD×EW . This approach is designed to approximate the leading EW corrections at higher orders. In the case of tt production these are NNLO EW contributions of order O(α 3 s α). The multiplicative approach is motivated by the fact that soft QCD and EW Sudakov logarithms factorise, with the latter typically leading to large negative corrections for boosted kinematics. Thus, when dominant NLO EW and NLO QCD corrections are at the same time induced by these two effects, the desired fixed order can be very well approximated via rescaling NLO EW corrections with NLO QCD K-factors. 6 Otherwise, if one is in a kinematical regime for which the dominant NLO EW or NLO QCD corrections are of different origin (i.e. not Sudakov or soft), the difference between the multiplicative and additive approaches given by the term Σ mixed in eq. (A.8) can be considered as an indication of theory uncertainty in that kinematics. It must be stressed that the perturbative orders involved in the additive approach are included exactly also in the multiplicative approach; the only addition the multiplicative approach introduces on top of the additive one is the approximated O(α 3 s α) contribution. One of the advantages of the multiplicative approach is the stabilisation of scale dependence. As we saw in sec. 3.1, when QCD and EW corrections are combined in the additive approach, the scale dependence at large p T,avt can exceed that of the NNLO QCD prediction. On the other hand, the large p T,avt limit is precisely the kinematic regime where the multiplicative approach is a good approximation and can be trusted: at large p T,avt the NLO EW and NLO QCD corrections are mainly induced by Sudakov logarithms and soft emissions, respectively, and as we just pointed out these two contributions factorise. The presence of large Sudakov logarithms in the NLO EW result at large p T,avt is easy to see since for Born kinematics large p T,t implies large p T,t which, in turn, implies largê s,t andû Mandelstam variables. That NLO QCD corrections at large p T,avt are mainly of soft origin can be shown with an explicit NLO calculation; by applying appropriate cuts on the jet, the top and/or the antitop one can easily see that the differential cross-section is dominated by kinematic configurations containing almost back-to-back hard top and antitop and a jet with small p T . Plots demonstrating this can be found in footnote 4.
In the following, for all observables Σ considered in this work, we present predictions in the multiplicative approach denoted as Σ QCD×EW . As a further check of the stability of the multiplicative approach we display also the quantity Σ QCD 2 ×EW , whose precise definition can be found in appendix A. Σ QCD 2 ×EW is defined analogously to Σ QCD×EW , but by rescaling NLO EW corrections via NNLO QCD K-factors. By comparing Σ QCD×EW and Σ QCD 2 ×EW one can further estimate the uncertainty due to mixed QCD-EW higher orders. Figure 4 shows the p T,avt and m(tt) distributions, while fig. 5 refers to y avt and y(tt). As in sec. 3.1, the plots on the left are produced using the LUXQED PDF set, while those on the right using the NNPDF3.0QED PDF set. We next describe the format of the plots.
Each plot consists of five insets, which all show ratios of different quantities over the central value of Σ QCD . In the first inset we compare the central-scale results for the three alternative predictions: Σ QCD+EW /Σ QCD (red line), Σ QCD×EW /Σ QCD (green line) and Σ QCD 2 ×EW /Σ QCD (violet line). These quantities are further displayed in the second, third and fourth inset, respectively, where not only the central value but also the scale dependence of the numerator is shown. In all cases we calculate the scale-uncertainty band as a scale-by-scale combination and subsequent variation in the 7-point approach. Scale variation bands have the same colour as the corresponding central-value line. For comparison we also display (grey band) the relative scale uncertainty of Σ QCD . Thus, the second inset is exactly the same as the second inset in the corresponding plots in sec. 3.1.  The last inset shows a comparison of the ratio Σ QCD+EW /Σ QCD including (red line) or not (orange line) the contribution Σ res , where "res" stands for residual and denotes the fact that Σ res are contributions to Σ EW that are expected to be small, regardless of the PDF set used (see eq. (A.6)).
As expected, the multiplicative approach shows much smaller dependence on the scale variation. This is particularly relevant for the tail of the p T,avt distribution, where the scale uncertainty of Σ EW alone is comparable in size with the one of Σ QCD ; with this reduction of the scale uncertainty the Σ QCD×EW and Σ QCD uncertainty bands do not overlap when LUXQED is used. In the case of m(tt) and y avt distributions, the Σ QCD×EW central-value predictions are typically larger in absolute value than those of Σ QCD+EW , while they are all almost of the same size for the y(tt) distribution. In the case of y avt the difference between the additive and multiplicative approaches is completely negligible compared to their scale uncertainty. Therefore, besides the kinematic region where Sudakov effects are the dom-  Figure 5. As in fig. 4 but for the y avt and y(tt) differential distributions.
inant contribution, the multiplicative and additive approaches are equivalent. Moreover, the difference between Σ QCD×EW and Σ QCD 2 ×EW is in general small; a sizeable difference between their scale dependences can be noted only in the tail of the p T,avt distribution. For all the reasons mentioned above we believe that the multiplicative approach should be preferred over the additive one and, indeed, it has been used for the calculation of our best predictions in sec. 2. As can be seen from figs. 4 and 5 and their threshold-zoomed-in versions in footnote 4 the difference between Σ QCD+EW and Σ QCD×EW for non-boosted kinematics is much smaller than the total theory uncertainty (scale+PDF) shown in fig. 1. Thus, the difference between the two approaches can be safely ignored in the estimation of the theory uncertainty. One should bear in mind that this conclusion depends on the choice of scale, which in our case, as explained in ref. [8], is based on the principle of fastest convergence. A different scale choice with larger K factors would likely artificially enhance the difference between Σ QCD+EW and Σ QCD×EW .
In the last inset in figs. 4 and 5 we compare the quantities Σ EW and Σ EW − Σ res , where the Σ res contribution is exactly included in both the additive and multiplicative approaches. As expected, one can see that the Σ res contribution is typically flat and very small. The only exception is the m(tt) distribution where a visible difference between the two curves (Σ EW and Σ EW − Σ res ) is present, especially in the tail. The Σ res contribution includes the squared EW tree-level diagrams, the O(α 2 ) contribution denoted as Σ

Impact of the photon PDF
In this section we quantify the impact on tt differential distributions of the difference between the photon densities provided by the LUXQED and NNPDF3.0QED PDF sets.
In other words, we repeat the study performed in ref. [34] for these two PDF sets since they were not considered in that work. We compare the size of the electroweak corrections with and without the photon PDF for both PDF sets. In each plot of fig. 6 we show the relative impact induced by the electroweak corrections (the ratio Σ EW /Σ QCD ; see definitions in Appendix A) for four cases: NNPDF3.0QED setting the photon PDF equal to zero (red) or not (green), and LUXQED setting the photon PDF equal to zero (violet) or not (blue). For the cases including the photon PDF, we also show the PDF-uncertainty band of Σ EW . The impact of the photon-induced contribution can be evaluated via the difference between the green and red lines in the case of NNPDF3.0QED and the difference between the blue and violet lines in the case of LUXQED. As can be seen in fig. 6, the impact of the photon PDF on the p T,avt , m(tt), y avt and y(tt) distributions is negligible in the case of LUXQED, while it is large and with very large uncertainties for the case of NNPDF3.0QED, as already pointed out in ref. [34] for NNPDF2.3QED. At very large p T,avt and m(tt) also LUXQED show a non-negligible relative PDF uncertainty, which is not induced by the photon but from the PDFs of the coloured partons at large x. We checked that a similar behaviour is exhibited also by NNPDF3.0QED when its photon PDF is set to zero.

Contributions from heavy boson radiation
In the calculation of EW corrections to QCD processes the inclusion of real emissions of massive gauge bosons (heavy boson radiation or HBR) is not mandatory since, due to the finite mass of the gauge bosons, real and virtual weak corrections are separately finite (albeit the virtual corrections are enhanced by large Sudakov logarithms). Furthermore, such emissions are typically resolved in experimental analyses and are generally considered  Figure 6. Impact of the photon PDF on the p T,avt , m(tt), y avt and y(tt) differential distributions at 13 TeV. The format of the plots is described in the text.
as a different process ttV (+X) with V = H, W ± , Z. For these reasons, the results in sec. 2 do not include HBR contributions. It is, nonetheless, interesting to estimate the contribution of HBR to inclusive tt production. Our motivation is threefold: First, resolved or not, HBR is a legitimate contribution to the tt(+X) final state considered in this work. Secondly, it is clear that one cannot guarantee that HBR is resolved with 100% efficiency. Therefore, it is mandatory to have a prior estimate for the size of the effect. Finally, we are unaware of prior works where the HBR contribution has been estimated in inclusive tt production. Recently, refs. [44,53] have provided estimates for HBR in the processes ttV (+X), with V = H, Z, W .
We have investigated the impact of HBR on all four distributions considered in this work: p T,avt , m(tt), y avt and y(tt). Our results are shown in fig. 7, where we plot the effect of HBR on the central scale normalised to the QCD prediction. We show separately the LO HBR effect of order O(α 2 s α) as well as the NLO QCD HBR prediction which includes terms of order O(α 3 s α). As a reference we also show the EW corrections for tt. In our calculations we include HBR due to H, W and Z. We are fully inclusive in HBR, i.e., no cuts on the emitted heavy bosons are applied. Clearly, any realistic experimental analysis will require an estimate of HBR subject to experimental cuts, but such an investigation would be well outside the scope of the present work.
From fig. 7 we conclude that the effect of HBR is generally much smaller than the EW corrections. In particular, higher-order QCD corrections to HBR are completely negligible, i.e. HBR is well described in LO for all the tt inclusive distributions and for the full kinematic ranges considered here. The absolute effect of HBR on the p T,avt distribution is positive and small; it never exceeds 2-3% (relative to the tt prediction at NNLO QCD accuracy) and is always much smaller than the EW correction. The only distribution where the HBR contribution is not negligible compared to the EW one is m(tt) computed with LUXQED. For this distribution the HBR correction is positive and only about half the absolute size of the (negative) EW correction. Still, the absolute size of the HBR, relative to the prediction at NNLO QCD accuracy, is within 1% and so its phenomenological relevance is unclear. The impact of HBR on the two rapidity distributions is tiny, typically within 3 permil of the NNLO QCD prediction.

Conclusions
In this work we derive for the first time predictions for all main top-quark pair differential distributions 7 with stable top quarks at the LHC at NNLO QCD accuracy and including the following EW corrections: the NLO EW effects of O(α 2 s α), all subleading NLO terms of order O(α s α 2 ) and O(α 3 ) as well as the LO contributions of order O(α s α) and O(α 2 ).
We present a detailed analysis of top-pair production at the LHC at 13 TeV and we find that the effect of EW corrections on differential distributions with stable top quarks is in general within the current total (scale+PDF) theory uncertainty. A notable exception is the p T,avt distribution in the boosted regime where the effect of EW corrections is significant with respect to the current total theory error. We have checked that similar conclusions apply also for LHC at 8 TeV. All results derived in this work in the multiplica-tive approach, for both 8 and 13 TeV, are available in electronic form 4 as well as with the ArXiv submission of this paper.
Providing phenomenological predictions for the LHC is only one of the motivations for the present study. In this work we also quantify the impact of the photon PDF on top-pair differential distributions and study the difference between the additive and multiplicative approaches for combining QCD and EW corrections. Moreover, we analyse the contribution from inclusive Heavy Boson Radiation on inclusive top-pair differential distributions.
In order to quantify the impact of the photon PDF, we use two recent PDF sets whose photon components are constructed within very different approaches. The first set, LUXQED, is based on the PDF4LHC15 set [54] and adds to it a photon contribution that is derived from the structure function approach of ref. [37]. The second set, NNPDF3.0QED, is based on the NNPDF3.0 family of PDFs and adds a photon component that is extracted from a fit to collider data. NNPDF3.0QED photon density has both a much larger central value and PDF uncertainty than those of LUXQED. On the other hand, the two sets are compatible within PDF errors and they both include QED effects in the DGLAP evolution on top of the usual NNLO QCD evolution.
We confirm the observations already made in ref. [34], namely, the way the photon PDF is included impacts all differential distributions. The size of this impact is different for the various distributions; the most significant impact can be observed in the p T,avt distribution at moderate and large p T where the net effect from EW corrections based on NNPDF3.0QED is rather small and with large PDF uncertainties, while using LUXQED it is negative, with small PDF uncertainties and comparable to the size of the NNLO QCD scale error. The m(tt) distribution displays even larger effects, but only at extremely high m(tt). The y(tt) distribution is also affected at large y(tt) values. 8 It seems to us that a consensus is emerging around the structure-function approach of ref. [37]. Given its appealing predictiveness, this approach will likely be utilised in the future in other PDF sets. Therefore, at present, it seems to us that as far as the photon PDF is concerned predictions based on the LUXQED set should be preferred.
Our best predictions in this work are based on the so-called multiplicative approach for combining QCD and EW corrections. We have also presented predictions based on the standard additive approach. In general, we find that the difference between the two approaches is small and well within the scale uncertainty band. The difference between the two approaches is more pronounced for the m(tt) and p T,avt distributions. Nevertheless, both approaches agree within the scale variation. The scale uncertainty is smaller within the multiplicative approach and, especially in the case of the p T,avt distribution, does not overlap with the NNLO QCD uncertainty band. We stress that these features may be sensitive to the choice of factorisation and renormalisation scales.
Since we are unaware of a past study of Heavy Boson Radiation (i.e. H, W ± and Z) in inclusive tt production, for completeness, we have also presented the impact of inclusive HBR on the inclusive top-pair differential spectrum. While it is often assumed that additional HBR emissions can be removed in the measurements, it is nevertheless instructive to consider the contribution of such final states. We find that, typically, the HBR contribution is negligible, except for the m(tt) distribution, where it tends to partially offset the EW correction (when computed with LUXQED). We have also checked that NLO QCD corrections to the LO HBR result are negligible for all inclusive tt distributions considered by us.
The LO (m + n = 2), NLO (m + n = 3) and NNLO (m + n = 4) contributions read Σ tt LO (α s , α) = α 2 s Σ 2,0 + α s αΣ 2,1 + α 2 Σ 2,2 ≡ Σ LO,1 + Σ LO,2 + Σ LO,3 , In order to simplify the notation, we further define the following purely QCD quantities and those involving EW corrections Throughout this work with the term "EW corrections" we refer to the quantity Σ EW , while the term "NLO EW corrections" will only refer to Σ NLO EW . In the additive approach, which is presented in section 3.1, QCD and electroweak corrections are combined through the linear combination Σ QCD+EW ≡ Σ QCD + Σ EW . (A.7) The so called "multiplicative approach", which has been discussed in sec. 3.1, is precisely defined in the following. The purpose of the multiplicative approach is to estimate the size of Σ NNLO,2 , which for convenience we rename Σ mixed and assuming complete factorisation of NLO QCD and NLO EW effects we estimate as In the regime where NLO QCD corrections are dominated by soft interactions and NLO EW by Sudakov logarithms, eq. (A.8) is a very good approximation, since the two effects factorise and are dominant. In other regimes Σ mixed can be used as an estimate of the leading missing mixed QCD-EW higher orders. The advantage of the inclusion of Σ mixed is the stabilisation of the scale dependence of the term Σ NLO EW , which in tt production has almost 9 the same functional form of Σ LO QCD . To this end we define the multiplicative 9 We say "almost" because this order receives also QCD corrections to the ΣLO EW contributions from the gγ and bb initial states. Besides these effects ΣNLO EW(µ2) = ΣNLO EW(µ1) In order to test the stability of the multiplicative approach under even higher mixed QCD-EW orders, we combine NNLO QCD corrections and NLO EW corrections in order to estimate, besides the Σ mixed term, also NNNLO contributions of order α 4 s α. Finally, we briefly describe how the dependence on the photon PDF enters the different perturbative orders. At LO and NLO accuracy, all contributions, with the exception of Σ LO QCD and Σ NLO QCD , depend on the photon PDF. The dominant photon-induced process is gγ → tt, which contributes to Σ LO EW and, via QCD corrections to this order, to Σ NLO EW . In addition, Σ NLO EW , but also Σ NLO,3 and Σ NLO,4 , receive contributions from the qγ → ttq andqγ → ttq processes. Moreover, in the case of Σ LO,3 and Σ NLO,4 , also the γγ initial state contributes. As already discussed in ref. [34], almost all of the photoninduced contribution arises form Σ LO EW . In this work, at variance with ref. [34], we also include the term Σ res in our calculations. However, since the size of Σ res is in general small, the previous argument still applies. The numerical impact of Σ res is discussed in sec. 3.2.
Given the structure of the photon-induced contributions described before, it is also important to note that, with LUXQED, the multiplicative approach is a better approximation of Σ mixed than in the case of NNPDF3.0QED. Indeed, the order Σ NLO EW contains also terms that can be seen as "QCD corrections" to the gγ contributions in Σ LO,2 (negligible only with the LUXQED), but are not taken into account in the multiplicative approach.