CERN-TH-2017-235

Nikhef/2017-064

Illuminating the photon content

[0.2cm] of the proton within a global PDF analysis

The NNPDF Collaboration:

[0.1cm] Valerio Bertone, Stefano Carrazza, Nathan P. Hartland, and Juan Rojo.

Department of Physics and Astronomy, VU University, NL-1081 HV Amsterdam,

and Nikhef Theory Group, Science Park 105, 1098 XG Amsterdam, The Netherlands

[0.1cm] Theoretical Physics Department, CERN, CH-1211 Geneva, Switzerland

Abstract

Precision phenomenology at the LHC requires accounting for both higher-order QCD and electroweak corrections as well as for photon-initiated subprocesses. Building upon the recent NNPDF3.1 fit, in this work the photon content of the proton is determined within a global analysis supplemented by the LUXqed constraint relating the photon PDF to lepton-proton scattering structure functions: NNPDF3.1luxQED. The uncertainties on the resulting photon PDF are at the level of a few percent, with photons carrying up to of the proton’s momentum. We study the phenomenological implications of NNPDF3.1luxQED at the LHC for Drell-Yan, vector boson pair, top quark pair, and Higgs plus vector boson production. We find that photon-initiated contributions can be significant for many processes, leading to corrections of up to . Our results represent a state-of-the-art determination of the partonic structure of the proton including its photon component.

###### Contents

## 1 Introduction

Recent progress in the computation of higher-order QCD corrections to LHC processes is such that the current state-of-the-art accuracy is NNLO, with even NLO calculations available in some relevant cases (see Ref. [1] for a review). At this level of theoretical precision, the inclusion of electroweak (EW) corrections becomes phenomenologically relevant. With this motivation, NLO EW corrections to hard-scattering matrix elements have been computed for many LHC processes, including single and double vector boson, inclusive jets and dijets, and top quark pair production, among others [2, 3, 4, 5, 6, 7, 8, 9, 10]. Alongside progress made in process-specific calculations, the automation of NLO EW calculations [11, 12, 13] has also advanced significantly.

In order to make the most of these developments in the calculation of higher-order QCD and EW corrections, equivalent progress in the determination of the parton distribution functions (PDFs) of the proton [14] is vital. In this respect, most of the recent global PDF analyses [15, 16, 17, 18, 19] are indeed based on NNLO QCD theory. On the other hand, PDF analyses that include QED and weak effects and a determination of the photon PDF are scarcer [20, 21, 22, 23]. Such QED PDF sets are required by consistency once EW corrections to matrix elements are included, as well as to account for the effects of photon-initiated (PI) subprocesses.

Indeed, the inclusion of QED and weak effects into a global PDF fit requires two main modifications. Firstly, hard-scattering matrix elements have to be corrected for EW effects where relevant. This also implies taking into account the contributions from photon-induced subprocesses. This leads to the second main modification, which is the introduction of an additional parton distribution quantifying the photon content of the proton. In turn, this requires generalising the DGLAP evolution equations to account for QED corrections. This generalisation is made possible thanks to the computation of the splitting functions up to and [24, 25], with the resulting QED-corrected evolution equations implemented in public PDF evolution codes such as APFEL [26], HOPPET [27], and QEDEVOL [28].

Until recently, two distinct strategies were adopted for the determination the photon PDF: model calculations and data-driven approaches. In the first case, the photon PDF is computed on the base of a theoretically motivated model ansatz. The original realisation of this strategy was the MRST04QED set [20], where the photon PDF was generated at some low scale by one-photon collinear emission off a model for the valence quarks. This led to a simple relation between the photon PDF and the up and down valence distributions, which was then evolved upwards in using the DGLAP equations corrected for contributions.

A similar approach was adopted in the CT14QED analysis [23]. However, differently from MRST04QED, the CT14QED determination includes also an elastic component, based on the equivalent photon approximation [29] that accounts for the fact that the proton can emit photons while remaining intact. This elastic component of the photon PDF can be determined by a QED calculation in a model-independent way [30, 31, 32, 33, 34] in terms of the electric and magnetic form factors of the proton. In addition, an estimate of the uncertainty in their model for the inelastic component of photon PDF was provided based on HERA isolated photon production data [35].

In the second strategy, first advocated by the NNPDF Collaboration, the photon PDF is treated on the same footing as the quark and gluon PDFs. Within this approach, the photon PDF is parametrised in a model-independent way using an artificial neural network and then constrained by LHC Drell-Yan measurements. This procedure was adopted in the NNPDF2.3/3.0QED determinations [36, 22, 37, 38]. The limited sensitivity of existing LHC data to PI contributions combined with the use of a flexible parametrisation resulted in large uncertainties on the photon distribution. A similar strategy was adopted in the recent analysis of Ref. [21] in which the ATLAS 8 TeV high-mass Drell-Yan data [39] was employed to constrain the photon PDF. Despite this dataset is particularly sensitive to the PI contribution, the resulting photon was still affected by large uncertainties.

Overcoming the limitations both two strategies, the LUXqed formalism presented in Refs. [40, 41] represented a breakthrough for the determination of the photon PDF. In this framework, both the elastic and the inelastic components of the photon PDF can be expressed in terms of the electromagnetic inclusive structure functions and from lepton-proton scattering by means of an exact QED calculation. This is very advantageous because these structure functions are known rather accurately both experimentally and theoretically. Accounting for the LUXqed constraints then leads to a reduction of the uncertainty of the photon PDF by more than an order of magnitude as compared to the NNPDF3.0QED data-driven determination.

Building upon the recent NNPDF3.1
fit [15], the goal of this paper is to perform a
global PDF analysis including QED corrections where the LUXqed
calculation is used to constrain the photon PDF.
The resulting PDF set, NNPDF3.1luxQED,
represents a state-of-the-art determination of the
partonic content of the proton including its photon
component.
The uncertainties on the photon PDF are now at the level
of a few percent, with
photons carrying up to of the total proton’s
momentum.^{1}^{1}1See also [42] for related studies in the MMHT framework.

We also take a first look at the phenomenological implications of NNPDF3.1luxQED for photon-initiated processes at the LHC. Previous studies based on NNPDF2.3/3.0QED indicated that PI contributions were potentially large, particularly at large invariant masses or transverse momenta, for processes such as Drell-Yan, pair, and top-quark pair production [43, 44, 45, 46, 47, 48]. Indeed, PI effects represented in some cases the dominant source of theoretical uncertainty. We now find that photon-initiated corrections computed with NNPDF3.1luxQED can be significant for many processes, leading to corrections of up to , and are fully consistent with previous estimates based on NNPDF3.0QED within uncertainties.

The outline of this paper is as follows. In Sect. 2 we present the settings of the global NNPDF3.1luxQED analysis in terms of input experimental data, theoretical calculations, and fitting strategy. Then in Sect. 3 we present the NNPDF3.1luxQED set, including a discussion of the momentum fraction of the proton carried by the photon, and in Sect. 4 we discuss some of its phenomenological implications for PI processes at the LHC. In Sect. 5 we summarise and discuss how our results and the code used to produce them are made publicly available. The full breakdown of the values in NNPDF3.1luxQED and its comparison with those in NNPDF3.1 is collected in Appendix A.

## 2 Fit settings

In this section we describe the fit configuration of the NNPDF3.1luxQED global analysis. We begin with a review of the input experimental dataset. We then discuss the theoretical framework, including a short summary of the relevant aspects of the LUXqed formalism along with the treatment of QED effects in the DGLAP evolution and the DIS structure functions. Finally, we present the strategy adopted to include the photon PDF in the global fit accounting for the LUXqed theoretical constraints.

### 2.1 Experimental data

The NNPDF3.1luxQED analysis is based on the same dataset as the recent NNPDF3.1 global fit [15]. This dataset includes fixed-target [49, 50, 51, 52, 53, 54, 55, 56] and HERA [19] inclusive DIS measurements; charm and bottom cross-sections from HERA [57]; fixed-target Drell-Yan production [58, 59, 60, 61]; Tevatron gauge boson and inclusive jet production [62, 63, 64, 65, 66]; along with electroweak boson production, inclusive jet, and cross-sections from ATLAS [67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81], CMS [82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93] and LHCb [94, 95, 96, 97, 98]. We refer to Ref. [15] for details about the implementation of each experiment.

For consistency, in this study we use exactly the same dataset as in NNPDF3.1, and in particular the same choice of kinematic cuts. Note that a number of those cuts were determined with the aim of minimising the potential effects from EW corrections and PI contributions. This choice implies that the kinematic regions more sensitive to PI effects are deliberately cut away. In addition, we do not include some recent measurements with known sensitivity to the photon PDF, such as the ATLAS high-mass Drell-Yan measurement at 8 TeV [39], since these were not part of the NNPDF3.1 dataset.

### 2.2 The LUXqed formalism

We briefly review the LUXqed formalism for the determination of the photon PDF, focusing on those features relevant to its implementation in a global analysis. For a comprehensive discussion we refer the reader to Refs. [40, 41]. In the LUXqed procedure, the photon PDF can be expressed in terms of the lepton-proton scattering inclusive structure functions and by means of an exact QED calculation as follows:

(2.1) |

where is the proton mass, is the factorisation scale, the running QED coupling, and the photon-quark splitting function. The lower integration limit in the integral is given by .

Note that the integral in in Eq. (2.1) extends up to . Therefore the LUXqed photon has an explicit dependence upon the elastic component of the structure functions, proportional to . This component can be expressed in terms of the electric and magnetic form factors and . In [40, 41], the elastic component of is determined using the form factors extracted from a fit to world data by the A1 collaboration [99] for GeV. The dipole model is then used to extrapolate the form factors to larger values of . A corresponding uncertainty due to the treatment of the large- extrapolation region is included in the evaluation of the photon PDF.

Furthermore, we observe that the integral in Eq. (2.1) requires an understanding of the structure functions down to potentially very low scales, well outside the region where perturbative QCD is applicable. Following the prescription of Refs. [40, 41], the integration in of the inelastic component () of and is achieved by combining parameterisations of experimental data with the perturbative computation in terms of PDFs where appropriate. Specifically, contributions to the inelastic structure functions come from two regions separated by . In the resonance region, defined GeV, the fit of the CLAS collaboration is used [100]. In order to assess the uncertainty due to this choice, the parametrisation of Ref. [101] is also considered.

The continuum region, defined as GeV, is further subdivided into two regions according the value of . For , with GeV, the GD11-P fit by HERMES [102, 103] is employed. For , structure functions are computed in terms of PDFs by means of their factorised expressions. The value of can be varied to estimate the uncertainty associated with this particular choice.

As a part of the present work, the LUXqed formalism has been implemented in an open-source public library, FiatLux [104], which has been used to produce the NNPDF3.1luxQED fits. The results obtained with FiatLux have been benchmarked with the original implementation used to produce the results of Refs. [40, 41], finding excellent agreement.

### 2.3 Theoretical calculations

The QCD calculations used in the present analysis are identical to those used in NNPDF3.1 [15]. DGLAP evolution and DIS structure functions are computed at NLO and NNLO accuracy in QCD using APFEL [26]. Heavy-quark mass effects in the structure functions are included using the FONLL general-mass variable-flavour-number scheme [105], and the charm PDF is fitted to data on an equal footing as the light quark and gluon PDFs [106, 107, 108]. Hadronic observables are computed at NLO using fast interpolation tables in the APPLgrid [109] and fastNLO [110] formats and combined to the DGLAP evolution kernels using APFELgrid [111]. NNLO corrections to hadronic processes are included by means of point-by-point NNLO/NLO -factors.

In the NNPDF3.1luxQED fit, QCD corrections are supplemented with QED effects. Concerning the evolution of PDFs, on top of the corrections, also the and splitting functions computed in Refs. [24, 25] are included in the DGLAP evolution equations. Additionally, corrections to the DIS coefficient functions are included. This introduces an additional (but mild) sensitivity to the photon PDF. The implementation in APFEL of the aforementioned QED corrections to the DIS structure functions and to the DGLAP evolution equations, together with the corresponding benchmarking, were presented in Ref. [21]. On the other hand, as in NNPDF3.1, pure weak corrections to hadronic observables are not accounted for.

Using NNPDF3.1luxQED as an input, we find that the cumulative effect on the photon-photon luminosity of the and corrections to the DGLAP splitting functions ranges between at low invariant masses and for high . Since these effects are larger than the typical uncertainties on the photon PDF determined through the LUXqed approach, it is important to take them into account. Concerning the DIS structure functions, using NNPDF3.1luxQED, one finds that the overall impact of the QED effects is at the permille level, except at large values of where they can be up to 2%.

In data-driven determinations of the photon PDF, it is in general necessary to include constraints from observables sensitive to PI processes. For instance, in the analysis of Ref. [21] PI contributions to the ATLAS 8 TeV high-mass Drell-Yan data [39] were computed via the aMCfast interface [112] to the Monte Carlo generator MadGraph5_aMC@NLO [113] and used to constrain the photon PDF. In the present study, the photon PDF is determined from a global analysis of hard-scattering data supplemented by the LUXqed theoretical constraint of Eq. (2.1). Given that the NNPDF3.1 dataset and the associated kinematic cuts were specifically designed to minimise the effects of PI contributions, one expects the impact of PI processes to hadronic cross-sections in NNPDF3.1luxQED to be minimal, and thus they are not included here. We have explicitly verified for some of the NNPDF3.1 datasets that this is an excellent approximation, see also the comparisons of Sect. 4. Nevertheless, the approach presented in this work is fully general, and future NNPDF analyses with QED corrections will include collider measurements characterised by sizeable PI contributions.

### 2.4 Fitting strategy

The determination of the NNPDF3.1luxQED set is performed by means of an iterative procedure. The starting point is a prior set of quark and gluon PDFs, in this case NNPDF3.1. From this PDF set, the high- inelastic component of the photon PDF is computed at GeV using Eq. (2.1), while for the other components the same inputs as in Ref. [41] are adopted. The resulting photon PDF is then evolved down to the parametrisation scale GeV and used as a fixed input in a refit of quark and gluon PDFs.

In this refit, the DGLAP evolution equations consistently include QED effects, the PI contribution to the DIS structure functions is taken into account, and the momentum sum rule reads

(2.2) |

This procedure, schematically illustrated in Fig. 2.1, is repeated until convergence is reached, in the sense that a stable photon PDF, and thus stable quark and gluon PDFs, are obtained. We consider stable the results of two consecutive iterations where the central value of the photon PDF varies by less than a fraction of its uncertainty.

As usual in the NNPDF approach, PDF uncertainties are represented by means of an ensemble of Monte Carlo replicas. Each replica is required to meet a set of quality criteria, discussed in Ref. [37], with fits failing these criteria being discarded. As the present study involves an iterative procedure, one must start with a sample of replicas large enough such that once all iterations have been completed a significant number of replicas still survives. To this end, here we use a prior sample of replicas.

Each replica will lead to different high- DIS structure functions and therefore, by virtue of Eq. (2.1), to a different photon PDF that is used as an external constraint in the following fit. The resulting quark and gluon PDFs are then used as an input to the following iteration of the fitting procedure, until convergence is reached. As the NNPDF3.1 dataset is (by construction) relatively insensitive to the photon PDF, the convergence is rapid and results are stable already after the second iteration. In future analyses, when hadronic measurements sensitive to PI contributions will be included, convergence is likely to be slower.

There are two main differences between our strategy and the direct application of Eq. (2.1) to NNPDF3.1. Firstly, the influence of the photon PDF in the DGLAP evolution equations and in the DIS structure functions is consistently taken into account during the fit of the quark and gluon PDFs. Secondly, the contribution of the photon PDF to the total momentum fraction is properly treated by imposing the momentum sum rule Eq. (2.2) during the fits. While these effects are likely to be small in this specific analysis, our framework is fully general and allows for the consistent inclusion of hadronic observables sensitive to the photon-initiated contributions.

In order to illustrate the convergence of the procedure, in Fig. 2.2 we show a comparison of the photon, gluon, up quark, and down quark PDFs at GeV between the first and the second iteration (labelled ITE1 and ITE2, respectively). We have verified that additional iterations leave the photon PDF unchanged, demonstrating that stability has been reached. For completeness, in Fig. 2.2 we also show the third and final iteration of the procedure (ITE3), where the additional LUXqed17 systematic variations are added to the photon PDF (see Sect. 2.5). As expected, these have the largest impact in the region , where the elastic contribution to the photon PDF is most important.

### 2.5 The uncertainties on the photon PDF

As mentioned above, the calculation of the photon PDF in terms of structure functions involves several contributions: the elastic component, the inelastic resonance component, and the inelastic low- and high- continuum components. Only the last component can be factorised in terms of PDFs and perturbative coefficient functions. Therefore, the ensemble of Monte Carlo replicas of the photon PDF accounts only for a part of the uncertainty, namely the one associated to the inelastic high- component. For a comprehensive estimate of the uncertainty one must also account for a number of additional sources of error.

The following sources of uncertainty are considered [41]: the elastic contribution from the A1 world proton form factor fits [99]; the parametrisation of the DIS structure functions in the resonance region [101, 100, 102]; the parametrisation of [103, 114, 115]; the scale at which low- and high- inelastic structure functions are matched; a twist-4 modification of the longitudinal structure function [116]; and finally an estimate of the missing higher-order corrections in the calculation of the DIS structure functions at high .

In the NNPDF3.1luxQED analysis, these uncertainties are introduced at the last iteration of the procedure. Once the quark and gluon PDFs from the -th iteration have been determined, they are used to construct . Then, for each photon PDF replica of this last iteration, , extra uncertainties are included as statistical fluctuations upon the photon PDF at GeV with correlations in , namely:

(2.3) |

where is an univariate Gaussian random number and is the normalised eigenvector for the -th systematic uncertainty in LUXqed17. Specifically, is obtained through the diagonalisation of the covariance matrix for the extra LUXqed17 uncertainties defined on a grid of points, using a similar method as that of Refs. [117, 118]. We have verified that this approach is numerically equivalent to using the corresponding Hessian eigenvectors of LUXqed17.

The photon PDF defined in Eq. (2.3) is finally used as an input for the -th fit iteration, to determine the final set of quark and gluons of NNPDF3.1luxQED.

## 3 The NNPDF3.1luxQED set

In this section we discuss the results of the iterative procedure outlined in Sect. 2.4: the NNPDF3.1luxQED fit. We focus on the NNLO case and comment where appropriate on any differences with respect to NLO. The overall fit quality of NNPDF3.1luxQED is at NLO and at NNLO. While there is some variation dataset-to-dataset, the global fit quality is identical to the corresponding NNPDF3.1 results at NNLO. See Appendix A for a full breakdown of the data description at NNLO and its comparison with NNPDF3.1.

### 3.1 The photon PDF

Here we compare our results for the photon PDF with those of the NNPDF3.0QED and LUXqed16/17 PDF sets. Comparisons with the latter are performed always at GeV, as the LUXqed16/17 sets are not defined below this scale. In the left panel of Fig. 3.1 we compare the NNPDF3.1luxQED photon PDF at GeV with the corresponding results from LUXqed16 and LUXqed17, normalised to the central value of the latter. The three determinations agree well across the full range, with central values always compatible within uncertainties. In addition, for the total uncertainties on the photon PDF from NNPDF3.1luxQED and LUXqed16/17 are identical. This feature is explained by the fact that in this region the uncertainties due to the elastic and low- inelastic structure functions dominate.

At medium- and small-, the NNPDF3.1luxQED photon exhibits somewhat smaller uncertainties. This is due to the use of a different set of quark and gluon PDFs determining the high- inelastic component, specifically NNPDF3.1 rather than the PDF4LHC15 set [119] used in LUXqed16/17. The contribution from the different error sources is further illustrated in the right panel of Fig. 3.1 where the same comparison including only the uncertainties due to the high- inelastic component is shown. The plot shows how at medium- to small- the contribution from the high- inelastic structure functions dominates the overall uncertainty.

In order to gauge the stability of the photon PDF with respect to the perturbative order of the QCD calculations used in the fit, in Fig. 3.2 we compare the photon PDFs from the NNPDF3.1luxQED NLO and NNLO fits, normalised to the central value of the former. The photon distributions are consistent within uncertainties, demonstrating good perturbative stability. Indeed, the shift due to the change in perturbative order is outside the PDF error bands only in the small- region, where the photon is sensitive to the prior PDF used for the computation of the high- inelastic component. In addition, we find that the photon PDF uncertainties are unaffected by the variation of the perturbative order.

In order to quantify the differences between photon PDFs determined from global analyses with and without imposing the LUXqed theoretical constraint, we compare NNPDF3.1luxQED with NNPDF3.0QED. In the following, the PDF uncertainties of NNPDF3.0QED are computed as 68% confidence-level (CL) intervals, with the central value taken to be the midpoint of the interval. In Fig. 3.3 we show the photon distributions from these two sets at GeV (left plot) and GeV (right plot). We find that both at low and high scales, in the region the two determinations agree within uncertainties. For instead, the NNPDF3.0QED photon undershoots NNPDF3.1luxQED by up to 40% at GeV. At high scales, PDF uncertainties in NNPDF3.0QED are at the level of a few percent at small but become as large as almost 100% at large . The uncertainties in NNPDF3.1luxQED are instead at the level of a few percent over the entire range in (see also Fig. 3.1).

It is also interesting to examine the relative size of the photon PDF with respect to the total quark singlet and gluon PDFs. This allows us to estimate where PI contribution to hadron-collider processes becomes sizeable as compared to quark- and gluon-initiated subprocesses. In Fig. 3.4 we compare the predictions of NNPDF3.0QED and NNPDF3.1luxQED for the (left) and (right) ratios at GeV. From the right panel of Fig. 3.4 we observe that the ratio is around over the entire range of . On the other hand, from the right panel of Fig. 3.4 we find that for the ratio becomes larger than the one. In fact, the ratio is as large as 10% for .

Therefore, we find that that complementing a data-driven determination of the photon PDF with the LUXqed theoretical constraints allows for a precise determination of the photon PDF in the entire range of relevant for applications at the LHC.

### 3.2 QED effects on the quark and gluon PDFs

In this section we study the quark and gluon PDFs in NNPDF3.1luxQED as compared to their corresponding QCD-only counterparts in NNPDF3.1. This comparison gauges the impact on quarks and gluons of three different QED effects: the modification of the momentum sum rule, the QED splitting functions in the DGLAP evolution equations, and the QED corrections to the DIS coefficient functions.

In Fig. 3.5 we show the singlet and gluon PDFs of the NNPDF3.1 and NNPDF3.1luxQED sets at GeV normalised to the central value of the former. While differences at the level of the singlet are small, differences for the gluon PDF are somewhat larger. Indeed, the NNPDF3.1luxQED gluon is smaller than its QCD counterpart by about 1% at and enhanced by about for . In both cases, the shift in the central values is at the edge of the corresponding PDF uncertainty band. The effect on the gluon PDF can be explained by observing that, as we will discuss in Sect. 3.4, the photon PDF can carry up to 0.5% of the proton momentum. This fraction is effectively subtracted from the singlet and gluon distributions by means of the sum rule, Eq. (2.2). However, the sum rule mostly affects the gluon PDF because the normalisation of the quark singlet is more tightly constrained from the DIS inclusive structure function data. We conclude that the back-reaction of QED effects onto the quark and gluon PDFs is small but not negligible, particularly for the latter.

For completeness, in Fig. 3.6 we show the same comparison as in Fig. 3.5 but now between NNPDF3.1luxQED and LUXqed17. Note that the quark and gluon PDFs of LUXqed17 correspond closely to those of the PDF4LHC15 set, differing only by a rescaling of the gluon PDF and by the QED contributions to the DGLAP evolution. Since the PDF4LHC15 set is built as a combination of three different PDF sets, namely CT14, MMHT14, and NNPDF3.0, it exhibits larger uncertainties that the individual sets. We find good compatibility between the singlet and gluon of NNPDF3.1luxQED and LUXqed17, with the PDF errors of the former being rather smaller. These reduced uncertainties are particularly noticeable for the medium and large- gluon PDF, due to the several gluon-sensitive experiments included in NNPDF3.1 such as top-quark pair distributions [120] and the boson [121].

### 3.3 Partonic luminosities

Next we compare partonic luminosities integrated over rapidity as a function of the final-state invariant mass for photon-photon and photon-quark initial states (see [44] for the definitions used). In Fig. 3.7 we compare the luminosity obtained with the NNPDF3.0QED, NNPDF3.1luxQED, and LUXqed17 sets for a centre-of-mass energy of TeV. From the left panel of Fig. 3.7 we observe good agreement between NNPDF3.0QED and NNPDF3.1luxQED. The two sets agree within uncertainties over the entire mass range considered, except for at GeV where NNPDF3.1luxQED overshoots NNPDF3.0QED. Evident once again is the effect of the LUXqed constraints on the photon PDF uncertainty, with the errors of just a few percent as compared to the determination that does not account for them. From Fig. 3.7 we also see that there is good agreement between LUXqed17 and NNPDF3.1luxQED both at the level of central values and ofuncertainties. Only at TeV NNPDF3.1luxQED tends to be a few percent larger than LUXqed17. Similar considerations hold for the luminosities.

Following the PDF-level comparisons presented in Sect. 3.1, it is instructive to also consider the ratios of the photon-photon luminosity over gluon-gluon and over quark-antiquark luminosities, and . These ratios are interesting since they provide an estimate of the relative importance of the PI contribution over quark- and gluon-initiated contributions as a function of . These ratios are shown in Fig. 3.8, where we compare the results from NNPDF3.0QED, NNPDF3.1luxQED, and LUXqed17. The ratio varies very mildly with , with a value . On the other hand, the ratio increases steeply with the final state invariant mass , beginning at at low invariant masses and growing up to for TeV. For large invariant masses, we therefore find that the two ratios and take similar values. For both cases, the results of the three PDF sets are consistent within uncertainties, with NNPDF3.1luxQED and LUXqed17 exhibiting significantly reduced errors as compared to NNPDF3.0QED.

### 3.4 The momentum fraction carried by photons

As the photon PDF carries a non-zero amount of the total proton momentum, it therefore contributes to the momentum sum rule, Eq. (2.2). Here we examine the momentum fraction carried by the photon PDF

(3.1) |

In Fig. 3.9 we show the value of in the NNPDF3.1luxQED, NNPDF3.0QED, and LUXqed17 sets as a function of the scale . In the right plot of Fig. 3.9 we also show the corresponding percent PDF uncertainties. For LUXqed17 we restrict the comparison to the validity region of this set, i.e. GeV.

From Fig. 3.9 we observe that, while the NNPDF3.0QED determination is affected by large PDF uncertainties, the other two sets lead to a compatible prediction for the photon momentum fraction for all values of . As expected from the PDF-level comparisons, there is a significant reduction in the uncertainty on the value of once the LUXqed theoretical constraints are accounted for. Indeed, while in NNPDF3.0QED the uncertainties in the photon momentum fraction range from around 50% to 100%, in NNPDF3.1luxQED the contribution of the photon PDF to the momentum of the proton is known with an accuracy better than 1% over the entire range in . Nevertheless, the central value of in NNPDF3.0QED turns out to be rather close to that of NNPDF3.1luxQED, highlighting the consistency between the two approaches.

NNPDF3.0QED | % | % |
---|---|---|

NNPDF3.1luxQED | % | % |

LUXqed17 |

In Tab. 3.1 we report the photon momentum fraction Eq. (3.1) both at the initial parametrisation scale GeV and at for the three PDF sets including the associated uncertainties. While in NNPDF3.0QED the photon momentum fraction at the initial scale is consistent with zero, in NNPDF3.1luxQED one finds a non-zero photon momentum fraction with very high statistical significance. In particular, the photon momentum fraction in NNPDF3.1luxQED increases from at low scales to at , with small uncertainties in both cases. For , the results of NNPDF3.1luxQED are fully consistent with those of LUXqed17, as also shown in Fig. 3.9.

## 4 Photon-initiated processes at the LHC

We shall now explore some of the implications of NNPDF3.1luxQED for LHC phenomenology. Specifically, we shall investigate the application of this new set to the study of Drell-Yan, vector-boson pair production, top-quark pair production, and the associated production of a Higgs boson with a boson. Representative PI diagrams contributing to these processes at the Born level are shown in Fig. 4.1. Our aim is to assess the relative size of the PI contributions with respect to quark- and gluon-initiated subprocesses at the TeV LHC. See also Refs. [43, 23, 44, 33, 45, 46, 47] for recent studies.

The results presented in this section have been obtained at LO in both the QCD and QED couplings using MadGraph5_aMC@NLO interfaced to APPLgrid through aMCfast. We will compare the predictions of NNPDF3.1luxQED to those of NNPDF3.0QED and LUXqed17. In all cases we will use the NNLO PDF sets, though the photon PDF depends only mildly on the perturbative order (see Fig. 3.2). PDF uncertainties for the NNPDF sets are defined as the 68% confidence level interval and the central value as the midpoint of this range. This is particularly relevant for NNPDF3.0QED for which, due to non-Gaussianity in the replica distribution, PDF errors computed as standard deviations can differ by up to one order of magnitude as compared to the 68% CL intervals.

### 4.1 Drell-Yan production

We begin by examining the role of PI contributions in neutral-current Drell-Yan
production.
We will study this process in three different kinematic regions of the
outgoing lepton pair: around the peak, at low invariant masses,
and at high invariant masses.
We start with the peak region, defined as
GeV, where is the lepton-pair invariant mass, and focus
on the central rapidity region , relevant for
ATLAS and CMS.^{2}^{2}2We have verified that similar results
hold for the forward rapidity region, ,
relevant for LHCb.
This region provides the bulk of the Drell-Yan measurements included
in modern PDF fits and therefore assessing the impact of PI
contributions is particularly important here.

In Fig. 4.2 we show the ratio of the PI contributions to the corresponding quark- and gluon-initiated contributions for Drell-Yan production as a function of at TeV in the peak region. We compare the predictions of NNPDF3.0QED, LUXqed17, and NNPDF3.1luxQED, with the PI contributions normalised to the central value of NNPDF3.1luxQED. For reference we also show the value of the PDF uncertainties in NNPDF3.1luxQED.

We find that PI effects for this process are at the permille level for but they become larger as we move away from the peak, reaching 3% for GeV. At the lower edge of the region the contribution of the PI channel exceeds the level of PDF uncertainty, highlighting the sensitivity of this distribution to the photon PDF. We find that NNPDF3.1luxQED and LUXqed17 lead to a larger PI contribution as compared to NNPDF3.0QED at low . As the PI contribution is only significant away from the -peak, where the bulk of the cross-section lies, these effects may be reasonably neglected in the integrated cross-sections.

Fig. 4.2 demonstrates that the PI contributions in NNPDF3.1luxQED and LUXqed17 lead to very similar results for Drell-Yan production around the peak. We have verified that this similarity holds also for the low and high mass kinematic regions, as well as for the rest of processes studied in this section. In the following discussion we will therefore restrict ourselves to comparisons between NNPDF3.0 and NNPDF3.1luxQED.

We now move to study the low- and high-mass regions, defined as GeV and GeV respectively. Drell-Yan low-mass measurements have been presented by ATLAS, CMS, and LHCb [80, 122, 123], with the two-fold motivation of providing input for PDF fits and to study QCD in complementary kinematic regimes. The high-mass region is relevant for BSM searches that exploit lepton-pair final states, such as those expected in the presences of new heavy gauge bosons or [124, 125].

In Fig. 4.3 we show the same comparison as in the right panel of Fig. 4.2 for the low- and high-mass regions. In the low-mass case, the PI effects are more significant than in the -peak region, being between 3% and 4% for most of the range, consistently larger that the PDF uncertainty. We find that PI effects in NNPDF3.1luxQED can be up to a factor three larger than in the NNPDF3.0QED case due to the corresponding differences in the photon PDF at small . In the case of the high-mass region, we observe that the effect of the PI contribution computed with NNPDF3.1luxQED is comparable to the PDF uncertainties for TeV, eventually becoming as large as of the QCD cross-section. These effects are markedly smaller than in NNPDF3.0QED, where shifts in the cross-section up to due to PI contributions were allowed within uncertainties.

To conclude this discussion on Drell-Yan at the LHC, we have evaluated the ratio of the LO PI contributions to the NLO QCD cross-sections for the kinematics of the ATLAS high-mass Drell-Yan measurements at 8 TeV [39]. Both the Bayesian reweighting study of the ATLAS paper [39] and the analysis of Ref. [21] indicate that this dataset has a considerable sensitivity to PI contributions if NNPDF3.0QED is used as a prior. Here we revisit this process to assess how the picture changes when using NNPDF3.1luxQED.

In Fig. 4.4 we show the ratio of PI over QCD contributions for the lepton-pair rapidity distributions in Drell-Yan at 8 TeV for two invariant mass bins, GeV GeV and GeV GeV. As can be seen, with NNPDF3.0QED the effects of the PI contribution at large invariant masses can be as large as 25% of the QCD cross-section. This shift is larger than the corresponding experimental uncertainties, which are typically at the percent level, explaining the sensitivity of NNPDF3.0QED to this dataset. From Fig. 4.4 we observe that the PI contribution becomes smaller when using NNPDF3.1luxQED, though its effect is still comparable to the experimental uncertainties. Indeed, the PI contribution in the highest invariant mass bin ranges between 6% of the QCD cross-section in the central rapidity region and 1% for . This comparison confirms that the PI contribution is important for the quantitative description of the Drell-Yan process above the peak.

### 4.2 Vector-boson pair production

The production of vector-boson pairs is of particular interest for the LHC physics program. Firstly, they probe the electroweak sector of the SM and provide bounds on possible anomalous couplings. Secondly, this final state appears in several BSM scenarios and therefore many searches involve the detection of pairs of weak vector bosons. When the invariant mass of the vector-boson pair is large, the PI contributions, arising already at the Born level (see Fig. 4.1), are known to be significant [126]. Here we will examine the case of opposite-sign production at the LHC 13 TeV.

To assess the size of the PI contribution to this process, in Fig. 4.5 we show the invariant mass distribution and the transverse momentum distributions of the bosons . The bosons are taken to be stable and required to be in the central rapidity region, . From the comparisons shown in Fig. 4.5 we find that for the distribution the PI contribution is larger than the PDF uncertainties over the entire range considered. In particular, when using NNPDF3.1luxQED, the size of the PI contribution with respect to the total cross-section increases from 1% at GeV up to 35% at TeV. As expected, the trend is similar using NNPDF3.0QED but with much larger uncertainties. An upwards shift of the cross-section by a factor two or larger would be allowed within PDF uncertainties in this case.

The picture is rather different for the case of the transverse momentum distribution of the individual bosons . Here we find that the PI contribution using NNPDF3.1luxQED are small, around the level, over the entire range in considered. Additionally, the effect of using NNPDF3.0QED instead is not so dramatic, with an increase in the cross-section of a few percent at most. The differences between the two distributions arise from the fact that the PI contribution to boson pair production is kinematically enhanced only in the large limit, irrespective of the value of . The results of Fig. 4.5 suggest that current measurements of this process from ATLAS and CMS [127, 128] might already be sensitive to the photon PDF.

### 4.3 Top-quark pair differential distributions

Next we turn to study the impact of the PI contributions on differential distributions in top-quark pair production (see also Refs. [5, 47]). In Fig. 4.6 we show the ratio of the PI contribution over the QCD cross-section for the invariant mass distribution of top-quark pairs and the transverse momentum of the single top quarks at the 13 TeV LHC. Unlike in the case of high-mass Drell-Yan production, we find that the PI contribution to top-quark pair production is negligible even for the highest values of and accessible at the LHC. Indeed, in the case of NNPDF3.1luxQED the size of the PI contribution is at the permille level at most. Therefore, in theoretical calculations of top-quark pair production with electroweak corrections, the PI contribution can be safely neglected.

From the comparisons in Fig. 4.6 we also see that the PI correction is somewhat larger in NNPDF3.0QED but with larger associated uncertainties. Even in this case, at the highest invariant masses the upper edge of the 68% CL interval indicates that corrections due to the PI contribution are at most 0.1%. We have also verified that the PI contribution to production is phenomenologically negligible also for other distributions such as the rapidity distribution of top quarks and top-antitop pairs, and , respectively.

### 4.4 Higgs production in association with a vector bosons

The last process that we consider in this section is Higgs production in association with a vector boson, see Fig. 4.1. PI corrections to this process are known to be significant. In the Higgs Cross Section Working Group prediction for the total cross-sections, where NNPDF3.0QED is used, the uncertainty due to the PI contribution is the dominant source of theory error [129]. To investigate this, we have generated and then combined exclusive samples for and with and without the PI contribution. Both the Higgs boson and the boson are required to be in the central rapidity region, and . No other kinematic cuts are applied.

In Fig. 4.7 we show the same comparison as in Fig. 4.2 for the Higgs transverse momentum and rapidity distributions. In the case of the distribution, we find that PI effects can be up to when using NNPDF3.1luxQED, with the largest effects localised at intermediate values of GeV. We also note that the shift induced by the PI contribution is bigger than the PDF uncertainties. Concerning the rapidity distribution, the PI contribution can be in the central rapidity region when using NNPDF3.1luxQED, while it becomes smaller as one moves to the forward region.

The comparisons of Fig. 4.7 illustrate that PI contributions are relevant also for Higgs boson physics, including the measurements of its couplings and branching fractions.

## 5 Summary

Parton distributions with QED effects and a photon PDF are an essential component in high-precision calculations of many LHC processes. Previous NNPDF QED sets adopted a data-driven strategy to determine the photon PDF, independently parametrising and then fitting it using constraints from Drell-Yan measurements at the LHC. While this strategy minimised the theoretical bias due to model assumptions, the lack of a precise experimental handle to constrain the photon PDF led to large uncertainties.

With the development of the LUXqed framework, it is now possible to constrain the photon PDF in terms of the accurately known inclusive structure functions in lepton-hadron scattering. In this work we have presented the NNPDF3.1luxQED set, where the photon content of the proton is determined by means of a global PDF analysis supplemented by the LUXqed theoretical constraint. As a result, the uncertainty upon the photon PDF is considerably reduced as compared to our previous NNPDF3.0QED determination, down now to the level of a few percent. We find that photons carry up to 0.5% of the total momentum of the proton, and that the overall impact of the various types of QED effects included in NNPDF3.1luxQED induce small but non-negligible modifications in the quark and gluon PDFs.

We have then presented a first exploration of the implications of NNPDF3.1luxQED for photon-initiated processes at the LHC. We determine that the impact of PI contributions is consistent within uncertainties with respect to previous estimates based on NNPDF3.0QED, and that they can be significant for many processes. For instance, we find corrections up to for high-mass Drell-Yan and up to for production. In many cases, PI processes can be either comparable with or larger than PDF uncertainties. The uncertainty associated with these PI effects is in itself at the level of a few percent, so their overall effect is a shift of the cross-sections as compared to the QCD-only calculation.

The NNPDF3.1luxQED set represents a state-of-the-art determination of the PDFs of the proton including its photon component, accounting for all relevant theoretical and experimental constraints. This set is therefore well suited for precision calculations of LHC processes. The NNPDF3.1luxQED sets are available via the LHAPDF6 interface [130]:

NNPDF31_nlo_as_0118_luxqed

NNPDF31_nnlo_as_0118_luxqed

while the FiatLux library, an open-source C++ implementation of the LUXqed formalism, can be obtained from:

https://github.com/scarrazza/fiatlux

together with the corresponding documentation.

### Acknowledgments

We are grateful to our colleagues in the NNPDF Collaboration for many illuminating discussions on the photon PDF and fits with QED corrections, and specially to Stefano Forte for comments on this manuscript. We are grateful to Aneesh Manohar, Paolo Nason, Gavin Salam, and Giulia Zanderighi for discussions about the LUXqed formalism. J. R. is grateful to Lucian Harland-Lang for discussions about PDF fits with QED effects.

V. B., N. H., and J. R. are supported by an European Research Council Starting Grant “PDF4BSM”. J. R. is also partially supported by the Netherlands Organization for Scientific Research (NWO). S. C. is supported by the HICCUP ERC Consolidator grant (614577) and by the European Research Council under the European Union’s Horizon 2020 research and innovation Programme (grant agreement n 740006).

## Appendix A NNPDF3.1luxQED fit quality

In this Appendix we collect the results of the fit quality in the NNPDF3.1luxQED analysis and compare it to those from its QCD-only counterpart, NNPDF3.1. As customary in NNPDF analyses, these values are computed using the experimental definition of the covariance matrix, while the definition [131] was instead used during the fits in order to avoid the D’Agostini bias.

In Table A.1 we list the values of from the NNPDF3.1 and NNPDF3.1luxQED NNLO fits for all the experiments included in the global analysis. We see that at the total dataset level, the fit quality is essentially identical in the two cases, yielding . Also at the level of individual experiments there is good agreement between the two sets, with some small differences that are consistent with statistical fluctuations. Recall that these two fits are statistically independent, and therefore one expects fluctuations of the order in the values of individual experiments.

NNPDF3.1 | NNPDF3.1luxQED | |

NMC | 1.30 | 1.31 |

SLAC | 0.75 | 0.71 |

BCDMS | 1.21 | 1.21 |

CHORUS | 1.11 | 1.11 |

NuTeV dimuon | 0.82 | 0.75 |

HERA I+II incl. | 1.16 | 1.16 |

HERA | 1.45 | 1.48 |

HERA | 1.11 | 1.11 |

DY E866 | 0.41 | 0.39 |

DY E886 | 1.43 | 1.43 |

DY E605 | 1.21 | 1.20 |

CDF rap | 1.48 | 1.48 |

CDF Run II jets | 0.87 | 0.88 |

D0 rap | 0.60 | 0.60 |

D0 asy | 2.70 | 2.68 |

D0 asy | 1.56 | 1.57 |

ATLAS total | 1.09 | 1.07 |

ATLAS 7 TeV 2010 | 0.96 | 0.96 |

ATLAS HM DY 7 TeV | 1.54 | 1.57 |

ATLAS low-mass DY 7 TeV | 0.90 | 0.88 |

ATLAS 7 TeV 2011 | 2.14 | 2.18 |

ATLAS jets 2010 7 TeV | 0.94 | 0.91 |

ATLAS jets 2.76 TeV | 1.03 | 1.02 |

ATLAS jets 2011 7 TeV | 1.07 | 1.06 |

ATLAS 8 TeV | 0.93 | 0.93 |

ATLAS 8 TeV | 0.94 | 0.91 |

ATLAS | 0.86 | 0.89 |

ATLAS rap | 1.45 | 1.22 |

CMS total | 1.06 | 1.05 |

CMS W asy 840 pb | 0.78 | 0.78 |

CMS W asy 4.7 fb | 1.75 | 1.74 |

CMS Drell-Yan 2D 2011 | 1.23 | 1.27 |

CMS W rap 8 TeV | 1.01 | 1.00 |

CMS jets 7 TeV 2011 | 0.84 | 0.81 |

CMS jets 2.76 TeV | 1.03 | 1.02 |

CMS 8 TeV | 1.32 | 1.32 |

CMS | 0.20 | 0.16 |

CMS rap | 0.94 | 0.99 |

LHCb total | 1.47 | 1.47 |

LHCb 940 pb | 1.49 | 1.48 |

LHCb 2 fb | 1.14 | 1.11 |

LHCb 7 TeV | 1.76 | 1.79 |

LHCb 8 TeV | 1.37 | 1.36 |

Total | 1.148 | 1.148 |

From the comparison shown in Table A.1 we can conclude that while QED effects lead to small differences at the quark and gluon PDF level (see Fig. 3.5), the fit quality is still very similar to the QCD-only fit. This implies that these small QED effects can be reabsorbed into the PDFs, leading to the same quantitative description of the datasets included in the analysis. Note that the situation is likely to be rather different once measurements directly sensitive to the photon content of the proton, such as the ATLAS high-mass Drell-Yan at 8 TeV, are included into the fit.

## References

- [1] G. Heinrich, QCD calculations for the LHC: status and prospects, in 5th Large Hadron Collider Physics Conference (LHCP 2017) Shanghai, China, May 15-20, 2017, 2017. arXiv:1710.04998.
- [2] S. Dittmaier, A. Huss, and C. Speckner, Weak radiative corrections to dijet production at hadron colliders, JHEP 1211 (2012) 095, [arXiv:1210.0438].
- [3] T. Becher and X. Garcia i Tormo, Electroweak Sudakov effects in and production at large transverse momentum, Phys. Rev. D88 (2013), no. 1 013009, [arXiv:1305.4202].
- [4] J. M. Campbell, R. K. Ellis, and C. Williams, Direct photon production at next-to-next-to-leading order, arXiv:1612.04333.
- [5] M. Czakon, D. Heymes, A. Mitov, D. Pagani, I. Tsinikos, and M. Zaro, Top-pair production at the LHC through NNLO QCD and NLO EW, arXiv:1705.04105.
- [6] C. Carloni Calame, G. Montagna, O. Nicrosini, and A. Vicini, Precision electroweak calculation of the production of a high transverse-momentum lepton pair at hadron colliders, JHEP 0710 (2007) 109, [arXiv:0710.1722].
- [7] Y. Li and F. Petriello, Combining QCD and electroweak corrections to dilepton production in FEWZ, Phys.Rev. D86 (2012) 094034, [arXiv:1208.5967].
- [8] R. Gavin, Y. Li, F. Petriello, and S. Quackenbush, W Physics at the LHC with FEWZ 2.1, Comput.Phys.Commun. 184 (2013) 208–214, [arXiv:1201.5896].
- [9] K. Mishra et al., Electroweak Corrections at High Energies, in Proceedings, 2013 Community Summer Study on the Future of U.S. Particle Physics: Snowmass on the Mississippi (CSS2013): Minneapolis, MN, USA, July 29-August 6, 2013, 2013. arXiv:1308.1430.
- [10] J. M. Campbell, D. Wackeroth, and J. Zhou, Study of weak corrections to Drell-Yan, top-quark pair, and dijet production at high energies with MCFM, Phys. Rev. D94 (2016), no. 9 093009, [arXiv:1608.03356].
- [11] S. Frixione, V. Hirschi, D. Pagani, H. S. Shao, and M. Zaro, Electroweak and QCD corrections to top-pair hadroproduction in association with heavy bosons, JHEP 06 (2015) 184, [arXiv:1504.03446].
- [12] S. Kallweit, J. M. Lindert, P. MaierhÃ¶fer, S. Pozzorini, and M. SchÃ¶nherr, NLO electroweak automation and precise predictions for W+multijet production at the LHC, JHEP 04 (2015) 012, [arXiv:1412.5157].
- [13] B. Biedermann, S. BrÃ¤uer, A. Denner, M. Pellen, S. Schumann, and J. M. Thompson, Automation of NLO QCD and EW corrections with Sherpa and Recola, Eur. Phys. J. C77 (2017) 492, [arXiv:1704.05783].
- [14] J. Gao, L. Harland-Lang, and J. Rojo, The Structure of the Proton in the LHC Precision Era, arXiv:1709.04922.
- [15] NNPDF Collaboration, R. D. Ball et al., Parton distributions from high-precision collider data, Eur. Phys. J. C77 (2017), no. 10 663, [arXiv:1706.00428].
- [16] S. Dulat, T.-J. Hou, J. Gao, M. Guzzi, J. Huston, P. Nadolsky, J. Pumplin, C. Schmidt, D. Stump, and C. P. Yuan, New parton distribution functions from a global analysis of quantum chromodynamics, Phys. Rev. D93 (2016), no. 3 033006, [arXiv:1506.07443].
- [17] L. A. Harland-Lang, A. D. Martin, P. Motylinski, and R. S. Thorne, Parton distributions in the LHC era: MMHT 2014 PDFs, Eur. Phys. J. C75 (2015) 204, [arXiv:1412.3989].
- [18] S. Alekhin, J. Blumlein, S. Moch, and R. Placakyte, Parton Distribution Functions, and Heavy-Quark Masses for LHC Run II, arXiv:1701.05838.
- [19] ZEUS, H1 Collaboration, H. Abramowicz et al., Combination of measurements of inclusive deep inelastic scattering cross sections and QCD analysis of HERA data, Eur. Phys. J. C75 (2015), no. 12 580, [arXiv:1506.06042].
- [20] A. D. Martin, R. G. Roberts, W. J. Stirling, and R. S. Thorne, Parton distributions incorporating QED contributions, Eur. Phys. J. C39 (2005) 155, [hep-ph/0411040].
- [21] xFitter Developers’ Team Collaboration, F. Giuli et al., The photon PDF from high-mass DrellâYan data at the LHC, Eur. Phys. J. C77 (2017), no. 6 400, [arXiv:1701.08553].
- [22] NNPDF Collaboration, R. D. Ball et al., Parton distributions with QED corrections, Nucl.Phys. B877 (2013) 290–320, [arXiv:1308.0598].
- [23] C. Schmidt, J. Pumplin, D. Stump, and C. P. Yuan, CT14QED parton distribution functions from isolated photon production in deep inelastic scattering, Phys. Rev. D93 (2016), no. 11 114015, [arXiv:1509.02905].
- [24] D. de Florian, G. F. R. Sborlini, and G. Rodrigo, QED corrections to the AltarelliâParisi splitting functions, Eur. Phys. J. C76 (2016), no. 5 282, [arXiv:1512.00612].
- [25] D. de Florian, G. F. R. Sborlini, and G. Rodrigo, Two-loop QED corrections to the Altarelli-Parisi splitting functions, JHEP 10 (2016) 056, [arXiv:1606.02887].
- [26] V. Bertone, S. Carrazza, and J. Rojo, APFEL: A PDF Evolution Library with QED corrections, Comput.Phys.Commun. 185 (2014) 1647, [arXiv:1310.1394].
- [27] G. P. Salam and J. Rojo, A Higher Order Perturbative Parton Evolution Toolkit (HOPPET), Comput. Phys. Commun. 180 (2009) 120–156, [arXiv:0804.3755].
- [28] R. Sadykov, Impact of QED radiative corrections on Parton Distribution Functions. , arXiv:1401.1133.
- [29] V. M. Budnev, I. F. Ginzburg, G. V. Meledin, and V. G. Serbo, The Two photon particle production mechanism. Physical problems. Applications. Equivalent photon approximation, Phys. Rept. 15 (1975) 181–281.
- [30] A. De Rujula and W. Vogelsang, On the photon constituency of protons, Phys.Lett. B451 (1999) 437, [hep-ph/9812231].
- [31] M. Gluck, C. Pisano, and E. Reya, The Polarized and unpolarized photon content of the nucleon, Phys. Lett. B540 (2002) 75–80, [hep-ph/0206126].
- [32] A. Martin and M. Ryskin, The photon PDF of the proton, Eur.Phys.J. C74 (2014) 3040, [arXiv:1406.2118].
- [33] L. A. Harland-Lang, V. A. Khoze, and M. G. Ryskin, The production of a diphoton resonance via photon-photon fusion, JHEP 03 (2016) 182, [arXiv:1601.07187].
- [34] L. A. Harland-Lang, V. A. Khoze, and M. G. Ryskin, The photon PDF in events with rapidity gaps, Eur. Phys. J. C76 (2016), no. 5 255, [arXiv:1601.03772].
- [35] ZEUS Collaboration, S. Chekanov et al., Measurement of isolated photon production in deep inelastic ep scattering, Phys.Lett. B687 (2010) 16–25, [arXiv:0909.4223].
- [36] R. D. Ball et al., Parton distributions with LHC data, Nucl.Phys. B867 (2013) 244, [arXiv:1207.1303].
- [37] NNPDF Collaboration, R. D. Ball et al., Parton distributions for the LHC Run II, JHEP 04 (2015) 040, [arXiv:1410.8849].
- [38] V. Bertone and S. Carrazza, Combining NNPDF3.0 and NNPDF2.3QED through the APFEL evolution code, PoS DIS2016 (2016) 031, [arXiv:1606.07130].
- [39] ATLAS Collaboration, G. Aad et al., Measurement of the double-differential high-mass Drell-Yan cross section in pp collisions at TeV with the ATLAS detector, JHEP 08 (2016) 009, [arXiv:1606.01736].
- [40] A. Manohar, P. Nason, G. P. Salam, and G. Zanderighi, How bright is the proton? A precise determination of the photon parton distribution function, Phys. Rev. Lett. 117 (2016), no. 24 242002, [arXiv:1607.04266].
- [41] A. V. Manohar, P. Nason, G. P. Salam, and G. Zanderighi, The Photon Content of the Proton. , arXiv:1708.01256.
- [42] L. A. Harland-Lang, R. Nathvani, R. S. Thorne, and A. D. Martin, MMHT PDFs: updates and outlook, Acta Phys. Polon. B48 (2017) 1011–1024, [arXiv:1704.00162].
- [43] V. Bertone, S. Carrazza, D. Pagani, and M. Zaro, On the Impact of Lepton PDFs, JHEP 11 (2015) 194, [arXiv:1508.07002].
- [44] M. L. Mangano et al., Physics at a 100 TeV pp collider: Standard Model processes, arXiv:1607.01831.
- [45] D. Bourilkov, Photon-induced Background for Dilepton Searches and Measurements in pp Collisions at 13 TeV. , arXiv:1606.00523.
- [46] E. Accomando, J. Fiaschi, F. Hautmann, S. Moretti, and C. H. Shepherd-Themistocleous, Photon-initiated production of a dilepton final state at the LHC: Cross section versus forward-backward asymmetry studies, Phys. Rev. D95 (2017), no. 3 035014, [arXiv:1606.06646].
- [47] D. Pagani, I. Tsinikos, and M. Zaro, The impact of the photon PDF and electroweak corrections on distributions, Eur. Phys. J. C76 (2016), no. 9 479, [arXiv:1606.01915].
- [48] M. Czakon, D. Heymes, A. Mitov, D. Pagani, I. Tsinikos, and M. Zaro, Top-quark pair production at NNLO QCD + NLO EW accuracy: Tevatron results, 2017. arXiv:1712.04842.
- [49] New Muon Collaboration, M. Arneodo et al., Accurate measurement of and , Nucl. Phys. B487 (1997) 3–26, [hep-ex/9611022].
- [50] New Muon Collaboration, M. Arneodo et al., Measurement of the proton and deuteron structure functions, and , and of the ratio , Nucl. Phys. B483 (1997) 3–43, [hep-ph/9610231].
- [51] BCDMS Collaboration, A. C. Benvenuti et al., A high statistics measurement of the proton structure functions and from deep inelastic muon scattering at high , Phys. Lett. B223 (1989) 485.
- [52] BCDMS Collaboration, A. C. Benvenuti et al., A high statistics measurement of the deuteron structure functions and from deep inelastic muon scattering at high , Phys. Lett. B237 (1990) 592.
- [53] L. W. Whitlow, E. M. Riordan, S. Dasu, S. Rock, and A. Bodek, Precise measurements of the proton and deuteron structure functions from a global analysis of the SLAC deep inelastic electron scattering cross-sections, Phys. Lett. B282 (1992) 475–482.
- [54] CHORUS Collaboration, G. Onengut et al., Measurement of nucleon structure functions in neutrino scattering, Phys. Lett. B632 (2006) 65–75.
- [55] NuTeV Collaboration, M. Goncharov et al., Precise measurement of dimuon production cross-sections in Fe and Fe deep inelastic scattering at the Tevatron, Phys. Rev. D64 (2001) 112006, [hep-ex/0102049].
- [56] D. A. Mason, Measurement of the strange - antistrange asymmetry at NLO in QCD from NuTeV dimuon data, . FERMILAB-THESIS-2006-01.
- [57] H1 , ZEUS Collaboration, H. Abramowicz et al., Combination and QCD Analysis of Charm Production Cross Section Measurements in Deep-Inelastic ep Scattering at HERA, Eur.Phys.J. C73 (2013) 2311, [arXiv:1211.1182].
- [58] NuSea Collaboration, J. C. Webb et al., Absolute Drell-Yan dimuon cross sections in 800-GeV/c p p and p d collisions, hep-ex/0302019.
- [59] J. C. Webb, Measurement of continuum dimuon production in 800-GeV/c proton nucleon collisions, hep-ex/0301031.
- [60] FNAL E866/NuSea Collaboration, R. S. Towell et al., Improved measurement of the anti-d/anti-u asymmetry in the nucleon sea, Phys. Rev. D64 (2001) 052002, [hep-ex/0103030].
- [61] G. Moreno et al., Dimuon production in proton - copper collisions at = 38.8-GeV, Phys. Rev. D43 (1991) 2815–2836.
- [62] CDF Collaboration, T. A. Aaltonen et al., Measurement of of Drell-Yan pairs in the Mass Region from Collisions at TeV, Phys. Lett. B692 (2010) 232–239, [arXiv:0908.3914].
- [63] D0 Collaboration, V. M. Abazov et al., Measurement of the shape of the boson rapidity distribution for + events produced at =1.96-TeV, Phys. Rev. D76 (2007) 012003, [hep-ex/0702025].
- [64] CDF Collaboration, T. Aaltonen et al., Measurement of the Inclusive Jet Cross Section at the Fermilab Tevatron p-pbar Collider Using a Cone-Based Jet Algorithm, Phys. Rev. D78 (2008) 052006, [arXiv:0807.2204].
- [65] D0 Collaboration, V. M. Abazov et al., Measurement of the muon charge asymmetry in W+X + X events at =1.96 TeV, Phys.Rev. D88 (2013) 091102, [arXiv:1309.2591].
- [66] D0 Collaboration, V. M. Abazov et al., Measurement of the electron charge asymmetry in decays in collisions at TeV, Phys. Rev. D91 (2015), no. 3 032007, [arXiv:1412.2862]. [Erratum: Phys. Rev.D91,no.7,079901(2015)].
- [67] ATLAS Collaboration, G. Aad et al., Measurement of the inclusive and cross sections in the electron and muon decay channels in pp collisions at = 7 TeV with the ATLAS detector, Phys.Rev. D85 (2012) 072004, [arXiv:1109.5141].
- [68] ATLAS Collaboration, G. Aad et al., Measurement of the high-mass Drell–Yan differential cross-section in pp collisions at =7 TeV with the ATLAS detector, Phys.Lett. B725 (2013) 223, [arXiv:1305.4192].
- [69] ATLAS Collaboration, G. Aad et al., Measurement of the Transverse Momentum Distribution of Bosons in Collisions at TeV with the ATLAS Detector, Phys.Rev. D85 (2012) 012005, [arXiv:1108.6308].
- [70] ATLAS Collaboration, G. Aad et al., Measurement of inclusive jet and dijet production in pp collisions at = 7 TeV using the ATLAS detector, Phys. Rev. D86 (2012) 014022, [arXiv:1112.6297].
- [71] ATLAS Collaboration, G. Aad et al., Measurement of the inclusive jet cross section in pp collisions at =2.76 TeV and comparison to the inclusive jet cross section at =7 TeV using the ATLAS detector, Eur.Phys.J. C73 (2013) 2509, [arXiv:1304.4739].
- [72] ATLAS Collaboration, G. Aad et al., Measurement of the cross section for top-quark pair production in collisions at TeV with the ATLAS detector using final states with two high-pt leptons, JHEP 1205 (2012) 059, [arXiv:1202.4892].
- [73] ATLAS Collaboration, G. Aad et al., Measurement of the production cross-section in pp collisions at sqrts = 7 TeV using kinematic information of lepton+jets events, ATLAS-CONF-2011-121, ATLAS-COM-CONF-2011-132.
- [74] ATLAS Collaboration, G. Aad et al., Measurement of the production cross-section in collisions at TeV using events with -tagged jets, ATLAS-CONF-2013-097, ATLAS-COM-CONF-2013-112.
- [75] ATLAS Collaboration, G. Aad et al., Measurement of the transverse momentum and distributions of DrellâYan lepton pairs in protonâproton collisions at TeV with the ATLAS detector, Eur. Phys. J. C76 (2016), no. 5 291, [arXiv:1512.02192].
- [76] ATLAS Collaboration, M. Aaboud et al., Precision measurement and interpretation of inclusive , and production cross sections with the ATLAS detector, arXiv:1612.03016.
- [77] ATLAS Collaboration, G. Aad et al., Measurement of the production cross-section using events with b-tagged jets in pp collisions at = 7 and 8 with the ATLAS detector, Eur. Phys. J. C74 (2014), no. 10 3109, [arXiv:1406.5375]. [Addendum: Eur. Phys. J.C76,no.11,642(2016)].
- [78] ATLAS Collaboration, M. Aaboud et al., Measurement of the production cross-section using events with b-tagged jets in pp collisions at =13 TeV with the ATLAS detector, Phys. Lett. B761 (2016) 136–157, [arXiv:1606.02699].
- [79] ATLAS Collaboration, G. Aad et al., Measurements of top-quark pair differential cross-sections in the lepton+jets channel in collisions at TeV using the ATLAS detector, Eur. Phys. J. C76 (2016), no. 10 538, [arXiv:1511.04716].
- [80] ATLAS Collaboration, G. Aad et al., Measurement of the low-mass Drell-Yan differential cross section at = 7 TeV using the ATLAS detector, JHEP 06 (2014) 112, [arXiv:1404.1212].
- [81] ATLAS Collaboration, G. Aad et al., Measurement of the boson transverse momentum distribution in collisions at = 7 TeV with the ATLAS detector, JHEP 09 (2014) 145, [arXiv:1406.3660].
- [82] CMS Collaboration, S. Chatrchyan et al., Measurement of the electron charge asymmetry in inclusive W production in pp collisions at = 7 TeV, Phys.Rev.Lett. 109 (2012) 111806, [arXiv:1206.2598].
- [83] CMS Collaboration, S. Chatrchyan et al., Measurement of the muon charge asymmetry in inclusive pp to WX production at = 7 TeV and an improved determination of light parton distribution functions, Phys.Rev. D90 (2014) 032004, [arXiv:1312.6283].
- [84] CMS Collaboration, S. Chatrchyan et al., Measurement of the differential and double-differential Drell-Yan cross sections in proton-proton collisions at 7 TeV, JHEP 1312 (2013) 030, [arXiv:1310.7291].
- [85] CMS Collaboration, S. Chatrchyan et al., Measurement of associated W + charm production in pp collisions at = 7 TeV, JHEP 02 (2014) 013, [arXiv:1310.1138].
- [86] CMS Collaboration, S. Chatrchyan et al., Measurement of the production cross section in the dilepton channel in pp collisions at = 8 TeV, JHEP 1402 (2014) 024, [arXiv:1312.7582].
- [87] CMS Collaboration, S. Chatrchyan et al., Measurement of the production cross section in the dilepton channel in collisions at TeV, JHEP 1211 (2012) 067, [arXiv:1208.2671].
- [88] CMS Collaboration, S. Chatrchyan et al., Measurement of the production cross section in collisions at TeV with lepton + jets final states, Phys.Lett. B720 (2013) 83–104, [arXiv:1212.6682].
- [89] CMS Collaboration, V. Khachatryan et al., Measurement of the differential cross section and charge asymmetry for inclusive production at TeV, Eur. Phys. J. C76 (2016), no. 8 469, [arXiv:1603.01803].
- [90] CMS Collaboration, V. Khachatryan et al., Measurement of the inclusive jet cross section in pp collisions at , Eur. Phys. J. C76 (2016), no. 5 265, [arXiv:1512.06212].
- [91] CMS Collaboration, V. Khachatryan et al., Measurement of the t-tbar production cross section in the e-mu channel in proton-proton collisions at sqrt(s) = 7 and 8 TeV, JHEP 08 (2016) 029, [arXiv:1603.02303].
- [92] CMS Collaboration, V. Khachatryan et al., Measurement of the differential cross section for top quark pair production in pp collisions at , Eur. Phys. J. C75 (2015), no. 11 542, [arXiv:1505.04480].
- [93] CMS Collaboration, V. Khachatryan et al., Measurement of the Z boson differential cross section in transverse momentum and rapidity in protonâproton collisions at 8 TeV, Phys. Lett. B749 (2015) 187–209, [arXiv:1504.03511].
- [94] LHCb Collaboration, R. Aaij et al., Inclusive and production in the forward region at TeV, JHEP 1206 (2012) 058, [arXiv:1204.1620].
- [95] LHCb Collaboration, R. Aaij et al., Measurement of the cross-section for production in collisions at TeV, JHEP 1302 (2013) 106, [arXiv:1212.4620].
- [96] CMS Collaboration, S. Chatrchyan et al., Measurements of differential jet cross sections in proton-proton collisions at TeV with the CMS detector, Phys.Rev. D87 (2013) 112002, [arXiv:1212.6660].
- [97] LHCb Collaboration, R. Aaij et al., Measurement of the forward boson production cross-section in collisions at TeV, JHEP 08 (2015) 039, [arXiv:1505.07024].
- [98] LHCb Collaboration, R. Aaij et al., Measurement of forward W and Z boson production in collisions at TeV, JHEP 01 (2016) 155, [arXiv:1511.08039].
- [99] A1 Collaboration, J. C. Bernauer et al., Electric and magnetic form factors of the proton, Phys. Rev. C90 (2014), no. 1 015206, [arXiv:1307.6227].
- [100] CLAS Collaboration, M. Osipenko et al., A Kinematically complete measurement of the proton structure function F(2) in the resonance region and evaluation of its moments, Phys. Rev. D67 (2003) 092001, [hep-ph/0301204].
- [101] M. E. Christy and P. E. Bosted, Empirical fit to precision inclusive electron-proton cross- sections in the resonance region, Phys. Rev. C81 (2010) 055213, [arXiv:0712.3731].
- [102] HERMES Collaboration, A. Airapetian et al., Inclusive Measurements of Inelastic Electron and Positron Scattering from Unpolarized Hydrogen and Deuterium Targets, JHEP 05 (2011) 126, [arXiv:1103.5704].
- [103] H. Abramowicz, E. M. Levin, A. Levy, and U. Maor, A Parametrization of sigma-T (gamma* p) above the resonance region Q**2 ¿= 0, Phys. Lett. B269 (1991) 465–476.
- [104] S. Carrazza, scarrazza/fiatlux: fiatlux v0.0.2, Dec., 2017.
- [105] S. Forte, E. Laenen, P. Nason, and J. Rojo, Heavy quarks in deep-inelastic scattering, Nucl. Phys. B834 (2010) 116–162, [arXiv:1001.2312].
- [106] R. D. Ball, V. Bertone, M. Bonvini, S. Forte, P. Groth Merrild, J. Rojo, and L. Rottoli, Intrinsic charm in a matched general-mass scheme, Phys. Lett. B754 (2016) 49–58, [arXiv:1510.00009].
- [107] R. D. Ball, M. Bonvini, and L. Rottoli, Charm in Deep-Inelastic Scattering, JHEP 11 (2015) 122, [arXiv:1510.02491].
- [108] NNPDF Collaboration, R. D. Ball, V. Bertone, M. Bonvini, S. Carrazza, S. Forte, A. Guffanti, N. P. Hartland, J. Rojo, and L. Rottoli, A Determination of the Charm Content of the Proton, Eur. Phys. J. C76 (2016), no. 11 647, [arXiv:1605.06515].
- [109] T. Carli et al., A posteriori inclusion of parton density functions in NLO QCD final-state calculations at hadron colliders: The APPLGRID Project, Eur.Phys.J. C66 (2010) 503, [arXiv:0911.2985].
- [110] fastNLO Collaboration, M. Wobisch, D. Britzger, T. Kluge, K. Rabbertz, and F. Stober, Theory-Data Comparisons for Jet Measurements in Hadron-Induced Processes, arXiv:1109.1310.
- [111] V. Bertone, S. Carrazza, and N. P. Hartland, APFELgrid: a high performance tool for parton density determinations, Comput. Phys. Commun. 212 (2017) 205–209, [arXiv:1605.02070].
- [112] V. Bertone, R. Frederix, S. Frixione, J. Rojo, and M. Sutton, aMCfast: automation of fast NLO computations for PDF fits, JHEP 1408 (2014) 166, [arXiv:1406.7693].
- [113] J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, et al., The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations, JHEP 1407 (2014) 079, [arXiv:1405.0301].
- [114] Jefferson Lab Hall C E94-110 Collaboration, Y. Liang et al., Measurement of R = sigma(L) / sigma(T) and the separated longitudinal and transverse structure functions in the nucleon resonance region, nucl-ex/0410027.
- [115] E143 Collaboration, K. Abe et al., Measurements of R = (L) / (T) for 0.03 ¡ x ¡ 0.1 and fit to world data, Phys. Lett. B452 (1999) 194–200, [hep-ex/9808028].
- [116] A. M. Cooper-Sarkar, I. Abt, B. Foster, M. Wing, V. Myronenko, and K. Wichmann, Study of HERA data at Low and Low , PoS DIS2016 (2016) 013, [arXiv:1605.08577].
- [117] S. Carrazza, S. Forte, Z. Kassabov, J. I. Latorre, and J. Rojo, An Unbiased Hessian Representation for Monte Carlo PDFs, Eur. Phys. J. C75 (2015), no. 8 369, [arXiv:1505.06736].
- [118] S. Carrazza, S. Forte, Z. Kassabov, and J. Rojo, Specialized minimal PDFs for optimized LHC calculations, Eur. Phys. J. C76 (2016), no. 4 205, [arXiv:1602.00005].
- [119] J. Butterworth et al., PDF4LHC recommendations for LHC Run II, J. Phys. G43 (2016) 023001, [arXiv:1510.03865].
- [120] M. Czakon, N. P. Hartland, A. Mitov, E. R. Nocera, and J. Rojo, Pinning down the large-x gluon with NNLO top-quark pair differential distributions, JHEP 04 (2017) 044, [arXiv:1611.08609].
- [121] R. Boughezal, A. Guffanti, F. Petriello, and M. Ubiali, The impact of the LHC Z-boson transverse momentum data on PDF determinations, JHEP 07 (2017) 130, [arXiv:1705.00343].
- [122] CMS Collaboration, V. Khachatryan et al., Measurements of differential and double-differential Drell-Yan cross sections in proton-proton collisions at 8 TeV, Eur. Phys. J. C75 (2015), no. 4 147, [arXiv:1412.1115].
- [123] N. Chiapolini, Electroweak physics, QCD and jets in the forward region, Tech. Rep. LHCb-PROC-2013-011. CERN-LHCb-PROC-2013-011, CERN, Geneva, Mar, 2013.
- [124] ATLAS Collaboration, M. Aaboud et al., Search for new high-mass phenomena in the dilepton final state using 36 fb of proton-proton collision data at = 13 TeV with the ATLAS detector, JHEP 10 (2017) 182, [arXiv:1707.02424].
- [125] CMS Collaboration, V. Khachatryan et al., Search for narrow resonances in dilepton mass spectra in proton-proton collisions at = 13 TeV and combination with 8 TeV data, Phys. Lett. B768 (2017) 57–80, [arXiv:1609.05391].
- [126] A. Bierweiler, T. Kasprzik, H. Kuhn, and S. Uccirati, Electroweak corrections to W-boson pair production at the LHC, JHEP 1211 (2012) 093, [arXiv:1208.3147].
- [127] ATLAS Collaboration, M. Aaboud et al., Measurement of the production cross section in collisions at a centre-of-mass energy of = 13 TeV with the ATLAS experiment, Phys. Lett. B773 (2017) 354–374, [arXiv:1702.04519].
- [128] CMS Collaboration, V. Khachatryan et al., Measurement of the cross section in pp collisions at 8 TeV and limits on anomalous gauge couplings, Eur. Phys. J. C76 (2016), no. 7 401, [arXiv:1507.03268].
- [129] LHC Higgs Cross Section Working Group Collaboration, D. de Florian et al., Handbook of LHC Higgs Cross Sections: 4. Deciphering the Nature of the Higgs Sector, arXiv:1610.07922.
- [130] A. Buckley, J. Ferrando, S. Lloyd, K. NordstrÃ¶m, B. Page, et al., LHAPDF6: parton density access in the LHC precision era, Eur.Phys.J. C75 (2015) 132, [arXiv:1412.7420].
- [131] The NNPDF Collaboration, R. D. Ball et al., Fitting Parton Distribution Data with Multiplicative Normalization Uncertainties, JHEP 05 (2010) 075, [arXiv:0912.2276].