Precise predictions for V+jets dark matter backgrounds

High-energy jets recoiling against missing transverse energy (MET) are powerful probes of dark matter at the LHC. Searches based on large MET signatures require a precise control of the $Z(\nu\bar\nu)+$jet background in the signal region. This can be achieved by taking accurate data in control regions dominated by $Z(\ell^+\ell^-)+$jet, $W(\ell\nu)+$jet and $\gamma+$jet production, and extrapolating to the $Z(\nu\bar\nu)+$jet background by means of precise theoretical predictions. In this context, recent advances in perturbative calculations open the door to significant sensitivity improvements in dark matter searches. In this spirit, we present a combination of state-of-the art calculations for all relevant $V+$jets processes, including throughout NNLO QCD corrections and NLO electroweak corrections supplemented by Sudakov logarithms at two loops. Predictions at parton level are provided together with detailed recommendations for their usage in experimental analyses based on the reweighting of Monte Carlo samples. Particular attention is devoted to the estimate of theoretical uncertainties in the framework of dark matter searches, where subtle aspects such as correlations across different $V+$jet processes play a key role. The anticipated theoretical uncertainty in the $Z(\nu\bar\nu)+$jet background is at the few percent level up to the TeV range.

Abstract High-energy jets recoiling against missing transverse energy (MET) are powerful probes of dark matter at the LHC. Searches based on large MET signatures require a precise control of the Z(νν)+ jet background in the signal region. This can be achieved by taking accurate data in control regions dominated by Z( + − )+ jet, W ( ν)+ jet and γ+ jet production, and extrapolating to the Z(νν)+ jet background by means of precise theoretical predictions. In this context, recent advances in perturbative calculations open the door to significant sensitivity improvements in dark matter searches. In this spirit, we present a combination of state-of-the-art calculations for all relevant V + jets processes, including throughout NNLO QCD corrections and NLO electroweak corrections supplemented by Sudakov logarithms at two loops. Predictions at parton level are provided together with detailed recommendations for their usage in experimental analyses based on the reweighting of Monte Carlo samples. Particular attention is devoted to the estimate of theoretical uncertainties in the framework of dark matter searches, where subtle aspects such as correlations across different V + jet processes play a key role. The anticipated theoretical uncertainty in the Z(νν)+ jet background is at the few percent level up to the TeV range.

Introduction
The signature of missing transverse energy (MET) is one of the most powerful tools in the interpretation of data from hadron colliders. In the Standard Model (SM), MET arises from the neutrinos from the decay of W and Z bosons, and it can be used in their identification and study, as well as in the identification and study of Higgs bosons, top quarks and other SM particles whose decay products include W or Z bosons. But MET is also an almost omnipresent feature of theories beyond the SM (BSM), where it can be associated to the decay of new particles to W and Z bosons, or directly to the production of new stable, neutral and weakly interacting particles. Typical examples are theories with dark matter (DM) candidates, or Kaluza-Klein theories with large extra dimensions. Depending on the details, MET is accompanied by other modeldiscriminating features, such as the presence of a small or large multiplicity of hard jets, or of specific SM particles. The experimental search for these extensions of the SM relies on a proper modeling of the SM backgrounds to the MET signature. The determination of these backgrounds is ideally done by using data control samples, but theoretical input is often helpful, or even necessary, to extend the experimental information from the control to the signal regions, or to extend the application range of the background predictions and to improve their precision [1][2][3].
In this paper we focus on the theoretical modeling of the SM V + jet backgrounds to inclusive production of large MET recoiling against one or more hadronic jets. These final states address a broad set of BSM models, where the production of an otherwise invisible final state is revealed by the emission of one or more high-p T jets from initial state radiation, where p T is the momentum in the transverse plane. 1 Recent publications by ATLAS [5] and CMS [6,7], relative to LHC data collected at √ s = 13 TeV, document in detail the current experimental approaches to the background evaluation. The leading background is Z(νν)+ jet production, followed by W ( ν)+ jet (in particular for = τ or when the lepton is outside of the detector). 2 The experimental constraints on Z(νν)+ jet production at large MET can be obtained from accurate measurements of V + jet production processes with visible vector-boson signatures. It is quite obvious, for example, that the measurement of Z( + − )+ jets with = e, µ is the most direct and reliable proxy for Z(νν)+ jets. This control sample, however, is statistics limited, due to the smaller branching ratio of Z bosons to charged leptons relative to neutrinos. To extrapolate the shape of the Z spectrum to the largest p T values, therefore, requires a theoretical prediction. The larger statistics of W ( ν)+ jets and γ+ jets events makes it possible to directly access the relevant p T range, but the relation between their spectra and the Z spectrum needs, once again, theoretical guidance.
To put things into a concrete perspective, Figure 1 shows the expected event rates, and the relative statistical uncertainty, for 300 fb −1 of integrated luminosity at 13 TeV. The extrapolation to the O(100 fb −1 ) and O(3000 fb −1 ) expected from the full run 2 and at the end of the full LHC programme, respectively, is straightforward. The Z( + − )+jets data allow for a direct estimate of the Z(νν)+jets rate with a statistical precision below 1% for p T up to about 600 GeV. Using the W ( ν)+jets or γ+jets data could in principle extend this range up to about 900 GeV. Beyond this value, the statistical precision of the W ( ν)+jets and γ+jets events remains a factor of two better than that of the Z(νν)+jets signal. In order to ensure that the theoretical systematics in the extrapolation from the W +jets and γ+jets rates to the Z+jets rates remains negligible with respect to the statistical uncertainty, the former should be kept at the level of a few percent up to p T ∼ 2 TeV, and around 10% up to 1 For a recent comprehensive review of DM models leading to this class of signatures, see e.g. [4]. 2 Other backgrounds (such as QCD multijets, tt or pairs of gauge bosons) are suppressed, and their contribution to the overall uncertainty is well below the percent level. L = 300 fb −1 Z(νν)+ jet Z(ℓ + ℓ − )+ jet W(ℓν)+ jet γ+ jet 1 p T ∼ 2.5 − 3 TeV, which is the ultimate kinematic reach for the Z(νν)+jets signal at the end of LHC data taking.
The main result of this work is to prove that, thanks to the recent theoretical advances, these goals can be met. This proof requires the analysis of a series of possible effects. On the one hand, the theoretical extrapolation to larger p T of the very precise Z( + − )+jets data requires firm control over the shape of the distribution. Several effects, from the choice of parton distribution functions (PDFs) to the choices made for the renormalisation and factorisation scales used in the calculations, can influence the extrapolation. On the other hand, the level of correlation between the W , γ and Z spectra must be kept under control. At large p T , in particular, large and process-dependent corrections arise due to the growth of the electroweak (EW) corrections, and these may spoil the correlation induced by pure QCD effects. For our analysis we shall use the most up-to-date theoretical predictions available today for the description of vector boson production at large p T . On the QCD side, we rely on the next-to-next-to-leading or-der (NNLO) calculations, which appeared recently for Z+jet [8][9][10][11][12], W +jet [13,14] and γ+jet [15,16] production. On the EW side, we apply full NLO calculations for Z+jet [17][18][19], W +jet [19,20] and γ+jet [21] production with off-shell decays of the Z and W bosons. Given the strong enhancement of EW Sudakov effects in the TeV region, we also include 2-loop logarithmic terms at next-to-leading logarithmic (NLL) accuracy for all V + jet processes [22][23][24][25]. An extensive assessment and discussion of the estimates of missing higherorder terms, and of the relative systematics, is given in the main body of this paper. In particular, in order to address non-trivial issues that arise in the context of dark matter searches, we introduce a global framework for the estimate of theoretical uncertainties in all V + jet processes, taking into account correlation effects across different processes and p T regions. Also the uncertainties associated with the combination of QCD and EW corrections are discussed in detail.
From the experimental perspective, the determination of the background composition in signal and control regions, and the modeling of other key aspects of experimental analyses (e.g. lepton identification and reconstruction, missing energy, etc.) require a theoretical description of the various V + jets processes at the particle level. Typically, this is provided by Monte Carlo (MC) samples based on multi-jet merging at LO or NLO QCD, and improvements based on higher-order theoretical calculations can be implemented through reweighting of MC events. For the fit of MC predictions to data, ATLAS and CMS analyses rely on the profile likelihood approach, where experimental and theoretical uncertainties are described in terms of nuisance parameters with Gaussian distributions. In this context, the correlations of theoretical uncertainties across p T bins (shape uncertainties) and across different V + jets processes play a key role for searches at large MET.
For the implementation of higher-order QCD and EW corrections and for the estimate of theoretical uncertainties in the experimental analysis framework, we propose a procedure based on a one-dimensional reweighting of MC samples. The proposed framework should enable the experiments to carry out their profile likelihood approach, quantifying the impact of the theoretical systematics in their analyses, and validating directly with data the reliability and robustness of the theoretical inputs. In this respect, we would like to stress that, independently of the application to BSM searches, the results in this paper provide a framework for incisive validations of the theoretical calculations. Furthermore, these results might allow for further constraints on PDFs [3,26].
If the experimental analyses of the MET+jets channel should confirm the usefulness of the approach we propose, the same framework could be adapted to more complex or exclusive final states, in which for example MET is accompanied by a large number of (hard) jets or by specific objects (photons, heavy quarks, Higgs, etc). These extensions are left for future studies.
The structure of this paper is as follows: In Section 2 we introduce the reweighting technique, to incorporate in a MC analysis the effect of higher-order corrections and of their systematic uncertainties including correlations. Section 3 describes details of the setup for our numerical calculations, the employed tools and methods, as well as the detailed definition of physics objects and observables to be used in the context of MC reweighting. In Section 4 we discuss higher-order QCD and EW corrections, including the contribution of photoninitiated processes and real vector boson emission. We present here our approach to the estimate of the various systematics, covering QCD scale, shape and processdependent uncertainties, as well as uncertainties arising from higher-order EW and mixed QCD-EW corrections. Section 5 contains our summary and conclusions. As detailed in Appendix A, results for all V +jets processes are available in form of one-dimensional histograms in the vector-boson p T covering central predictions and all mentioned uncertainties. Technical plots on the individual sources of QCD and EW uncertainties are documented in Appendix B.

Reweighting of Monte Carlo samples
The reweighting of MC samples is an approximate, but straightforward and easy to implement method of combining (N)LO MC simulations with (N)NLO QCD+NLO EW perturbative calculations and to account for the respective uncertainties in a systematic way. The following formula describes the one-dimensional reweighting of MC samples for V + jet production (V = γ, Z, W ± ) in a generic variable x, .
In the case at hand, i.e. V + jet production, the onedimensional parameter x should be understood as the vector-boson transverse momentum, x = p (V ) T , while y generically denotes the remaining variables of the fully differential kinematic dependence of the accompanying QCD and QED activity, including both extra jet and photon radiation, as well as leptons and neutrinos from hadron decays. It is implicitly understood that d dx d dy σ depends on x and y, while in d dx σ the variables y are integrated out.
The labels MC and TH in Eq. (1) refer to Monte Carlo and higher-order theoretical predictions, respectively, and the related uncertainties are parametrised through nuisance parameters ε TH , ε MC . Our recommendations for theory uncertainties in Section 4 are formulated in terms of intervals for the related nuisance parameters, which pragmatically should be understood as the 1σ range of Gaussian uncertainties. Monte Carlo uncertainties, described by ε MC , must be correlated in the numerator and denominator on the r.h.s. of Eq. (1), while they can be kept uncorrelated across different processes (apart from Z(νν) + jet and Z( + − ) + jet).
We note that, as opposed to an approach based only on ratios of p T distributions, where theory is used for extrapolations across different processes at fixed p T , MC reweighting is more powerful as it supports all possible extrapolations across different processes and p T regions. In particular, it makes it possible to exploit V + jet precision measurements at moderate p T in order to constrain Z(νν) + jet production in the TeV region.
A further advantage of the reweighting approach (1) lies in the fact that the three terms on the r.h.s. of Eq. (1) do not need to be computed with the same numerical setup (parameters, cuts, observables, etc.). More precisely, only the definition of the variable x and the binning of its distribution need to be the same in all three terms. Scale choices, QCD and EW input parameters and PDFs should be the same only in the numerator and denominator of but can be chosen in a different way in d dx σ TH , provided that QCD and EW corrections themselves are computed using the same settings. Vice versa, possible cuts must be identical only in the numerator and denominator of For an optimal combination of higher-order calculations and MC predictions, two conditions should be fulfilled. On the one hand, theory calculations should describe the distribution in the reweighting variable with higher (or at least equal) precision as compared to the MC sample, On the other hand, the MC sample should be more accurate than TH calculations in describing the correlation between x and all other variables y, More precisely, condition (6) needs to be fulfilled only for those aspects of V + jet events that are relevant for the actual experimental analysis.
As concerns the first condition, we note that, depending on the choice of the observable x, using stateof-the-art theory calculations that involve higher-order QCD and EW corrections may not guarantee that Eq. (5) is fulfilled. In fact, there are a number of aspects, i.e. resolved multi-jet emissions, the resummation of soft logarithms in the region of small vector-boson p T , soft QCD radiation of non-perturbative origin, multiple photon radiation, or neutrinos and charged leptons resulting from hadron decays, for which fixed-order perturbative calculations of pp → V + jet are less accurate than MC simulations.
Thus, the reweighting variable x should be defined such as to have minimal sensitivity to the abovementioned aspects. In this respect, due to its reduced sensitivity to multiple jet emissions, the vector-boson p T is a natural choice. However, in order to fulfil Eq. (5), the region p (V ) T M V should be excluded from the reweighting procedure, unless QCD Sudakov logarithms are resummed to all orders in the theoretical calculations. Moreover, in order to simultaneously fulfil conditions (5) and (6), any aspect of the reconstructed vector-boson p T that is better described at MC level should be excluded from the definition of x and included in y. This applies, as discussed in Section 3, to multiple photon emissions off leptons, and to possible isolation prescriptions for the soft QCD radiation that surrounds leptons or photons. In general, purely nonperturbative aspects of MC simulations, i.e. MPI, UE, hadronisation and hadron decays, should be systematically excluded from the definition of the reweighting variable x. Thus, impact and uncertainties related to this non-perturbative modelling will remain as in the original MC samples.
It should be stressed that the above considerations are meant for dark-matter searches based on the inclusive MET distribution, while more exclusive searches that exploit additional information on hard jets may involve additional subtleties. In particular, for analyses that are sensitive to multi-jet emissions, using the inclusive vector-boson p T as the reweighting variable would still fulfil Eq. (5), but the lack of QCD and EW corrections to V + 2 jets production in MC simulations could lead to a violation of Eq. (6). In analyses that are sensitive to the tails of inclusive jet-p T and H T distributions this issue is very serious, and QCD+EW corrections should be directly implemented at MC level using multi-jet merging [19].
In general, as a sanity check of the reweighting procedure, we recommend verifying that, for reasonable choices of input parameters and QCD scales, (N)NLO QCD calculations and (N)LO merged MC predictions for vector-boson p T distributions are in reasonably good agreement within the respective uncertainties. Otherwise, in case of significant MC mismodelling of the d dx σ (V ) distribution, one should check the reliability of the MC in extrapolating TH predictions from the reweighting distribution to other relevant observables.
In general, one could check whether the onedimensional reweighting via the variable x in Eq. (1) can in fact reproduce the dependence of the corrections in other kinematic variables that are relevant for the experimental analysis. To this end, distributions of σ (V ) w.r.t. another kinematic variable x should be calculated upon integrating Eq. (1). Switching on and off the corrections on the r.h.s. of Eq. (1) in σ (V ) TH and taking the ratio of the obtained differential cross sections σ (V ) , produces the relative correction to the x distribution that could be compared to the corresponding result directly calculated from σ (V ) TH . 3 Finally, it is crucial to check that state-of-the-art predictions for absolute dσ/dp T distributions agree with data for the various visible final states.

Setup for numerical predictions
In this section we specify the physics objects (Section 3.1), acceptance cuts and observables (Section 3.2), input parameters (Section 3.3) and tools (Section 3.4) used in the theoretical calculations for pp → W ± /Z/γ+ jet.
The definitions of physics objects, cuts and observables-which specify the setup for the reweighting procedure discussed in Section 2-should be adopted both for theoretical calculations and for their Monte Carlo counterpart in the reweighting factor (3). The details of the reweighting setup are designed such as to avoid any possible deficit in the perturbative predictions (e.g. due to lack of resummation at small p T ) and any bias due to non-perturbative aspects of Monte Carlo simulations (e.g. leptons and missing energy from hadron decays). Let us also recall that this setup is completely independent of the physics objects, cuts and observables employed in the experimental analyses.
As concerns input parameters and PDFs, the recommendation of Section 3.3 should be applied to all QCD and EW higher-order calculations. In particular, it is mandatory to compute (N)NLO QCD and EW corrections in the same EW input scheme, otherwise NLO EW accuracy would be spoiled. Instead, Monte Carlo simulations and the corresponding d dx σ MC contributions to the reweighting factor (3) do not need to be based on the same input parameters and PDFs used for theory predictions.
We recommend handling W/Z + jet production and decay on the Monte Carlo side as the full processes pp → / ν/νν + jet, i.e. with a consistent treatment of off-shell effects, as is done on the theory side.

Definition of physics objects
In the following we define the various physics objects relevant for higher-order perturbative calculations and for the reweighting in the Monte Carlo counterparts in Eq. (3).

Neutrinos
In parton-level calculations of pp → / ν/νν +jet, neutrinos originate only from vector-boson decays, while in Monte Carlo samples they can arise also from hadron decays. In order to avoid any bias in the reweighting procedure, only neutrinos arising from Z and W decays at Monte Carlo truth level should be considered.

Charged leptons
Distributions in the lepton p T and other leptonic observables are known to be highly sensitive to QED radiative corrections, and the differences in the treatment of QED radiation on Monte Carlo and theory side can lead to a bias in the reweighting procedure. This should be avoided by using dressed leptons, i.e. recombining all leptons with nearly collinear photons that lie within a cone For the radius of the recombination cone we employ the standard value R rec = 0.1, which allows one to capture the bulk of the collinear final-state radiation, while keeping contamination from large-angle photon radiation at a negligible level. All lepton observables as well as the kinematics of reconstructed W and Z bosons are defined in terms of dressed leptons, and, in accordance with standard experimental practice, both muons and electrons should be dressed. In this way differences between electrons and muons, = e, µ, become negligible, and the reweighting function needs to be computed only once for a generic lepton flavour . Similarly as for neutrinos, only charged leptons that arise from Z and W decays at Monte Carlo truth level should be considered. Concerning QCD radiation in the vicinity of leptons, no lepton isolation requirement should be imposed in the context of the reweighting procedure. Instead, in the experimental analysis lepton isolation cuts can be applied in the usual manner.

Z and W bosons
The off-shell four-momenta of W and Z bosons are defined as where the leptons and neutrinos that result from Z and W decays are defined as discussed above.

Photons
At higher orders in QCD, photon production involves final-state q → qγ splittings that lead to collinear singularities when QCD radiation is emitted in the direction of the photon momentum. Since such singularities are of QED type, they are not cancelled by corresponding virtual QCD singularities. Thus, in order to obtain finite predictions in perturbation theory, the definition of the pp → γ+ jet cross section requires a photonisolation prescription that vetoes collinear q → qγ radiation while preserving the cancellation of QCD infrared singularities.
To this end, in this study we adopt Frixione's isolation prescription [27], which limits the hadronic transverse energy within a smooth cone around the photon by requiring where the sum runs over all quarks/gluons and hadrons at parton level and Monte Carlo level, respectively, while p T,i and p T,γ denote the transverse momenta of partons/hadrons and photons. The p T -fraction ε 0 , the cone size R 0 , and the exponent n are free parameters that allow one to control the amount of allowed QCD radiation in the vicinity of the photon. The photon-isolation prescription is applicable to QCD as well as to EW higher-order corrections. At NLO EW, γ+ jet production involves bremsstrahlung contributions with two final-state photons. In this case, at least one isolated photon is required. The other photon might become soft, guaranteeing the cancellation of related soft and collinear singularities in the virtual EW corrections. In case of two isolated photons in the final state, the hardest photon is considered. In particular, an explicit photon isolation prescription is mandatory at NLO EW in order to prevent uncancelled singularities from q → qγ splittings in the O(α 2 α S ) mixed EW-QCD contributions from qq → qqγ and crossing-related channels.
As a consequence of q → qγ collinear singularities and the need to apply a photon isolation prescription, QCD corrections to pp → γ+ jet behave differently as compared to Z/W + jet production. Such differences can be important even at the TeV scale, where one might naively expect that massive and massless vector bosons behave in a universal way from the viewpoint of QCD dynamics. Instead, the presence of collinear q → qV singularities at (N)NLO QCD implies a logarithmic sensitivity to the vector-boson masses, which results, respectively, in ln(R 0 ) and ln(p T,V /M V ) terms for the case of massless and massive vector bosons at p T,V M W,Z . A quantitative understanding of these differences and their implications on the correlation of QCD uncertainties between γ+ jet and Z+ jet production is crucial for the extrapolation of γ+ jet measurements to Z+ jet dark-matter backgrounds. To this end, as discussed in Section 4, we propose a systematic approach based on the idea that, at large p T,V , the pp → γ+ jet process can be split into a dominant part with universal QCD dynamics (in the sense that QCD effects in γ+ jet and Z/W + jet production are strongly correlated) and a remnant contribution that has to be handled as uncorrelated in the treatment of QCD uncertainties. To achieve this, we introduce a modified photon isolation prescription, which is designed such as to render the QCD dynamics of γ+ jet and Z/W + jet production as similar as possible at high p T . To this end we define a dynamic cone radius which is chosen in such a way that the invariant mass of a photon-jet pair with R γj = R dyn and p T,j = ε 0 p T,γ corresponds to the Z-boson mass, i.e.
where the first identity is valid in the small-R approximation. In this way, using a smooth isolation with R 0 = R dyn (p T,γ , ε 0 ) mimics the role of the Zand Wboson masses as regulators of collinear singularities in Z/W +jet production at high p T , while using a fixed cone radius R 0 would correspond to an effective M γj cut well beyond M Z,W , resulting in a more pronounced suppression of QCD radiation in γ+ jet production as compared to Z/W + jet. Specifically, as default photon selection for the theoretical predictions 4 in this study we use the dynamic cone isolation defined through Eq. (10) and Eq. (11), with parameters ε 0,dyn = 0.1, n dyn = 1, Note that, in order to prevent that the veto against collinear QCD radiation is applied to an excessively large region of phase space, the dynamic cone radius in Eq. (13) is limited to R dyn ≤ 1.0. As a result of this upper bound, for p T,γ < M Z ε −1/2 0,dyn 290 GeV the cone radius is kept fixed, and the impact of collinear QCD radiation starts to be significantly enhanced as compared to the case of Z/W + jet production. Vice versa, for p T,γ > M Z ε −1/2 0,dyn , thanks to the dynamic isolation cone (13), QCD effects in γ+ jet and Z/W + jet production become closely related, and the degree of correlation between QCD uncertainties across all V + jet processes can be described with the prescription of Eqs. (37)- (38).
For a realistic assessment of theoretical uncertainties, one should also consider the fact that photon isolation prescriptions used in experimental analyses differ in a significant way from the dynamic prescription of Eq. (13). To this end, we recommend to repeat the reweighting procedure using theory predictions for γ+ jet based on a standard Frixione isolation (10) with fixed cone radius and parameters that mimic typical experimental selections at particle level [28], The difference between γ+ jet MC samples reweighted in the dynamic-and fixed-cone setup should be taken as an additional uncertainty for pp → γ+ jet. As ingredients for this uncertainty estimate we provide higherorder QCD predictions (without uncertainties) with fixedcone isolation (14) besides the full set of pp → γ+ jet 4 The same isolation prescription used for theory predictions should be applied also to their MC counterparts dσ MC /dx in the context of the reweighting procedure.
predictions and uncertainties with dynamic photon isolation (see Appendix A). In the EW corrections, differences between the two photon isolation prescriptions are well below the percent level. Thus predictions for γ+ jet at (n)NLO EW are provided only with the dynamic cone prescription of Eq. (13). In Figure 2 we present a comparison of the NLO QCD K-factors for W/Z+ jet and γ+ jet production with dynamic and fixed cone isolation. For p T,γ < 290 GeV, where both isolation prescriptions correspond to a fixed cone radius, the QCD corrections to pp → γ+ jet grow rapidly with decreasing p T . At low p T , due to the smaller cone size, fixed isolation (R 0 = 0.4) leads to more pronounced corrections as compared to dynamic isolation (R 0 = 1.0), but the slopes of the corresponding γ+ jet K-factors are quite similar to each other and very different as compared to the ones for pp → W/Z+ jet. In the case of fixed isolation, this difference persists also in the high-p T regime (apart form the accidental agreement of K-factors at p T,V ≈ 800 GeV). Instead, in the case of dynamic photon isolation, at large p T the QCD corrections to γ+ jet and W/Z+ jet production turn out to be remarkably similar, both in shape and size. As expected, the onset of this universal behaviour is located close to p T,γ = 290 GeV, where the isolation radius R 0,dyn starts varying with p T in a way that rejects QCD radiation with M γj < ∼ M W,Z . The differences between γ+ jet and W/Z+ jet K-factors remain as small as a few percent up to the TeV scale.

QCD partons and photons inside jets
In order to avoid any bias due to the different modelling of jets in MC simulations and perturbative calculations, theory calculations and reweighting should be performed at the level of inclusive vector-boson p T distributions, without imposing any requirement on the recoiling jet(s). Predictions presented in this study are thus independent of specific jet definitions or jet cuts.
Concerning the composition of the recoil, we observe that, at NLO EW, q → qγ splittings can transfer an arbitrary fraction of the recoiling momentum from QCD partons to photons. In particular, in pp → V γj contributions of O(α 2 α S ), the photon can carry up to 100% of the recoil momentum. Such contributions involve soft QCD singularities that are cancelled by including also virtual QCD corrections to pp → V γ. In order to minimise double counting with diboson production, 5 V γ production at LO is not included in the   for W + jet, Z+ jet, and γ+ jet production with dynamic photon isolation (13) and standard fixed-cone isolation (14). On the right corresponding ratios of K-factors are shown, and the dotted lines indicate K-factor variations of ±0.05. EW corrections to pp → V j. In practice, as demonstrated in Figure 3, the relative weight of pp is well below the percent level. Thus the impact of O(α 2 α S ) contributions from hard V γ production, which are included in this study, should be completely negligible.

Cuts and observables
Theoretical calculations and the reweighting of MC samples should be performed in a fully inclusive V + jet setup, imposing a single cut with p T,W ± and p T,Z defined as in Section 3.1. The cut (15) is crucial in order to avoid the region where perturbative predictions suffer form the lack of QCD resummation. 6 For leptons and MET we do not apply any p T or rapidity cuts. Moreover, we do not impose any restrictions on QCD radiation in the vicinity of leptons and MET. Also QCD radiation is handled in a fully inclusive way, i.e. the presence of a recoiling jet is not explicitly required, and, as discussed in Section 3.1, at NLO EW the recoil can be entirely carried by a photon. Here we want to stress again that of course the particle-level analysis of the reweighted Monte Carlo samples can (and will) involve a more exclusive event selection than used for the reweighting itself.
The differential distributions to be used for the reweighting of the various pp → V + jet processes and process-specific selection cuts to be applied in addition to Eq. (15) are summarised in Table 1. In the case of pp → νν+ jet all three neutrino species are added, while for all other Z and W decays only a single lepton generation is considered. For pp → + − + jet an extra invariant-mass cut is applied in order to avoid far off-shell contributions, especially from γ * → + − at low invariant mass. The relatively low value of the lower cut, m > 30 GeV, is intended to minimise cross section loss due to photon radiation that shifts events from the Z-peak region down to lower invariant mass (see Figure 4). This choice guarantees a reduced sensitivity with respect to the modelling of QED radiation. process extra cuts observable comments Table 1: Extra selection cuts, in addition to Eq. (15), and observables for the various V + jet processes. Alternative predictions for γ+jet production are provided also for the case of a standard Frixione isolation with parameters (14).

Input parameters, PDFs and QCD scales
Input parameters and PDFs employed for theoretical predictions in this study are specified in the following. Let us recall that, as discussed in Section 2, Monte Carlo samples used in the experimental analyses do not need to be generated with the same input parameters and PDFs used for higher-order theoretical predictions.
In the calculation of pp → / ν/νν/γ + jet we use the gauge-boson masses [30] M Z = 91.1876 GeV, M W = 80.385 GeV, (17) and the corresponding widths, The latter are obtained from state-of-the-art theoretical calculations. For the top-quark [30] and Higgs-boson [31] masses and widths we use All unstable particles are treated in the complexmass scheme [32], where width effects are absorbed into the complex-valued renormalised masses For W +jet and Z+jet production processes the EW couplings are derived from the gauge-boson masses and the Fermi constant, G µ = 1.16637 × 10 −5 GeV −2 , using while for γ+jet production the EW coupling is chosen to be [30] α = α(0) = 1/137.035999074 .
In both schemes the weak mixing angle θ w is determined by and becomes complex-valued. The G µ -scheme guarantees an optimal description of pure SU(2) interactions at the EW scale. It is the scheme of choice for W + jet production, and it provides a very good description of Z + jet production as well. The α(0) scheme to be used for γ+jet, on the other hand, expresses the fact that onshell photons effectively couple at a scale Q 2 =0 . The CKM matrix is assumed to be diagonal and we checked at LO and NLO QCD that for W +jet production the difference with respect to a non-diagonal CKM matrix is always well below 1%. For the choice of renormalisation and factorisation scales and variations thereof we refer to Section 4.1.
For the calculation of hadron-level cross sections at (N)NLO QCD + (n)NLO EW we employ the LUXqed_ plus_PDF4LHC15_nnlo_100 PDF set, which is based on PDF4LHC NNLO PDFs [33][34][35][36][37][38] supplemented with QED effects [39]. The same PDF set, and the related α S value, is used throughout, i.e. also in the relevant LO and NLO ingredients used in the estimate of theoretical uncertainties. At the level of precision discussed in this study also the uncertainty on the value of α S becomes relevant. Given 1% uncertainty on the measured value of α S this results in an overall 1-2% normalisation uncertainty on the differential p T distributions. However, one should keep in mind that in the process ratios this uncertainty cancels completely and thus it is irrelevant for background estimates in DM searches at high-MET. Consistently with the five-flavour number scheme employed in the PDFs, b-quarks are treated as massless partons, and channels with initial-state b-quarks are taken into account. All light quarks, including bottom quarks, are treated as massless particles, and top-quark loops are included up to NLO throughout. Matrix elements at (N)NLO are evaluated using the five-flavour running of the strong coupling supported by the PDFs and, for consistency, top-quark loops are renormalised in the decoupling scheme. For the NNLO QCD coefficient no top-quark loops are considered.
For the assessment of PDF uncertainties the PDF4LHC prescription [33] is adopted. In addition to standard PDF variations, also additional LUXqed variations for the photon PDF are applied. For more details see more details in Sections 4.3-4.4.

Computational frameworks
The theoretical predictions presented in Section 4 include corrections up to NNLO QCD and NLO EW, as well as Sudakov EW effects at O(α 2 ). They have been obtained by means of a variety of methods and tools, as detailed in the following.
The NLO QCD and NLO EW calculations for all pp → V + jet processes have been performed with Mu-nich+OpenLoops and/or Sherpa+OpenLoops. In these automated frameworks [19,40,41] virtual amplitudes are provided by the OpenLoops program [42,43], combined with the Collier tensor reduction library [44] or with CutTools [45]. The remaining tasks are supported by the two independent and fully automated Monte Carlo generators Munich [46] and Sherpa [47][48][49][50]. Additionally, we carefully validated the NLO EW predictions against the results of Refs. [17,18,20]. The NLO EW calculations for pp → V + 2 jets performed to test the factorisation of QCD and EW corrections have been checked against the one of Ref. [51] for pp → Z + 2 jets in Ref. [21]. The NLO EW amplitudes for all V +jet processes in OpenLoops have been supplemented with the one-and two-loop analytical Sudakov logarithms of Refs. [22][23][24][25]52].
The NNLO QCD predictions for Z+jet production have been obtained with the parton-level event generator NNLOjet, which provides the necessary infrastructure to perform fully differential calculations at NNLO using the antenna subtraction formalism [53][54][55][56][57][58][59][60][61]. The computation of pp → W +jet through NNLO is based on the N -jettiness subtraction scheme for NNLO calculations [13]. The above-cut contribution within the Njettiness subtraction was obtained using Munich+OpenLoops. The NNLO QCD prediction for the pp → γ+jet process is based on the calculations of Refs. [15,16] and has been obtained using MCFM [62]. In order to ensure the correctness of the numerical implementation of cuts and other parameters in the NNLO codes, a detailed comparison has been performed at the level of the NLO QCD results as described above.

Higher-order QCD and EW predictions
Precise theory predictions for V + jet production require QCD and EW higher-order corrections, mixed QCD-EW contributions, as well as photon-induced contributions, In this section we present theoretical predictions that include corrections up to NNLO QCD and NLO EW supplemented by EW Sudakov logarithms at two loops. Moreover, we introduce a coherent theoretical framework for the combination of EW and QCD calculations for the various V + jet production processes and for the assessment of the corresponding remaining sources of theoretical uncertainty. State-of-the art QCD and EW predictions and the related theoretical uncertainties are discussed in Section 4.1 and 4.2 respectively. Section 4.3 is devoted to photon-induced channels and Section 4.4 to PDF uncertainties, while in Section 4.5 we discuss the real emission of vector bosons, and mixed corrections of O(αα S ) are addressed in Section 4.6 by means of a factorised combination of QCD and EW corrections.
To illustrate the effect of higher-order corrections and uncertainties we present a series of numerical results for pp → V + jet at a centre-of-mass energy of 13 TeV in the setup specified in Section 3. In particular, pp → γ+ jet predictions are based on the dynamic photon isolation (13). As anticipated in Section 3.1, this prescription provides a very convenient basis for the systematic modelling of the correlation of QCD uncertainties between the various V + jet production processes (see Section 4.1).
Vector-boson p T spectra are plotted starting at 80 GeV, but for the sake of a complete documentation data sets are provided above 30 GeV (see Appendix A). However, we note that in the region of p T < ∼ 100 GeV there are potential sources of systematics that we are not controlling or even discussing, as they would require a separate study. These arise from the resummation of QCD Sudakov logarithms or from non-perturbative effects (e.g. an order Λ QCD average shift of the vector boson p T associated with the asymmetry of colour flow in the final state). Furthermore, as shown later, a reliable correlation between the Z/W spectra and the photon spectrum requires p T to be large enough so that fragmentation contributions in γ+jet production become small. We also expect that in the p T regions up to a few hundred GeV the statistics are sufficient to guarantee that experimental analyses of missing-E T backgrounds can entirely rely on the direct measurement of the Z spectrum measured via Z → + − . As a result, we believe that our conclusions on the systematic uncertainties are most reliable and useful for experimental applications in the region of p T larger than 100-200 GeV.

Higher-order QCD predictions
For perturbative QCD predictions at LO, NLO and NNLO we use the generic notation with k = 0, 1 or 2. Wherever possible, nominal predictions are provided at NNLO QCD, i.e. including terms up to 8 O(αα 3 S ). However, as ingredients for the assessment of some theory uncertainties, also LO and NLO QCD contributions will be used.
For convenience, results at N k LO QCD are systematically expressed in terms of LO predictions and relative correction factors defined through We calculate all N k LO and LO cross sections with one and the same set of NNLO PDFs as discussed in Section 3.3. The dependence on the renormalisation and factorisation scales, µ = (µ R , µ F ), is absorbed into the K-factors, while LO predictions on the r.h.s. of Eq. (27) are taken at the central scale, µ 0 = (µ R,0 , µ F,0 ). For the central scale we adopt the commonly used choice where the total transverse energy,Ĥ T , is defined as the scalar sum of the transverse energy of all parton-level final-state objects, Also quarks (q), gluons (g) and photons that are radiated in the (N)NLO QCD or EW corrections are included inĤ T , and the vector-boson transverse energy, E T,V , is computed using the total (off-shell) fourmomentum of the corresponding decay products, i.e.
In order to guarantee infrared safety at NLO EW, the scale (29) must be insensitive to collinear photon emissions off charged fermions. To this end, the vector-boson transverse energies defined in Eq. (30) should be computed in terms of dressed leptons as specified in Section 3.1, while |p T,γ | contributions to Eq. (29) should involve only photons that have not been recombined with charged leptons. It is worth to note that µ 0 ≈ p T,V at large p T,V .

Pure QCD uncertainties
The uncertainty associated with the truncation of the perturbative expansion in α S is estimated by means of factorisation and renormalisation scale variations. We consider standard seven-point variations applying, respectively, factor-two rescalings, i.e.
where i = 0, . . . 6. Nominal predictions and related uncertainties are defined as the centre and the half-width of the band resulting from the above variations. In terms of K-factors this corresponds to Since the shift resulting form the symmetrisation of scale variations in Eq. (32) is encoded in the K-factors, also the LO K-factor differs from one. Constant scale variations mainly affect the overall normalisation of p T -distributions and tend to underestimate shape uncertainties, which play an important role in the extrapolation of low-p T measurements to high p T . Thus, for a reasonably conservative estimate of shape uncertainties, we introduce an additional variation, where the standard scale uncertainty (33) is supplemented by a shape distortion ω shape (x), with |ω shape (x)| ≤ 1 and ω shape (x) → ±1 at high and small transverse momentum, respectively. The function ω shape is defined as and as reference transverse momentum we choose the value p T,0 = 650 GeV, which corresponds (in logarithmic scale) to the middle of the range of interest, 0.2-2 TeV. As illustrated in Figure 5, the function ω shape (x) induces asymmetric variations that cover ±75% of the standard scale variation band for p T ∈ [250, 1750] GeV. Note that, in the combination of the uncertainties (33) and (35), our choice to have an additional shape variation augments the standard scale uncertainty by a fac- Besides shape uncertainties, also the correlation of QCD uncertainties across V + jet processes plays a key role in fits of the Z(νν)+ jet dark matter background, and the quantitative understanding of such process correlations belongs to the most important theoretical aspects in dark matter searches. From the viewpoint of QCD interactions, the processes pp → W + jet and pp → Z+ jet are quite similar to each other at p T,V M W,Z . Thus, the respective QCD uncertainties are expected to be strongly correlated. However, due to the presence of q → qγ collinear singularities and the need to suppress them with an appropriate photon-isolation prescription, higher-order QCD contributions to γ+ jet production can behave in a significantly different way as compared to the case of pp → W/Z+ jet. In order to reduce such differences, we adopt the dynamic photon isolation approach defined in Eq. (13). As discussed in Section 3.1, this prescription renders the QCD dynamics of pp → γ+ jet and pp → Z/W + jet processes almost universal. As a result, QCD K-factors K depend only very weakly 9 on V at high p T , and in this situation the small process-dependent part of QCD K-factors can be used as an estimator of the degree of correlation across processes. To this end we consider the highest available term in the perturbative expansion, and as estimate of unknown process correlation effects we take the difference of the known QCD K-factors with respect to Z+ jet production, This process correlation uncertainty can be assessed using the central scale (28) throughout. Applying it to nominal predictions, i.e. replacing K N k LO , amounts to doubling or removing K-factor differences between processes. The choice of Z+ jet production as reference process in Eq. (38) is arbitrary, but changing the reference process has very little impact on process correlations since the resulting overall shift cancels to a large extent in ratios of V + jet cross sections.
The above prescription should be regarded as conservative, since parts of the available K-factors are downgraded from the status of known higher-order corrections to uncertainties. However, thanks to the fact that the V + jet K-factors of the same order k are strongly correlated, N k LO , the resulting losses of accuracy in the nominal N k LO predictions for individual processes are rather small.
For the application to experimental analyses, it is important to keep in mind that the above modelling of process correlations assumes a close similarity of QCD effects between all pp → V + jet processes, which is achieved, in the present study, by means of the dynamic photon isolation (13). Thus, as discussed in Section 3.1, experimental analyses that employ a different photon isolation approach require an additional γ+ jet specific uncertainty.
The above uncertainties can be parametrised through a set of independent nuisance parameters, ε QCD , and combined using LO QCD (µ 0 ). The nuisance parameters ε QCD,1 , ε QCD,2 and ε QCD,3 should be Gaussian distributed with one standard deviation corresponding to the range ε QCD,i ∈ [−1, +1]. These parameters should be kept uncorrelated, but each ε QCD,i -variation should be applied in a correlated way across p T bins and processes, since correlation effects are consistently implemented in the δ (i) K

Numerical results
Predictions for V +jet distributions and their ratios at LO, NLO and NNLO QCD are presented in Figures 6-8 as well as in Figures 18-19 (see Appendix B). In Figures 7 and 18-19, scale uncertainties (33), shape uncertainties (35), and process-correlation uncertainties (38) are shown separately, while in Figures 6 and 8 the three QCD uncertainties are combined in quadrature. Here and in the following W denotes W + and W − combined.
At high transverse momentum, we find that QCD corrections and uncertainties for the various V + jet production processes behave in a very similar way. At NLO the corrections amount to 40-60% with residual uncertainties around 10-20%, while NNLO corrections increase the cross section by 5-10% and reduce the combined uncertainty to 3-10%. Scale variations δ (1) K N k LO and shape variations δ (2) K N k LO are the dominant sources of uncertainty in p T -distributions. Their contributions are very similar across V + jet processes. Thus in the ratios scale and shape variations largely cancel, and the process-correlation uncertainty δ (3) K N k LO tends to dominate.
The ratio plots ( Figure 8) allow one to appreciate small differences in the QCD dynamics of the various V + jet processes. As reflected in the Z/W ratio, the NLO and NNLO corrections for the corresponding proceses are almost identical, with differences below 1-2% up to one TeV. Only at very large p T the NLO and also NNLO corrections to W +jet grow faster than in the case of Z+jet. This results in an increase of the process-correlation uncertainty δ (3) K NLO up to about 5% beyond p T = 2 TeV.
As can be seen in the Z/γ and W/γ ratios, the higher-order QCD corrections to γ+jet production behave very similarly as for Z+ jet and W + jet production at large p T . This is the result of the dynamic photon isolation (13), which guarantees that the differences in the NLO and NNLO corrections remain below 3-4% for p T > 200 GeV. Instead, at lower p T the behaviour of γ+ jet production changes drastically due to mass effects, which results in sizeable process-correlation uncertainties. 10 Note that for p T ≈ 300 GeV the NLO process-correlation uncertainty in pp → γ+jet is accidentally very small (see Figure 18) yielding a pinch in the total QCD uncertainty for the Z/γ and the W/γ ratios (see also Figure 19). However, one should keep in mind that an additional analysis-dependent photonisolation uncertainty (see Section 3.1) has to be considered for these ratios.
In general, comparing QCD predictions at different orders we observe a good convergence of the perturbative expansion, and the fact that process ratios receive very small corrections both at NLO and NNLO provides strong evidence for the universality of QCD dynamics is all V + jet processes. Results at NNLO provide also a crucial test of the goodness of the proposed approach for the estimate of QCD uncertainties and their correlations. In particular, the remarkable consistency between NNLO and NLO predictions in Figure 8 confirms that QCD uncertainties for process ratios are as small as 1-2%.

Electroweak corrections
For EW higher-order corrections we use the notation, NLO EW denotes exact O(α 2 α S ) contributions, and 'NNLO Sud' stands for O(α 3 α S ) EW Sudakov logarithms in NLL approximation (see below). Their combination is dubbed nNLO EW as it accounts for the dominant EW effects at NNLO. While our power counting   does not consider the extra factor α associated with vector-boson decays, all predictions for pp → W/Z+ jet at (N)NLO QCD + NLO EW are at the level of the full processes, pp → ν/ /νν+ jet, including off-shell effects and NLO EW corrections in decays. Since EW Sudakov logarithms do not enter W and Z decays, they are applied only at the level of pp → V + jet production, including off-shell decays at LO. The EW corrections, similarly as for the QCD ones, are also expressed in terms of correction factors with respect to LO QCD, where EW stands for NLO EW or nNLO EW. At variance with Eq. (27), here the EW κ-factors are defined by taking the factorized LO cross section at the same QCD scales, µ = (µ R , µ F ), as in the higher-order EW prediction. In this way, since QCD scale variations at LO QCD and (n)NLO EW have almost identical impact, the relative EW correction is essentially indepen-dent of µ. Thus, in practice, κ EW can be computed at the fixed reference scale, while the scale dependence of σ LO QCD (µ) in Eq. (41). Moreover, the EW correction factor κ (V ) EW is rather insensitive to the choice of PDF set as long as it is derived from cross sections that are based on the same PDFs. Analogously to Eq. (40), nNLO EW correction factors are split into a full NLO part and an NNLO Sudakov part, At NLO EW, all relevant contributions of O(α 2 α S ) are included. In the qq channel, and in all crossingrelated channels, they comprise the following types of corrections: (a.1) virtual EW corrections to qq → V g; (a.2) qq → V gγ photon bremsstrahlung; correction to the qq → V g channel and NLO QCD correction to the γq → V q channel. As discussed in Section 4.3, given the relatively small impact of γq → V q processes at O(α 2 ), photon-induced contributions of O(α S α 2 ) will not be included in the present study; (a.6) real-boson emission, i.e. pp → V V j, contributes at O(α 2 α S ). As discussed in Section 4.5, in order to avoid double counting with diboson production, such contributions should be treated as separate background samples and not as part of the EW corrections to pp → V j.
At very high transverse momentum, EW corrections are strongly enhanced by Sudakov effects, and the inclusion of higher-order Sudakov logarithms becomes manda-tory in order to achieve few-percent level accuracy. In the high-p T regime, where all energy scales are far above the weak-boson mass scale, higher-order virtual EW corrections to hard scattering cross sections can be described by means of resummation formulas of the type 12 [63,64] dσ EW = exp where γ, χ and ξ are anomalous dimensions depending on the EW quantum numbers of the scattering particles. The hard cross section has the form and the correction factors δ yields exp . . . = 1 + α π δ At NLL level, which is the logarithmic accuracy at which NNLO Sudakov effects are known for V + jet production [22][23][24][25]52], the following types of logarithms are available, where M = M W ∼ M Z , Q 2 ij = |(p i ±p j ) 2 | are the various Mandelstam invariants built from the hard momentap i of the V + jet production process and Q 2 = Q 2 12 =ŝ. In this work we will employ the explicit NLL Sudakov results of Refs. [22][23][24][25]52], which have been implemented, in addition to exact NLO QCD+NLO EW amplitudes, in the OpenLoops matrix-element generator [19,40]. Let us recall that the results of Refs. [22][23][24][25]52] are based on the high-energy limit of virtual one-and two-loop corrections regularised with a fictitious photon mass of order M W . This generates logarithms of the form α n ln k (ŝ/M 2 W ) that correspond to the combination of virtual one-and two-loop EW corrections plus corresponding photon radiation contributions up to an effective cut-off scale of order M W . In the case of V + jet production, for physical observables that are inclusive with respect to photon radiation, this approximation is accurate at the one-percent level [21,22,25].
In this work we will employ full EW results at NLO and NLL Sudakov logarithms at NNLO. In the notation of Eqs. (41)-(43), for fully-differential partonic cross sections, this implies

Pure EW uncertainties
Assuming that the NLL Sudakov approximation at NNLO is comparably accurate as at NLO, we can consider unknown Sudakov logarithms beyond NNLO as the dominant source of EW uncertainty at high p T . Such Sudakov terms of relative O(α 3 ) can be easily estimated via naive exponentiation, which implies the following relations between NLO, NNLO and NNNLO terms, Sud .
Based on these relations, we estimate the uncertainty due to unknown high-p T EW effects beyond NNLO as which is an approximate implementation of Eq. (51), obtained by neglecting effects from angular integration, replacing δ Sud by the full NLO EW correction, and multiplying the term δ Sud by a factor two, in order to be conservative.
Besides Sudakov exponentiation effects, we introduce a second source of uncertainty, defined, at nNLO EW level, as 5% of the absolute full NLO EW correction, In general, the non-Sudakov factor κ NLO hard = ( α π )δ (1) hard can amount to several percent, e.g. due to photonbremsstrahlung effects in highly exclusive observables. However, for the boson-p T distributions considered in this paper, where dressed leptons are used, the quality of the Sudakov approximation observed in Figure 9 indicates that κ NLO hard is very small. Nevertheless, to be conservative, in Eq. (53) we choose a prefactor that allows for effects as large as κ NLO hard = 5%. As a second motivation, the uncertainty (53) which corresponds to δ hard . This limit should be conservative enough to hold also in situations where the NLO hard correction is accidentally small with respect to its NNLO counterpart.
In order to account for the limitations of the Sudakov approximation at nNLO in a sufficiently conservative way, we introduce an additional source of uncertainty defined as the difference between the rigorous NLL Sudakov approximation (50) and a naive exponentiation of the full NLO EW correction, This expression provides an estimate of the typical size of terms of type δ Sud . In correspondence to the nNLO uncertainties of Eqs. (52), (53) and (56), at NLO EW we introduce uncertainties δ (i) κ NLO EW , defined as Here the first term is the direct transposition of Eq. (52) to NLO. It accounts for the unknown O(α 2 ) Sudakov terms δ Sud in Eq. (51) supplemented with an extra factor of two. As explained in the following, the second uncertainty in Eq. (57) is the NLO counterpart of the nNLO EW uncertainty (53). The latter accounts for unknown O(α 2 ) terms of type (54) and (55), which correspond to the intrinsic uncertainty of the employed Sudakov approximation at nNLO. At NLO EW the situation is different, since the calculations are exact, i.e. there are no unknown terms of O(α). Thus, we assume an uncertainty δ (2) hard . We do not consider additional uncertainties of type α π 2 δ (2) Sud since they are already covered by the first term in Eq. (57). As estimate of the size of the unknown δ (2) hard coefficient, following the discussion of Eq. (55), we impose a very generous upper bound to the ratio between δ (2) hard and δ (1) hard . To be conservative, at NLO EW we adopt a ten times looser bound as compared to nNLO EW, i.e. we require δ Similarly as for QCD uncertainties, the EW uncertainties in Eqs. (52), (53), (56) and Eq. (57), can be parametrised in terms of nuisance parameters ε EW and combined via where EW stands for NLO EW or nNLO EW. The nuisance parameters ε (V ) EW,i should be Gaussian distributed with one standard deviation corresponding to the range ε (V ) EW,i ∈ [−1, +1], and their variations should be applied in a correlated way across p T -bins. Since the first uncertainty (52) reflects the universal exponentiation properties of Sudakov EW corrections, which permits to predict the magnitude and size of the dominant higherorder corrections for each individual processes, this variation should be correlated across processes, i.e. a single nuisance parameter should be used, In contrast, the remaining EW uncertainties (53) and (56) describe subleading NNLO effects whose sign, magnitude and process dependence are unknown. Thus these uncertainties should be treated as uncorrelated, i.e. independent nuisance parameters ε

Numerical results
Predictions for V +jet distributions and their ratios at LO, NLO EW and nNLO EW are presented in Contrary to the case of QCD corrections, higherorder EW effects have a significant impact on the shapes of p T distributions as well as a pronounced dependence on the scattering process. This behaviour is mainly due to the p T dependence of EW Sudakov logarithms and their dependence on the SU(2) charges of the produced vector bosons.
As can be seen in Figure 9, the vector-boson p T spectra receive negative EW corrections that grow with p T and become very sizable in the tails. At the TeV scale, NLO EW effects reach 20-50% for Z+jet and W +jet production, and 10-15% for γ+jet production. As expected from exponentiation, NNLO Sudakov logarithms have positive sign. Thus they compensate in part the impact of NLO EW corrections.
In Figure 9 exact NLO EW results are also compared to the NLL Sudakov approximation at the same order, denoted as nLO EW. The observed agreement indicates that the Sudakov approximation at NLO works very well, thereby supporting the usage of EW Su-  dakov logarithms at NNLO. Moreover, the fact that nNLO EW results are well consistent with NLO predictions supplemented by the corresponding uncertainties (57) provides an important confirmation of the goodness of the proposed approach for the estimate of EW uncertainties.
The importance of NLO and nNLO EW corrections for different processes and the role of individual uncertainties is shown in more detail in Figure 10. Regarding the size of EW uncertainties we observe that the inclusion of nNLO EW corrections is crucial in order to achieve few-percent accuracy in the tails, while uncertainties at NLO EW can be as large as 10% or beyond.
As shown in Figure 11, the various ratios of p T distributions and their shape receive significant EW corrections, with the largest effects observed in the Z( + − )/γ and W/γ ratios. In these ratios the remain-ing combined EW uncertainties are at the level of few percent in the TeV range, reaching about 5% for p T,V 2 TeV. Interestingly, also the Z( + − )/Z(νν) and W − /W + ratios receive non-negligible EW corrections. In the case of the W − /W + ratio this is due to the behaviour of mixed QCD-EW interference contributions at high p T , which yield relevant (negative) contributions in W + +jet production but less in W − +jet production. As for the Z( + − )/Z(νν) ratio, the observed EW effects can be attributed to p T -migration effects induced by QED radiation off leptons. At moderate p T,Z , the invariant mass of photon-lepton pairs that lie inside the recombination cone ∆R γ < 0.1 is well below M Z . Thus a significant fraction of the Z → + − γ phase space does not undergo photon-lepton recombination, and photon radiation results in a negative mass and momentum shift for the + − system. The Z-mass shift is  typically not sufficient to push Z → + − γ events outside the inclusive m window defined in Section 3.2. However, the reduction of the reconstructed p T, results in a negative correction to the Z( + − )/Z(νν) ratio. Vice versa, for p T,Z > ∼ 1 TeV the recombination cone ∆R γ < 0.1 covers photon-lepton invariant masses up to p T,Z ∆ γ > M Z , i.e. beyond the Z → + − γ phase space. As a result, p T, starts capturing a nonnegligible amount of ISR QED radiation, which results in a positive shift of p T, and thus in a positive correction to the Z( + − )/Z(νν) ratio. Note, that the quantitative impact of such corrections depends on the choice of the m mass window. Thus, for a consistent implementation of the predictions presented in this study it is crucial to reweight MC samples using the m window defined in Section 3.2. Moreover, in order to guarantee a consistent extrapolation of QED radiative effects to the m window employed in experimental analyses, it is mandatory to employ MC samples that account for QED radiation off leptons.

Photon-induced production and QED effects on PDFs
Higher-order QCD and EW calculations for pp → V + jet require PDFs at a corresponding accuracy level, i.e. including also QED corrections. The effect of QED interactions on parton densities is twofold. Firstly they introduce a photon parton distribution and so open up partonic channels such as γq → V q . Secondly they modify the quark (and even gluon) PDFs both through QED effects in the initial conditions and especially in the DGLAP evolution.
Photon-induced V + jet production is accounted for by the term d dx σ (V ) γ−ind. in eq. (25). It might become relevant in the TeV range, especially in the case of W + jet production [19,20], where the initial-state photon directly couples to a virtual W boson in the t-channel. Such contributions are suppressed by a relative factor α/α S and can be treated at LO, which corresponds to γq → V q at O(α 2 ) or, if necessary, at NLO QCD, i.e. up to order O(α 2 α S ). This order comprises: The latter can also be understood as photon-induced quark-bremsstrahlung NLO EW contribution to the dominant qq channel. See the contributions of type (a.5) in Section 4.2. Figure 12 illustrates the impact of photon-induced V + jet production at LO according to three recent PDF sets that implement QED corrections. Effects of the order of 5-10% for W +jet can be observed in the TeV region if CT14qed_inc [65] or LUXqed PDFs [39] are used. Much larger effects are found with NNPDF30qed [66,67]. The impact of photon-induced production to Z+jet (and also γ+jet) processes on the other hand is negligible [17,18].
For the description of PDFs and their uncertainties we will use the LUXqed PDFs and their intrinsic uncertainties, given that this set of parton distributions implements a model-independent, data-driven determination of the photon distribution. From Figure 12 one sees that the LUXqed uncertainties for γp → V + jet are small. Using the CT14qed_inc PDFs, based on a non-perturbative model with limited data-based constraints for the inelastic contribution, would result in fairly similar photon-induced cross sections but somewhat larger uncertainties (not shown) as compared to LUXqed PDFs. The NNPDF30qed parton distributions are model independent and data driven, but are based on a different approach from LUXqed for deducing the photon distribution from data, which results in large uncertainties in the photon-induced component, of the order of 100% for pp → + ν + jet at p T, = 1 TeV [20].
We have verified that the NLO QCD corrections to photon-induced production have an impact at the percent level relative to O(α 2 ) and can safely be omitted. This implies that γp → V + jet can be regarded as independent processes. Thus photon-induced V + jet production can be either included through the partonlevel predictions provided in this study or handled as separate background processes through dedicated MC simulations.
Concerning the size of the QED effects on the QCD partons, Figure 13 examines the two main parton luminosities that contribute to the Z+jet process, i.e. gΣ = 2 i (L gqi + L gqi ) (which dominates) and qq = 2 i L qiqi (which accounts for the remaining 15%−30%). It shows the ratio of these luminosities in LUXqed_plus_ PDF4LHC15_nnlo relative to the PDF4LHC15_nnlo set on which it is based. The ratio is given as a function of half the partonic invariant mass, M/2, which is commensurate with the p T of the Z.
Most of the difference between the LUXqed set and PDF4LHC15_nnlo results in Figure 13 comes from the QED effects in the DGLAP evolution [70], with photon emission during the evolution reducing the momentum in the quarks. This effect reaches about 2% at 2 TeV for the gΣ luminosity. There is also a part of the correction associated with the impact of QED effects on the initial partons. In the LUXqed set this has been approximated by absorbing the photon momentum from the gluon distribution in PDF4LHC15_nnlo and keeping   the quarks unchanged at a scale of 10 GeV. This is an ad-hoc procedure, however, insofar as the photon carries only 0.3% of the proton momentum (at a scale of 10 GeV), the uncertainty associated with the arbitrariness of this choice should be below 1%.

PDF uncertainties
The role of PDF uncertainties can be significant especially at high-p T , where PDFs tend to be less precise.
In Figure 14 we illustrate the effect of PDF uncertainties within LUXqed (for the quark and gluon uncertainties based on PDF4LHC15_nnlo_100) for the different V +jets processes and process ratios at NLO QCD. Up to about 800 GeV the PDF uncertainties on the nominal p T distributions remain below 2%. In the tails of the distributions the PDF uncertainties significantly increase. They grow beyond 5% for p T 1.5 TeV. In the Z/W ratio the PDF uncertainties cancel almost completely and remain below 0.5(2)% up to p T ≈ 800(1500) GeV. In the Z/γ and W/γ ratios the PDF uncertainties are at the level of 1 − 2% up to p T ≈ 1300 GeV, while the W − /W + ratio is subject to PDF uncertainties beyond 5% already for p T 1 TeV, driven by uncertainties on the u/d-ratio at large Bjorken-x [3].
To keep track of PDF uncertainties in the combination of QCD and EW corrections we introduce a generalised set of QCD nuisance parameters, ε QCD = (ε QCD,1 , ε QCD,2 , ε QCD,3 , ε PDF,1 , ε PDF,2 , . . . ),  Fig. 13: Impact of QED effects on the two partonic luminosities (gΣ and qq) that contribute dominantly to the Z+jet cross section. The luminosity for producing a system of mass M from two flavours a and b is defined as xs , M 2 ) and the gΣ luminosity corresponds to 2 i (L gqi + L gqi ), while the qq luminosity corresponds to 2 i L qiqi , where i runs over quark flavours. The solid red lines correspond to the ratio of luminosities obtained with the LUXqed_plus_PDF4LHC15_nnlo_100 [39] and PDF4LHC15_nnlo_100 [33] sets, where a given M/2 value corresponds roughly to the same p T,Z . The bands represent the PDF4LHC15_nnlo_100 uncertainty, shown for comparison. which comprises QCD scale and shape variations, as well as process-correlation and PDF uncertainties. To this end we extent Eq. (39), introducing a sum over the 107 independent Hessian PDF replicas provided by the PDF set LUXqed_plus_ PDF4LHC15_nnlo. Such a combination corresponds to the PDF4LHC recommendation as detailed in Eq. (20) of Ref. [33]. These PDF variations should be applied in a fully correlated way across processes and p T bins. As specified in more detail in Section 4.6, the various uncertainties parametrised through Eq. (60) should be applied at the level of QCD calculations and treated on the same footing in the combination of QCD and EW corrections.

Real-boson emission
Inclusive diboson production (in particular pp → V V +jets) can be understood as the real-emission counterpart to NLO EW corrections to pp → V + jet.
Both contributions are separately finite and well defined if V = W, Z. Although they are expected to cancel against each other to a certain (typically small) extent, in practice one should only make sure that both types of processes, pp → V + jet and pp → V V (+jets) with leptonic and hadronic decays of the V , are included in the analysis, and, in order to avoid double counting, contributions of type V V (+jets) should be included in separate diboson MC samples and not as EW correction effects in V + jets samples. Unless a very strong cancellation is observed (which is typically not the case), there is no reason to worry about the possible correlation of uncertainties in V + jets and V V (+jets) production, i.e. one can treat the respective uncertainties as uncorrelated.
As concerns the accuracy of MC simulations of pp → V V (+jets), it is important to notice that a large diboson background to inclusive vector-boson production at high p T is expected to arise from pp → V V j topologies with a hard back-to-back V j system accompanied by a relatively soft extra vector boson. This calls for a reliable description of V V + jet including QCD (and possibly EW) corrections. Thus we recommend the use of merged diboson samples that include at least one extra jet at matrix-element level. At the TeV scale, the EW corrections to pp → V V + jet can become quite large [71,72] and should ultimately be included, together with the corresponding QCD corrections [73][74][75][76][77][78][79][80].

Combination of QCD and electroweak corrections
The combination (25) of higher-order predictions presented in the previous sections can be cast in the form, LO (x, µ) corresponds to the standard additive combination of QCD and EW corrections as defined, respectively, in Eq. (27) and Eqs. (41)- (43). Note that the scaledependent LO QCD K-factor in Eq. (63) is due to the fact that QCD and EW correction factors are normalised to σ LO QCD (µ), respectively. Mixed QCD-EW corrections of relative O(αα S ) are not known to date. However, it is possible to obtain an improved prediction that partially includes such mixed effects by combining higher-order EW and QCD cor-rections through a factorised prescription 13 , The higher-order terms induced by this factorised formula can be written as where κ (V ) N k LO denotes the pure higher-order contribution to the QCD K-factor, i.e.
in analogy with the definition of the κ EW correction factor (41). The prescription (64) is motivated by the factorisation of QCD corrections from the large Sudakovenhanced EW corrections at high energies [64] and by the observation that in cases where the multiplicative and additive approach are far apart from each other, such as in the presence of giant K-factors [19,83], the former turns out to be much more reliable. In general, when QCD and EW corrections are simultaneously enhanced, the O(αα S ) mixed terms that are controlled by the multiplicative prescription can become quite significant. We also note that, thanks to the fact that the relative EW correction factors κ (V ) EW (x) are essentially insensitive to QCD scale variations, the scale dependence of the multiplicative combination (64) is similar as for pure N k LO QCD predictions. In contrast, the additive approach (63) can suffer from sizable scale uncertainties when EW corrections become large.
In order to estimate the typical size of higher-order effects that are not captured by the factorised prescription (64), we cast mixed QCD-EW corrections of O(αα S ) in the form (67) and to model the non-factorising term we use the simple Ansatz 14 The expectation that the bulk of QCD and EW corrections factorise implies that the absolute value of the free process-dependent factors ξ (V ) should be well below one. Note that Eq. (68) is equivalent to i.e. we assume that non-factorising EW-QCD mixed terms are proportional to the difference between the additive and multiplicative combination of QCD and EW corrections. The NLO EW corrections to pp → V +2 jets [19,51], which represent a real-virtual contribution to the unknown mixed EW-QCD NNLO corrections to V + jet production, can provide useful insights into the typical size of the ξ (V ) factors and the goodness of the Ansatz (67)- (68). In particular, starting from the O(αα S ) contributions to Eq. (67), it is possible to establish a relation between nonfactorising NNLO mixed corrections and the differences 14 As discussed below, the goodness of this naive Ansatz will be justified by fitting it to a realistic estimator of δκ between NLO EW K-factors for V + 2 jet and V + 1 jet production. To this end, we consider the identity which is obtained by multiplying both sides of Eq. (70) by the LO QCD cross section for pp → V + 1 jet and restricting the phase space to real-virtual contributions with V + 2 jet final states. This restriction is implemented by means of an N -jettiness [84] resolution parameter τ cut , as described in more detail below, and the above equation should be understood as definition of δκ NNLO mix (x, τ cut ), which will be used as estimator of δκ NNLO mix (x) in Eq. (70). In Eq. (71) we use the notation NLO EW (x), and we keep the µ-dependence as implicitly understood, since the term δκ (V ) NNLO mix (x, τ cut ) is expected to be quite stable with respect to scale variations. Instead, the τ cut parameter plays an important role since it acts as a cutoff of infrared QCD singularities in the regions where the second jet becomes soft or collinear. Based on the universal behaviour of IR QCD effects, such singularities are expected to factorise into identical singular factors on the left-and the right-hand side of Eq. (71). Thus, while the δκ (V ) NNLO mix (x, τ cut ) term on the right-hand side depends on τ cut , this dependence is expected to be free from large τ cut -logarithms and thus reasonably mild.
As anticipated above, solving for δκ NNLO mix we obtain the relation which allows us to estimate non-factorising mixed effects in terms of the difference between the V + 2-jet and V + 1-jet EW κ-factors. To this end, we will match the estimator (72) to the Ansatz (68). More precisely, we will fix the free coefficients ξ (V ) in Eq. (68) in such a way that for the whole x-spectrum and within an appropriately chosen τ cut range. Thanks to the cancellation of IR QCD singularities in Eq. (72), the resulting ξ (V ) coefficients should be reasonably stable with respect to the choice of the resolution parameter. Thus, τ cut can be varied in a rather wide range. In principle one could even consider the τ cut → 0 limit of Eq. (73). However, given that two-loop mixed EW-QCD contributions are not taken into account, this limit does not converge towards the full NNLO result corresponding to τ cut = 0. Moreover, for very small values of τ cut the numerator and denominator of κ V +2 jets NLO EW (x, τ cut ) are dominated by universal τ cut -logarithms that should cancel against virtual two-loop terms, and since such logarithms factorise, their dominance can result in an underestimation of non-factorising effects. Vice versa, excessively large values of τ cut can lead to an overestimation of non-factorising effects. This is due to the fact that increasing τ cut enhances the difference between EW κfactors in Eq. (73) but also suppresses the cross section of the V + 2-jet subprocess, rendering it a less and less significant estimator of the behaviour of mixed corrections for inclusive V + jet production. Thus, excessively small or large values of τ cut should be avoided. Based on the above considerations, for the fit of the ξ (V ) coefficients we require that Eq. (73) is fullfilled in a wide τ cut -range while keeping the σ V +2 jet /σ V +1 jet ratio at order one, in such a way that the V + 2 jet cross section is neither too suppressed nor too enhanced. This procedure is implemented using an N -jettiness cut parameter [84]. More precisely, we use the dimensionless one-jettiness parameter where the p i are light-like vectors for each of the initial beams and the hardest final-state jet, and the Q i characterise their respective hardness, which we set as Q i = 2E i . The hardest final-state jet is defined by applying an anti-k T algorithm with R=1 to all final-state partons. 15 The q k denote the four-momenta of any such final-state parton, and √ŝ is the partonic centre-of-mass energy. All quantities are defined in the hadronic centreof-mass system.
To isolate two-jet configurations against one-jet configurations we require τ 1 > τ cut , and the cut is varied in the range 0.001 ≤ τ cut ≤ 0.04. As demonstrated in Figure 15, this choice keeps the σ V +2 jet /σ V +1 jet ratio around order one, as desired. Moreover, we observe that the estimator (73) remains quite stable with respect to τ cut variations (see the solid lines in the right plot). Non-factorising effects turn out to be generally very small. They exceed the percent level only in the TeV tails of the distributions. As illustrated by the gray band in Figure 15 (right), setting guarantees an acceptable matching of the Ansatz (68) to the estimator (73). More precisely, for W + jet production the shape of the Ansatz (68) tends to overestimate the uncertainty in the p T range between one and 15 In order to guarantee a proper cancellation of QCD and EW singularities, the jet algorithm is applied to all QCD partons and photons, excluding photons that are recombined with leptons, as well as the leading identified photon in case of the γ+jets process.
two TeV. However, we have checked that the Ansatz becomes much less adequate if the full EW correction in Eq. (67) is replaced by its non-Sudakov part.
The rather small values of the ξ (V ) coefficients confirm that the bulk of the EW and QCD corrections factorise. However, in the case of W + jet and γ+ jet production, the relative size of non-factorising corrections appears to be rather significant. This is due to the behaviour of the EW κ-factors in the multi-TeV region, where the difference between the EW κ-factors for pp → V + 1 jet and pp → V + 2 jet is enhanced by the presence of mixed EW-QCD interference contributions in channels of type qq → qqV (see the contributions of type a.5 in Section 4.2). More precisely, EW-QCD interference effects of O(α S α 2 ) enhance the EW corrections to pp → V + 1 jet as a result of the opening of the qq channel at NLO EW, while in pp → V + 2 jet the EW K-factor is not enhanced since the qq channel is already open at LO. Based on this observation, and also due to the fact that the main effect of the opening of the qq channel is already reflected in the NLO QCD K-factor for V + 1 jet production, the above mentioned EW-QCD interference effects could be excluded from the factorisation prescription (64) and treated as a separate contribution. As illustrated by the dashed curves in Figure 15, this approach would lead to a drastic reduction of non-factorising effects, especially for γ+ jet production. Nevertheless, given that the effects observed in Figure 15 are subdominant with respect to current PDF and statistical uncertainties, in the present study we refrain from implementing such a splitting.

Combination of QCD and EW corrections with related uncertainties
Based on the above analysis, we recommend to combine QCD and EW corrections according to the multiplicative prescription (67), treating the non-factorising term (68) as uncertainty and using the estimated ξ (V ) factors given in Eq. (75). Including QCD and EW uncertainties as specified in Eq. (39) and Eq. (58), this leads to  the combination formula where the uncertainty associated with non-factorising mixed EW-QCD terms reads The related nuisance parameter, ε mix , should be Gaussian distributed with one standard deviation corresponding to the range ε mix ∈ [−1, +1]. Given that mixed uncertainties have been estimated using a proxy of the full NNLO QCD-EW calculation, it would be reasonable to assume some degree of correlation across different V + jet processes. However, for simplicity in this study we keep ε mix variations fully uncorrelated, bearing in mind that this approach is probably too conservative. In Figure 16 we compare the additive and multiplicative combinations of QCD and EW corrections showing also the corresponding uncertainty estimate (77) for various V +jet processes.

Summary and conclusions
The precise control of SM backgrounds, and notably of pp → Z(νν)+ jets, is crucial in order to maximise the potential of MET+jets searches at the LHC. Such backgrounds can be predicted directly using QCD and EW calculations. Alternatively, QCD and EW calculations can be used to relate them to experimental data for similar V + jet production processes, i.e. pp → γ+ jets, pp → W ( ν)+ jets and pp → Z( + − )+ jets.
In this article we have presented predictions for inclusive vector-boson p T distributions based on the most advanced calculations available today, bringing together results from a number of groups so as to have perturbative QCD to NNLO accuracy, EW corrections to NLO accuracy and additionally the inclusion of 2-loop EW Sudakov logarithms.  Fig. 16: Comparison of additive (green) and multiplicative (red) combination of (N)NLO QCD and nNLO EW corrections for various pp → V +jet processes at 13 TeV. The red band corresponds to the mixed QCD-EW uncertainty (77). The (N)NLO QCD result without EW corrections is shown in blue. The combination at NLO QCD is shown on the left and at NNLO QCD on the right.
A substantial part of our study concerned uncertainty estimates. In particular we proposed and applied various new approaches for uncertainty estimates and correlations across processes and p T regions.
We defined the uncertainties due to normal QCD scale variations in a way that gives a strong correlation across different p T regions, Eq. (33). We then supplemented it with a shape uncertainty that is anticorrelated across p T , Eqs. (35)- (36). To address the long-standing problem of evaluating the correlations between uncertainties for different processes, we separated the uncertainty into process-independent and process-dependent components. The universal component was taken to be composed of the overall scale and shape uncertainties for the reference Z + jet process. The process-dependent component, which is generally small, was determined by considering the difference between suitably normalised K-factors for the different processes, Eq. (38). This amounts to a conservative choice of taking the uncertainty on ratios as the difference between the best available prediction and the one at one order lower.
Special attention was devoted to the correlation of Z/W + jet and γ+ jet production. In that case a substantial non-universal contribution is associated with the masslessness of the photon and the need to control collinear divergent q → qγ radiation through a photonisolation prescription. We introduced a novel photonisolation prescription with a dynamically chosen isolation radius, Eq. (11), designed to suppress q → γq radiative effects in a way that is similar to the effect of the masses of the Z and W bosons in the case of q → V q splittings at large p T . Such a dynamic isolation allows one to split γ+jet production into a quasiuniversal part, which can be treated on the same footing as Z + jet and W + jet production, and a non-universal part which is kept uncorrelated. The non-universal part is given by the difference between the cross sections with conventional and dynamic photon isolation prescriptions.
For pure EW corrections we considered three uncertainty sources for unknown higher order contributions. These address unknown Sudakov logarithms beyond NNLO and/or NLL accuracy, as well as unknown hard (non-Sudakov) EW corrections beyond NLO and process-correlation effects.
One potentially large source of uncertainty arises from mixed QCD and EW corrections, given that both O(α S ) and O(α) NLO corrections can be large and that the O(αα S ) NNLO corrections are not currently known. We chose a multiplicative scheme for combining EW and QCD corrections. To obtain an estimate of unknown O(αα S ) corrections not captured by this factorised ansatz, we studied the NLO EW corrections to V + 2 jet production, which represent the real-virtual part of a full O(αα S ) calculation for V + jet production. Based on this analysis, we concluded that it is reasonable to assume that the multiplicative combination of QCD and EW corrections describes the full O(αα S ) correction with a relative uncertainty that varies between 10-20% for pp → W/Z+ jet and 40% for pp → γ+ jet.
Overall, QCD corrections are substantial, a few tens of percent at NLO, and up to 10% at NNLO. The NNLO results are consistent with the NLO predictions within our prescription for the uncertainty bands of the latter. This is true not just for absolute cross sections and their shapes, but also for ratios of cross sections. These ratios are remarkably stable across LO, NLO and NNLO QCD corrections, see Figure 8. Using dynamic photon isolation, this statement holds true also for the γ + jet process at p T 300 GeV.
The EW corrections to V +jet cross sections amount to a few tens of percent in the TeV region, see Figure 9. In the ratios they cancel only in part, due to the sensitivity of EW effects to the SU(2) charges of the produced vector bosons. At the TeV scale, the NNLO Sudakov logarithms can reach the several percent level and their systematic inclusion is an important ingredient in order to achieve percent precision at very high p T .
In Figure 17 we summarize our uncertainty estimates for the different V +jet processes and process ratios. Here we combine in quadrature all sources of perturbative uncertainties at N(N)LO QCD ⊗ nNLO EW and we overlay the remaining PDF uncertainties. For convenience, PDF variations have been assessed using NNLO PDFs in combination with NLO QCD calculations, but can be safely applied to the NNLO QCD results.
The nominal p T distributions at N(N)LO QCD ⊗ nNLO EW are constrained at the 10(5)% level up to about 1 TeV and at the 20(10)% level up to about 2 TeV. In the process ratios these uncertainties cancel to a large extent. In particular, in the Z/W ratio remaining uncertainties are at the level of only 1-2% up to 1 TeV and below 5% up to 2 TeV. Similarly, the Z/γ ratio is constrained at the 5% level up to 2 TeV. Noteworthy, including the NNLO QCD corrections the process ratios remain very stable and in particular within the uncertainty estimates based on NLO QCD. This reflects the fact that QCD uncertainties are very well under control: taking at face value the NNLO QCD systematics we are at the level of few percent all the way up to the multi-TeV scale (see Figure 8), and at large p T we are dominated by EW and PDF uncertainties. The latter are below the perturbative uncertainties in all nominal distributions and all but the W − /W + ratio, where a precise measurement at high p T could help to improve PDF fits. In this respect, we note that the theoretical uncertainty for the W − /W + ratio is entirely dominated by mixed QCD-EW effects and is most likely overestimated due to our conservative assumption of keeping such uncertainties uncorrelated across processes (see Section 4.6).
We also discussed photon-induced contributions and QED corrections to PDFs. In this context, for a precise prediction of the γ-PDF we have advocated the use of the LUXqed_plus_PDF4LHC15_nnlo PDFs, which implement a data driven determination of the γ-PDF. For a consistency treatment of O(α) effects in the PDFs, the LUXqed_plus_PDF4LHC15_nnlo distributions should be used in all photon-, quark-, and gluon-induced channels. 16 Photon-induced effects are negligible in Z+ jet and γ+ jet production, but their impact on pp → W + jet, and thus on the W/Z and W/γ ratios, can reach the 5% level at the TeV scale 17 (see Figure 12).
Our predictions are provided in the form of tables for the central predictions and for the different uncertainty sources. Each uncertainty source is to be treated as a 1-standard deviation uncertainty and pragmatically associated with a Gaussian-distributed nuisance parameter.
The predictions are given at parton level as distributions of the vector boson p T , with loose cuts and inclusively over other radiation. They are intended to be propagated to an experimental analysis using Monte Carlo parton shower samples whose inclusive vectorboson p T distribution has been reweighted to agree with our parton-level predictions. The impact of additional 16 This is automatically achieved by reweighting MC samples generated with arbitrary PDFs with our complete N k LO QCD × nNLO EW predictions based on LUXqed_plus_PDF4LHC15_nnlo PDFs. Vice versa, restricting the reweighting to pure EW corrections and using MC samples based on different PDFs can lead to inconsistencies at O(α). 17 Note that photon-induced contributions are not included in the summary plots of Figure 17.  Some caution is needed in implementing the results of this paper: for example the uncertainty prescriptions are tied to the use of the central values that we provide.
If an experiment relies on central values that differ, e.g. through the use of MC samples that are not reweighted to our nominal predictions, then the uncertainty scheme that we provide may no longer be directly applicable. Furthermore, for searches that rely on features of the event other than missing transverse momentum, one should be aware that our approach might need to be extended. This would be the case notably for any observable that relies directly on jet observables, whether related to the recoiling jet or vetoes on additional jets.
Overall, it is possible to obtain precise theoretical control both for vector-boson p T distributions, and for their ratios, at the level of a few percent. We expect this precision, across a wide range of p T , to be of significant benefit in MET+jets searches, notably enabling reliable identification or exclusion of substantially smaller BSM signals than was possible so far. In fact, since the release of the first version of this paper, the background estimates we propose here have been adopted in analyses by ATLAS [85] and CMS [86]. Table 2: List of processes, highest available QCD and EW order, and process labels used in data files (see Table 3). Predictions for pp → ν ν + jet are available only at NLO QCD, but corresponding NNLO QCD corrections and uncertainties can be taken from pp → + − + jet. normalisation according to Eqs. (62) and (76). For instance, at NNLO QCD ⊗ nNLO EW the relative impact of QCD and EW uncertainties should be δ (i) K NNLO /K NNLO and δ (i) K EW /(1+κ EW ), respectively.
Concerning QCD contributions, predictions at NNLO QCD should be combined with uncertainties at the same order. However, before higher-order QCD calculations are thoroughly validated against high-statistics measurements at moderate transverse momenta, theory uncertainties should be assessed in a more conservative way. To this end, we advocate the usage of NNLO QCD nominal predictions in combination with NLO QCD uncertainties, while keeping all EW effects at nNLO EW level.
All predictions and uncertainties for pp → γ+ jet are based on the dynamic photon isolation prescription introduced in Section 3.1. As explained therein, this requires an extra γ+ jet specific uncertainty, which needs to be evaluated by means of a separate reweighting in a standard Frixione isolation setup with fixed cone. Corresponding theoretical predictions at NLO QCD are denoted as K (γ,fix) NLO (x) in Table 3.

Appendix B: QCD and EW uncertainties
In this appendix we present a series of technical plots that illustrate the relative importance of the various sources of QCD and EW uncertainties discussed in Sections 4.1-4.2. The impact of individual QCD uncertainties, δ (i) K N k LO , in p T spectra and ratios is illustrated in Figures 18-19. Similar plots for the three types of EW uncertainties, δ (i) κ    Fig. 19: Ratios of p T -distributions for various pp → V +jet processes at LO and NLO QCD. The related scale uncertainties (33), shape uncertainties (35) and process-correlation uncertainties (38) are correlated amongst all processes as discussed in Section 4.1.    57), is plotted.