The impact of the photon PDF and electroweak corrections on $t \bar t$ distributions
Abstract
We discuss the impact of EW corrections on differential distributions in topquark pair production at the LHC and future hadron colliders, focusing on the effects of initialstate photons. Performing a calculation at NexttoLeading Order QCD+EW accuracy, we investigate in detail the impact of photoninitiated channels on central values as well as PDF and scale uncertainties, both at order and . 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 toppair production, in particular topquark and 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 photoninduced contributions for a correct determination of both the central value and the uncertainties of theoretical predictions.
∎
CP31628
1 Introduction
The production of a pair of top quarks () 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 production can feature leptons, missing transverse energy and ()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 production gives also access to many of the topquark properties and can be exploited for the determination of the PartonDistributionFunctions (PDF) of the proton Beneke et al. (2012a); Czakon et al. (2013a); Guzzi et al. (2015). Furthermore, since the top quark is the heaviest particle in the SM, it is often a proxy to new physics in BeyondtheStandardModel (BSM) scenarios. Thus, 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 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 put in place in order to improve the accuracy of theoretical predictions: inclusive and differential crosssections have been calculated up to NexttoNexttoLeadingOrder (NNLO) accuracy in QCD Czakon et al. (2013b, 2016), and including soft gluon effects up to NexttoNexttoLeadingLogarithmic (NNLL) accuracy Moch and Uwer (2008); Kidonakis (2009); Czakon et al. (2009); Ahrens et al. (2010); Beneke et al. (2011); Kidonakis (2010); Ahrens et al. (2011a, b); Beneke et al. (2012b); Cacciari et al. (2012); Czakon et al. (2013a); Pecjak et al. (2016). Notably, higherorder 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 () may not be respected. On the one hand, when large scales are probed (), the socalled 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 nonnegligible size, leading to possible compensations of the corrections induced by Sudakov logarithms.
While weak Beenakker et al. (1994); Kuhn et al. (2006); Bernreuther et al. (2006a); Kuhn et al. (2007); Bernreuther et al. (2006b, 2008); Manohar and Trott (2012); Kühn et al. (2015); Campbell et al. (2015), QED Hollik and Kollar (2008) and EW (weak+QED) Bernreuther and Si (2010); Hollik and Pagani (2011); Kuhn and Rodrigo (2012); Bernreuther and Si (2012) corrections to 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 photoninduced subprocesses has been addressed only in Hollik and Kollar (2008) for only the channel at LO and using the only PDF set with the photon distribution available at that time, i.e., MRST2004QED Martin et al. (2005).
In this work we (re)evaluate, besides the NLO QCD corrections, the complete set of NLO EW corrections to topquark pair production and carefully assess their impact on differential distributions for proton–proton collisions at 8, 13 and 100 TeV centreofmass energy. For the first time, photoninduced contributions are taken into account both at , from the channel at treelevel, and at , from the as well as the 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 Alwall et al. (2014), thanks to an extension of the code that allows to automatically calculate NLO QCD and EW corrections Frixione et al. (2014, 2015).
In order to have a reliable estimate of the photoninduced contribution and of its uncertainties, we evaluate the impact of electroweak corrections with both the NNPDF2.3QED Ball et al. (2013) and CT14QED Schmidt et al. (2016) 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
We show that, besides higherorder QCD corrections, it is also necessary to take into account EW corrections, including photoninduced 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 on NLOaccurate event generators and the compatibility with experimental data for the distribution at the LHC 8TeV measurements Collaboration (2015) and in view of the measurements at 13 TeV. Furthermore, we show that differential observables in toppair production, in particular topquark and rapidities, can be used to improve the determination of the photon PDF within the NNPDF approach, while in the CT14QED approach production is not sensitive to the photoninduced contributions.
The structure of our paper is the following: 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 photoninduced 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 QCDonly partons and DGLAP evolution. We give our conclusions and 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 Frixione et al. (2014, 2015), where the 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 Frixione et al. (2014) and precisely define the quantities included in our calculation, specifying those that depend on the photon PDF.
In the case of process a generic observable can be expanded simultaneously in the QCD and EW coupling constants as:
(2.1) 
Following this notation, one can separate the LO () and NLO () contributions as follows:
(2.2) 
In our results we include the , , and terms. We checked that the remaining terms are subleading as expected, giving results of the order or below 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
(2.3)  
(2.4) 
(2.5)  
(2.6)  
(2.7) 
In the following text with the term “EW corrections” we will in general refer to the quantity , while we will use “NLO EW corrections” for . At variance with refs. Frixione et al. (2014, 2015), in our calculation we do not consider the effect due to the HeavyBosonRadiation (HBR). Although the LO cross sections of , processes may in principle contribute at the same perturbative order of to the inclusive production, in this work we assume that and final states can always be distinguished.
The quantities that depend on the photon PDF are and . The process contributes to and to . In addition, receives contributions from the and processes at the tree level, which feature initialstate singularities that have to be subtracted. As in the case of 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 de Florian et al. (2016a, b).
Our calculation is performed using the following input parameters
(2.8) 
and setting all the other fermion masses to zero. All masses are renormalised onshell and all decay widths are set to zero. The renormalisation of is performed in the fiveflavour scheme, while EW parameters are chosen in the scheme, with
(2.9) 
Since NLO EW corrections of to hadroproduction do not involve the renormalisation of , 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 () and factorisation () scales defined as
(2.10) 
where the sum of the transverse masses runs over all the finalstate particles. In all cases theoretical uncertainties due to missing higher orders are estimated via independent variations of and in the interval .
It is worth to note that the NNPDF2.3QED set is in the variableflavour scheme with six active flavours, which for is equivalent to the sixflavour scheme. On the contrary, in CT14QED the active flavours are five, leading to the fiveflavour scheme also for . As we said, we renormalise in the fiveflavour 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 be easily performed by following the recipe described in Andersen et al. (2016) and based on Cacciari et al. (1998), which, at NLO QCD + NLO EW accuracy, has a direct effect only on the induced contribution to . 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 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 details the different modelling of the photon PDF in the NNPDF2.3QED, CT14QED
Although for all the three PDF sets the DGLAP evolution is performed at NLO QCD + LO QED accuracy
In fig. 1 we show the photon PDF for the different sets at the scales GeV. As it can be seen, these assumptions affect the dependence on and for both the central value and the uncertainty band

Consistently with the approach pursued for coloured partons, in NNPDF2.3QED no functional form is specified for the photon PDF at the initial scale, . The photon PDF is only constrained to be positive. In a first step, PDF replicas for all the partons are fit together from DISdata only. Afterwards, they are further constrained by DrellYan data form the LHC RunI at 7 TeV. At variance with DIS, neutralcurrent DrellYan is sensitive to the photon PDF already at LO, and it can put stronger constraints on . Because of the positivity requirement for the photon PDF, the replicas distribute in a very nonGaussian 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 , in this region uncertainties are very large and the central value alone can be misleading, leading even to an unphysical peak at very large , which can be seen in fig. 1.

The CT14QED and MRST2004QED sets are based on a completely different assumption for at the initial scale . Since the amount of data for constraining the shape of the photon PDF is limited, is chosen to be described by an ansatz; the photon PDF at is assumed to be completely determined by the valence quark distributions. Specifically, in CT14QED and MRST2004QED the photon parameterisation at the initial scale GeV reads
(3.1) In eq. (3.1) corresponds to the convolution of the splitting function with the socalled “primordial” quark distributions . In the case of MRST2004QED, are valencelike model distributions, i.e., they are not those fitted within the global MRST2004 Martin et al. (2004) analysis. In the case of CT14QED, correspond to the initial up and down valence distributions from CT14NLO Dulat et al. (2016).
More importantly, in CT14QED is set equal to in order to obtain a dependence on a single parameter. Conversely, in MRST2004QED, the coefficients are given by , where are the “Current Mass” (CM) of the quarks ( MeV and 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 CM for the up and the down quark. The CT14QED set leaves freedom for the value of in the coefficients, which can be alternatively parameterised by the momentum fraction carried by the photon at the initial scale. The constraints on and consequently on the photon PDF are then obtained by fitting the ZEUS data Chekanov et al. (2010) for the production of isolated photons in DeepInelasticScattering (DIS), , a process which is sensitive to the photon PDF already at LO. In this way, the constraint is set at 90% CL for the photon PDF, at the initial scale of GeV. For this reason, the uncertainty band in CT14QED corresponds to the area between the and 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 CM 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 QCD and QED evolution in CT14QED and MRST2004QED, 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 the central values of the three sets are different, although the CT14QED and NNPDF2.3QED photons are well compatible within the uncertainty. At large , while CT14QED and MRST2004QED converge to similar values, the prediction from NNPDF2.3QED is different. This trend is particularly visible at small and for large scales, where the photon PDF mostly stems from the quark densities via the DGLAP evolution. However, it is worth to remind that a new correct implementation of the coupled QCD+QED DGLAP equations is now also available in APFEL Bertone et al. (2014) and has been described in Bertone et al. (2015) in the context of the determination of new sets of PDFs including also lepton members.
^{6} 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 Bertone et al. (2015). This set is the same of NNPDF2.3QED at the initial scale , 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 , while for small and large it converges to the values of CT14QED and MRST2004QED. This proves that the discrepancy between the different sets at small and large is completely driven by the different treatment of the DGLAP evolution.
In fig. 2 we show the photongluon luminosity as function of the invariant mass for the three PDF sets considered at 8, 13 and 100 TeV, with the factorisation scale set equal to . In the first inset we show for every luminosity the ratio of their values at and , while in the second(third) inset we show at () the ratio of every luminosity over the corresponding central value with NNPDF2.3QED. For the entire range in considered at the three different protonprotoncollision energies, the CT14QED and NNPDF2.3QED photongluon luminosities are barely compatible or not compatible at all. The slightly better compatibility with 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 and large , 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 luminosity is very close to the NNPDF2.3QED one and inside its uncertainty band. The insensitivity of our study on the different treatment of QCD+QED DGLAP running is in principle expected for most of the predictions at the LHC, since large factorisation scales are typically not correlated with small values of . Moreover, at small values of , where small values of can be probed, the PDFs of quarks and gluons are much larger than the photon PDF, leading to a relative suppression of photoninitiated contributions.
In order to determine at which values of the photon–gluon luminosity is actually expected to be relevant in hadroproduction,
in fig. 3 we show the photongluon, gluongluon and quarkantiquark luminosities for NNPDF2.3QED sets at 8, 13 and 100 TeV. These plots will also become useful to understand the dependence of on the scale definition.
The solid lines refer to the luminosities with , with the uncertainty band shown only for the photongluon case. The samecolour dashed lines are the corresponding quantities with . The ratios between predictions with and with is shown in the first inset. The luminosities at are in general of the same order or smaller than at , with the exception of the 100 TeV case for . In particular at 8 and 13 TeV, in the range of the luminosities of and strongly reduce at with respect to , while the photongluon luminosity is less sensitive to the value of . Since in the TeV range the negative contributions from Sudakov logarithms in the and channels and the positive photongluon contribution tend to cancel each other, larger scales lead to larger values of . This effect will be discussed both in sect. 4 for predictions at 13 TeV and in sect. 5, where the 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 and luminosities over the luminosity at () in order to better identify at which scales the contribution is expected to be relevant. From these insets, the hierarchy between and channels can also be easily derived.
In conclusion, the different results for photoninduced 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 at the initial scale 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 toppair invariant mass , the transverse momentum of the top quark , and the rapidities of the top quark and pair . Predictions for the LHC at 13 TeV are shown in figs. 49, while those for a Future Circular Collider (FCC) at 100 TeV are shown in figs. 1015. Unless differently specified, results are obtained with no cut imposed on the finalstate particles.
In each figure we show two plots for the same observable, displaying in the left plot, denoted as “no ”, 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 to gauge the impact of the photon PDF on both the central value and 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 by the photoninduced 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 is expected to be stable under future PDF determination.
In each plot the main panel includes the distributions for , and predictions as defined in equations (2.3)(2.7). The four insets below the main panel display ratios of several quantities over the central value. In the first inset we plot this ratio for the central values of , and , i.e., the other perturbative orders which we consider. The second and third inset respectively present the ratio including the scale and PDF uncertainties for the numerator. In the last inset we present the 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 the right, including
effects due to the photon PDF, also show (in the last inset) the ratio computed with the CT14QED PDF set.
The first general feature that can be noticed in all the distributions for the LHC at 13 TeV is the different behaviour 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.3QED.
Indeed, by comparing the last inset in equivalent plots with (right) and without (left) photoninduced 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 photoninduced processes originates from the contribution of the channel. As anticipated in sect. 3, the contribution is negligible
As can be seen in figs. 4 and 5, at 13 TeV the photoninitiated contributions can compensate the Sudakov suppression in the toppair invariant mass and topquark transverse momentum distributions. At very large scales ( or ) 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 with both quantities evaluated at the scale . This quantity indicates that the cancellation between Sudakov logarithms strongly depends also on the scale choice, thus the scale dependence of 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 term alone are directly compared with the PDF uncertainties, the grey band. In the tail of the and distributions on the right, PDF uncertainties for separately and are of the same size. This is consistent with the difference between the PDF uncertainty for the total prediction with and without effects form 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 . Just for reference, we decided to show in the second inset the ratio with scale uncertainties, which are reduced w.r.t. the 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 . However, one should consider that the scale uncertainty is further reduced to very few percents by NNLO corrections, which have already been calculated Czakon et al. (2013b, 2016). Moreover, when normalised distributions () 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 will also investigate the impact of EW corrections and of the photon PDF.
While distributions do not change by requiring , , i.e. top and antitop produced in the central region of the detector, these cuts have an effect on predictions for high . This is shown in fig. 6, which is equivalent to fig. 5, but with cuts. At large invariant masses, tops are preferably produced in the forward or backward region, due to the  and channel diagrams in the channel, which is much less suppressed than the channels featuring only channel diagrams at LO QCD. Rapidity cuts suppress the contribution, as well as contributions, but also enhance the typical value of the partonic Mandelstam variables and . Consequently, with those cuts, we observe larger values for the Sudakov logarithms (plot on the left) and a similar behaviour for the photoninduced contributions. Moreover, we can notice that, as in the distributions in fig. 4, also in the tails of the plots in fig. 6 the ratio 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 section originating from the initial state, which is decreasing in the tail of distribution and, only with rapidity cuts applied, of the distribution.
In figs. 7 and 8 we show predictions for and distributions. EW corrections are quite flat in the plots on the left, while in the plots on the right they are enhanced by the initial state in the peripheral region, especially for the distribution. In this region photoninitiated contributions and their uncertainties become relevant since the photon PDF is sampled at rather large , where the suppression w.r.t the quark and gluon PDFs is reduced. At variance with the tail of or , where the 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, largerapidity 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 production at the LHC 13 TeV, in fig. 9 we also show the same plots of fig. 7 with a cut applied. EW corrections in the left plot are larger in the central region, where, due to large values of the and Mandelstam variables, the effect of Sudakov logarithms is enhanced. Again, this effect can be compensated by photoninduced 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 are negligible, while they are sizeable in the plots on the right. Also in this case the trend of the ratio displayed in the first inset is correlated with the fraction of the cross section originating from the 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 distributions, the transverse momentum of the antitop, the distribution is different from when terms are included, because of the charge asymmetry Kuhn and Rodrigo (1999); Hollik and Pagani (2011); Kuhn and Rodrigo (2012); Bernreuther and Si (2012). However, we observed the same qualitative behaviour for the photon PDF and the EW corrections in and distributions. Moreover, the channel cannot give a contribution to the numerator of the charge asymmetry (see also ref. Pagani (2011)), and all the other contributions to the asymmetry have already been investigated in ref. Bernreuther and Si (2012).
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.
In order to find a large contribution form the 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 and with up to 12 TeV and up to 30 TeV. First of all, we remind that at a FCC it would be possible with integrated luminosity to probe such small cross sections and large hardscales. Thus, the following discussion is not purely academic. Sudakov logarithms can induce up to corrections for and for , 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 final state; it is a general feature for all the production processes (see also the EW section in Mangano et al. (2016)). At such large scales, also the 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 distribution by applying the cuts, i.e., mimicking the expected coverage of the future detector. Sudakov logarithms are enhanced, photon contributions are not visibly changed, and the 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 analyses already published both by ATLAS and CMS for the 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 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 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 in to 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. Khachatryan et al. (2015), which are based on the dilepton and lepton+jets events collected with an integrated luminosity of . In particular, we focus on the distributions of the transverse momentum and rapidity of the top quark,
For each differential measurement of ATLAS and CMS, we show two plots with the corresponding theoretical predictions obtained with the MMHT2014 PDF set HarlandLang et al. (2015) (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 neither LO QED contribution in the DGLAP running nor the photon PDF, as typically done in the experimental analyses and in QCDonly 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 it 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 , e.g. the 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 term, and 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., (blue), (as in sect. 4, red) and (green). In the main panel of each plot we present the predictions, without theory uncertainties, compared to data. In the first inset we display the ratio of the 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 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 predictions with 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 Khachatryan et al. (2015) are shown in fig. 16, while ATLAS distributions taken from ref. Aad et al. (2016) and ref. Aad et al. (2015) are shown in fig. 17 and fig. 18, respectively.
Before discussing each plot, it is important to note that two different kinds of distributions are shown here, namely, differential distributions () and normalised distributions ().
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 on the predictions at NLO QCD accuracy via the ratio .
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 distributions.
However, it is reasonable to expect that in the RunII 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 it is also discussed in ref. Kuhn et al. (2007) and supported by the first CMS analyses CMS (2016a, b) on 13 TeV data collected with only integrated luminosity. Nevertheless, it is important to note how results obtained with MMHT2014 and NNPDF2.3QED are qualitatively different; the cancellation between photoninduced contributions and Sudakov logarithms at large also strongly depends on the factorisation scale used. This is particularly important because the dependence of the 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.
The case of the and distributions in the plots in the second and third lines of both figs. 16 and 17 clearly shows the relevance of EW corrections and photoninitiated 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 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 photoninduced processes, EW corrections are flat and thus their impact cancels in the 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 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 photoninduced contributions, a priori QED effects cannot be neglected.
The situation for the distribution in the first line of fig. 17 is very similar to the one for the 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 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 and distributions with the NNPDF2.3QED set, the dependence of the EW corrections on the factorisation scale is also particularly 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 different effects. First of all, unlike QCD corrections, EW corrections to and initial states do not involve neither renormalisation nor the PDF counter terms that stabilise the dependence. Thus, without photoninitiated contributions, both the and dependence
of and is the same; it is different for . This can be clearly seen in the left plot of fig. 18, where in addition to the ratio in the third inset we included also the
ratio
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 photoninduced contributions on topquark differential distributions 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 accuracy: the CT14QED and NNPDF2.3QED PDF sets. While contributions due to the photon PDF are negligible with CT14QED, this is not the case for NNPDF2.3QED, where such contributions at 13 TeV are sizeable and are affected by large PDF uncertainties. At high and , the photoninduced contributions can compensate the negative contributions of Sudakov logarithms or even change the sign of the EW corrections. Furthermore, we have shown that such a compensation strongly depends on the scale choice. In rapidity distributions, the impact of the initial state is sizeable in the 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, photoninduced channels become less relevant for a fixed value of or , since smaller values of are probed and consequently the quark and gluon PDFs are much larger than the photon PDF. At 100 TeV, photoninduced channels are important only in the very boosted regime ( or ), where Sudakov logarithms are negative and above the level.
For the same reason, at the LHC photoninduced 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 and .
In conclusion, our analyses showed two important points. First, differential observables in toppair production, in particular and 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 photoninduced 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 used in CT14QED. The second point suggests also that a priori not only EW corrections but also photoninduced contributions may be in general important for processes at the LHC and must not be neglected. For the case of production, given the experimental precision that has been reached at the LHC and especially will be achieved in the RunII, a combination of higherorder QCD corrections, beyond the NLO, and EW corrections, including photoninduced 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 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 SklodovskaCurie grant agreement No 660171 and in part by the ILP LABEX (ANR10LABX63), in turn supported by French state funds managed by the ANR within the “Investissements d’Avenir” programme under reference ANR11IDEX000402.
Footnotes
 journal: Eur. Phys. J. C
 We acknowledge that two new studies on the determination of the photon PDF Manohar et al. (2016); HarlandLang et al. (2016) have appeared during the publication of our work.
 CT14QED provides two kinds of sets, one with only the incoherent component of the photon PDF and another one with both the coherent and 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 including both components are very similar, in the and ranges relevant for our study, to those from the incoherentonly photon PDF with .
 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.
 The red dotted curve labelled as APFEL_NN23 will be explained after in the text.
 During the writing of this article also a new set NNPDF3.0QED has been released Bertone and Carrazza (2016) and it includes the correct DGLAP running at (N)NLO QCD + LO QED accuracy.
 In our calculation there are also contribution from quarkphoton initial states, however, as it 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.
 Only for this ratio, also the quantity has been evaluated with the CT14QED PDF set.
 This statement could be done only a posteriori, after having performed the calculation and cannot be generalised to other processes. See for instance ref. Baglio et al. (2013).
 We provided a few representative results also in Mangano et al. (2016).
 In the plots of this section, the 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 distribution is the average of the distribution of the rapidity of the top quark and antiquark. However, differences with real and distributions are in practice invisible.
 Consistently with what is done in the experimental analyses, we use as the sum of the values of the bins of the distribution considered, without including overflows.
 In the case of normalised distributions, we actually display the ratio of with both numerator and the denominator evaluated at QCD+EW accuracy over the same quantity at QCD accuracy , in order to correctly identify the effect of the electroweak corrections on the observable considered.
 The size of NNLO scale uncertainties quoted in this section are taken from ref. Czakon et al. (2016). In particular the authors of ref. Czakon et al. (2016) explicitly considered the same and the CMS distributions of fig. 16 that are discussed in detail in the text.
 It is worth to note that without photon PDF this ratio corresponds to the NLO EW factor, , which can be also seen as the ratio , where is the multiplicative combination of results at NLO QCD and EW accuracy and we have defined as usual . Thus the second and third insets in the left plot can be seen as a comparison of and . On the other hand, the inclusion of channels in , such as in the right plot, invalidates the previous argument introducing a nonnegligible contribution to .
References
 M. Beneke, P. Falgari, S. Klein, J. Piclum, C. Schwinn, M. Ubiali, and F. Yan, JHEP 07, 194 (2012a), arXiv:1206.2454 [hepph] .
 M. Czakon, M. L. Mangano, A. Mitov, and J. Rojo, JHEP 07, 167 (2013a), arXiv:1303.7215 [hepph] .
 M. Guzzi, K. Lipka, and S.O. Moch, JHEP 01, 082 (2015), arXiv:1406.0386 [hepph] .
 M. Czakon, P. Fiedler, and A. Mitov, Phys. Rev. Lett. 110, 252004 (2013b), arXiv:1303.6254 [hepph] .
 M. Czakon, D. Heymes, and A. Mitov, Phys. Rev. Lett. 116, 082003 (2016), arXiv:1511.00549 [hepph] .
 S. Moch and P. Uwer, Phys. Rev. D78, 034003 (2008), arXiv:0804.1476 [hepph] .
 N. Kidonakis, Phys. Rev. Lett. 102, 232003 (2009), arXiv:0903.2561 [hepph] .
 M. Czakon, A. Mitov, and G. F. Sterman, Phys. Rev. D80, 074017 (2009), arXiv:0907.1790 [hepph] .
 V. Ahrens, A. Ferroglia, M. Neubert, B. D. Pecjak, and L. L. Yang, JHEP 09, 097 (2010), arXiv:1003.5827 [hepph] .
 M. Beneke, P. Falgari, and C. Schwinn, Nucl. Phys. B842, 414 (2011), arXiv:1007.5414 [hepph] .
 N. Kidonakis, Phys. Rev. D82, 114030 (2010), arXiv:1009.4935 [hepph] .
 V. Ahrens, A. Ferroglia, M. Neubert, B. D. Pecjak, and L.L. Yang, JHEP 09, 070 (2011a), arXiv:1103.0550 [hepph] .
 V. Ahrens, A. Ferroglia, M. Neubert, B. D. Pecjak, and L. L. Yang, Phys. Lett. B703, 135 (2011b), arXiv:1105.5824 [hepph] .
 M. Beneke, P. Falgari, S. Klein, and C. Schwinn, Nucl. Phys. B855, 695 (2012b), arXiv:1109.1536 [hepph] .
 M. Cacciari, M. Czakon, M. Mangano, A. Mitov, and P. Nason, Phys. Lett. B710, 612 (2012), arXiv:1111.5869 [hepph] .
 B. D. Pecjak, D. J. Scott, X. Wang, and L. L. Yang, Phys. Rev. Lett. 116, 202001 (2016), arXiv:1601.07020 [hepph] .
 W. Beenakker, A. Denner, W. Hollik, R. Mertig, T. Sack, and D. Wackeroth, Nucl. Phys. B411, 343 (1994).
 J. H. Kuhn, A. Scharf, and P. Uwer, Eur. Phys. J. C45, 139 (2006), arXiv:hepph/0508092 [hepph] .
 W. Bernreuther, M. Fuecker, and Z. G. Si, Phys. Lett. B633, 54 (2006a), arXiv:hepph/0508091 [hepph] .
 J. H. Kuhn, A. Scharf, and P. Uwer, Eur. Phys. J. C51, 37 (2007), arXiv:hepph/0610335 [hepph] .
 W. Bernreuther, M. Fuecker, and Z.G. Si, Phys. Rev. D74, 113005 (2006b), arXiv:hepph/0610334 [hepph] .
 W. Bernreuther, M. Fucker, and Z.G. Si, Phys. Rev. D78, 017503 (2008), arXiv:0804.1237 [hepph] .
 A. V. Manohar and M. Trott, Phys. Lett. B711, 313 (2012), arXiv:1201.3926 [hepph] .
 J. H. Kühn, A. Scharf, and P. Uwer, Phys. Rev. D91, 014020 (2015), arXiv:1305.5773 [hepph] .
 J. M. Campbell, D. Wackeroth, and J. Zhou, Proceedings, 23rd International Workshop on DeepInelastic Scattering and Related Subjects (DIS 2015), PoS DIS2015, 130 (2015), arXiv:1508.06247 [hepph] .
 W. Hollik and M. Kollar, Phys. Rev. D77, 014008 (2008), arXiv:0708.1697 [hepph] .
 W. Bernreuther and Z.G. Si, Nucl. Phys. B837, 90 (2010), arXiv:1003.3926 [hepph] .
 W. Hollik and D. Pagani, Phys. Rev. D84, 093003 (2011), arXiv:1107.2606 [hepph] .
 J. H. Kuhn and G. Rodrigo, JHEP 01, 063 (2012), arXiv:1109.6830 [hepph] .
 W. Bernreuther and Z.G. Si, Phys. Rev. D86, 034026 (2012), arXiv:1205.6580 [hepph] .
 A. D. Martin, R. G. Roberts, W. J. Stirling, and R. S. Thorne, Eur. Phys. J. C39, 155 (2005), arXiv:hepph/0411040 [hepph] .
 J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, O. Mattelaer, H. S. Shao, T. Stelzer, P. Torrielli, and M. Zaro, JHEP 07, 079 (2014), arXiv:1405.0301 [hepph] .
 S. Frixione, V. Hirschi, D. Pagani, H. S. Shao, and M. Zaro, JHEP 09, 065 (2014), arXiv:1407.0823 [hepph] .
 S. Frixione, V. Hirschi, D. Pagani, H. S. Shao, and M. Zaro, JHEP 06, 184 (2015), arXiv:1504.03446 [hepph] .
 R. D. Ball, V. Bertone, S. Carrazza, L. Del Debbio, S. Forte, A. Guffanti, N. P. Hartland, and J. Rojo (NNPDF), Nucl. Phys. B877, 290 (2013), arXiv:1308.0598 [hepph] .
 C. Schmidt, J. Pumplin, D. Stump, and C. P. Yuan, Phys. Rev. D93, 114015 (2016), arXiv:1509.02905 [hepph] .
 A. Manohar, P. Nason, G. P. Salam, and G. Zanderighi, (2016), arXiv:1607.04266 [hepph] .
 L. A. HarlandLang, V. A. Khoze, and M. G. Ryskin, (2016), arXiv:1607.04635 [hepph] .
 C. Collaboration (CMS), (2015).
 D. de Florian, G. F. R. Sborlini, and G. Rodrigo, Eur. Phys. J. C76, 282 (2016a), arXiv:1512.00612 [hepph] .
 D. de Florian, G. F. R. Sborlini, and G. Rodrigo, (2016b), arXiv:1606.02887 [hepph] .
 J. R. Andersen et al. (2016) arXiv:1605.04692 [hepph] .
 M. Cacciari, M. Greco, and P. Nason, JHEP 05, 007 (1998), arXiv:hepph/9803400 [hepph] .
 A. D. Martin, R. G. Roberts, W. J. Stirling, and R. S. Thorne, Phys. Lett. B604, 61 (2004), arXiv:hepph/0410230 [hepph] .
 S. Dulat, T.J. Hou, J. Gao, M. Guzzi, J. Huston, P. Nadolsky, J. Pumplin, C. Schmidt, D. Stump, and C. P. Yuan, Phys. Rev. D93, 033006 (2016), arXiv:1506.07443 [hepph] .
 S. Chekanov et al. (ZEUS), Phys. Lett. B687, 16 (2010), arXiv:0909.4223 [hepex] .
 V. Bertone, S. Carrazza, and J. Rojo, Comput. Phys. Commun. 185, 1647 (2014), arXiv:1310.1394 [hepph] .
 V. Bertone, S. Carrazza, D. Pagani, and M. Zaro, JHEP 11, 194 (2015), arXiv:1508.07002 [hepph] .
 V. Bertone and S. Carrazza, (2016), arXiv:1606.07130 [hepph] .
 J. Baglio, L. D. Ninh, and M. M. Weber, Phys. Rev. D88, 113005 (2013), arXiv:1307.4331 .
 J. H. Kuhn and G. Rodrigo, Phys. Rev. D59, 054017 (1999), arXiv:hepph/9807420 [hepph] .
 D. Pagani, Proceedings, 11th Hellenic School and Workshops on Elementary Particle Physics and Gravity (CORFU2011), PoS CORFU2011, 027 (2011), arXiv:1205.6182 [hepph] .
 M. L. Mangano et al., (2016), arXiv:1607.01831 [hepph] .
 V. Khachatryan et al. (CMS), Eur. Phys. J. C75, 542 (2015), arXiv:1505.04480 [hepex] .
 G. Aad et al. (ATLAS), Phys. Rev. D93, 032009 (2016), arXiv:1510.03818 [hepex] .
 G. Aad et al. (ATLAS), (2015), arXiv:1511.04716 [hepex] .
 M. Aliev, H. Lacker, U. Langenfeld, S. Moch, P. Uwer, and M. Wiedermann, Comput. Phys. Commun. 182, 1034 (2011), arXiv:1007.1327 [hepph] .
 L. A. HarlandLang, A. D. Martin, P. Motylinski, and R. S. Thorne, Eur. Phys. J. C75, 204 (2015), arXiv:1412.3989 [hepph] .
 Measurement of the inclusive and differential tt production cross sections in lepton + jets final states at 13 TeV, Tech. Rep. CMSPASTOP16008 (CERN, Geneva, 2016).
 Measurement of the differential cross section for production in the dilepton final state at , Tech. Rep. CMSPASTOP16011 (CERN, Geneva, 2016).