Photons in the proton: implications for the LHC

We construct a set of parton distribution functions (PDFs), based on the recent NNPDF4.0 PDF set, that also include a photon PDF. The photon PDF is constructed using the LuxQED formalism, while QED evolution accounting for O(alpha), O(alpha alphas) and O(alpha^2) corrections is implemented and benchmarked by means of the EKO code. We investigate the impact of QED effects on NNPDF4.0, and compare our results both to our previous NNPDF3.1QED PDF set and to other recent PDF sets that include the photon. We assess the impact of photon-initiated processes and electroweak corrections on a variety of representative LHC processes, and find that they can reach the 5% level in vector boson pair production at large invariant mass.


Introduction
High precision physics at the LHC, especially in its high-luminosity era [1][2][3], will demand the inclusion of electroweak corrections in the computation of theoretical predictions for hard processes [4].This requires an extension of the set of proton parton distribution functions (PDFs).In particular a photon PDF has to be provided, and evolution equations need to be supplemented with QED splitting functions.The photon PDF enables the inclusion of photon-initiated processes, which typically are enhanced in the high-mass and large transverse-momentum tails of the distributions.In principle, at high-enough scales, proton PDFs should also include PDFs for leptons [5], gauge bosons [6] and indeed for the full set of 52 standard model fields [7], or even new hypothetical particles such as dark photons [8].Even at the LHC lepton PDFs are needed in searches for leptoquarks [9] or exotic resonances that couple to leptons [10].However, in practice, the main requirement of current precision physics at the LHC is the availability of PDF sets that also include a photon PDF, that mixes upon combined QED×QCD evolution with the standard quark and gluon PDFs.We will henceforth refer to such PDF sets as "QED PDFs".
Initial attempts at the construction of QED PDF sets relied on models for the photon PDF at the initial evolution scale [11,12].A first data-driven determination of QED PDFs, based on a fit using the NNPDF2.3methodology [13,14], resulted in a photon PDF with large uncertainties.The fact that determining the photon PDF from the data yields a result affected by large uncertainties was more recently confirmed in determinations based on fitting to the data PDFs with a fixed functional forms within the xFitter methodology [15] applied to high-mass ATLAS Drell-Yan distributions [16].
A breakthrough in the determination of QED PDFs was achieved in 2016 in Refs.[17,18] (see also related results in Refs.[19,20]), where it was shown that the photon PDF can be computed perturbatively in QED, given as input the proton structure functions at all scales, from the elastic (Q 2 → 0) to the deep-inelastic (Q 2 → ∞) regimes -the so-called LuxQED method.Since then, this LuxQED framework has been the basis of all QED PDF sets [21][22][23][24][25].
While the effect of the inclusion of the photon PDF on other PDFs is small, it is not negligible within current uncertainties.For instance, the photon typically carries a fraction of the proton momentum which is about two orders of magnitude smaller than that of the gluon, so the corresponding depletion of the gluon momentum fraction is relevant at the percent level, which is comparable to the size of the current uncertainty on the gluon PDF.Precision calculations of LHC processes with electroweak effects thus require a consistent global QED PDF determination.
With this motivation, we construct a QED PDF set based on the recent NNPDF4.0PDF determination [26,27], including QED corrections to parton evolution up to O (αα s ) and O α 2 .We follow closely the methodology developed for the NNPDF3.1QEDPDFs in Ref. [21], based on an iterative procedure.Namely, we evaluate the photon PDF at some fixed scale using the LuxQED formula and the structure function computed from an existing PDF set, now NNPDF4.0.We evolve it together with the other PDFs to the initial parametrization scale Q 0 = 1.65 GeV.We then re-determine the quark and gluon PDFs while also including this photon PDF as a boundary condition to the QCD×QED evolution.We compute the photon PDF again using LuxQED, and we iterate until convergence.
All results presented in this paper are obtained using a new theory implementation, now adopted by default in the NNPDF public code [27].In particular, the APFEL [28] and APFELgrid [29] codes have been replaced by a suite of newly developed tools including the evolution equations solver EKO [30], the DIS structure functions calculator YADISM [31], and the interpolator of hard-scattering cross-sections PineAPPL [32,33].Taken together, they provide a pipeline for the efficient automatization of the computation of theory predictions [34].An integral component of this theory pipeline is the use of interpolation grids in the format provided by PineAPPL [32,33], which can be interfaced to Monte Carlo generators.Among these is mg5 aMC@NLO [35], which automates NLO QCD+EW calculations for a wide range of LHC processes where knowledge of QED PDFs is crucial.
The outline of this paper is as follows.First, Sect. 2 reviews the theoretical framework underlying the NNPDF4.0QEDdetermination and in particular the implementation of QED evolution in EKO.Then the NNPDF4.0QEDPDFs are presented in Sect.3, where they are compared to the previous NNPDF3.1QEDPDF set, and to other recent QED PDF sets.Implications for LHC phenomenology are studied in Sect. 4 by means of the PineAPPL interface to mg5 aMC@NLO.Conclusions and an outline of future developments are finally presented in Sect. 5.Because this paper is, as mentioned, the first to make use of a new theory pipeline, and specifically its implementation of combined QCD×QED evolution, a set of benchmarks is collected in two appendices.Specifically, in Appendix A we benchmark NNPDF4.0PDFs determined using the old and new pipeline and show that they are indistinguishable; in Appendix B the new implementation of joint QCD and QED evolution of PDFs in the EKO code is discussed and benchmarked against the previous implementation in the APFEL code.

Evolution and determination of QED PDFs
In this section, we discuss the structure of combined QED×QCD evolution and we briefly review the methodology used to determine QED PDFs.This methodology is based on the LuxQED formalism [17,18], and was developed for the construction of the NNPDF3.1QEDPDF set [21].
The scale dependence of the strong and electromagnetic couplings is in turn determined by coupled renormalization group equations of the form in terms of the coefficients of the corresponding QCD and QED beta functions.Consistently with the treatment of evolution equations, we include pure QCD contributions up to NNLO, namely β (2,0) QCD , β (3,0)  QCD , and β (4,0)  QCD ; pure QED up to NLO, namely β (0,2) QED and β (0,3) QED ; and the leading mixed terms β (2,1)  QCD and β (1,2)  QED [40].We adopt the MS scheme.Schemes in which the electroweak coupling does not run, such as the G µ scheme, are commonly used in the computation of electroweak corrections, but MS is more convenient when considering combined QCD and QED corrections [41,42].Equations (2.1-2.4) must then be simultaneously solved with a common scale µ.
The solution is most efficiently obtained in a maximally decoupled basis in quark flavor space.This requires adopting a suitable combination of quark and antiquark flavors such that the sum over j in Eq. (2.1) contains the smallest possible number of entries.In the case of QCD-only evolution, this is achieved in the so-called evolution basis, in which one separates off the singlet combination Σ = i (q i + qi ), which mixes with the gluon, and then one constructs nonsinglet combinations of individual C-even (sea-like) and C-odd (valence-like) q ± i = q i ± qi quark and antiquark flavors, each of which evolves independently.Because the photon couples differently to up-like u k = {u, c, t} and down-like d k = {d, s, b} quarks, a different basis choice is necessary when also including QED.To this end, given n f active quark flavors, we split them into n u up-like and n d down-like flavors, such that n f = n u + n d , and define the four combinations Also, we separate off the QCD contributions to the anomalous dimensions, by rewriting the perturbative expansion of the anomalous dimensions, Eq. (2.2), as where γ ij (a s ) contains the pure QCD contributions and γij contains both the pure QED and the mixed QCD×QED corrections.The maximally decoupled evolution equations are then constructed as follows.The nonsinglet combinations evolve independently according to nonsinglet evolution equations of the form ) (2.10) The valence sum and difference combinations, defined as satisfy coupled evolution equations in terms of the linear combinations of anomalous dimensions Finally, the gluon and photon satisfy coupled evolution equations together with the quark singlet sum and difference combinations, defined as These coupled evolution equations read with Γ a 4 × 4 anomalous dimension matrix of the form where γps qq ′ = e 2 q e 2 q ′ γps [38], and the combinations ⟨γ gq ⟩, ⟨γ qg ⟩, γg∆q , γ∆qg , ⟨γ γq ⟩, ⟨γ qγ ⟩, γγ∆q , and γ∆qγ are constructed analogously to those listed in Eq. (2.13).
The basis defined by Eqs.(2.7), (2.8), (2.11) and (2.14) is denoted as the unified evolution basis.Compared to the basis used in APFEL [28] to solve the QCD×QED evolution equations, our definitions of Σ ∆ and V ∆ differ due to the prefactors n d /n u that make the basis fully orthogonal.In the presence of scale-independent intrinsic heavy quarks (e.g. a charm-quark PDF in a three-flavor scheme [43]), a further decomposition needs to be applied [30].
The unified flavor basis has been implemented in the EKO code [30], which, as discussed in the introduction, is now used to solve evolution equations as part of the new theory pipeline.The basic ingredient of the EKO code is the construction of evolution kernel operators where f (Q 2 ) is a vector whose components are all PDF flavors, including all active quark and antiquark flavors in a suitable basis, the photon and the gluon.Formally, the EKOs are given by where γ is the full matrix of anomalous dimensions, and P denotes path ordering.In the unified evolution basis the evolution kernel is a block diagonal matrix, with individual diagonal entries for the nonsinglet combinations Eqs.(2.7,2.8)(without path-ordering), a 2 × 2 block for the valence combinations Eq. (2.11), and a 4 × 4 block in the singlet sector.
A variety of implementations of the solution of the evolution equations, each of those corresponding to a determination of the EKO, and which differ by higher-order terms, are available in the EKO code.These are discussed in Ref. [30].At NNLO in QCD, the solutions are essentially indistinguishable in the data region.The solution adopted here is the iterated-exact (EXA) of Ref. [30], while the truncated (TRN) solution was adopted for the NNPDF3.1QED[21] and NNPDF4.0[26] PDF sets, using the APFEL [28] implementation.These different implementations (at NNLO) are benchmarked and shown to be equivalent in Appendix B, where we also discuss the motivations for this choice.

Construction of the photon PDF
As mentioned in Sect. 1, the photon PDF is determined using the LuxQED [17,18] formalism, implemented in a PDF fit using the same methodology as the one used for the NNPDF3.1QEDPDF set [21], but now with the new theory pipeline.The LuxQED result amounts to proving that the photon PDF is perturbatively determined in QED by knowledge of the proton inclusive structure functions F 2 and F L : where m p is the mass of the proton and P γq is the photon-quark splitting function.Equation (2.18) holds in the MS scheme, including terms of order a em and a 2 em ln µ 2 /m 2 p , times the accuracy of the QCD determination of the structure functions.
The integration over the scale Q 2 and Bjorken-x dependence of the structure functions F i (x, Q 2 ) in Eq. (2.18) includes four different regions and corresponding contributions to the structure functions: an elastic contribution at x = 1, a resonance contribution when x is large and close to x = 1, an inelastic nonperturbative contribution at low Q 2 and intermediate x and an inelastic perturbative region for intermediate x and large Q 2 .The first three contributions must be determined by fits to lepton-proton scattering data, but the latter contribution may be determined by expressing the structure functions through perturbative QCD factorization in terms of quark and gluon PDFs.In a global PDF determination, one may choose to determine the photon PDF using the LuxQED formula Eq. (2.18) at a single chosen scale µ 2 = Q 2 γ , and then evolve jointly the photon and all other PDFs through the QCD×QED evolution equations discussed in Sect.2.1.An alternative option is to determine the photon PDF using Eq.(2.18) at all scales.The two choices are equivalent because the LuxQED photon Eq. (2.18) satisfies the joint QED×QCD evolution equations to the accuracy of the LuxQED formula itself, though they differ by higher-order corrections and also due to the fact that the LuxQED photon is partly determined from a parametrization of data.The former choice was made for the construction of our previous NNPDF3.1QEDset [21], which we follow here, also in the choice of parametrization of the structure functions in the non-perturbative region.The same choice was made in Refs.[21][22][23][24][25], though in Refs.[24,25] a PDF set in which the photon PDF is determined using Eq.(2.18) at all scales was also presented.
The relative size of the elastic, inelastic non-perturbative (including resonances) and elastic perturbative contributions to γ(x, µ 2 ) is displayed in Fig. 2.1 as a function of x for four different choices of the scale µ = Q γ at which the LuxQED formula Eq. (2.18) is used to determine the photon (see also Fig. 18 in Ref. [18]).If the LuxQED formula is used at the scale Q γ = 1.65 GeV at which the NNPDF4.0PDFs are parametrized, the photon is entirely determined by the elastic and non-perturbative contributions, but as the scale is increased, an increasingly large contribution comes from the perturbative region.By the time a value Q γ ≳ 100 GeV is reached, the photon at small x ∼ < 10 −3 is almost entirely determined by the perturbative contribution, while the elastic and non-perturbative contributions only remain dominant at large x ∼ > 0.1 where however the photon PDF is tiny.Choosing a large value of Q γ has the dual advantage that the LuxQED result is more accurate at high scale, because it includes contributions of order a em and a 2 em ln µ 2 /m 2 p , and also that one is then mostly relying on the accurate perturbative determination of the structure functions, that exploits global information on quarks and gluon PDFs and not just the leptonproton scattering data.As in Ref. [21], we choose µ = Q γ = 100 GeV, as also advocated in Ref. [18]. 1 We will discuss the dependence of our results on the choice of scale Q γ in Sect.3. Equation (2.18) must be viewed as a constraint on the set of photon, quark and gluon PDFs that are simultaneously determined.Because of the small impact of the photon PDF on the other PDFs, it was suggested in Refs.[18,21] that the constraint can be implemented iteratively.Namely, we first determine the photon PDF using Eq.(2.18) by means of structure functions that are determined from an existing PDF set at a given scale µ = Q γ .This photon PDF is then evolved, by solving joint QCD×QED evolution equations with the unchanged given set, to a chosen PDF parametrization scale where it is taken as fixed.All other PDFs are then re-determined, with the constraint that the momentum sum rule now also includes a contribution from the given (fixed) photon, i.e.
The photon PDF at a scale µ = Q γ is then determined again from Eq. (2.18), in which structure functions are obtained from the new fit.The procedure is iterated until convergence, which was achieved in two iterations in the NNPDF3.1QEDdetermination [21].
Here we follow the same procedure, starting with a re-determination of the NNPDF4.0PDF set with the new pipeline, which is compared and shown to be equivalent to the published NNPDF4.0 in Appendix A.

The NNPDF4.0QED parton distributions
We present here the NNPDF4.0QEDPDFs: we first summarize our procedure, then discuss the effect on fit quality of the inclusion of the photon PDF, examine the photon PDF itself, also in comparison to other determinations, and finally study the photon momentum fraction.

Construction of the NNPDF4.0QED parton set
As mentioned in Sect. 2 all the methodological aspects and settings of the PDF determination are the same as used for the underlying pure QCD NNPDF4.0PDF [26], but now using a new theory pipeline.Even if in principle theory predictions should be independent of implementation details, in practice differences may arise, e.g.due to issues of numerical accuracy.Also, in the process of transitioning to this new pipeline, a few minor bugs in data implementation were uncovered and fixed.Finally, the new pipeline includes a new implementation of heavy-quark mass effect in deep-inelastic structure functions that differs by subleading terms from the previous one.Benchmarks showing the equivalence of the old and new pipeline are briefly presented in Appendix A. Because of this equivalence, NNPDF4.0PDF replicas produced with the new pipeline should be considered equivalent to the published ones, and indeed for phenomenological applications (specifically those presented in Sect.4) we will compare the results obtained using NNPDF4.0QEDPDFs to pure QCD results obtained using the published NNPDF4.0 replicas.Nevertheless, in this section only, for all comparisons between pure QCD and QCD×QED, we will use NNPDF4.0pure QCD replicas generated using the new pipeline, in order to avoid even small confounding effects.Note however that, as discussed in the end of Sect.2.1, the QED PDFs are based on the EXA solution of the evolution equations, that differs by subleading terms from the TRN solution used in the pure QCD fit.The effect of this difference is assessed in Appendix B and is very small at NNLO.However, all comparisons between pure QCD and QED PDFs also include the effect of this change.
The NNPDF4.0QEDPDFs are determined at NLO and NNLO by supplementing with a photon PDF the pure QCD PDF set, according to the methodology outlined in the previous section.However, all theory predictions are obtained as in the pure QCD determination: hence in particular no photon-induced contributions are included, and thus the only effect of the inclusion of a photon PDF is through its mixing with other PDFs.This is justified because the NNPDF4.0dataset was constructed including cuts that remove all datapoints for which the effect of electroweak corrections is larger than the experimental uncertainties (see Ref. [26], Sect.4.1).The inclusion of the electroweak corrections, which will allow relaxing these cuts, is left for future work.As in Ref. [21], the final PDF set is obtained after two iterations, with a third iteration providing a check of convergence.

Fit quality
Table 3.1 displays the statistical estimators obtained using a set of 100 NNPDF4.0QEDNLO and NNLO PDF replicas, compared to their QCD-only counterparts, generated using the new theory pipeline.Specifically, we show the χ2 , for both the full dataset and for datasets grouped by process; for the full dataset we also show the average over replicas of the training and validation figures of merit ⟨E tr ⟩ rep and ⟨E val ⟩ rep , and the average χ 2 over replicas χ 2 rep , all as defined in Table 9 of Ref. [45].Note that χ 2 and χ 2 rep are computed using the experimental covariance matrix, while, as in all NNPDF determinations, the figure of merit used for minimization is computed using the t 0 covariance matrix [46].
It is clear from Table 3.1 that all estimators corresponding are essentially unchanged by the inclusion of QED corrections, with slightly larger differences seen at NLO than at NNLO: the impact of QED effects on fit quality is negligible, both globally and for individual processes.Specifically, the training and validation figures of merit in the pure QCD and QCD×QED determinations differ by less than one sigma.As mentioned, these differences also include the effect of switching from the TRN solution of evolution equations in the pure QCD fit to the EXA solution.The impact of this change is yet smaller than that of the QED corrections, but it slightly reduces it.

The photon PDF
In Fig. 3.1, the NNLO NNPDF4.0QEDphoton PDF at Q = 100 GeV is compared to its counterpart in our previous NNPDF3.1QED[21] set, and in the recent QED PDF sets MSHT20QED [23], and CT18QED [25]. 2  Here and elsewhere in this paper all uncertainties correspond to 1σ. Results agree at the percent level, despite the fact that quark and gluon PDFs in these sets can display much larger differences.This is a consequence of the fact that in all these PDF sets the photon PDF is determined with the LuxQED formalism in terms of the proton structure function, and that the latter, in turn, is well constrained by experimental data both at high and low scale.In fact, we have checked that the dominant contribution to the difference between the NNPDF3.1QEDand NNPDF4.0QEDphoton PDFs seen in Fig. 3.1 is due to the difference between the TRN and EXA solutions of evolution equations (see Section 3.1 and Appendix B), i.e. to higher-order corrections, with the residual difference being due to the change in the PDFs used in order to compute   Top left: the photon PDF at Q = 100 GeV in NNPDF4.0QEDNNLO compared to its NNPDF3.1QEDcounterpart as a ratio to the central value of the former.Top right: the relative PDF uncertainties on the photon PDF in these two determinations.Bottom: same as top panels now comparing NNPDF4.0QED with MSHT20QED [23] and CT18QED [24] (all NNLO).Bands correspond to 1σ uncertainties.  the structure function.The fact that the 3.1 and 4.0 PDFs are compatible within uncertainties is thus a consequence of the fact that the uncertainties due to higher-order corrections and to PDFs are correctly accounted for by the LuxQED construction, while the NNPDF4.0PDFs are backward-compatible with the NNPDF3.1 PDFs [26].
The uncertainty on the photon PDF is completely dominated by theoretical uncertainties on the LuxQED procedure, which include [18] missing higher order corrections, uncertainties on the experimentally measured low-scale structure function and so on.In our uncertainty determination we follow Ref. [18]; MSHT also mostly follows this reference, with an extra higher-twist contribution due to the low choice of scale Q γ , [22] and indeed finds a very similar uncertainty.A somewhat more conservative uncertainty estimate is provided by CT18QED [25], which also adopts a somewhat different determination of the elastic contribution.
In Fig. 3.2 the central value of the photon PDF in all these sets is compared at different scales, using the native scale dependence of each PDF set.Uncertainties are not shown in order not to clutter the plot.Note that differences in the scale dependence of the PDFs compared in the Figure also arise due to somewhat different treatments of the QCD×QED evolution equations.Specifically, the NNPDF3.1QEDPDF set adopts a numerical implementation of the TRN solution (see Appendix B).However, one would expect differences to be mostly driven by the mixing with the quark and gluon PDFs.The fact that differences grow at low scales suggests that this is indeed the case.Even so, all photon PDFs agree within about 3% for all x ≳ 10 −3 , even at the lowest scale Q = 1.65 GeV.
In Fig. 3.2 we also show the central photon PDF which is found by repeating the NNPDF4.0QEDdetermination with the photon PDF determined at a scale Q γ = 10 GeV instead of the default Q γ = 100 GeV.If a low value of Q γ is adopted, the upper limit of integration in Q 2 , Eq. ( 2  to the LuxQED formula may then become relevant.We see that the shift in central photon PDF which is found by making this choice instead of the default one is at most of the order of the uncertainty on the photon PDF.Again, this shows that the uncertainty on the LuxQED procedure is correctly estimated. We finally assess the impact of the inclusion of a photon PDF on the other PDFs, by comparing the NNLO NNPDF4.0QED and NNPDF4.0(pure QCD) PDFs.In Fig. 3.3 we present this comparison at Q = 100 GeV.It is clear that the impact of the inclusion of the photon is moderate, with the NNPDF4.0(pure QCD) and NNPDF4.0QEDPDFs generally differing by less than one sigma and always in agreement within uncertainties.The largest effects are seen in the gluon PDF, which is suppressed at the percent level due to the momentum fraction transferred from the gluon to the photon.

The photon momentum fraction
The main impact of the photon on other PDFs is through its contribution to the momentum sum rule, Eq. (2.19).We quantify this contribution by evaluating the photon momentum fraction In Fig. 3.4 the momentum fractions carried by the photon and by the gluon PDFs in the NNPDF4.0QEDset are shown (in percentage) as a function of scale.For the photon the result is compared to that of NNPDF3.1QED, and for the gluon to that of NNPDF4.0 (pure QCD).The photon momentum fraction in the NNPDF3.1QEDand NNPDF4.0QEDPDF sets is essentially the same: the photon carries around 0.2% of the proton momentum at a low (Q∼1 GeV) scale, growing logarithmically with Q up to around 0.6% at the multi-TeV scale.The momentum fraction carried by the gluon is reduced by a comparable amount upon inclusion of the photon.

Implications for LHC phenomenology
We now study the phenomenological implications of the NNPDF4.0QEDPDF set.First we compare parton luminosities to those computed using other QED PDF sets.Then we assess the impact of QED corrections on selected processes, by comparing to the pure QCD case calculations that include photon-induced contributions, and also by directly comparing results obtained using NNPDF4.0(pure QCD) and NNPDF4.0QEDPDFs.

Luminosities
The phenomenological effect on the parton luminosity of the inclusion of a photon PDF is both direct, through the presence of photon-induced partonic channels, and indirect, through the effect of the photon on other PDFs, mostly through the depletion of the gluon that is necessary in order to preserve the momentum sum rule, as discussed in Sect.The luminosities for all other parton channels are shown in Fig. 4.2, where we compare NNPDF4.0QEDboth to the previous set NNPDF3.1QED and to the pure QCD NNPDF4.0.Because, as shown in Sect.3, the effect of the inclusion of the photon on the other PDFs is moderate, the comparison between the two QED sets is very similar to the comparison between the pure QCD NNPDF4.0 and NNPDF3.1 shown in Fig. 9.1 of Ref. [26], where it was shown that NNPDF4.0 are backward compatible, i.e. generally agree with NNPDF3.1 within uncertainties.The effect of the inclusion of QED corrections is, as expected, mostly seen in the gluon-gluon channel, where it leads to a suppression of a few percent in the m X ∼ 100 GeV region in order to account for the transfer of a small amount of momentum to the photon.The comparison to other PDF sets is dominated by the differences in the quark and gluon PDFs, which are rather more significant than the difference in the photon PDF, and thus very similar to the corresponding comparison of pure QCD luminosities shown in Fig. 9.3 of Ref. [26].

Physics processes
We now study the impact of the photon PDF on a few representative processes: Drell-Yan production (neutral-and charged-current), Higgs production in gluon-gluon fusion, in vector boson fusion, and in associated production with weak bosons, diboson production, and top-quark pair production, all at the LHC with center-of-mass energy √ s = 14 TeV.We have computed theory predictions for these processes exploiting the PineAPPL [32,33] interface to the automated QCD and EW calculations provided by mg5 aMC@NLO [35].PineAPPL produces interpolation grids, accurate to NLO in both the strong and electroweak couplings, that are independent of PDFs.They therefore make it easy to vary the input PDF set, since the same grid can be used for all PDF sets considered.As we are interested in assessing differences between PDF sets, rather than in doing precision phenomenology, we do not include NNLO QCD corrections, and we only show PDF uncertainties.For a detailed list of parameters and cuts used in the calculation of these grids, see Sect. 9 of Ref. [26].
For all processes, we display predictions in figures below with a standardized format, as follows.We show the absolute distributions (top panels) and the ratio to the central value obtained using NNPDF4.0QEDPDFs and including photon-induced channels, which we call NNPDF4.0QED(bottom panels).In the left plots we compare NNPDF4.0QEDto: NNPDF4.0QEDbut with no photon-initiated channels; NNPDF4.0 pure QCD; NNPDF3.1QED.In the right plots we compare NNPDF4.0QED to MSHT20QED and CT18QED.The left plots allow for assessing the overall size of the QED corrections (by comparison of the QED and pure QCD results), and disentangling the size of the photon-initiated contributions (by comparison of predictions with the photon-initiated channels switched on and off) and the impact of the changes in the quark and gluon PDFs due to QED effects (by comparison of NNPDF4.0 with NNPDF4.0QED with the photon-initiated channels switched off).
In Fig. 4.4 we show results for inclusive Drell-Yan production both in the neutral-current and chargedcurrent channels.For neutral current we show results for the invariant mass distribution of the dilepton pair, while for charged-current we display the rapidity distribution of the lepton.Predictions for the Higgs rapidity distribution are shown in Fig. 4.5 for gluon fusion, associated production with a W + boson, and vector boson fusion.Finally, in Fig. 4.6 we show predictions for weak-boson pair production (W + W − and W + Z), as a function of the dilepton transverse momentum, and for top-quark pair production as a function of the invariant mass of the top quark pair.
For charged-current Drell-Yan, the QED corrections have essentially no effect.In all the other cases that we consider, the effects of the QED correction fall in one of two categories.Either upon inclusion of QED effects we see an enhancement of the cross-section, which is only present when the photon-induced contribution is included, while the NNPDF4.0QEDresult without photon-induced contribution is very close to the pure QCD result.Or else we see a suppression of the cross-section, but with the NNPDF4.0QEDresult with and without the photon-induced contribution very close to each other.Of course, the former case can be explained with the presence of a sizable positive photon-induced contribution, while the latter case is explained by the suppression of the gluon luminosity due to the transfer of momentum fraction from the gluon to the photon seen in Fig. 4.2.
The first situation -enhancement due to the photon-induced contribution -is observed in neutralcurrent dilepton production, where the enhancement increases with invariant mass and can reach up to 5% at the TeV scale.Note that the photon-induced contribution to the dilepton final state can proceed also through t-and u-channel leading-order diagrams.The fact that the QED enhancement is absent at the Z peak suggests that this non-resonant contribution provides the dominant part of the photon-induced contribution.A similar situation occurs in W + W − and ZW + , where the enhancement increases with p T and reaches 5% in the former case and 2% in the latter case for transverse momenta in the TeV range.Finally, the enhancement is also observed in associate Higgs production with W + and in vector boson fusion.In the former case the enhancement is largest at forward rapidity, where it reaches 4%, and it decreases to 2% for the largest rapidity y H = 2.5.In the latter case the enhancement is very moderate, around 1%, and almost independent of the rapidity (though slightly decreasing as the rapidity increases).
The second situation -suppression, independent of the photon-induced contribution -is found in processes that proceed through gluon fusion.The effect is clearly seen in Higgs production in gluon fusion, where the suppression is weakly dependent on rapidity, varying between 2% at central rapidity and about 1% at the largest rapidity.A similar, but more moderate effect, is also seen in top pair production, where the suppression is of order 1%, essentially independent of the invariant mass of the top pair.
When comparing results obtained using different QED PDF sets, be they NNPDF3.1QED,MSHT20QED or CT18QED, we observe that differences are essentially driven by the difference in quark and gluon PDFs.This is a direct consequence of the similarity of the photon PDF in all sets.Indeed, all comparisons are quite similar to those shown in Sect.9.3 of Ref. [26], where the same processes were studied in pure QCD.From top to bottom: neutral-current dilepton production as a function of the dilepton invariant mass m ℓ l; W + production as a function of the antilepton pseudo-rapidity ηl; W − production as a function of the antilepton pseudorapidity η l .For each process, we display the absolute distributions (top panels) and the ratio to the central value obtained with NNPDF4.0QEDPDFs, including photon-initiated channels (bottom panels).In the left panels the full NNPDF4.0result is compared to NNPDF4.0 (QCD only), NNPDF3.1QED,NNPDF4.0QED(no photon-initiated); in the right panels it is compared to MSHT20QED [23] and CT18QED [25].Note that the experimentally measurable quantity is shown, so for dilepton production the t-and u-channel photon-induced contribution is also included.
All in all, we conclude that the inclusion of QED corrections is important for precision phenomenology at the percent level, even in cases in which the photon-induced contribution is negligible, such as Higgs production in gluon fusion.Here neglecting the indirect effect of including the photon PDF results in an overestimation of the peak cross-section (and thus the total cross-section) by about 2%, thus biasing the prediction by an amount that is of the same order as the PDF uncertainty, but not included in it.

Summary and outlook
We have presented a new determination of QED PDFs based on the NNPDF4.0set of parton distributions, using the methodology previously adopted for the construction of the NNPDF3.1QEDPDFs [21].This methodology implements the LuxQED procedure [17,18] to determine the photon PDF.Results are consistent with previous studies: specifically, we find that QED effects have a small but non-negligible impact mostly on the gluon PDF, and that photon-initiated contributions are most important for high-mass process such as neutral-current Drell-Yan production.This PDF determination is based on a new NNPDF pipeline for producing theory predictions [34].Thanks to the integration of QED evolution in this pipeline, and in particular thanks to the use of the EKO evolution code [30], the production of QED variants of NNPDF determinations is now essentially automated, and will become the default in future releases.Indeed, there is in general no reason to switch off the photon PDF, even if it is a small correction, nor to neglect its effect on the momentum fraction carried by the gluon.The NNPDF4.0QEDPDF sets are made publicly available via the LHAPDF6 interface, http://lhapdf.hepforge.org/ .
All sets are delivered as sets of N rep = 100 Monte Carlo replicas.Specifically, we provide the NLO and NNLO global fits constructed with the settings defined in Sect. 2 and denoted as NNPDF40 nlo as 01180 qed NNPDF40 nnlo as 01180 qed These sets are also made available via the NNPDF collaboration website https://nnpdf.mi.infn.it/nnpdf4-0-qed/ .
They should be considered the QED PDF counterparts of the published NNPDF4.0QCD-only PDF sets [26].
We also make available the set of 100 NNPDF4.0(QCD only) replicas used to produce the comparisons shown in Sect.3.These differ from the published NNPDF4.0PDFs because they have been produced using the new theory pipeline, and include some minor bug corrections in the implementation of the NNPDF4.0dataset.These are called: NNPDF40 nlo as 01180 qcd NNPDF40 nnlo as 01180 qcd .
The equivalence of these replicas to the published NNPDF4.0 replicas is demonstrated in Appendix A. They are made available for completeness, on the NNPDF website only.
The NNPDF4.0QED determination is part of a family of developments based on the NNPDF4.0PDF set, and aimed at increasing its accuracy.These will also include a determination of the theory uncertainty on NNPDF4.0PDFs based on the methodology of Refs.[48,49], and a first PDF determination based on NNPDF methodology at approximate N 3 LO [50].All of these, as well as their combination, will become part of the default NNPDF methodology in future releases.
A natural development of QED PDFs is the full inclusion of electroweak corrections in theory predictions used for PDF determination, which enables a widening of both the set of processes and the kinematic range that may be used for PDF determination, which are currently constrained by the need of ensuring that electroweak corrections are small.This is especially important as electroweak effects become more relevant in regions of phase space sensitive to the large-x PDFs, which are in turn relevant for new physics searches [51,52].This development will be greatly facilitated by the availability, within the new pipeline, of the PineAPPL interface to automated Monte Carlo generators, such as mg5 aMC@NLO.All these developments will help achieve PDF determination with percent or sub-percent accuracy.

A The new NNPDF theory pipeline
As mentioned in Sect.3, this work is based on a new pipeline for the calculation of theoretical predictions.This new theory pipeline is described in Ref. [34]; it supersedes and replaces the one used for the NNPDF4.0determination, which was based on a combination of different pieces of code of various origin, specifically APFEL [28] for PDF evolution and DIS structure function computation, and APFELgrid [29] for the generation of interpolation grids of NLO partonic matrix elements.This last piece of code made use of APPLgrid [53] and FastNLO [54] interpolators.The main benefit of the new pipeline is its unified, yet modular and flexible, structure.
This theory pipeline has been extensively benchmarked for numerical accuracy, including QCD evolution, the computation of structure functions, the interpolation of grids, and the interfacing to mg5 aMC@NLO.
Using this new pipeline, interpolation grids for a number of processes have been recomputed, as discussed below.Therefore Tables 2.1-2.5 in Ref. [26] should be updated accordingly.
Specifically, in this code heavy quark mass effects are accounted for using a new implementation of the FONLL prescription [56].
• The following datasets of Table 2.4 and Table 2.5 in Ref. [26] have been recomputed using mg5 aMC@NLO (all grids, in the PineAPPL format, are available at https://github.com/NNPDF/pineapplgrids): Grids for all remaining datasets have been converted to the PineAPPL format without recomputing them.Of course, none of these changes except the new FONLL implementation should make any difference to the extent that the previous and new theory implementations are both numerically accurate.The new FONLL implementation differs from the previous one by subleading corrections and thus it does introduce NNLO differences in the NLO fit and N 3 LO differences in the NNLO fit; these differences are confined to the charm and bottom mass corrections to deep-inelastic structure functions and thus only affect a small number of datapoints, at NNLO at the sub-percent level.
Finally, in the process of transitioning to the new pipeline, a few bugs were discovered in the implementation of a few datapoints, such as incorrect normalization, incorrect scale assignment or incorrect bin size.These corrections in practice have a negligible impact, as they involve a handful of points out of more than 4500.
We have generated variants of the published NNPDF4.0NLO and NNLO fits using the new theory pipeline.In order to illustrate the equivalence between these fits and the NNPDF4.0fits, obtained with the previous theory pipeline, we compare each pair of fits.All fits are made of 100 replicas.The corresponding statistical estimators are reported in Table A.1.In addition to the fit quality estimators, we also show the ϕ estimator, defined in Eq. (4.6) of Ref. [58], which measures the ratio of the average (correlated) PDF uncertainty on datapoints over the experimental uncertainty, and is thus a fairly sensitive measure of fit quality.All estimators are seen to coincide at NNLO (note that the statistical uncertainty on the χ 2 with the given number of datapoints is 0.02), demonstrating explicitly the equivalence between the two replica sets.All estimators are also very close at NLO: differences are compatible with the statistical uncertainty on the χ 2 .[45]) between central values and uncertainties computed using the two replica sets at Q = 100 GeV.For a set of 100 replicas, d ∼ 1 indicates statistical equivalence, while d ∼ 10 indicates one-sigma differences.At NNLO, distances fluctuate about 1, indicating that the two replica sets are drawn by the same underlying distribution.At NLO, distances are somewhat larger, around the half-sigma level for most PDFs and reaching 1σ in some cases; this is mostly due to the different FONLL implementation [56]  respectively at NLO and NNLO.The complete agreement between the two pairs of NNLO replica sets, and the small differences at NLO, are manifest.

B Solution of evolution equations
We summarize results and benchmarks concerning the solution to evolution equations.First, we comparatively review solutions to evolution equations, for which a different choice is made in this work in comparison to the previous NNPDF3.1QEDand NNPDF4.0PDF determinations: namely the use of an exact (EXA) solution instead of the truncated (TRN) solution [30].We then study the impact of this different choice on  PDFs.Having established the equivalence of (pure QCD) PDFs obtained using either solution, we perform an inversion test of the exact QED×QCD solution implemented in EKO, which we adopt, in order to verify its accuracy.We finally benchmark this EKO QED×QCD solution against the implementation in the APFEL code, which was used in the previous NNPDF3.1QEDPDF determination.

B.1 Exact vs. expanded solutions: formal aspects
Perturbative solutions of QCD evolution equations of the form are based on the observation that if the equation is solved exactly at leading order, with the leading-order term in the beta function Eq. (2.3), then the solution f LO (Q 2 ) in terms of a boundary condition f LO (Q 2 0 ) is a pure leading log (LL) function, i.e. it is a function of a s L only, where L = ln(Q 2 /Q 2 0 ), rather than depending on a s and L separately.However, if the next-to leading order contributions γ (2,0) ij and β 1 are also included, then the solution is next-to-leading log (NLL) accurate, but the exact solution Eqs. (2.16) and (2.17) includes terms to all logarithmic orders.A solution that only includes NLL terms, namely with the structure where both f LO and f NLO are pure LL functions, can be constructed by linearization of the exact solution, namely by expanding the exact solution and by truncating the expansion.Various intermediate options are also possible.The argument can be repeated at any logarithmic accuracy.Of course all these solutions are equivalent up to subleading logarithmic terms.
A truncated solution to the pure QCD evolution equation of the form Eq. (B.2) can be determined in closed form to all orders in Mellin space [59], by diagonalizing order by order the anomalous dimension matrix.This corresponds to the solution referred to as truncated in Ref. [30].However, this strategy fails for combined QED×QCD evolution, because the QED and QCD anomalous dimension matrices do not commute, and thus cannot be diagonalized simultaneously.Because a s and a em depend on scale in different ways, this implies that the anomalous dimension matrices γ ij evaluated at different scales do not commute: This means that the LO solution, constructed only including the LO contributions to the anomalous dimen- sions, takes the form of a path-ordered exponential that cannot be written in closed form.Note that the commutator terms are not subleading: for example, the quadratic term is proportional to a s a em L 2 , so it is a LL contribution.
The problem persists to all orders, and requires a truncated solution to be also given as a path-ordered exponential.For instance, including NLO QCD and LO QED contributions, the evolution equation is The perturbative NLO QCD solution then has the form Eq. (B.2), where in the general case f LO and f NLO are LO accurate in QCD and QED, but may also include subleading contributions, i.e. f Substituting the perturbative expansion Eq. (B.2) in the NLO Eq. (B.5), and using the LO solution, Eq. (B.4), leads to Again, due to the non-commutativity of γ (1,0) ij and γ (0,1) ij , the solution for f NLO is a path-ordered exponential and cannot be given in closed analytic form.This continues to be the case at higher orders.
A truncated solution can be constructed numerically, by starting with the exact path-ordered solution Eqs. (2.16-2.17),and then expanding numerically.Such an approach, however, involves approximating  A.1, now comparing a set of 100 NNPDF4.0replicas obtained using the default truncated solution of evolution equations with a set obtained using the exact solution (pure NNLO QCD theory).Both PDF sets are produced using the new theory pipeline.higher-order derivatives by finite differences [28], which may lead to numerical instabilities.This method was adopted in APFEL: APFEL is an x-space evolution code, so the analytic truncated solution of Ref. [59] cannot be used, and the expansion is performed numerically anyway, even in the case of pure QCD.On the other hand, in the EKO evolution code the pure QCD truncated solution is implemented through the analytic solution of Ref. [59], while the numerical path-ordered solution only needs to be used for the exact solution Eqs. (2.16) and (2.17).Therefore, in order to ensure greater accuracy, we adopt the exact solution for the QED×QCD case.The impact of this choice is benchmarked in the next subsection.

B.2 Exact vs. expanded solution: impact on PDFs
We wish to quantify the impact of the choice of solution: exact, Eqs.(2.16-2.17), vs. truncated, based on an expansion of the form of Eq. (B.2).To this goal, we have determined a set of NNPDF4.0NNLO (pure QCD) PDF replicas, but now using the exact solution, instead of the truncated solution used for the published NNPDF4.0PDF sets.We have used the new theory pipeline.
The statistical indicators for these two sets of replicas are compared in Table B.1, and are seen to be indistinguishable.The corresponding PDFs are compared in Fig. B.1, where we display the distances (defined in Ref. [45]) between the central values and uncertainties of all PDFs in either set, and in Fig. B.2, where we specifically compare the gluon and antiup PDFs at Q = 100 GeV.It is clear that for x ∼ > 10 −3 the PDFs are statistically indistinguishable.At small x they start differing at the half-σ level, with differences at most reaching the one-σ level for the gluon.This is consistent with the fact that the exact and truncated solutions differ by higher-order perturbative corrections that go beyond the NNLO accuracy of the computation.Indeed, it is well known that at small x the perturbative convergence starts deteriorating because of high-energy logarithms that need resummation in order for PDF determination to be accurate [60].The qualitative behavior of the PDFs shown in Fig. B.2 agrees with this explanation: the exact solution exponentiates a set of subleading small-x logarithms which are linearized in the truncated solution.It follows that in the small-x region, where there is no data, these contributions lead to a stronger rise of the gluon, which then feeds back onto the quark-antiquark sea.

B.3 Exact solution: invertibility
The unexpanded EKO Eq. (2.17) manifestly satisfies the exact inversion property Hence, checking that evolving a PDF back and forth between two scales gives back the starting PDF provides a stringent test of the implementation of evolution equations and its accuracy.
We have performed this check, by starting with the NNPDF3.1QEDNNLO PDFs at Q a = 100 GeV, then evolving down to Q b = 1.65 GeV and back to Q a = 100 GeV.Note that this evolution crosses back and forth the bottom quark threshold, so this also checks the accuracy of the inversion of matching conditions between the N f = 4 and N f = 5 flavor schemes.Results are displayed in Fig. B.3, where we compare the initial and final PDFs.Differences are at most at the permille level, except at very large x where they can reach the percentage level, but PDFs are becoming rapidly very small.We conclude that exact evolution as implemented in the EKO code has (at least) this permille accuracy.

B.4 Benchmarking of QCD×QED evolution: EKO vs. APFEL
We finally compare the implementation of QCD×QED evolution in the EKO code used in this paper, with the APFEL code used for the NNPDF3.1QEDPDF determination.The comparison is performed by taking as input the NNPDF3.1QEDNNLO PDF set at the initial parametrization scale of Q 0 = 1.65 GeV, evolving to Q = 100 GeV and determining the percentage difference for all PDFs.
Results are shown in Fig. B.4 for three pairwise comparison.First, we compare APFEL exact vs. truncated evolution (green curve); the APFEL truncated result is the published NNPDF3.1 PDF set, as given by public LHAPDF grids.Then we compare EKO vs. APFEL evolution using the exact solution in both cases, and the default settings of either code as respectively used in this work and for the NNPDF3.1QEDPDF set (blue curve).These settings differ in the running of the couplings: in the APFEL settings the coefficients β (2,1)  QCD , β (0,3) QED and β (1,2)  QED in Eqs.(2.3) and (2.4) are neglected, i.e., the two equations are decoupled and a em runs at leading order.Finally, we compare EKO vs. APFEL evolution using the same solution and the same settings (red curve), namely exact solution and running of the couplings as in APFEL.Note that the scale on the y axis is logarithmic, and that it shows percentage difference (so 10 −3 denotes a relative difference of 10 −5 ).

Figure 3
Figure 3.1.Top left: the photon PDF at Q = 100 GeV in NNPDF4.0QEDNNLO compared to its NNPDF3.1QEDcounterpart as a ratio to the central value of the former.Top right: the relative PDF uncertainties on the photon PDF in these two determinations.Bottom: same as top panels now comparing NNPDF4.0QED with MSHT20QED[23] and CT18QED[24] (all NNLO).Bands correspond to 1σ uncertainties.

Figure 3 . 3 .
Figure 3.3.Comparison of PDFs in the NNPDF4.0QEDand the NNPDF4.0(pure QCD) sets, shown as a ratio to the former at Q = 100 GeV.Bands correspond to 1σ uncertainties.From left to right and from top to bottom, the up, anti-up, down, anti-down, strange, and gluon PDFs are shown.

Figure 3 . 4 .
Figure 3.4.Left: the percentage momentum fraction carried by the photon PDF γ(x, Q 2 ) in NNPDF4.0QEDand in NNPDF3.1QED as a function of the scale Q, where the bands indicate 1σ uncertainties.Right: same for the momentum fraction carried by the gluon PDF in NNPDF4.0QEDand in NNPDF4.0 (pure QCD).

3 . 3 - 3 . 4 .
The photon-induced contributions for NNPDF4.0QEDNNLO are compared in Fig.4.1 to their counterpart in NNPDF3.1QED(top) and in MSHT20QED and CT18QED (bottom).The level of agreement is high and directly follows from that seen between the respective photon PDFs in Figs.3.1

Figure 4 . 1 .
Figure 4.1.The photon-induced contributions to the luminosity at the LHC with √ s = 14 TeV as a function of the invariant mass m X for NNPDF4.0QEDNNLO compared to its NNPDF3.1QEDcounterpart (top) and compared to MSHT20QED and CT18QED, all shown as a ratio to NNPDF4.0QED.

Figure 4 . 2 .
Figure 4.2.Comparison of parton luminosities in the NNPDF4.0QEDand the NNPDF4.0(pure QCD) sets, shown as a ratio to the former for the LHC with √ s = 14 TeV.From left to right and from top to bottom the gluon-gluon, gluon-quark, quark-antiquark and quark-quark luminosities are shown.

Figure 4 . 4 .
Figure 4.4.Predictions for inclusive Drell-Yan production at the LHC with center-of-mass energy √ s = 14 TeV, computed at NLO accuracy in the QCD and electroweak couplings.All uncertainties shown are PDF uncertainties only.From top to bottom: neutral-current dilepton production as a function of the dilepton invariant mass m ℓ l; W + production as a function of the antilepton pseudo-rapidity ηl; W − production as a function of the antilepton pseudorapidity η l .For each process, we display the absolute distributions (top panels) and the ratio to the central value obtained with NNPDF4.0QEDPDFs, including photon-initiated channels (bottom panels).In the left panels the full NNPDF4.0result is compared to NNPDF4.0 (QCD only), NNPDF3.1QED,NNPDF4.0QED(no photon-initiated); in the right panels it is compared to MSHT20QED[23] and CT18QED[25].Note that the experimentally measurable quantity is shown, so for dilepton production the t-and u-channel photon-induced contribution is also included.

Figure 4 . 5 .
Figure 4.5.Same as Fig. 4.4 but for the rapidity distribution of the Higgs for production in gluon-gluon fusion (top panel), in association with a W + boson (middle panel) and in vector-boson fusion (bottom panel).

Figure 4 . 6 .
Figure 4.6.Same as Fig. 4.4 but for the dilepton transverse momentum distribution for weak boson pair production (W + W − and W + Z) and for the invariant mass distribution for top-quark pair production.No acceptance cuts on the decay productions of the W, Z bosons and top quark have been imposed.

Figure A. 1 .
Figure A.1.Distances at Q = 100 GeV between the central values (left) and uncertainties (right) of the 100 NNPDF4.0PDF replicas at NLO (top) and NNLO (bottom) whose statistical indicators are compared in Table A.1.
. The PDFs themselves are shown in Figs.A.2-A.3, and the parton luminosities for the LHC at √ s = 14 TeV are shown in Figs.A.4-A.5,

Figure A. 2 .
Figure A.2.The NLO PDFs at Q = 100 GeV from the 100 replica sets compared in Table A.1 and Fig. A.1.Results are shown normalized to the central value of the published set.Bands correspond to 1σ uncertainties.From left to right and from top to bottom, we show the up, anti-up, down, anti-down, strange and gluon PDFs.

Figure A. 4 .
Figure A.4. Same as in Fig. A.2 for the partonic luminosities at the LHC with √ s = 14 TeV as a function of the invariant mass m X .The gg (left) and qg (right) luminosities are shown in the top row, the q q (left), and qq (right) luminosities in the bottom row.

Figure A. 5 .
Figure A.5. Same as in Fig. A.4 at NNLO.

Figure B. 1 .
Figure B.1.Statistical distances between the central values (left) and uncertainties (right) of the 100 NNPDF4.0PDF replicas whose statistical indicators are compared in Table B.1.

Figure B. 2 .
Figure B.2.The gluon (left) and anti-up quark (right) PDFs at Q = 100 GeV from the PDF sets compared in Table B.1 and Fig. B.1.Results are shown normalized to the central value of the PDFs obtained with truncated evolution (published NNPDF4.0).

Table 3 .
1. Statistical estimators for NNPDF4.0QEDNLO and NNLO, compared to NNPDF4.0 pure QCD.From top to bottom: total χ 2 per number of data points, average over replicas of the training and validation figures of merit ⟨E tr ⟩ rep and ⟨E val ⟩ rep , average χ 2 over replicas χ 2 rep , χ 2 for datasets grouped by process.The total number of data points is 4424 (4616) at NLO (NNLO).
Table B.1.Same as Table