# The impact of the photon PDF and electroweak corrections on $t^-t$ distributions

D Pagani, I Tsinikos, M Zaro

► To cite this version:

D Pagani, I Tsinikos, M Zaro. The impact of the photon PDF and electroweak corrections on $t^-t$ distributions. European Physical Journal C: Particles and Fields, 2016, 76 (9), pp.479. 10.1140/epjc/s10052-016-4318-z. hal-01360893

HAL Id: hal-01360893

https://hal.sorbonne-universite.fr/hal-01360893

Submitted on 6 Sep 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# The impact of the photon PDF and electroweak corrections on $t\bar{t}$ distributions

D. Pagani<fnref n="1" />, I. Tsinikos<fnref n="1" />, M. Zaro<fnref n="2" />,<fnref n="3" />,a

Received: 16 June 2016 / Accepted: 12 August 2016

© The Author(s) 2016. This article is published with open access at Springerlink.com

**Abstract** We discuss the impact of EW corrections on differential distributions in top-quark pair production at the LHC and future hadron colliders, focussing on the effects of initial-state photons. Performing a calculation at Next-to-Leading Order QCD+EW accuracy, we investigate in detail the impact of photon-initiated channels on central values as well as PDF and scale uncertainties, both at order $\alpha_s\alpha$ and $\alpha_s^2\alpha$. We present predictions at 13 and 100 TeV, and provide results for the 8 TeV differential measurements performed by ATLAS and CMS. A thorough comparison of results obtained with the NNPDF2.3QED and CT14QED PDF sets is performed. While contributions due to the photon PDF are negligible with CT14QED, this is not the case for NNPDF2.3QED, where such contributions are sizeable and show large PDF uncertainties. On the one hand, we show that differential observables in top-pair production, in particular top-quark and $t\bar{t}$ rapidities, can be used to improve the determination of the photon PDF within the NNPDF approach. On the other hand, with current PDF sets, we demonstrate the necessity of including EW corrections and photon-induced contributions for a correct determination of both the central value and the uncertainties of theoretical predictions.

## 1 Introduction

The production of a pair of top quarks ($t\bar{t}$) is one of the key processes in the physics program at the LHC. On the one hand, this process is an ideal testing ground for strong interactions, which can be probed in the perturbative regime, and more in general for the standard-model (SM) interactions. On the other hand, because the experimental signatures emerging from $t\bar{t}$ production can feature leptons, miss-

ing transverse energy and (b)-jets, this process enters as a background in several experimental analyses and searches, notably those related to the Higgs boson and new physics. The study of $t\bar{t}$ production gives also access to many of the top-quark properties and can be exploited for the determination of the parton-distribution functions (PDFs) of the proton [1–3]. Furthermore, since the top quark is the heaviest particle in the SM, it is often a proxy to new physics in beyond-the-standard-model (BSM) scenarios. Thus, $t\bar{t}$ production can be exploited to precisely test SM predictions as well as to search for or put constraints on BSM interactions.

For these reasons, a proper modelling of the $t\bar{t}$ signal and a sound understanding of all effects that can lead to distortions in the distributions are mandatory. This explains the great efforts which have been made in order to improve the accuracy of theoretical predictions: inclusive and differential cross sections have been calculated up to next-to-next-to-leading-order (NNLO) accuracy in QCD [4,5], and including soft gluon effects up to next-to-next-to-leading-logarithmic (NN-LL) accuracy [2,6–16]. Notably, higher-order QCD corrections lead to the reduction of scale uncertainties down to very few percents. At this level of accuracy effects due to electroweak (EW) corrections cannot be neglected, especially for differential distributions, where the expected hierarchy between the strong and electroweak interactions ($\alpha_s \gg \alpha$) may not be respected. On the one hand, when large scales are probed ($Q \gg m_W$), the so-called Sudakov logarithms render EW corrections large and negative. On the other hand, EW corrections receive also contributions from channels that feature photons in the initial state, thus they depend on the PDF of the photon. As we will show in this work, the impact of photons in the initial state strongly depends on the PDF set employed and can give a positive contribution of non-negligible size, leading to possible compensations of the corrections induced by Sudakov logarithms.

<footnote n="1">¹Centre for Cosmology, Particle Physics and Phenomenology (CP3), Université Catholique de Louvain, Chemin du Cyclotron 2, 1348 Louvain-la-Neuve, Belgium</footnote>

<footnote n="2">²Sorbonne Universités, UPMC University Paris 06, UMR 589, LPTHE, 75005 Paris, France</footnote>

<footnote n="3">³CNRS, UMR 7589, LPTHE, 75005 Paris, France</footnote>

<footnote>a e-mail: marco.zaro@cern.ch</footnote>
While weak [17–25], QED [26] and EW (weak+QED) [27–30] corrections to $t\bar{t}$ production at hadron colliders have been known for quite some time and the effect of Sudakov logarithms has been quantified and studied, the impact of photon-induced subprocesses has been addressed only in [26] for only the $g\gamma \to t\bar{t}$ channel at LO and using the only PDF set with the photon distribution available at that time, i.e., MRST2004QED [31].

In this work we (re-)evaluate, besides the NLO QCD corrections, the complete set of NLO EW corrections to top-quark pair production and carefully assess their impact on differential distributions for proton-proton collisions at 8, 13 and 100 TeV centre-of-mass energy. For the first time, photon-induced contributions are taken into account both at $O(\alpha_s\alpha)$, from the $g\gamma$-channel at tree level, and at $O(\alpha_s^2\alpha)$, from the $g\gamma$ as well as the $q\gamma(\bar{q}\gamma)$ initial states arising in the NLO EW corrections. For all the results presented here, the calculation has been performed within the MADGRAPH5_AMC@NLO framework [32], thanks to an extension of the code that allows one to automatically calculate NLO QCD and EW corrections [33,34]. In order to have a reliable estimate of the photon-induced contribution and of its uncertainties, we evaluate the impact of electroweak corrections with both the NNPDF2.3QED [35] and the CT14QED [36] PDF sets, the only modern sets that include LO QED contribution in the DGLAP evolution and, with very different assumptions and strategies in the extraction from data, the photon distribution.<fnref n="1" />

We show that, besides higher-order QCD corrections, it is also necessary to take into account EW corrections, including photon-induced contributions, for a correct determination of both the central value and the uncertainty of theoretical predictions. This is particularly important in the context of the ongoing discussion of NLO-accurate event generators and the compatibility with experimental data for the $p_T(t)$ distribution at the LHC 8-TeV measurements [39] and in view of the measurements at 13 TeV. Furthermore, we show that differential observables in top-pair production, in particular top-quark and $t\bar{t}$ rapidities, can be used to improve the determination of the photon PDF within the NNPDF approach, while in the CT14QED approach $t\bar{t}$ production is not sensitive to the photon-induced contributions.

The structure of our paper is as follows: in Sect. 2 we present the framework employed to perform our calculation and discuss relevant input parameters. In Sect. 3 we discuss differences among existing PDF sets which provide a photon distribution (NNPDF2.3QED, CT14QED and the older MRST2004QED set). We describe the different theoretical approaches employed in the sets, and we compare central values and uncertainties for the photon PDF and the parton

luminosities relevant for our calculation. In Sect. 4 we show predictions at 13 and 100 TeV, and we compare results with and without the contribution of photon-induced processes and using the NNPDF2.3QED or the CT14QED PDF set. In Sect. 5 we analyse the impact of EW corrections and the photon PDF for specific measurements performed by ATLAS and CMS at 8 TeV. We compare the usage of NNPDF2.3QED with a standard modern set with QCD-only partons and DGLAP evolution. We give our conclusions and an outlook in Sect. 6.

## 2 Calculation setup and input parameters

The calculation has been performed in a completely automated way and we do not describe here the technical details; they can be found in [33,34], where the $t\bar{t}H$ process has been calculated in the same framework. In the following, we only want to match the notation of this paper to the one introduced in [33] and precisely define the quantities included in our calculation, specifying those that depend on the photon PDF.

In the case of $pp \to t\bar{t}$ process a generic observable $\Sigma^{t\bar{t}}$ can be expanded simultaneously in the QCD and EW coupling constants as

$$ \Sigma^{t\bar{t}}(\alpha_s, \alpha) = \sum_{m+n \ge 2} \alpha_s^m \alpha^n \Sigma_{m+n,n}. \qquad (2.1) $$

Following this notation, one can separate the LO ($m+n=2$) and NLO ($m+n=3$) contributions as follows:

$$ \begin{aligned} \Sigma_{\text{LO}}^{t\bar{t}}(\alpha_s, \alpha) &= \alpha_s^2 \Sigma_{2,0} + \alpha_s \alpha \Sigma_{2,1} + \alpha^2 \Sigma_{2,2} \\ &\equiv \Sigma_{\text{LO},1} + \Sigma_{\text{LO},2} + \Sigma_{\text{LO},3}, \\ \Sigma_{\text{NLO}}^{t\bar{t}}(\alpha_s, \alpha) &= \alpha_s^3 \Sigma_{3,0} + \alpha_s^2 \alpha \Sigma_{3,1} + \alpha_s \alpha^2 \Sigma_{3,2} + \alpha^3 \Sigma_{3,3} \\ &\equiv \Sigma_{\text{NLO},1} + \Sigma_{\text{NLO},2} + \Sigma_{\text{NLO},3} + \Sigma_{\text{NLO},4}. \end{aligned} \qquad (2.2) $$

In our results we include the $\Sigma_{\text{LO},1}$, $\Sigma_{\text{LO},2}$, $\Sigma_{\text{NLO},1}$ and $\Sigma_{\text{NLO},2}$ terms. We checked that the remaining terms are subleading as expected, giving results of the order or below 1% of the LO,1 contribution both in the total cross section and at the differential level. In order to help the reader we further define the quantities

$$ \Sigma_{\text{LO QCD}} \equiv \Sigma_{\text{LO},1}, \qquad \Sigma_{\text{NLO QCD}} \equiv \Sigma_{\text{NLO},1}, \qquad (2.3) $$

$$ \Sigma_{\text{LO EW}} \equiv \Sigma_{\text{LO},2}, \qquad \Sigma_{\text{NLO EW}} \equiv \Sigma_{\text{NLO},2}, \qquad (2.4) $$

$$ \Sigma_{\text{QCD}} \equiv \Sigma_{\text{LO QCD}} + \Sigma_{\text{NLO QCD}}, \qquad (2.5) $$

$$ \Sigma_{\text{EW}} \equiv \Sigma_{\text{LO EW}} + \Sigma_{\text{NLO EW}}, \qquad (2.6) $$

$$ \Sigma_{\text{QCD+EW}} \equiv \Sigma_{\text{QCD}} + \Sigma_{\text{EW}}. \qquad (2.7) $$

In the following text by the term “EW corrections” we will in general refer to the quantity $\Sigma_{\text{EW}}$, while we will

<footnote n="1">¹ We acknowledge that two new studies of the determination of the photon PDF [37,38] have appeared during the publication of our work.</footnote>
use “NLO EW corrections” for $\Sigma_{\text{NLO EW}}$. At variance with Refs. [33,34], in our calculation we do not consider the effect due to the heavy-Boson-radiation (HBR). Although the LO cross sections of $pp \to t\bar{t}V$, $V = W^{\pm}, Z, H$ processes may in principle contribute at the same perturbative order of $\Sigma_{\text{NLO EW}}$ to the inclusive $t\bar{t}$ production, in this work we assume that $t\bar{t}$ and $t\bar{t}V$ final states can always be distinguished.

The quantities that depend on the photon PDF are $\Sigma_{\text{LO EW}}$ and $\Sigma_{\text{NLO EW}}$. The $g\gamma \to t\bar{t}$ process contributes to $\Sigma_{\text{LO EW}}$ and to $\Sigma_{\text{NLO EW}}$. In addition, $\Sigma_{\text{NLO EW}}$ receives contributions from the $q\gamma \to t\bar{t}q$ and $\bar{q}\gamma \to t\bar{t}\bar{q}$ processes at the tree level, which feature initial-state singularities that have to be subtracted. As in the case of $q\bar{q} \to t\bar{t}\gamma$ processes, the subtracted QED singularities are taken into account in the DGLAP evolution in MRST2004QED, NNPDF2.3QED and CT14QED. Thus, for theoretical consistency, these three PDF sets should be preferred whenever NLO EW corrections are computed. However, since the QED accuracy of DGLAP evolution is only LO in all the three PDF sets, one could not technically claim at the moment NLO QED and thus NLO EW accuracy for hadronic predictions. Even before considering the numerical results in Sects. 4 and 5, this issue points to the necessity of a better determination of the photon PDF by means of fits which include NLO QED effects in the DGLAP evolution [40,41].

Our calculation is performed using the following input parameters:

$$ \begin{aligned} m_t &= 173.3 \text{ GeV}, & m_H &= 125.09 \text{ GeV}, \\ m_W &= 80.385 \text{ GeV}, & m_Z &= 91.1876 \text{ GeV}, \end{aligned} \qquad (2.8) $$

and setting all the other fermion masses to zero. All masses are renormalised on-shell and all decay widths are set to zero. The renormalisation of $\alpha_s$ is performed in the five-flavour scheme, while EW parameters are chosen in the $G_\mu$-scheme, with

$$ G_{\mu} = 1.1663787 \cdot 10^{-5} \text{ GeV}^{-2}. \qquad (2.9) $$

Since NLO EW corrections of $O(\alpha_s^2\alpha)$ to $t\bar{t}$ hadroproduction do not involve the renormalisation of $\alpha$, the choice of a different EW scheme will not change our results in a visible way. The CKM matrix is taken as the identity.

Unless differently specified, we use a dynamical reference scale for the central values of the renormalisation ($\mu_r$) and factorisation ($\mu_f$) scales defined as

$$ \mu = \frac{H_T}{2} = \frac{1}{2} \sum_i m_{T,i}, \qquad (2.10) $$

where the sum of the transverse masses runs over all the final-state particles. In all cases theoretical uncertainties due

to missing higher orders are estimated via independent variations of $\mu_r$ and $\mu_f$ in the interval $[\mu/2, 2\mu]$.

It is worth to note that the NNPDF2.3QED set is in the variable-flavour scheme with six active flavours, which for $\mu > m_t$ is equivalent to the six-flavour scheme. On the contrary, in CT14QED the active flavours are five, leading to the five-flavour scheme also for $\mu > m_t$. As we said, we renormalise $\alpha_s$ in the five-flavour scheme for all the predictions; for the comparison between the NNPDF2.3QED and CT14QED results we simply change the PDF set without modifying the calculation framework. The change of scheme can easily be performed by following the recipe described in [42] and based on [43], which, at NLO QCD + NLO EW accuracy, has a direct effect only on the $q\bar{q}$-induced contribution to $\Sigma_{\text{NLO QCD}}$. We explicitly verified that the numerical impact of such a change of renormalisation scheme is always much smaller of the scale uncertainty and furthermore cancels in any ratio involving $\Sigma_{\text{NLO QCD}}$ both at the numerator and the denominator. Thus, it has not any influence on the discussion presented in this work.

### 3 Photon PDF and parton luminosities

In this section we discuss in some detail the different modelling of the photon PDF in the NNPDF2.3QED, CT14QED<fnref n="2" /> and MRST2004QED sets.

Although for all the three PDF sets the DGLAP evolution is performed at NLO QCD + LO QED accuracy,<fnref n="3" /> very different and crucial assumptions underlie the determination of the photon PDF $\gamma(x, Q)$. These differences mainly concern the following three aspects. First, the ansatz for $\gamma(x, Q)$ at the initial scale $Q = Q_0$. Second, the different datasets which are used in the fit. Third, the practical implementation of the DGLAP evolution from the initial scale $Q_0$ to the scale $Q$.

In Fig. 1 we show the photon PDF for the different sets at the scales $Q = (3, 173, 5000)$ GeV. As can be seen, these assumptions affect the dependence on $x$ and $Q$ for both the central value and the uncertainty band.<fnref n="4" /> The main reasons

<footnote n="2">² CT14QED provides two kinds of sets, one with only the incoherent component of the photon PDF and another one with both the coherent and the incoherent components. In our work we have used the first set. We have checked that the inclusion of the coherent component in the photon PDF does not significantly alter our findings. The predictions obtained with the photon PDF with momentum fraction $p_0^\gamma = 0.00\%$ including both components are very similar, in the $x$ and $Q$ ranges relevant for our study, to those from the incoherent-only photon PDF with $p_0^\gamma = 0.14\%$.</footnote>

<footnote n="3">³ In the case of NNPDF2.3QED, PDFs at NNLO QCD + LO QED accuracy as well as at LO QCD + LO QED accuracy are also available. However, we considered here only the NLO QCD + LO QED case, consistently with the other PDF sets discussed.</footnote>

<footnote n="4">⁴ The red dotted curve labelled as APFEL_NN23 will be explained after in the text.</footnote>
<figure><img src="image_3.png" /><figcaption>Fig. 1 Comparison of the photon PDF for MRST2004QED, CT14QED, NNPDF2.3QED and APFEL_NN23 sets at $Q = 3$ GeV, $Q = 173$ GeV and $Q = 5000$ GeV</figcaption></figure>

for the differences can be traced to the different assumption for the photon PDF at the initial scale and in differences in the QCD+QED evolution. In particular:

- Consistently with the approach pursued for coloured partons, in NNPDF2.3QED no functional form is specified for the photon PDF at the initial scale, $\gamma(x, Q_0)$. The photon PDF is only constrained to be positive. In a first step, PDF replicas for all the partons are fit together from DIS-data only. Afterwards, they are further constrained by Drell-Yan data from the LHC Run-I at 7 TeV. At variance with DIS, neutral-current Drell-Yan is sensitive to the photon PDF already at LO, and it can put stronger constraints on $\gamma(x, Q_0)$. Because of the positivity requirement for the photon PDF, the replicas distribute in a very non-Gaussian way around the central value. The prescription suggested in order to determine a 68% CL uncertainty band consists in the evaluation of the symmetric error including 68 of the 100 replicas around the central value. Since no model is assumed for the photon PDF and no data are present for large $x$, in this region uncertainties are very large and the central value alone can be misleading, leading even to an unphysical peak at very large $x$, which can be seen in Fig. 1.

- The CT14QED and MRST2004QED sets are based on a completely different assumption for $\gamma(x, Q)$ at the initial scale $Q = Q_0$. Since the amount of data for constraining the shape of the photon PDF is limited, $\gamma(x, Q_0)$ is chosen to be described by an ansatz; the photon PDF at $Q = Q_0$ is assumed to be completely determined by the valence quark distributions. Specifically, in CT14QED and MRST2004QED the photon parameterisation at the initial scale $Q_0$ GeV reads

$$f_{\gamma/p}(x, Q_0) = \frac{\alpha}{2\pi} \left( A_u e_u^2 \tilde{P}_{\gamma q} \otimes u^0(x) + A_d e_d^2 \tilde{P}_{\gamma q} \otimes d^0(x) \right). \quad (3.1)$$

In Eq. (3.1) $\tilde{P}_{\gamma q} \otimes f^0(x)$ corresponds to the convolution of the splitting function $\tilde{P}_{\gamma q}(x)$ with the so-called “primordial” quark distributions $f^0(x)$. In the case of MRST2004QED, $f^0(x)$ are valence-like model distributions, i.e., they are not those fitted within the global MRST2004 [44] analysis. In the case of CT14QED, $f^0(x)$ correspond to the initial up and down valence distributions from CT14NLO [45].

More importantly, in CT14QED $A_u$ is set equal to $A_d$ in order to obtain a dependence on a single parameter. Conversely, in MRST2004QED, the coefficients $A_i$ are given by $A_i = \ln(Q_0^2/m_i^2)$, where the $m_i$ are the “Current Masses” (CMs) of the quarks ($m_u = 6$ MeV and $m_d = 10$ MeV). Thus, the case of CT14QED is a generalisation of the original ansatz used in MRST2004QED, where the photon PDF is simply determined by the CMs for the up and the down quark. The CT14QED set leaves freedom for the value of $m_i$ in the $A_i$ coefficients, which can be alternatively parameterised by the momentum fraction $p_\gamma^0$ carried by the photon at the initial scale. The constraints on $p_\gamma^0$ and consequently on the photon PDF are then obtained by fitting the ZEUS data [46] for the production of isolated photons in deep-inelastic-scattering (DIS), $ep \rightarrow e\gamma + X$, a process which is sensitive to the photon PDF already at LO. In this way, the constraint $p_0^\gamma \lesssim 0.14\%$ is set at 90% CL for the photon PDF, at the initial scale of $Q_0 = 1.295$ GeV. For this reason, the uncertainty band in CT14QED corresponds to the area between the $p_0^\gamma = 0.14\%$ and $p_0^\gamma = 0\%$ predictions; no central, or preferred value is given in this approach. Since in MRST2004QED the photon PDF is simply determined from the quark PDF assuming CMs for the up and the down quark, no uncertainty band is provided. Moreover, for this set, no constraints from data enter in the determination of the photon PDF.

- All three PDF sets implement the DGLAP evolution at NLO QCD+LO QED accuracy. However, while the scale is evolved simultaneously with the same value for the
<figure><img src="image_4.png" /><figcaption>Fig. 2 Comparison of the $g\gamma$ luminosity for MRST2004QED, CT14QED, NNPDF2.3QED and APFEL_NN23 sets at 8, 13 and 100 TeV</figcaption></figure>

QCD and QED evolution in CT14QED and MRST2004-QED, in NNPDF2.3QED the two scales are run independently. The effect from the different evolution can be seen by comparing the plots in Fig. 1, which display the same quantities at different factorisation scales. At $Q = 3$ GeV the central values of the three sets are different, although the CT14QED and NNPDF2.3QED photons are well compatible within the uncertainty. At large $Q$, while CT14QED and MRST2004QED converge to similar values, the prediction from NNPDF2.3QED is different. This trend is particularly visible at small $x$ and for large scales, where the photon PDF mostly stems from the quark densities via the DGLAP evolution. However, it is worth to recall that a new correct implementation of the coupled QCD+QED DGLAP equations is now also available in APFEL [47] and has been described in [48] in the context of the determination of new sets of PDFs including also lepton members.<fnref n="5" /> In Fig. 1 we included also a red dashed line labelled APFEL_NN23, which corresponds to the central value for the set denoted as “C2” in [48]. This set is the same of NNPDF2.3QED at the initial scale $Q_0 = 1$ GeV, but it is different at larger scales due to the different treatment of the evolution. Indeed, the photon PDF in APFEL_NN23 is very close to the one of NNPDF2.3QED only at large $x$, while for small $x$ and large $Q$ it converges to the values of CT14QED and MRST-2004QED. This proves that the discrepancy between the

different sets at small $x$ and large $Q$ is completely driven by the different treatment of the DGLAP evolution.

In Fig. 2 we show the photon-gluon luminosity as a function of the invariant mass $M$ for the three PDF sets considered at 8, 13 and 100 TeV, with the factorisation scale $Q$ set equal to $M$. In the first inset we show for every luminosity the ratio of their values at $Q = M$ and $Q = m_t$, while in the second(third) inset we show at $Q = M(Q = m_t)$ the ratio of every luminosity over the corresponding central value with NNPDF2.3QED. For the entire range in $M$ considered at the three different proton-proton-collision energies, the CT14QED and NNPDF2.3QED photon-gluon luminosities are barely compatible or not compatible at all. The slightly better compatibility with $Q = m_t$ is mainly due to the larger uncertainty band in NNPDF2.3QED.

It is important to note that the behaviour of the photon PDF at small $x$ and large $Q$, the region where differences among the PDF sets are the largest (see Fig. 1) and are determined by the different QCD+QED DGLAP evolution, is not of particular interest for our study. Indeed, as can be seen in Fig. 2, the APFEL_NN23 $g\gamma$ luminosity is very close to the NNPDF2.3QED one and inside its uncertainty band. The insensitivity of our study to the different treatment of QCD+QED DGLAP running is in principle expected for most of the predictions at the LHC, since large factorisation scales $Q$ are typically not correlated with small values of $x$. Moreover, at small values of $M$, where small values of $x$ can be probed, the PDFs of quarks and gluons are much larger than the photon PDF, leading to a relative suppression of photon-initiated contributions.

<footnote n="5">⁵ During the writing of this article also a new set NNPDF3.0QED has been released [49] and it includes the correct DGLAP running at (N)NLO QCD + LO QED accuracy.</footnote>
<figure><img src="image_5.png" /><figcaption>Fig. 3 Comparison of the $g\gamma$, $gg$ and $q\bar{q}$ luminosities for the NNPDF2.3QED set at 8, 13 and 100 TeV</figcaption></figure>

In order to determine at which values of $M$ the photon-gluon luminosity is actually expected to be relevant in $t\bar{t}$ hadroproduction, in Fig. 3 we show the photon-gluon, gluon-gluon and quark-antiquark luminosities for NNPDF2.3QED sets at 8, 13 and 100 TeV. These plots will also become useful to understand the dependence of $\Sigma_{EW}/\Sigma_{LO}$ QCD on the definition of the scale.<fnref n="6" />

The solid lines refer to the luminosities with $Q = M$, with the uncertainty band shown only for the photon-gluon case. The same-colour dashed lines are the corresponding quantities with $Q = m_t$. The ratios between predictions with $Q = M$ and with $Q = m_t$ is shown in the first inset. The luminosities at $Q = M$ are in general of the same order or smaller than at $Q = m_t$, with the exception of the 100 TeV case for $M \lesssim 10$ TeV. In particular at 8 and 13 TeV, in the range of $M \lesssim 1$ TeV the luminosities of $gg$ and $q\bar{q}$ strongly reduce at $Q = M$ with respect to $Q = m_t$, while the photon-gluon luminosity is less sensitive to the value of $Q$. Since in the TeV range the negative contributions from Sudakov logarithms in the $gg$ and $q\bar{q}$ channels and the positive photon-gluon contribution tend to cancel each other, larger scales lead to larger values of $\Sigma_{EW}/\Sigma_{LO}$ QCD. This effect will be discussed both in Sect. 4 for predictions at 13 TeV and in Sect. 5, where the $\Sigma_{QCD+EW}/\Sigma_{QCD}$ ratio will be compared to the experimental accuracy reached at the LHC measurements at 8 TeV. In the second (third) insets we also show the ratios of the $gg$ and $q\bar{q}$ luminosities over the $g\gamma$ luminosity at $Q = M(Q = m_t)$ in order to better identify at which scales the $g\gamma$ contribution

is expected to be relevant. From these insets, the hierarchy between $gg$ and $q\bar{q}$ channels can also be easily derived.

In conclusion, the different results for photon-induced processes and consequently for EW corrections obtained with the different PDF sets discussed in this article have to be attributed mainly to the different assumptions that underlie the determination of $\gamma(x, Q)$ at the initial scale $Q = Q_0$ and the choice of the factorisation scale, but not to the different treatment of the DGLAP NLO QCD+LO QED evolution.

## 4 Predictions at 13 and 100 TeV

In this section we discuss the impact of the EW corrections and the photon PDF on several distributions at 13 and 100 TeV. In particular, we focus on the top-pair invariant mass $m(t\bar{t})$, the transverse momentum of the top quark $p_T(t)$, and the rapidities of the top quark $y(t)$ and $t\bar{t}$ pair $y(t\bar{t})$. Predictions for the LHC at 13 TeV are shown in Figs. 4, 5, 6, 7, 8, 9, while those for a future circular collider (FCC) at 100 TeV are shown in Figs. 10, 11, 12, 13, 14, 15. Unless stated otherwise, results are obtained with no cut imposed on the final-state particles.

In each figure we show two plots for the same observable, displaying in the left plot, denoted as “no $\gamma$”, predictions with the photon PDF artificially set equal to zero. The reason behind this choice is manifold. First of all, the comparison of left and right plots allows one to gauge the impact of the photon PDF on both the central value and the PDF uncertainties of the electroweak contributions. Furthermore, in the plots on the left it is possible to observe the impact of the Sudakov logarithms, which can be hidden by the compensation given

<footnote n="6">⁶ In our calculation there are also contribution from quark-photon initial states, however, as will be discussed in Sect. 4, their contribution is small so we decided to exclude the corresponding luminosities from the discussion in this section for the sake of clarity.</footnote>
<figure><img src="image_6.png" /><figcaption>Fig. 4 Differential distributions for the $p_T(t)$ at 13 TeV. The format of the plots is described in detail in the text</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Fig. 5 Differential distributions for the $m(t\bar{t})$ at 13 TeV. The format of the plots is described in detail in the text</figcaption></figure>

by the photon-induced processes included in the plots on the right. Last but not least, whereas the EW corrections in the right plots strongly depend on the PDF set used and may possibly change with future improved determination of the photon PDF, the plots on the left display the subset of the EW corrections that are expected to be stable under future PDF determination.

In each plot the main panel includes the distributions for $\Sigma_{LO} \text{QCD}$, $\Sigma_{QCD}$ and $\Sigma_{QCD+EW}$ predictions as defined in equations (2.3)–(2.7). The four insets below the main panel

display ratios of several quantities over the $\Sigma_{LO} \text{QCD}$ central value. In the first inset we plot this ratio for the central values of $\Sigma_{LO} \text{EW}$, $\Sigma_{NLO} \text{QCD}$ and $\Sigma_{NLO} \text{EW}$, i.e., the other perturbative orders which we consider. The second and third insets, respectively, present the $\Sigma_{QCD+EW}/\Sigma_{LO} \text{QCD}$ ratio including the scale and PDF uncertainties for the numerator. In the last inset we present the $\Sigma_{EW}/\Sigma_{LO} \text{QCD}$ ratio with the PDF uncertainties for the numerator. All the results in the plots are obtained with the NNPDF2.3QED PDF set, with the PDF uncertainty computed at the 68% CL. However, plots on
<figure><img src="image_8.png" /><figcaption>**Fig. 6** Differential distributions for the $m(t\bar{t})$ at 13 TeV with $|y(t)|$, $|y(\bar{t})|$ < 2.5 cuts applied. The format of the plots is described in detail in the text</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 7** Differential distributions for the $y(t)$ at 13 TeV. The format of the plots is described in detail in the text</figcaption></figure>

the right, including effects due to the photon PDF, also show
(in the last inset) the ratio $\Sigma_{EW} / \Sigma_{LO\,QCD}$ computed with the
CT14QED PDF set.<fnref n="7" /> For this particular set two predictions are
shown, which correspond to a momentum fraction carried by
the photon at the initial scale of 0 and 0.14%; as explained

in Sect. 3 this interval corresponds to the uncertainty band of
CT14QED.

The first general feature that can be noticed in all the dis-
tributions for the LHC at 13 TeV is the different behaviours
of the EW corrections with NNPDF2.3QED and CT14QED in
the plots on the right. The effects induced by the photon PDF
with CT14QED is much smaller w.r.t. the case of NNPDF2.3-
QED. Indeed, by comparing the last inset in equivalent plots

<footnote n="7">⁷ Only for this ratio, also the $\Sigma_{LO\,QCD}$ quantity has been evaluated with the CT14QED PDF set.</footnote>
<figure><img src="image_10.png" /><figcaption>**Fig. 8** Differential distributions for the $y(\bar{t}\bar{t})$ at 13 TeV. The format of the plots is described in detail in the text</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 9** Differential distributions for the $y(t)$ at 13 TeV with the cut $m(t\bar{t}) > 1$ TeV. The format of the plots is described in detail in the text</figcaption></figure>

with (right) and without (left) photon-induced contributions, it is easy to see that effects due to the CT14QED photon PDF are essentially invisible, regardless of the momentum fraction (in the 0–0.14% range). Furthermore, by comparing the first inset in the two plots of each figure, it is possible to notice how the dominant contribution from photon-induced processes originates from the $\Sigma_{LO\,EW}$ contribution of the $g\gamma$ channel. As anticipated in Sect. 3, the $q\gamma$ contribution is

negligible<fnref n="8" /> and will not be further discussed in the following. Moreover, when we will refer to the effect of photon PDF, unless differently specified, we will always understand the case of NNPDF2.3QED PDF set in order to be conservative.

<footnote n="8">⁸ This statement could be done only *a posteriori*, after having performed the calculation and cannot be generalised to other processes. See for instance Ref. [50].</footnote>
<figure><img src="image_12.png" /><figcaption>Fig. 10 Differential distributions for the $p_T(t)$ at 100 TeV. The format of the plots is described in detail in the text</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>Fig. 11 Differential distributions for the $m(t\bar{t})$ at 100 TeV. The format of the plots is described in detail in the text</figcaption></figure>

As can be seen in Figs. 4 and 5, at 13 TeV the photon-initiated contributions can compensate the Sudakov suppression in the top-pair invariant mass and top-quark transverse momentum distributions. At very large scales ($p_T(t) > 1.7$ TeV or $m(t\bar{t}) > 3$ TeV) their size is larger than the effect due to Sudakov logarithms. In all the figures at 13 TeV, in the last inset of the plots on the right, we also included the central value of the ratio $\Sigma_{EW}/\Sigma_{LO\,QCD}$ with both quantities evaluated at the scale $\mu = m_t$. This quantity indicates

that the cancellation between Sudakov logarithms strongly depends also on the scale choice, thus the scale dependence of $\Sigma_{EW}$ should also be taken into account as a source of uncertainties. The origin of this dependence will be discussed in detail in Sect. 5, in the context of the measurements at 8 TeV.

In the last inset, the PDF uncertainties of the $\Sigma_{EW}$ term alone are directly compared with the $\Sigma_{LO\,QCD}$ PDF uncertainties, the grey band. In the tail of the $p_T(t)$ and $m(t\bar{t})$
<figure><img src="image_14.png" /><figcaption>**Fig. 12** Differential distributions for the $y(t)$ at 100 TeV. The format of the plots is described in detail in the text</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Fig. 13** Differential distributions for the $y(t\bar{t})$ at 100 TeV. The format of the plots is described in detail in the text</figcaption></figure>

distributions on the right, PDF uncertainties for separately $\Sigma_{EW}$ and $\Sigma_{LO\,QCD}$ are of the same size. This is consistent with the difference between the PDF uncertainty for the total prediction $\Sigma_{QCD+EW}$ with and without effects from photon in the initial state, which is shown in the third inset and, in the plots on the right, is much larger w.r.t the case of $\Sigma_{LO\,QCD}$. Just for reference, we decided to show in the second inset the ratio $\Sigma_{QCD+EW}/\Sigma_{LO\,QCD}$ with scale uncertainties, which are reduced w.r.t. the $\Sigma_{LO\,QCD}$ case, the grey band. The effects

due to the EW corrections and the impact of the photon PDF is in general smaller than the scale uncertainty of $\Sigma_{QCD+EW}$. However, one should consider that the scale uncertainty is further reduced to very few percents by NNLO corrections, which have already been calculated [4,5]. Moreover, when normalised distributions ($d\sigma/\sigma$) are considered, even at NLO accuracy scale uncertainties are strongly reduced. In Sect. 5 we will comment both effects for specific distributions measured at the LHC at 8 TeV by ATLAS and CMS, where we
<figure><img src="image_16.png" /><figcaption>Fig. 14 Integrated distributions for the $p_T(t)$ at 100 TeV. The format of the plots is described in detail in the text</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>Fig. 15 Integrated distributions for the $m(t\bar{t})$ at 100 TeV. The format of the plots is described in detail in the text</figcaption></figure>

will also investigate the impact of EW corrections and of the photon PDF.

While $p_T(t)$ distributions do not change by requiring $|y(t)|, |y(\bar{t})| < 2.5$, i.e. top and anti-top produced in the central region of the detector, these cuts have an effect on predictions for high $m(t\bar{t})$. This is shown in Fig. 6, which is equivalent to Fig. 5, but with $|y(t)|, |y(\bar{t})| < 2.5$ cuts. At large invariant masses, tops are preferably produced in the forward or backward region, due to the $t$- and $u$-channel diagrams in the $gg$-channel, which is much less suppressed than

the $q\bar{q}$ channels featuring only $s$-channel diagrams at LO QCD. Rapidity cuts suppress the $gg$ contribution, as well as $g\gamma$ contributions, but also enhance the typical value of the partonic Mandelstam variables $t$ and $\bar{t}$. Consequently, with those cuts, we observe larger values for the Sudakov logarithms (plot on the left) and a similar behaviour for the photon-induced contributions. Moreover, we can notice that, as in the $p_T(t)$ distributions in Fig. 4, also in the tails of the plots in Fig. 6 the ratio $\Sigma_{NLO\,QCD}/\Sigma_{LO\,QCD}$ in the first insets decrease, at variance with Fig. 5, where cuts are not applied.
This trend is correlated with the fraction of the cross sec-
tion originating from the gg initial state, which is decreasing
in the tail of p<sub>T</sub>(t) distribution and, only with rapidity cuts
applied, of the m(t<span style="text-decoration: overline">t</span>) distribution.

In Figs. 7 and 8 we show predictions for y(t) and y(t̃) distributions. EW corrections are quite flat in the plots on the left, while in the plots on the right they are enhanced by the gγ initial state in the peripheral region, especially for the y(t̃) distribution. In this region photon-initiated contributions and their uncertainties become relevant since the photon PDF is sampled at rather large x, where the suppression w.r.t. the quark and gluon PDFs is reduced. At variance with the tail of p<sub>T</sub>(t) or m(t̃t), where the gγ contribution is also sizeable, the cross section in the peripheral region is not largely suppressed w.r.t. to the total cross section. For these reasons, large-rapidity bins can be exploited to set additional constraints on the photon PDF in a theoretical framework such as the one adopted in NNPDF2.3QED, while with an assumption à la CT14QED the sensitivity on the photon PDF remains very small. Given the large cross section in t̃t production at the LHC 13 TeV, in Fig. 9 we also show the same plots of Fig. 7 with a m(t̃t) > 1 TeV cut applied. EW corrections in the left plot are larger in the central region, where, due to large values of the ̃t and û Mandelstam variables, the effect of Sudakov logarithms is enhanced. Again, this effect can be compensated by photon-induced processes as shown in the plot on the right. As already observed for all the previous distributions, only in the plots on the left the PDF uncertainties related to Σ<sub>EW</sub> are negligible, while they are sizeable in the plots on the right. Also in this case the trend of the ratio Σ<sub>NLO QCD</sub>/Σ<sub>LO QCD</sub> displayed in the first inset is correlated with the fraction of the cross section originating from the q̃q initial state, as discussed for Figs. 4 and 6.

It is worth to notice that while plots in Fig. 4 would be
identical in the case of p<sub>T</sub>(t̃) distributions, the transverse
momentum of the antitop, the y(t̃) distribution is different
from y(t) when Σ<sub>QCD+EW</sub> terms are included, because of
the charge asymmetry [28–30,51]. However, we observed
the same qualitative behaviour for the photon PDF and the
EW corrections in y(t) and y(t̃) distributions. Moreover, the
gγ channel cannot give a contribution to the numerator of the
charge asymmetry (see also Ref. [52]), and all the other con-
tributions to the asymmetry have already been investigated
in Ref. [30].

The distributions that have been discussed for the LHC at
13 TeV are also presented for a FCC at 100 TeV, with larger
ranges for the abscissae.<fnref n="9" /> By comparing plots in Figs. 10,
11, 12, 13 with their corresponding ones at 13 TeV, it can be
noticed that the impact of the photon PDF is strongly reduced
at 100 TeV. In each figure, the plot on the right (with photons
in the initial state) does not exhibit any qualitatively different

behaviour w.r.t. the plot on the left. The smaller impact of the
photon-induced contributions at 100 TeV w.r.t. the 13 TeV
case is due to the different range of x spanned in the PDFs;
keeping the hardness of the process fixed, a larger energy
of the hadronic collisions corresponds to probing smaller
values of x, where parton luminosities involving photons are
suppressed with respect to those involving QCD partons, as
shown in Fig. 3. For the same reason, the impact of the photon
PDF at the LHC at 8 TeV is even larger than at 13 TeV, as will
be discussed in the next session. Moreover, at 100 TeV, for
a given value of p<sub>T</sub>(t) or m(t̃t), EW corrections are slightly
smaller than at 13 TeV also in the plots on the left, i.e., without
considering photon-induced contributions.

In order to find a large contribution form the gγ initial state at a FCC, it is necessary to probe very large scales. In Figs. 14 and 15 we respectively plot, using the same layout of the previous plots, the cumulative distributions σ(p<sub>T</sub>(t) > p<sub>T,cut</sub>) and σ(m(t̃t) > m<sub>cut</sub>) with p<sub>T,cut</sub> up to 12 TeV and m<sub>cut</sub> up to 30 TeV. First of all, we recall that at a FCC it would be possible with 10 ab⁻¹ integrated luminosity to probe such small cross sections and large hard-scales. Thus, the following discussion is not purely academic. Sudakov logarithms can induce up to ~ -40 % corrections for σ(p<sub>T</sub>(t) > 12 TeV) and ~ -30 % for σ(m(t̃t) > 30 TeV), rendering EW corrections not only important for precision physics, but essential for obtaining a sensible result. For precision physics, Sudakov logarithms have to be considered also beyond the NLO EW accuracy and possibly resummed. The central role of EW corrections at a FCC at 100 TeV is not peculiar for the t̃t final state; it is a general feature for all the production processes (see also the EW section in [53]). At such large scales, also the gγ luminosity is not negligible, and indeed the size and the PDF uncertainty of the EW corrections is very different in the left and right plots. Once again, the prediction obtained with CT14QED is similar to the case where the photon PDF has been set to zero. As in the case of 13 TeV, we observed similar changes in the tail of the m(t̃t) distribution by applying the |y(t)|, |y(t̃)| < 4 cuts, i.e., mimicking the expected coverage of the future detector. Sudakov logarithms are enhanced, photon contributions are not visibly changed, and the Σ<sub>NLO QCD</sub>/Σ<sub>LO QCD</sub> ratio in the first inset decreases in the tail as observed in Fig. 14.

5 Impact on ATLAS and CMS measurements at 8 TeV

In this section we discuss the relevance of the photon PDF
and EW corrections for experimental measurements at the
LHC. Specifically, we will consider their impact on the anal-
yses already published both by ATLAS and CMS for the t\bar{t}
differential distributions at 8 TeV.

We want to stress that our aim is not to perform a direct
comparison of the SM “best prediction” with experimental

<footnote n="9">⁹ We provided a few representative results also in [53].</footnote>
data. Indeed, at differential level, the best accuracy reached in perturbative QCD is NNLO or NLO+NNLL. Terms beyond NLO in QCD are not taken into account in this analysis, but they are crucial both for a correct determination of the central value and for reducing the scale uncertainty. Here, we want to show that a reliable comparison between SM predictions and experimental data in $t\bar{t}$ production cannot be based only on purely QCD corrections. Indeed, the contribution of EW corrections, especially the uncertainty due to the photon PDF and the scale choice, cannot be ignored and has to be taken into account. In particular cases, such as the normalised distributions discussed in the following, they can yield the largest effect, playing a primary role in the determination of precise SM predictions.

After the discussion and the results presented in the previous section, it is clear that EW corrections at 8 TeV are expected to show a larger dependence on the photon PDF w.r.t. the 13 TeV case. On the other hand, smaller cross sections and integrated luminosities may in practice not lead to enough events to probe this effect. In order to evaluate their real impact at the experimental level, we take into account in our analysis both the theory uncertainties (scale and PDF) and the errors from the experimental measurements. Moreover, this study serves also as a motivation to perform this kind of analysis for future 13 TeV data, where larger cross sections and luminosities should considerably decrease the experimental errors.

In the following, we consider CMS data from Ref. [54], which are based on the dilepton and lepton+jets events collected with an integrated luminosity of 19.7 fb⁻¹. In particular, we focus on the distributions of the transverse momentum and rapidity of the top quark,<fnref n="10" /> and of the rapidity of the $t\bar{t}$ pair. Moreover, we consider ATLAS data based on lepton+jets events collected with an integrated luminosity of 20.3 fb⁻¹. Distributions for the top-quark transverse momentum in the boosted regime are taken from Ref. [55], while those for the invariant mass of the $t\bar{t}$ pair, the absolute value of the rapidity of the top and of the $t\bar{t}$ pair are taken from Ref. [56]. In the ATLAS analyses, part of the electroweak corrections (the purely weak contribution calculated in [20]) has been taken into account via an updated version of HATHOR [57]. Thus, although the Sudakov enhancement has been taken into account, the QED corrections and especially the dependence on the photon PDF, which is a crucial aspect of our discussion, have not been considered yet.

For each differential measurement of ATLAS and CMS, we show two plots with the corresponding theoretical predictions

obtained with the MMHT2014 PDF set [58] (plot on the left), and with NNPDF2.3QED (plots on the right), respectively. At variance with what has been done in Sect. 4, we do not set the photon PDF to zero in the left plots, but we directly use a PDF set that does not include the LO QED contribution in the DGLAP running nor the photon PDF, as typically done in the experimental analyses and in QCD-only predictions. In the plots on the right we decided to use NNPDF2.3QED in order to be conservative, since the effects due to the photon PDF are the largest for this PDF set. We explicitly verified that, as could be expected from the discussion in Sect. 4, using the CT14QED PDFs or setting the photon PDF in NNPDF2.3QED to zero, the relative corrections induced by $\Sigma_{EW}$, e.g. the $\Sigma_{QCD+EW}/\Sigma_{QCD}$ ratio, are the same as those obtained with MMHT2014. Indeed the effect of the LO QED terms in the DGLAP running, which is present in CT14QED and NNPDF2.3QED but not in MMHT2014, cancels in the ratio. On the other hand, this effect can have an impact on the cross section, especially on the $\Sigma_{LO\,QCD}$ term, and it may be important for a direct comparison with data. This is another reason, besides the theoretical consistency, for which NNPDF2.3QED or CT14QED should be in general preferred when NLO EW calculations are performed.

In each plot we present our predictions for a specific measurement, using three different definitions for the renormalisation and factorisation scales, i.e., $\mu = m(t\bar{t})$ (blue), $\mu = H_T/2$ (as in Sect. 4, red) and $\mu = m_t$ (green). In the main panel of each plot we present the $\Sigma_{QCD+EW}$ predictions, without theory uncertainties, compared to data. In the first inset we display the ratio of the $\Sigma_{QCD+EW}$ over data, including scale (light band) and PDF (dark band) uncertainties linearly added. In the second inset we show the central value of the ratio $\Sigma_{QCD+EW}/\Sigma_{QCD}$ for the three different scales and also the relative error in the measurement for each bin, with the error bar centred around one. In the third inset we zoom the PDF uncertainty already included in the main panel. Specifically, we show the PDF uncertainty band of the $\Sigma_{QCD+EW}$ predictions with $\mu = H_T/2$ normalised to their central value and we compare them again with the relative error in the measurement for each bin. Distributions corresponding to CMS analyses [54] are shown in Fig. 16, while ATLAS distributions taken from Refs. [55] and [56] are shown in Figs. 17 and 18, respectively.

Before discussing each plot, it is important to note that two different kinds of distributions are shown here, namely, differential distributions ($d\sigma$) and normalised distributions ($d\sigma/\sigma$).<fnref n="11" /> The $m(t\bar{t})$ distribution in Fig. 17 and the boosted $p_T(t)$ distribution in Fig. 18 are differential distributions, while all the other plots actually display normalised distributions

<footnote n="10">¹⁰ In the plots of this section, the $p_T(t)$ distribution is actually the average of the distribution of the transverse momentum of the top quark and antiquark, consistently with what has been done in the experimental analyses. Similarly, the $y(t)$ distribution is the average of the distribution of the rapidity of the top quark and antiquark. However, differences with real $p(t)$ and $y(t)$ distributions are in practice invisible.</footnote>

<footnote n="11">¹¹ Consistently with what is done in the experimental analyses, we use as $\sigma$ the sum of the values of the bins of the distribution considered, without including overflows.</footnote>
<figure><img src="image_18.png" /><figcaption>**Fig. 16** Comparison of CMS data from Ref. [54] and NLO QCD+EW predictions computed with MMHT2104 (**left**) and NNPDF2.3QED (**right**) PDF sets</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>**Fig. 17** Comparison of ATLAS data from Ref. [56] and NLO QCD+EW predictions computed with MMHT2104 (**left**) and NNPDF2.3QED (**right**) PDF sets</figcaption></figure>
tions, which exhibit much smaller scale and PDF uncertainties. Indeed, for normalised distributions both types of uncertainties are calculated directly for the $d\sigma/\sigma$ ratio, yielding large cancellations between the numerator and the denominator. As we said, our aim is not to perform a direct comparison of data and SM predictions, which would need QCD corrections beyond the NLO. Nevertheless, it is interesting to note how the scale uncertainty shown in the main panel and in the first inset of the plots, which reflects the typical QCD accuracy (NLO) of available Monte Carlo event generators, is of the same order of the experimental precision for all the distributions considered here. Thus, although NNLO and NLO+NNLL have shown to reduce the dependence of the predictions on the renormalisation and factorisation scale, a few comments can be done about the scale choice for NLO simulations. The plots for the $m(t\bar{t})$ distribution in Fig. 17 suggests that $m(t\bar{t})$ is not the natural scale for this process and indeed it artificially leads to particularly large corrections also beyond NLO (see Ref. [16]). Similarly, the discrepancy between data and predictions observed in Ref. [5] at NLO QCD for large $p_T(t)$, but not at NNLO QCD with fixed scale $\mu = m_t$, can be simply reduced at NLO using the scale $\mu = H_T/2$, as shown in the plots in the first line of Fig. 16.

The most important information of the plots is, however, contained in the second and third insets. In the former we compare the experimental precision with the relative corrections induced by $\Sigma_{EW}$ on the predictions at NLO QCD accuracy via the ratio $\Sigma_{QCD+EW}/\Sigma_{QCD}$.<fnref n="12" /> In the latter we compare the experimental precision with the relative PDF uncertainties for the $\Sigma_{QCD+EW}$ predictions. In the following, we discuss these comparisons for every measurement analysed here.

In the plots of the first line in Fig. 16, we can see that the size of electroweak corrections is always smaller than the precision of the experimental measurement for $p_T(t)$ distributions. However, it is reasonable to expect that in the Run-II at 13 TeV, with larger cross sections and higher luminosities, it will be possible to reach an experimental precision such that this effect could be probed, as is also discussed in Ref. [20] and supported by the first CMS analyses [59,60] on 13 TeV data collected with only $2.3 \text{ fb}^{-1}$ integrated luminosity. Nevertheless, it is important to note how results obtained with MMHT2014 and NNPDF2.3QED are qualitatively different; the cancellation between photon-induced contributions and Sudakov logarithms at large $p_T(t)$ also strongly depends on the factorisation scale used. This is particularly important because the dependence of the $\Sigma_{QCD+EW}/\Sigma_{QCD}$ ratio on

the PDF set and on the central scale choice is of the same order, or even exceeds, the scale uncertainty of the NNLO QCD predictions, which is at the percent level.<fnref n="13" /> The PDF uncertainties for $\Sigma_{QCD+EW}$ predictions are smaller than the experimental precision both in the right and left plot, although slightly larger values appear in the NNPDF2.3QED case at large $p_T(t)$ due to the photon-induced contribution.

The case of the $y(t)$ and $y(t\bar{t})$ distributions in the plots in the second and third lines of both Figs. 16 and 17 clearly shows the relevance of EW corrections and photon-initiated processes for a trustworthy comparison of theory predictions and data measurements at the LHC. As can be seen in the first insets, scale uncertainties are smaller than the experimental errors and, for the $y(t)$ distribution, they are further reduced to 1% level or even less by NNLO QCD corrections, thus they are negligible. The theory uncertainties, as can be noted in the second insets, are completely dominated by PDF uncertainties. What is not negligible is the impact of the electroweak corrections. In the case of MMHT2014, without photon-induced processes, EW corrections are flat and thus their impact cancels in the $d\sigma/\sigma$ ratio. On the contrary, in the case of NNPDF2.3QED, for large rapidities their impact is comparable to the experimental error, with a negligible dependence on the factorisation scale. Also, the enhancement in the peripheral region induces a small negative effect in the central region because of the fixed normalisation. PDF uncertainties for the $\Sigma_{QCD+EW}$ are smaller or of the same order of the experimental precision in the case of MMHT2014, while they are larger in the case of NNPDF2.3QED for large rapidities. All these effects are more prominent in the case of the ATLAS measurement (Fig. 17), where experimental errors are smaller. It is also worth to note that for rapidity distributions EW corrections do not receive Sudakov enhancements and, at this level of accuracy, besides the photon-induced contributions, a priori QED effects cannot be neglected.

The situation for the $m(t\bar{t})$ distribution in the first line of Fig. 17 is very similar to the one for the $p_T(t)$ distribution, which has been already discussed for the CMS measurement (first line of Fig. 16) and it is also shown for the boosted regime in the ATLAS analysis in Fig. 18. However, these two differential distributions exhibit much larger scale uncertainties than the normalised $p_T(t)$ distribution in Fig. 16. The inclusion of NNLO QCD corrections is crucial in order to pin down the scale uncertainties below the 5% level, pointing to the necessity of a combination of EW and NNLO QCD corrections for a reliable comparison of theory predictions and data measurements at the LHC. In all the $p_T(t)$ and $m(t\bar{t})$ distributions with the NNPDF2.3QED set, the dependence of

<footnote n="12">¹² In the case of normalised distributions, we actually display the ratio of $d\sigma/\sigma$ with both numerator and the denominator evaluated at QCD+EW accuracy over the same quantity at QCD accuracy $(d\sigma/\sigma)|_{QCD+EW}/(d\sigma/\sigma)|_{QCD}$, in order to correctly identify the effect of the electroweak corrections on the observable considered.</footnote>

<footnote n="13">¹³ The size of NNLO scale uncertainties quoted in this section are taken from Ref. [5]. In particular the authors of Ref. [5] explicitly considered the same $p_T(t)$ and the $y(t)$ CMS distributions of Fig. 16 that are discussed in detail in the text.</footnote>
<figure><img src="image_20.png" /><figcaption>**Fig. 18** Comparison of ATLAS data from Ref. [55] and NLO QCD+EW predictions computed with MMHT2104 (**left**) and NNPDF2.3QED (**right**) PDF sets</figcaption></figure>

the EW corrections on the factorisation scale is also particu-
larly evident, and in the following we describe the reason of
this feature.

The dependence of the size (and possibly of the sign) of
EW corrections on the scale choice can be traced to differ-
ent effects. First of all, unlike QCD corrections, EW cor-
rections to $q\bar{q}$ and $gg$ initial states do not involve neither
$\alpha_s$ renormalisation nor the $O(\alpha_s)$ PDF counter terms that
stabilise the $\mu_f$ dependence. Thus, without photon-initiated
contributions, the $\mu_r$ and $\mu_f$ dependences of $\Sigma_{LO\,QCD}$ and
of $\Sigma_{NLO\,EW}$ are the same; it is different for $\Sigma_{QCD}$. This can
be clearly seen in the left plot of Fig. 18, where in addition
to the $\Sigma_{QCD+EW}/\Sigma_{QCD}$ ratio in the third inset we included
also the $(\Sigma_{LO\,QCD}+\Sigma_{EW})/\Sigma_{LO\,QCD}$ ratio,<fnref n="14" /> where the scale
dependence is the same in the numerator and the denom-
inator. On top of this, as anticipated in the last point in
Sect. 3, $gg$, $q\bar{q}$ and $g\gamma$ luminosities have a different depen-
dence on $\mu_f$. Moreover, at LO the $g\gamma$ channel contributes
only to $\Sigma_{LO\,EW}$, which features only one power of $\alpha_s$ and
thus a different dependence on $\mu_r$ w.r.t. $\Sigma_{LO\,QCD}$. While the
Sudakov logarithms are negative and proportional to the $gg$

and $q\bar{q}$ contributions to $\Sigma_{LO\,QCD}$, the $g\gamma$ contributions to $\Sigma_{EW}$ are positive. Consequently, the cancellation of this two classes of contributions both entering $\Sigma_{EW}$ depends on $\mu$. All these effects are present in any distribution on the right with NNPDF2.3QED. However, in the rapidity distributions there is no Sudakov enhancement and $m_t \sim H_T/2$, especially for large rapidities, so only a small difference for the case $\mu = m(\bar{t}t)$ is visible.

In conclusion, the impact of the photon PDF and also its
dependence on the scale choice should be in general taken
into account for the determination of the theory uncertainties.

**6 Conclusions and outlook**

In this work we studied the impact of EW corrections and
photon-induced contributions on top-quark differential di-
tributions at 8, 13 and 100 TeV. We compared predictions
with two different modern PDF sets including the photon
density and DGLAP evolution at NLO QCD + LO QED accu-
racy: the CT14QED and NNPDF2.3QED PDF sets. While contri-
butions due to the photon PDF are negligible with CT14QED,
this is not the case for NNPDF2.3QED, where such contribu-
tions at 13 TeV are sizeable and are affected by large PDF
uncertainties. At high $p_T(t)$ and $m(t\bar{t})$, the photon-induced
contributions can compensate the negative contributions of
Sudakov logarithms or even change the sign of the EW cor-
rections. Furthermore, we have shown that such a compen-
sation strongly depends on the scale choice. In rapidity dis-
tributions, the impact of the $g\gamma$ initial state is sizeable in the

<footnote n="14">¹⁴ It is worth to note that without photon PDF this ratio corresponds to the NLO EW K-factor, $K^{\text{EW}} = (\Sigma_{\text{LO QCD}} + \Sigma_{\text{NLO EW}})/\Sigma_{\text{LO QCD}}$, which can also be seen as the ratio $\Sigma_{\text{QCD}\times\text{EW}}/\Sigma_{\text{QCD}}$, where $\Sigma_{\text{QCD}\times\text{EW}} = \Sigma_{\text{LO QCD}} \times K^{\text{QCD}} \times K^{\text{EW}}$ is the multiplicative combination of results at NLO QCD and EW accuracy and we have defined as usual $K^{\text{QCD}} = \Sigma_{\text{QCD}}/\Sigma_{\text{LO QCD}}$. Thus the second and third insets in the left plot can be seen as a comparison of $\Sigma_{\text{QCD+EW}}$ and $\Sigma_{\text{QCD\timesEW}}$. On the other hand, the inclusion of $g\gamma$-channels in $\Sigma_{\text{LO EW}}$, such as in the right plot, invalidates the previous argument introducing a non-negligible contribution to $\Sigma_{\text{LO EW}}$.</footnote>
peripheral region and much larger than NLO EW corrections, which do not receive large Sudakov enhancements in these kinematical configurations.

Increasing the energy of the collider, photon-induced channels become less relevant for a fixed value of $p_T(t)$ or $m(t\bar{t})$, since smaller values of $x$ are probed and consequently the quark and gluon PDFs are much larger than the photon PDF. At 100 TeV, photon-induced channels are important only in the very boosted regime ($p_T(t) \gtrsim 5$ TeV or $m(t\bar{t}) \gtrsim 10$ TeV), where Sudakov logarithms are negative and above the 20% level.

For the same reason, at the LHC photon-induced contributions are relatively larger at 8 TeV than at 13 TeV. We computed their size for the same differential (and normalised) distributions already analysed by ATLAS and CMS, taking into account both experimental errors and theory uncertainties. Data from 8 TeV at large rapidities already appear to be sensitive to the photon PDF; with smaller experimental errors, as expected at 13 TeV, such a sensitivity may be reached also at large $p_T(t)$ and $m(t\bar{t})$.

In conclusion, our analyses showed two important points. First, differential observables in top-pair production, in particular $y(t)$ and $y(t\bar{t})$ rapidities, can be used to improve the determination of the photon PDF within the NNPDF approach. Second, given the status of the current available PDF sets, both EW corrections and photon-induced contribution have to be taken into account for a correct determination of both central values and uncertainties of theoretical predictions. The first point is particularly relevant for a future fit with LHC 13 TeV data at (N)NLO QCD as well NLO QED accuracy, which will presumably return a better determination of the photon PDF and consequently, with smaller uncertainties, a more solid comparison with the ansatz for $\gamma(x, Q_0)$ used in CT14QED. The second point suggests also that a priori not only EW corrections but also photon-induced contributions may be in general important for processes at the LHC and must not be neglected. For the case of $t\bar{t}$ production, given the experimental precision that has been reached at the LHC and especially the precision that will be achieved in Run-II, a combination of higher-order QCD corrections, beyond the NLO, and EW corrections, including photon-induced contributions, is mandatory for providing a precise and reliable SM predictions. This is relevant not only for top physics but also for all the analyses that feature $t\bar{t}$ production as background, notably those related to the Higgs boson and new physics.

**Acknowledgments** First and foremost, we want to thank Michelangelo Mangano for having encouraged us to perform this study and for his valuable comments. We also thank Valerio Bertone, Stefano Carrazza, Andrea Giammanco, Fabio Maltoni, Alexander Mitov and Carl Schmidt for interesting discussions and clarifications on different subjects treated in this paper. This work is done in the context of and supported in part (D.P.) by the ERC Grant 291377, LHCtheory: Theoretical predictions and analyses of LHC physics: advancing the precision frontier. I.T. is supported by the F.R.S.-FNRS "Fonds de la

Recherche Scientifique" (Belgium) and in part by the Belgian Federal Science Policy Office through the Interuniversity Attraction Pole P7/37. M.Z. is supported by the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie Grant Agreement No 660171 and in part by the ILP LABEX (ANR-10-LABX-63), in turn supported by French state funds managed by the ANR within the "Investissements d'Avenir" programme under Reference ANR-11-IDEX-0004-02.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP³.

## References

1. M. Beneke, P. Falgari, S. Klein, J. Piclum, C. Schwinn, M. Ubiali, F. Yan, JHEP 07, 194 (2012a). arXiv:1206.2454 [hep-ph]

2. M. Czakon, M.L. Mangano, A. Mitov, J. Rojo, JHEP 07, 167 (2013a). arXiv:1303.7215 [hep-ph]

3. M. Guzzi, K. Lipka, S.-O. Moch, JHEP 01, 082 (2015). arXiv:1406.0386 [hep-ph]

4. M. Czakon, P. Fiedler, A. Mitov, Phys. Rev. Lett. 110, 252004 (2013b). arXiv:1303.6254 [hep-ph]

5. M. Czakon, D. Heymes, A. Mitov, Phys. Rev. Lett. 116, 082003 (2016). arXiv:1511.00549 [hep-ph]

6. S. Moch, P. Uwer, Phys. Rev. D 78, 034003 (2008). arXiv:0804.1476 [hep-ph]

7. N. Kidonakis, Phys. Rev. Lett. 102, 232003 (2009). arXiv:0903.2561 [hep-ph]

8. M. Czakon, A. Mitov, G.F. Sterman, Phys. Rev. D 80, 074017 (2009). arXiv:0907.1790 [hep-ph]

9. V. Ahrens, A. Ferroglia, M. Neubert, B.D. Pecjak, L.L. Yang, JHEP 09, 097 (2010). arXiv:1003.5827 [hep-ph]

10. M. Beneke, P. Falgari, C. Schwinn, Nucl. Phys. B 842, 414 (2011). arXiv:1007.5414 [hep-ph]

11. N. Kidonakis, Phys. Rev. D 82, 114030 (2010). arXiv:1009.4935 [hep-ph]

12. V. Ahrens, A. Ferroglia, M. Neubert, B.D. Pecjak, L.-L. Yang, JHEP 09, 070 (2011a). arXiv:1103.0550 [hep-ph]

13. V. Ahrens, A. Ferroglia, M. Neubert, B.D. Pecjak, L.L. Yang, Phys. Lett. B 703, 135 (2011b). arXiv:1105.5824 [hep-ph]

14. M. Beneke, P. Falgari, S. Klein, C. Schwinn, Nucl. Phys. B 855, 695 (2012b). arXiv:1109.1536 [hep-ph]

15. M. Cacciari, M. Czakon, M. Mangano, A. Mitov, P. Nason, Phys. Lett. B 710, 612 (2012). arXiv:1111.5869 [hep-ph]

16. B.D. Pecjak, D.J. Scott, X. Wang, L.L. Yang, Phys. Rev. Lett. 116, 202001 (2016). arXiv:1601.07020 [hep-ph]

17. W. Beenakker, A. Denner, W. Hollik, R. Mertig, T. Sack, D. Wackeroth, Nucl. Phys. B 411, 343 (1994)

18. J.H. Kuhn, A. Scharf, P. Uwer, Eur. Phys. J. C 45, 139 (2006). arXiv:hep-ph/0508092 [hep-ph]

19. W. Bernreuther, M. Fuecker, Z.G. Si, Phys. Lett. B 633, 54 (2006a). arXiv:hep-ph/0508091 [hep-ph]

20. J.H. Kuhn, A. Scharf, P. Uwer, Eur. Phys. J. C 51, 37 (2007). arXiv:hep-ph/0610335 [hep-ph]

21. W. Bernreuther, M. Fuecker, Z.-G. Si, Phys. Rev. D 74, 113005 (2006b). arXiv:hep-ph/0610334 [hep-ph]

22. W. Bernreuther, M. Fucker, Z.-G. Si, Phys. Rev. D 78, 017503 (2008). arXiv:0804.1237 [hep-ph]
23. A.V. Manohar, M. Trott, Phys. Lett. B 711, 313 (2012). [arXiv:1201.3926 [hep-ph]]

24. J.H. Kühn, A. Scharf, P. Uwer, Phys. Rev. D 91, 014020 (2015). [arXiv:1305.5773 [hep-ph]]

25. J.M. Campbell, D. Wackeroth, J. Zhou, Proceedings, 23rd International Workshop on Deep-Inelastic Scattering and Related Subjects (DIS 2015), PoS DIS2015, 130 (2015). [arXiv:1508.06247 [hep-ph]]

26. W. Hollik, M. Kollar, Phys. Rev. D 77, 014008 (2008). [arXiv:0708.1697 [hep-ph]]

27. W. Bernreuther, Z.-G. Si, Nucl. Phys. B 837, 90 (2010). [arXiv:1003.3926 [hep-ph]]

28. W. Hollik, D. Pagani, Phys. Rev. D 84, 093003 (2011). [arXiv:1107.2606 [hep-ph]]

29. J.H. Kuhn, G. Rodrigo, JHEP 01, 063 (2012). [arXiv:1109.6830 [hep-ph]]

30. W. Bernreuther, Z.-G. Si, Phys. Rev. D 86, 034026 (2012). [arXiv:1205.6580 [hep-ph]]

31. A.D. Martin, R.G. Roberts, W.J. Stirling, R.S. Thorne, Eur. Phys. J. C 39, 155 (2005). [arXiv:hep-ph/0411040 [hep-ph]]

32. J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, O. Matelaer, H.S. Shao, T. Stelzer, P. Torrielli, M. Zaro, JHEP 07, 079 (2014). [arXiv:1405.0301 [hep-ph]]

33. S. Frixione, V. Hirschi, D. Pagani, H.S. Shao, M. Zaro, JHEP 09, 065 (2014). [arXiv:1407.0823 [hep-ph]]

34. S. Frixione, V. Hirschi, D. Pagani, H.S. Shao, M. Zaro, JHEP 06, 184 (2015). [arXiv:1504.03446 [hep-ph]]

35. R.D. Ball, V. Bertone, S. Carrazza, L. Del Debbio, S. Forte, A. Guffanti, N. P. Hartland, J. Rojo (NNPDF), Nucl. Phys. B 877, 290 (2013). [arXiv:1308.0598 [hep-ph]]

36. C. Schmidt, J. Pumplin, D. Stump, C.P. Yuan, Phys. Rev. D 93, 114015 (2016). [arXiv:1509.02905 [hep-ph]]

37. A. Manohar, P. Nason, G.P. Salam, G. Zanderighi, How bright is the proton? A precise determination of the photon PDF, CERN-TH-2016-155 (2016). [arXiv:1607.04266 [hep-ph]]

38. L.A. Harland-Lang, V.A. Khoze, M.G. Ryskin, Photon-initiated processes at high mass, IPPP-16-67 (2016). [arXiv:1607.04635 [hep-ph]]

39. C. Collaboration (CMS), (2015)

40. D. de Florian, G.F.R. Sborlini, G. Rodrigo, Eur. Phys. J. C 76, 282 (2016a). [arXiv:1512.00612 [hep-ph]]

41. D. de Florian, G.F.R. Sborlini, G. Rodrigo, (2016b). [arXiv:1606.02887 [hep-ph]]

42. J.R. Andersen et al. (2016). [arXiv:1605.04692 [hep-ph]]

43. M. Cacciari, M. Greco, P. Nason, JHEP 05, 007 (1998). [arXiv:hep-ph/9803400 [hep-ph]]

44. A.D. Martin, R.G. Roberts, W.J. Stirling, R.S. Thorne, Phys. Lett. B 604, 61 (2004). [arXiv:hep-ph/0410230 [hep-ph]]

45. S. Dulat, T.-J. Hou, J. Gao, M. Guzzi, J. Huston, P. Nadolsky, J. Pumplin, C. Schmidt, D. Stump, C.P. Yuan, Phys. Rev. D 93, 033006 (2016). [arXiv:1506.07443 [hep-ph]]

46. S. Chekanov et al. (ZEUS), Phys. Lett. B 687, 16 (2010). [arXiv:0909.4223 [hep-ex]]

47. V. Bertone, S. Carrazza, J. Rojo, Comput. Phys. Commun. 185, 1647 (2014). [arXiv:1310.1394 [hep-ph]]

48. V. Bertone, S. Carrazza, D. Pagani, M. Zaro, JHEP 11, 194 (2015). [arXiv:1508.07002 [hep-ph]]

49. V. Bertone, S. Carrazza, (2016). [arXiv:1606.07130 [hep-ph]]

50. J. Baglio, L.D. Ninh, M.M. Weber, Phys. Rev. D 88, 113005 (2013). [arXiv:1307.4331]

51. J.H. Kuhn, G. Rodrigo, Phys. Rev. D 59, 054017 (1999). [arXiv:hep-ph/9807420 [hep-ph]]

52. D. Pagani, Proceedings, 11th Hellenic School and Workshops on Elementary Particle Physics and Gravity (CORFU2011), PoS CORFU2011, 027 (2011). [arXiv:1205.6182 [hep-ph]]

53. M.L. Mangano et al., (2016). [arXiv:1607.01831 [hep-ph]]

54. V. Khachatryan et al. (CMS), Eur. Phys. J. C 75, 542 (2015). [arXiv:1505.04480 [hep-ex]]

55. G. Aad et al. (ATLAS), Phys. Rev. D 93, 032009 (2016). [arXiv:1510.03818 [hep-ex]]

56. G. Aad et al. (ATLAS), (2015). [arXiv:1511.04716 [hep-ex]]

57. M. Aliev, H. Lacker, U. Langenfeld, S. Moch, P. Uwer, M. Wiedermann, Comput. Phys. Commun. 182, 1034 (2011). [arXiv:1007.1327 [hep-ph]]

58. L.A. Harland-Lang, A.D. Martin, P. Motylinski, R.S. Thorne, Eur. Phys. J. C 75, 204 (2015). [arXiv:1412.3989 [hep-ph]]

59. Measurement of the inclusive and differential $t\bar{t}$ production cross sections in lepton + jets final states at 13 TeV, Tech. Rep. CMS-PAS-TOP-16-008 (CERN, Geneva, 2016)

60. Measurement of the differential cross section for $t\bar{t}$ production in the dilepton final state at $\sqrt{s} = 13$ TeV, Tech. Rep. CMS-PAS-TOP-16-011 (CERN, Geneva, 2016)