arxiv: v2 [hep-ph] 16 Oct 2017

Transcription

1 Eur. Phys. J. C manuscript No. (will be inserte by the eitor arxiv: v2 [hep-ph] 6 Oct 27 Precise preictions for V +jets ark matter backgrouns J. M. Linert, S. Pozzorini 2, R. Boughezal 3, J. M. Campbell 4, A. Denner 5, S. Dittmaier 6, A. Gehrmann-De Rier 2,7, T. Gehrmann 2, N. Glover, A. Huss 7, S. Kallweit 8, P. Maierhöfer 6, M. L. Mangano 8, T.A. Morgan, A. Mück 9, F. Petriello 3,, G. P. Salam a,8, M. Schönherr 2, C. Williams Institute for Particle Physics Phenomenology, Department of Physics, University of Durham, Durham, DH 3LE, UK 2 Physik-Institut, Universität Zürich, Winterthurerstrasse 9, CH-857 Zürich, Switzerlan 3 High Energy Physics Division, Argonne National Laboratory, Argonne, IL 6439, USA 4 Fermilab, P.O.Box 5, Batavia, IL 65, USA 5 Universität Würzburg, Institut für Theoretische Physik un Astrophysik, 9774 Würzburg, Germany 6 Albert-Luwigs-Universität Freiburg, Physikalisches Institut, 794 Freiburg, Germany 7 Institute for Theoretical Physics, ETH, CH-893 Zürich, Switzerlan 8 Theoretical Physics Department, CERN, CH-2 Geneva 23, Switzerlan 9 Institut für Theoretische Teilchenphysik un Kosmologie, RWTH Aachen University, D-5256 Aachen, Germany Department of Physics & Astronomy, Northwestern University, Evanston, IL 628, USA Department of Physics, University at Buffalo, The State University of New York, Buffalo 426 USA Receive: ate / Accepte: ate Abstract High-energy jets recoiling against missing transverse energy (MET are powerful probes of ark matter at the LHC. Searches base on large MET signatures require a precise control of the Z(ν ν+ jet backgroun in the signal region. This can be achieve by taking accurate ata in control regions ominate by Z(l + l + jet, W (lν+ jet an γ+ jet prouction, an extrapolating to the Z(ν ν+ jet backgroun by means of precise theoretical preictions. In this context, recent avances in perturbative calculations open the oor to significant sensitivity improvements in ark matter searches. In this spirit, we present a combination of state-of-the-art calculations for all relevant V + jets processes, incluing throughout NN QCD corrections an N electroweak corrections supplemente by Suakov logarithms at two loops. Preictions at parton level are provie together with etaile recommenations for their usage in experimental analyses base on the reweighting of Monte Carlo samples. Particular attention is evote to the estimate of theoretical uncertainties in the framework of ark matter searches, where subtle aspects such as correlations across ifferent V + jet processes play a key role. The anticipate theoretical uncertainty in the Z(ν ν+ jet backgroun is at the few percent level up to the TeV range. a on leave from CNRS, UMR 7589, LPTHE, F-755, Paris, France Introuction CERN-TH-27-2, CERN-LPCC-27-2, FERMILAB-PUB-7-52-T, IPPP/7/38, ZU TH 2/7 The signature of missing transverse energy (MET is one of the most powerful tools in the interpretation of ata from haron colliers. In the Stanar Moel (SM, MET arises from the neutrinos from the ecay of W an Z bosons, an it can be use in their ientification an stuy, as well as in the ientification an stuy of Higgs bosons, top quarks an other SM particles whose ecay proucts inclue W or Z bosons. But MET is also an almost omnipresent feature of theories beyon the SM (BSM, where it can be associate to the ecay of new particles to W an Z bosons, or irectly to the prouction of new stable, neutral an weakly interacting particles. Typical examples are theories with ark matter (DM caniates, or Kaluza- Klein theories with large extra imensions. Depening on the etails, MET is accompanie by other moeliscriminating features, such as the presence of a small or large multiplicity of har jets, or of specific SM particles. The experimental search for these extensions of the SM relies on a proper moeling of the SM backgrouns to the MET signature. The etermination of these backgrouns is ieally one by using ata control samples, but theoretical input is often helpful, or even necessary, to exten the experimental information from the control to the signal regions, or to exten the application range of the backgroun preictions an to improve their precision [ 3].

2 2 In this paper we focus on the theoretical moeling of the SM V + jet backgrouns to inclusive prouction of large MET recoiling against one or more haronic jets. These final states aress a broa set of BSM moels, where the prouction of an otherwise invisible final state is reveale by the emission of one or more high-p T jets from initial state raiation, where p T is the momentum in the transverse plane. Recent publications by ATLAS [5] an CMS [6, 7], relative to LHC ata collecte at s = 3 TeV, ocument in etail the current experimental approaches to the backgroun evaluation. The leaing backgroun is Z(ν ν+ jet prouction, followe by W (lν+ jet (in particular for l = τ or when the lepton is outsie of the etector. 2 The experimental constraints on Z(ν ν+ jet prouction at large MET can be obtaine from accurate measurements of V + jet prouction processes with visible vector-boson signatures. It is quite obvious, for example, that the measurement of Z(l + l + jets with l = e, µ is the most irect an reliable proxy for Z(ν ν+ jets. This control sample, however, is statistics limite, ue to the smaller branching ratio of Z bosons to charge leptons relative to neutrinos. To extrapolate the shape of the Z spectrum to the largest p T values, therefore, requires a theoretical preiction. The larger statistics of W (lν+ jets an γ+ jets events makes it possible to irectly access the relevant p T range, but the relation between their spectra an the Z spectrum nees, once again, theoretical guiance. To put things into a concrete perspective, Figure shows the expecte event rates, an the relative statistical uncertainty, for 3 fb of integrate luminosity at 3 TeV. The extrapolation to the O( fb an O(3 fb expecte from the full run 2 an at the en of the full LHC programme, respectively, is straightforwar. The Z(l + l +jets ata allow for a irect estimate of the Z(ν ν+jets rate with a statistical precision below % for p T up to about 6 GeV. Using the W (lν+jets or γ+jets ata coul in principle exten this range up to about 9 GeV. Beyon this value, the statistical precision of the W (lν+jets an γ+jets events remains a factor of two better than that of the Z(ν ν+jets signal. In orer to ensure that the theoretical systematics in the extrapolation from the W +jets an γ+jets rates to the Z+jets rates remains negligible with respect to the statistical uncertainty, the former shoul be kept at the level of a few percent up to p T 2 TeV, an aroun % up to For a recent comprehensive review of DM moels leaing to this class of signatures, see e.g. [4]. 2 Other backgrouns (such as QCD multijets, t t or pairs of gauge bosons are suppresse, an their contribution to the overall uncertainty is well below the percent level. N = σ(pt,v > p min T,V L / N 8 N QCD for 3 TeV Z(ν ν+ jet Z(l + l + jet W(lν+ jet γ+ jet L = 3 fb p min T,V [GeV] Fig. : Prouction rates for V + jet(s, for various ecay channels, as a function of the minimum p T of the vector boson. Decays into l ± = e ±, µ ± an ν e, ν ν, ν τ are inclue. The number of events, N, is normalize to 3 fb of LHC ata at s = 3 TeV, an inclues the basic selection cuts liste in the main boy of the paper. The log lower panel shows the statistical uncertainties, calculate as / N. The gray ban in the lower panel inicates the regime of % statistical uncertainty. p T TeV, which is the ultimate kinematic reach for the Z(ν ν+jets signal at the en of LHC ata taking. The main result of this work is to prove that, thanks to the recent theoretical avances, these goals can be met. This proof requires the analysis of a series of possible effects. On the one han, the theoretical extrapolation to larger p T of the very precise Z(l + l +jets ata requires firm control over the shape of the istribution. Several effects, from the choice of parton istribution functions (PDFs to the choices mae for the renormalisation an factorisation scales use in the calculations, can influence the extrapolation. On the other han, the level of correlation between the W, γ an Z spectra must be kept uner control. At large p T, in particular, large an process-epenent corrections arise ue to the growth of the electroweak (EW corrections, an these may spoil the correlation inuce by pure QCD effects. For our analysis we shall use the most up-to-ate theoretical preictions available toay for the escription of vector boson prouction at large p T. On the QCD sie, we rely on the next-to-next-to-leaing or-

3 3 er (NN calculations, which appeare recently for Z+jet [8 2], W +jet [3, 4] an γ+jet [5, 6] prouction. On the EW sie, we apply full N calculations for Z+jet [7 9], W +jet [9, 2] an γ+jet [2] prouction with off-shell ecays of the Z an W bosons. Given the strong enhancement of EW Suakov effects in the TeV region, we also inclue 2-loop logarithmic terms at next-to-leaing logarithmic (NLL accuracy for all V + jet processes [22 25]. An extensive assessment an iscussion of the estimates of missing higherorer terms, an of the relative systematics, is given in the main boy of this paper. In particular, in orer to aress non-trivial issues that arise in the context of ark matter searches, we introuce a global framework for the estimate of theoretical uncertainties in all V + jet processes, taking into account correlation effects across ifferent processes an p T regions. Also the uncertainties associate with the combination of QCD an EW corrections are iscusse in etail. From the experimental perspective, the etermination of the backgroun composition in signal an control regions, an the moeling of other key aspects of experimental analyses (e.g. lepton ientification an reconstruction, missing energy, etc. require a theoretical escription of the various V + jets processes at the particle level. Typically, this is provie by Monte Carlo (MC samples base on multi-jet merging at or N QCD, an improvements base on higher-orer theoretical calculations can be implemente through reweighting of MC events. For the fit of MC preictions to ata, ATLAS an CMS analyses rely on the profile likelihoo approach, where experimental an theoretical uncertainties are escribe in terms of nuisance parameters with Gaussian istributions. In this context, the correlations of theoretical uncertainties across p T bins (shape uncertainties an across ifferent V + jets processes play a key role for searches at large MET. For the implementation of higher-orer QCD an EW corrections an for the estimate of theoretical uncertainties in the experimental analysis framework, we propose a proceure base on a one-imensional reweighting of MC samples. The propose framework shoul enable the experiments to carry out their profile likelihoo approach, quantifying the impact of the theoretical systematics in their analyses, an valiating irectly with ata the reliability an robustness of the theoretical inputs. In this respect, we woul like to stress that, inepenently of the application to BSM searches, the results in this paper provie a framework for incisive valiations 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 shoul confirm the usefulness of the approach we propose, the same framework coul be aapte to more complex or exclusive final states, in which for example MET is accompanie by a large number of (har jets or by specific objects (photons, heavy quarks, Higgs, etc. These extensions are left for future stuies. The structure of this paper is as follows: In Section 2 we introuce the reweighting technique, to incorporate in a MC analysis the effect of higher-orer corrections an of their systematic uncertainties incluing correlations. Section 3 escribes etails of the setup for our numerical calculations, the employe tools an methos, as well as the etaile efinition of physics objects an observables to be use in the context of MC reweighting. In Section 4 we iscuss higher-orer QCD an EW corrections, incluing the contribution of photoninitiate processes an real vector boson emission. We present here our approach to the estimate of the various systematics, covering QCD scale, shape an processepenent uncertainties, as well as uncertainties arising from higher-orer EW an mixe QCD EW corrections. Section 5 contains our summary an conclusions. As etaile in Appenix A, results for all V +jets processes are available in form of one-imensional histograms in the vector-boson p T covering central preictions an all mentione uncertainties. Technical plots on the iniviual sources of QCD an EW uncertainties are ocumente in Appenix B. 2 Reweighting of Monte Carlo samples The reweighting of MC samples is an approximate, but straightforwar an easy to implement metho of combining (N MC simulations with (NN QCD+N EW perturbative calculations an to account for the respective uncertainties in a systematic way. The following formula escribes the one-imensional reweighting of MC samples for V + jet prouction (V = γ, Z, W ± in a generic variable x, x y σ(v (ε MC, ε TH = ( [ x y σ(v MC (ε x MC σ(v TH (ε ] TH. MC (ε MC In the case at han, i.e. V + jet prouction, the oneimensional parameter x shoul be unerstoo as the vector-boson transverse momentum, x = p (V T, while y generically enotes the remaining variables of the fully ifferential kinematic epenence of the accompanying QCD an QED activity, incluing both extra jet an photon raiation, as well as leptons an neutrinos from

4 4 haron ecays. It is implicitly unerstoo that x y σ epens on x an y, while in xσ the variables y are integrate out. The labels MC an TH in Eq. ( refer to Monte Carlo an higher-orer theoretical preictions, respectively, an the relate uncertainties are parametrise through nuisance parameters ε TH, ε MC. Our recommenations for theory uncertainties in Section 4 are formulate in terms of intervals for the relate nuisance parameters, < ε TH,k <, (2 which pragmatically shoul be unerstoo as the σ range of Gaussian uncertainties. Monte Carlo uncertainties, escribe by ε MC, must be correlate in the numerator an enominator on the r.h.s. of Eq. (, while they can be kept uncorrelate across ifferent processes (apart from Z(ν ν + jet an Z(l + l + jet. We note that, as oppose to an approach base only on ratios of p T istributions, where theory is use for extrapolations across ifferent processes at fixe p T, MC reweighting is more powerful as it supports all possible extrapolations across ifferent processes an p T regions. In particular, it makes it possible to exploit V + jet precision measurements at moerate p T in orer to constrain Z(ν ν + jet prouction in the TeV region. A further avantage of the reweighting approach ( lies in the fact that the three terms on the r.h.s. of Eq. ( o not nee to be compute with the same numerical setup (parameters, cuts, observables, etc.. More precisely, only the efinition of the variable x an the binning of its istribution nee to be the same in all three terms. Scale choices, QCD an EW input parameters an PDFs shoul be the same only in the numerator an enominator of R MC (x, y = x y σ(v MC MC, (3 but can be chosen in a ifferent way in TH, provie that QCD an EW corrections themselves are compute using the same settings. Vice versa, possible cuts must be ientical only in the numerator an enominator of R TH/MC (x = TH MC, (4 y σ(v while particle-level MC preictions, x MC, can be subject to more exclusive or inclusive cuts in the experimental analysis. For an optimal combination of higher-orer calculations an MC preictions, two conitions shoul be fulfille. On the one han, theory calculations shoul escribe the istribution in the reweighting variable with higher (or at least equal precision as compare to the MC sample, [ ] TH [ ] MC. (5 On the other han, the MC sample shoul be more accurate than TH calculations in escribing the correlation between x an all other variables y, [ x y σ(v MC MC ] [ x y σ(v TH TH ]. (6 More precisely, conition (6 nees to be fulfille only for those aspects of V + jet events that are relevant for the actual experimental analysis. As concerns the first conition, we note that, epening on the choice of the observable x, using stateof-the-art theory calculations that involve higher-orer QCD an EW corrections may not guarantee that Eq. (5 is fulfille. In fact, there are a number of aspects, i.e. resolve multi-jet emissions, the resummation of soft logarithms in the region of small vector-boson p T, soft QCD raiation of non-perturbative origin, multiple photon raiation, or neutrinos an charge leptons resulting from haron ecays, for which fixe-orer perturbative calculations of pp V + jet are less accurate than MC simulations. Thus, the reweighting variable x shoul be efine such as to have minimal sensitivity to the abovementione aspects. In this respect, ue to its reuce sensitivity to multiple jet emissions, the vector-boson p T is a natural choice. However, in orer to fulfil Eq. (5, the region p (V T M V shoul be exclue from the reweighting proceure, unless QCD Suakov logarithms are resumme to all orers in the theoretical calculations. Moreover, in orer to simultaneously fulfil conitions (5 an (6, any aspect of the reconstructe vector-boson p T that is better escribe at MC level shoul be exclue from the efinition of x an inclue in y. This applies, as iscusse in Section 3, to multiple photon emissions off leptons, an to possible isolation prescriptions for the soft QCD raiation that surrouns leptons or photons. In general, purely nonperturbative aspects of MC simulations, i.e. MPI, UE, haronisation an haron ecays, shoul be systematically exclue from the efinition of the reweighting variable x. Thus, impact an uncertainties relate to this non-perturbative moelling will remain as in the original MC samples.

5 5 It shoul be stresse that the above consierations are meant for ark-matter searches base on the inclusive MET istribution, while more exclusive searches that exploit aitional information on har jets may involve aitional subtleties. In particular, for analyses that are sensitive to multi-jet emissions, using the inclusive vector-boson p T as the reweighting variable woul still fulfil Eq. (5, but the lack of QCD an EW corrections to V + 2 jets prouction in MC simulations coul lea to a violation of Eq. (6. In analyses that are sensitive to the tails of inclusive jet-p T an H T istributions this issue is very serious, an QCD+EW corrections shoul be irectly implemente at MC level using multi-jet merging [9]. In general, as a sanity check of the reweighting proceure, we recommen verifying that, for reasonable choices of input parameters an QCD scales, (NN QCD calculations an (N merge MC preictions for vector-boson p T istributions are in reasonably goo agreement within the respective uncertainties. Otherwise, in case of significant MC mismoelling of the istribution, one shoul check the reliability of the MC in extrapolating TH preictions from the reweighting istribution to other relevant observables. In general, one coul check whether the oneimensional reweighting via the variable x in Eq. ( can in fact reprouce the epenence of the corrections in other kinematic variables that are relevant for the experimental analysis. To this en, istributions of σ (V w.r.t. another kinematic variable x shoul be calculate upon integrating Eq. (. Switching on an off the corrections on the r.h.s. of Eq. ( in σ (V TH an taking the ratio of the obtaine ifferential cross sections σ (V, prouces the relative correction to the x istribution that coul be compare to the corresponing result irectly calculate from σ (V TH.3 Finally, it is crucial to check that state-of-the-art preictions for absolute σ/p T istributions agree with ata for the various visible final states. 3 Setup for numerical preictions In this section we specify the physics objects (Section 3., acceptance cuts an observables (Section 3.2, input parameters (Section 3.3 an tools (Section 3.4 use in the theoretical calculations for pp W ± /Z/γ+ jet. The efinitions of physics objects, cuts an observables which specify the setup for the reweighting proceure iscusse in Section 2 shoul be aopte 3 This proceure shoul be restricte to variables x that can be escribe with goo accuracy both in perturbative calculations an in the MC simulations. both for theoretical calculations an for their Monte Carlo counterpart in the reweighting factor (3. The etails of the reweighting setup are esigne such as to avoi any possible eficit in the perturbative preictions (e.g. ue to lack of resummation at small p T an any bias ue to non-perturbative aspects of Monte Carlo simulations (e.g. leptons an missing energy from haron ecays. Let us also recall that this setup is completely inepenent of the physics objects, cuts an observables employe in the experimental analyses. As concerns input parameters an PDFs, the recommenation of Section 3.3 shoul be applie to all QCD an EW higher-orer calculations. In particular, it is manatory to compute (NN QCD an EW corrections in the same EW input scheme, otherwise N EW accuracy woul be spoile. Instea, Monte Carlo simulations an the corresponing contributions to the reweighting factor (3 o not nee to be base on the same input parameters an PDFs use for theory preictions. We recommen hanling W/Z + jet prouction an ecay on the Monte Carlo sie as the full processes pp ll/lν/νν + jet, i.e. with a consistent treatment of off-shell effects, as is one on the theory sie. 3. Definition of physics objects MC In the following we efine the various physics objects relevant for higher-orer perturbative calculations an for the reweighting in the Monte Carlo counterparts in Eq. (3. Neutrinos In parton-level calculations of pp ll/lν/νν+jet, neutrinos originate only from vector-boson ecays, while in Monte Carlo samples they can arise also from haron ecays. In orer to avoi any bias in the reweighting proceure, only neutrinos arising from Z an W ecays at Monte Carlo truth level shoul be consiere. Charge leptons Distributions in the lepton p T an other leptonic observables are known to be highly sensitive to QED raiative corrections, an the ifferences in the treatment of QED raiation on Monte Carlo an theory sie can lea to a bias in the reweighting proceure. This shoul be avoie by using resse leptons, i.e. recombining all leptons with nearly collinear photons that lie within a cone R lγ = φ 2 lγ + η2 lγ < R rec. (7

6 6 For the raius of the recombination cone we employ the stanar value R rec =., which allows one to capture the bulk of the collinear final-state raiation, while keeping contamination from large-angle photon raiation at a negligible level. All lepton observables as well as the kinematics of reconstructe W an Z bosons are efine in terms of resse leptons, an, in accorance with stanar experimental practice, both muons an electrons shoul be resse. In this way ifferences between electrons an muons, l = e, µ, become negligible, an the reweighting function nees to be compute only once for a generic lepton flavour l. Similarly as for neutrinos, only charge leptons that arise from Z an W ecays at Monte Carlo truth level shoul be consiere. Concerning QCD raiation in the vicinity of leptons, no lepton isolation requirement shoul be impose in the context of the reweighting proceure. Instea, in the experimental analysis lepton isolation cuts can be applie in the usual manner. Z an W bosons The off-shell four-momenta of W an Z bosons are efine as p µ W + = p µ l + + p µ ν l, p µ W = p µ l + p µ ν l, (8 p µ Z = pµ l + + p µ l, p µ Z = pµ ν l + p µ ν l, (9 where the leptons an neutrinos that result from Z an W ecays are efine as iscusse above. Photons At higher orers in QCD, photon prouction involves final-state q qγ splittings that lea to collinear singularities when QCD raiation is emitte in the irection of the photon momentum. Since such singularities are of QED type, they are not cancelle by corresponing virtual QCD singularities. Thus, in orer to obtain finite preictions in perturbation theory, the efinition of the pp γ+ jet cross section requires a photonisolation prescription that vetoes collinear q qγ raiation while preserving the cancellation of QCD infrare singularities. To this en, in this stuy we aopt Frixione s isolation prescription [27], which limits the haronic transverse energy within a smooth cone aroun the photon by requiring p T,i Θ(R R iγ ( i=partons/harons ( n cos R ɛ p T,γ R R, cos R where the sum runs over all quarks/gluons an harons at parton level an Monte Carlo level, respectively, while p T,i an p T,γ enote the transverse momenta of partons/harons an photons. The p T -fraction ε, the cone size R, an the exponent n are free parameters that allow one to control the amount of allowe QCD raiation in the vicinity of the photon. The photon-isolation prescription is applicable to QCD as well as to EW higher-orer corrections. At N EW, γ+ jet prouction involves bremsstrahlung contributions with two final-state photons. In this case, at least one isolate photon is require. The other photon might become soft, guaranteeing the cancellation of relate soft an collinear singularities in the virtual EW corrections. In case of two isolate photons in the final state, the harest photon is consiere. In particular, an explicit photon isolation prescription is manatory at N EW in orer to prevent uncancelle singularities from q qγ splittings in the O(α 2 α S mixe EW QCD contributions from qq qqγ an crossing-relate channels. As a consequence of q qγ collinear singularities an the nee to apply a photon isolation prescription, QCD corrections to pp γ+ jet behave ifferently as compare to Z/W + jet prouction. Such ifferences can be important even at the TeV scale, where one might naively expect that massive an massless vector bosons behave in a universal way from the viewpoint of QCD ynamics. Instea, the presence of collinear q qv singularities at (NN QCD implies a logarithmic sensitivity to the vector-boson masses, which results, respectively, in ln(r an ln(p T,V /M V terms for the case of massless an massive vector bosons at p T,V M W,Z. A quantitative unerstaning of these ifferences an their implications on the correlation of QCD uncertainties between γ+ jet an Z+ jet prouction is crucial for the extrapolation of γ+ jet measurements to Z+ jet ark-matter backgrouns. To this en, as iscusse in Section 4, we propose a systematic approach base on the iea that, at large p T,V, the pp γ+ jet process can be split into a ominant part with universal QCD ynamics (in the sense that QCD effects in γ+ jet an Z/W + jet prouction are strongly correlate an a remnant contribution that has to be hanle as uncorrelate in the treatment of QCD uncertainties. To achieve this, we introuce a moifie photon isolation prescription, which is esigne such as to rener the QCD ynamics of γ+ jet an Z/W + jet prouction as similar as possible at high p T. To this en we efine a ynamic cone raius R yn (p T,γ, ε = M Z p T,γ ε, (

7 7 which is chosen in such a way that the invariant mass of a photon-jet pair with R γj = R yn an p T,j = ε p T,γ correspons to the Z-boson mass, i.e. M 2 γj p T,γ p T,j R 2 γj = ε p 2 T,γR 2 yn = M 2 Z, (2 where the first ientity is vali in the small-r approximation. In this way, using a smooth isolation with R = R yn (p T,γ, ε mimics the role of the Z- an W - boson masses as regulators of collinear singularities in Z/W +jet prouction at high p T, while using a fixe cone raius R woul correspon to an effective M γj cut well beyon M Z,W, resulting in a more pronounce suppression of QCD raiation in γ+ jet prouction as compare to Z/W + jet. Specifically, as efault photon selection for the theoretical preictions 4 in this stuy we use the ynamic cone isolation efine through Eq. ( an Eq. (, with parameters ε,yn =., n yn =, R,yn = min {, R yn (p T,γ, ε yn, }. (3 Note that, in orer to prevent that the veto against collinear QCD raiation is applie to an excessively large region of phase space, the ynamic cone raius in Eq. (3 is limite to R yn. As a result of this upper boun, for p T,γ < M Z ε /2,yn 29 GeV the cone raius is kept fixe, an the impact of collinear QCD raiation starts to be significantly enhance as compare to the case of Z/W + jet prouction. Vice versa, for p T,γ > M Z ε /2,yn, thanks to the ynamic isolation cone (3, QCD effects in γ+ jet an Z/W + jet prouction become closely relate, an the egree of correlation between QCD uncertainties across all V + jet processes can be escribe with the prescription of Eqs. (37 (38. For a realistic assessment of theoretical uncertainties, one shoul also consier the fact that photon isolation prescriptions use in experimental analyses iffer in a significant way from the ynamic prescription of Eq. (3. To this en, we recommen to repeat the reweighting proceure using theory preictions for γ+ jet base on a stanar Frixione isolation ( with fixe cone raius an parameters that mimic typical experimental selections at particle level [28], ε,fix =.25, n fix = 2, R,fix =. (4 The ifference between γ+ jet MC samples reweighte in the ynamic- an fixe-cone setup shoul be taken as an aitional uncertainty for pp γ+ jet. As ingreients for this uncertainty estimate we provie higherorer QCD preictions (without uncertainties with fixecone isolation (4 besies the full set of pp γ+ jet 4 The same isolation prescription use for theory preictions shoul be applie also to their MC counterparts σ MC /x in the context of the reweighting proceure. preictions an uncertainties with ynamic photon isolation (see Appenix A. In the EW corrections, ifferences between the two photon isolation prescriptions are well below the percent level. Thus preictions for γ+ jet at (nn EW are provie only with the ynamic cone prescription of Eq. (3. In Figure 2 we present a comparison of the N QCD K-factors for W/Z+ jet an γ+ jet prouction with ynamic an fixe cone isolation. For p T,γ < 29 GeV, where both isolation prescriptions correspon to a fixe cone raius, the QCD corrections to pp γ+ jet grow rapily with ecreasing p T. At low p T, ue to the smaller cone size, fixe isolation (R = leas to more pronounce corrections as compare to ynamic isolation (R =, but the slopes of the corresponing γ+ jet K-factors are quite similar to each other an very ifferent as compare to the ones for pp W/Z+ jet. In the case of fixe isolation, this ifference persists also in the high-p T regime (apart form the acciental agreement of K-factors at p T,V 8 GeV. Instea, in the case of ynamic photon isolation, at large p T the QCD corrections to γ+ jet an W/Z+ jet prouction turn out to be remarkably similar, both in shape an size. As expecte, the onset of this universal behaviour is locate close to p T,γ = 29 GeV, where the isolation raius R,yn starts varying with p T in a way that rejects QCD raiation with M γj < M W,Z. The ifferences between γ+ jet an W/Z+ jet K-factors remain as small as a few percent up to the TeV scale. QCD partons an photons insie jets In orer to avoi any bias ue to the ifferent moelling of jets in MC simulations an perturbative calculations, theory calculations an reweighting shoul be performe at the level of inclusive vector-boson p T istributions, without imposing any requirement on the recoiling jet(s. Preictions presente in this stuy are thus inepenent of specific jet efinitions or jet cuts. Concerning the composition of the recoil, we observe that, at N 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 % of the recoil momentum. Such contributions involve soft QCD singularities that are cancelle by incluing also virtual QCD corrections to pp V γ. In orer to minimise ouble counting with iboson prouction, 5 V γ prouction at is not inclue in the 5 Diboson backgrouns, incluing pp V γ, can be inclue through separate Monte Carlo samples in the experimental analyses.

8 σ N QCD V+j /σ V+j Z+j W+j γ (fix +j γ (yn +j K QCD γ+j /K QCD V+j K QCD γ (fix +j /KQCD W+j K QCD γ (fix +j /KQCD Z+j K QCD γ (yn +j /KQCD W+j K QCD γ (yn +j /KQCD Z+j Fig. 2: Comparison of N QCD K-factors (left for W + jet, Z+ jet, an γ+ jet prouction with ynamic photon isolation (3 an stanar fixe-cone isolation (4. On the right corresponing ratios of K-factors are shown, an the otte lines inicate K-factor variations of ± σ Vγ /σ Vj γ W W + Z 2 3 p T [GeV] Fig. 3: Ratios of istributions in the vector-boson transverse momenta for pp V γ versus pp V j at with µ R,F = H T /2. The vector bosons V = W ±, Z, γ are on shell an s = 3 TeV. EW corrections to pp V j. In practice, as emonstrate in Figure 3, the relative weight of pp V γ at O(α 2 versus pp V j at O(αα S is well below the percent level. Thus the impact of O(α 2 α S contributions from har V γ prouction, which are inclue in this stuy, shoul be completely negligible. 3.2 Cuts an observables Theoretical calculations an the reweighting of MC samples shoul be performe in a fully inclusive V + jet setup, imposing a single cut perturbative preictions suffer form the lack of QCD resummation. 6 For leptons an MET we o not apply any p T or rapiity cuts. Moreover, we o not impose any restrictions on QCD raiation in the vicinity of leptons an MET. Also QCD raiation is hanle in a fully inclusive way, i.e. the presence of a recoiling jet is not explicitly require, an, as iscusse in Section 3., at N EW the recoil can be entirely carrie by a photon. Here we want to stress again that of course the particle-level analysis of the reweighte Monte Carlo samples can (an will involve a more exclusive event selection than use for the reweighting itself. The ifferential istributions to be use for the reweighting of the various pp V + jet processes an process-specific selection cuts to be applie in aition to Eq. (5 are summarise in Table. In the case of pp ν ν+ jet all three neutrino species are ae, while for all other Z an W ecays only a single lepton generation is consiere. For pp l + l + jet an extra invariant-mass cut is applie in orer to avoi far off-shell contributions, especially from γ l + l at low invariant mass. The relatively low value of the lower cut, m ll > 3 GeV, is intene to minimise cross section loss ue to photon raiation that shifts events from the Z-peak region own to lower invariant mass (see Figure 4. This choice guarantees a reuce sensitivity with respect to the moelling of QED raiation. p T,V > 3 GeV for V = W ±, Z, γ, (5 with p T,W ± an p T,Z efine as in Section 3.. The cut (5 is crucial in orer to avoi the region where 6 See e.g. the comparison of NNPS against fixe-orer preictions in Figure 3 of Ref. [29].

9 9 process extra cuts observable comments pp l + ν l + jet none p T,l+ ν l l = e or µ pp l ν l + jet none p T,l ν l l = e or µ pp ν l ν l + jet none p T,νl ν l l = e + µ + τ pp l + l + jet m ll > 3 GeV p T,l+ l l = e or µ pp γ+ jet ynamic isolation p T,γ ( (3 Table : Extra selection cuts, in aition to Eq. (5, an observables for the various V + jet processes. Alternative preictions for γ+jet prouction are provie also for the case of a stanar Frixione isolation with parameters (4. pp Z(l + l + 3 TeV σ/m l + l [pb/gev] N EW σ/σ m l + l [GeV] Fig. 4: Dilepton invariant-mass istribution in pp l + l +jet for m ll [3, 2] GeV comparing an N EW. Collinear lepton photon pairs with R γl <. are recombine. The following binning is aopte for istributions in the reconstructe vector-boson transverse momenta, p T [3, 4,..., 4, 5, 2, 25..., 95,, GeV, 2, 3, 4, 6..., 28, 3, 65]. 3.3 Input parameters, PDFs an QCD scales (6 Input parameters an PDFs employe for theoretical preictions in this stuy are specifie in the following. Let us recall that, as iscusse in Section 2, Monte Carlo samples use in the experimental analyses o not nee to be generate with the same input parameters an PDFs use for higher-orer theoretical preictions. In the calculation of pp ll/lν/νν/γ + jet we use the gauge-boson masses [3] M Z = GeV, M W = GeV, (7 an the corresponing withs, Γ Z = GeV, Γ W = GeV. (8 The latter are obtaine from state-of-the-art theoretical calculations. For the top-quark [3] an Higgs-boson [3] masses an withs we use M t = 73.2 GeV, M H = 25 GeV, (9 an 7 Γ t =.339 GeV, Γ H = GeV. (2 7 Besies loop iagrams with top quarks an Higgs bosons, the N EW corrections to pp W ± + jet receive O(α 2 α S bremsstrahlung contributions from qb q W ± b channels that involve s-channel top-quark propagators an thus require a finite top-quark with, for which we use the N QCD value Γ t =.339 GeV. However, at the perturbative orer consiere in this stuy, such topologies arise only in QCD EW interference terms that o not give rise to Breit Wigner resonances. The epenence of our results on Γ t is thus completely negligible.

10 All unstable particles are treate in the complexmass scheme [32], where with effects are absorbe into the complex-value renormalise masses µ 2 i = M 2 i iγ i M i for i = W, Z, t. (2 For W +jet an Z+jet prouction processes the EW couplings are erive from the gauge-boson masses an the Fermi constant, G µ = GeV 2, using 2 sin 2 θ w µ 2 W α = G µ π, (22 while for γ+jet prouction the EW coupling is chosen to be [3] α = α( = / (23 In both schemes the weak mixing angle θ w is etermine by sin 2 θ w = cos 2 θ w = µ2 W µ 2, (24 Z an becomes complex-value. The G µ -scheme guarantees an optimal escription of pure SU(2 interactions at the EW scale. It is the scheme of choice for W + jet prouction, an it provies a very goo escription of Z + jet prouction as well. The α( scheme to be use for γ+jet, on the other han, expresses the fact that onshell photons effectively couple at a scale Q 2 =. The CKM matrix is assume to be iagonal an we checke at an N QCD that for W +jet prouction the ifference with respect to a non-iagonal CKM matrix is always well below %. For the choice of renormalisation an factorisation scales an variations thereof we refer to Section 4.. For the calculation of haron-level cross sections at (NN QCD + (nn EW we employ the LUXqe_ plus_pdf4lhc5_nnlo_ PDF set, which is base on PDF4LHC NN PDFs [33 38] supplemente with QED effects [39]. The same PDF set, an the relate α S value, is use throughout, i.e. also in the relevant an N ingreients use in the estimate of theoretical uncertainties. At the level of precision iscusse in this stuy also the uncertainty on the value of α S becomes relevant. Given % uncertainty on the measure value of α S this results in an overall 2% normalisation uncertainty on the ifferential p T istributions. However, one shoul keep in min that in the process ratios this uncertainty cancels completely an thus it is irrelevant for backgroun estimates in DM searches at high-met. Consistently with the five-flavour number scheme employe in the PDFs, b-quarks are treate as massless partons, an channels with initial-state b-quarks are taken into account. All light quarks, incluing bottom quarks, are treate as massless particles, an top-quark loops are inclue up to N throughout. Matrix elements at (NN are evaluate using the five-flavour running of the strong coupling supporte by the PDFs an, for consistency, top-quark loops are renormalise in the ecoupling scheme. For the NN QCD coefficient no top-quark loops are consiere. For the assessment of PDF uncertainties the PDF4LHC prescription [33] is aopte. In aition to stanar PDF variations, also aitional LUXqe variations for the photon PDF are applie. For more etails see more etails in Sections Computational frameworks The theoretical preictions presente in Section 4 inclue corrections up to NN QCD an N EW, as well as Suakov EW effects at O(α 2. They have been obtaine by means of a variety of methos an tools, as etaile in the following. The N QCD an N EW calculations for all pp V + jet processes have been performe with Munich+OpenLoops an/or Sherpa+OpenLoops. In these automate frameworks [9, 4, 4] virtual amplitues are provie by the OpenLoops program [42, 43], combine with the Collier tensor reuction library [44] or with CutTools [45]. The remaining tasks are supporte by the two inepenent an fully automate Monte Carlo generators Munich [46] an Sherpa [47 5]. Aitionally, we carefully valiate the N EW preictions against the results of Refs. [7, 8, 2]. The N EW calculations for pp V + 2 jets performe to test the factorisation of QCD an EW corrections have been checke against the one of Ref. [5] for pp Z + 2 jets in Ref. [2]. The N EW amplitues for all V +jet processes in OpenLoops have been supplemente with the one- an two-loop analytical Suakov logarithms of Refs. [22 25, 52]. The NN QCD preictions for Z+jet prouction have been obtaine with the parton-level event generator NNjet, which provies the necessary infrastructure to perform fully ifferential calculations at NN using the antenna subtraction formalism [53 6]. The computation of pp W +jet through NN is base on the N-jettiness subtraction scheme for NN calculations [3]. The above-cut contribution within the N- jettiness subtraction was obtaine using Munich+OpenLoops. The NN QCD preiction for the pp γ+jet process is base on the calculations of Refs. [5, 6] an has been obtaine using MCFM [62]. In orer to ensure the correctness of the numerical implementation of cuts an other parameters in the NN

11 coes, a etaile comparison has been performe at the level of the N QCD results as escribe above. 4 Higher-orer QCD an EW preictions Precise theory preictions for V + jet prouction require QCD an EW higher-orer corrections, mixe QCD EW contributions, as well as photon-inuce contributions, TH = QCD + x σ(v EW + x σ(v mix + γ in.. (25 In this section we present theoretical preictions that inclue corrections up to NN QCD an N EW supplemente by EW Suakov logarithms at two loops. Moreover, we introuce a coherent theoretical framework for the combination of EW an QCD calculations for the various V + jet prouction processes an for the assessment of the corresponing remaining sources of theoretical uncertainty. State-of-the art QCD an EW preictions an the relate theoretical uncertainties are iscusse in Section 4. an 4.2 respectively. Section 4.3 is evote to photon-inuce channels an Section 4.4 to PDF uncertainties, while in Section 4.5 we iscuss the real emission of vector bosons, an mixe corrections of O(αα S are aresse in Section 4.6 by means of a factorise combination of QCD an EW corrections. To illustrate the effect of higher-orer corrections an uncertainties we present a series of numerical results for pp V + jet at a centre-of-mass energy of 3 TeV in the setup specifie in Section 3. In particular, pp γ+ jet preictions are base on the ynamic photon isolation (3. As anticipate in Section 3., this prescription provies a very convenient basis for the systematic moelling of the correlation of QCD uncertainties between the various V + jet prouction processes (see Section 4.. Vector-boson p T spectra are plotte starting at 8 GeV, but for the sake of a complete ocumentation ata sets are provie above 3 GeV (see Appenix A. However, we note that in the region of p T < GeV there are potential sources of systematics that we are not controlling or even iscussing, as they woul require a separate stuy. These arise from the resummation of QCD Suakov logarithms or from non-perturbative effects (e.g. an orer Λ QCD average shift of the vector boson p T associate with the asymmetry of colour flow in the final state. Furthermore, as shown later, a reliable correlation between the Z/W spectra an the photon spectrum requires p T to be large enough so that fragmentation contributions in γ+jet prouction become small. We also expect that in the p T regions up to a few hunre GeV the statistics are sufficient to guarantee that experimental analyses of missing-e T backgrouns can entirely rely on the irect measurement of the Z spectrum measure via Z l + l. As a result, we believe that our conclusions on the systematic uncertainties are most reliable an useful for experimental applications in the region of p T larger than 2 GeV. 4. Higher-orer QCD preictions For perturbative QCD preictions at, N an NN we use the generic notation QCD = N k QCD, (26 with k =, or 2. Wherever possible, nominal preictions are provie at NN QCD, i.e. incluing terms up to 8 O(ααS 3. However, as ingreients for the assessment of some theory uncertainties, also an N QCD contributions will be use. For convenience, results at N k QCD are systematically expresse in terms of preictions an relative correction factors efine through (µ = K(V (x, µ QCD (µ. (27 N k QCD N k We calculate all N k an cross sections with one an the same set of NN PDFs as iscusse in Section 3.3. The epenence on the renormalisation an factorisation scales, µ = (µ R, µ F, is absorbe into the K-factors, while preictions on the r.h.s. of Eq. (27 are taken at the central scale, µ = (µ R,, µ F,. For the central scale we aopt the commonly use choice µ R, = µ F, = µ = Ĥ T/2, (28 where the total transverse energy, Ĥ T, is efine as the scalar sum of the transverse energy of all parton-level final-state objects, Ĥ T = E T,V + p T,i. (29 i {q,g,γ} Also quarks (q, gluons (g an photons that are raiate in the (NN QCD or EW corrections are inclue in Ĥ T, an the vector-boson transverse energy, E T,V, is compute using the total (off-shell fourmomentum of the corresponing ecay proucts, i.e. E 2 T,Z = p 2 T,l + l + m2 l + l, E 2 T,W = p 2 T,lν + m 2 lν, E 2 T,γ = p 2 T,γ. (3 8 Here an in the following we aopt a power counting that oes not inclue the extra factor α associate with vector-boson ecays.

12 2 In orer to guarantee infrare safety at N EW, the scale (29 must be insensitive to collinear photon emissions off charge fermions. To this en, the vector-boson transverse energies efine in Eq. (3 shoul be compute in terms of resse leptons as specifie in Section 3., while p T,γ contributions to Eq. (29 shoul involve only photons that have not been recombine with charge leptons. It is worth to note that µ p T,V at large p T,V. Pure QCD uncertainties The uncertainty associate with the truncation of the perturbative expansion in α S is estimate by means of factorisation an renormalisation scale variations. We consier stanar seven-point variations applying, respectively, factor-two rescalings, i.e. µ i µ = (,, (2, 2, (.5,.5, (2,, (, 2, (,.5, (.5,, (3 where i =, Nominal preictions an relate uncertainties are efine as the centre an the half-with of the ban resulting from the above variations. In terms of K-factors this correspons to K (V N k (x = 2 δ ( K (V N k (x = 2 with [ [ K (V,max N k K (V,max N k K (V,max N k (x = max { K (V ] (x + K(V,min N k (x, (32 ] (x K(V,min N k (x, (33 } N k (x, µ i i 6, { K (V,min N k (x = min K (V N k (x, µ i i 6 }. (34 Since the shift resulting form the symmetrisation of scale variations in Eq. (32 is encoe in the K-factors, also the K-factor iffers from one. Constant scale variations mainly affect the overall normalisation of p T -istributions an ten to unerestimate 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 introuce an aitional variation, δ (2 K (V N k (x = ω shape(x δ ( K (V (x, (35 N k where the stanar scale uncertainty (33 is supplemente by a shape istortion ω shape (x, with ω shape (x an ω shape (x ± at high an small transverse momentum, respectively. The function ω shape is efine as ω shape (p T = tanh [ ( ] pt ln p T, = p2 T p2 T, p 2 T +, (36 p2 T, an as reference transverse momentum we choose the value p T, = 65 GeV, which correspons (in logarithmic scale to the mile of the range of interest,.2 2 TeV. As illustrate in Figure 5, the function ω shape (x inuces asymmetric variations that cover ±75% of the stanar scale variation ban for p T [25, 75] GeV. Note that, in the combination of the uncertainties (33 an (35, our choice to have an aitional shape variation augments the stanar scale uncertainty by a factor + ωshape 2 (p T 2. Besies shape uncertainties, also the correlation of QCD uncertainties across V + jet processes plays a key role in fits of the Z(ν ν+ jet ark matter backgroun, an the quantitative unerstaning of such process correlations belongs to the most important theoretical aspects in ark matter searches. From the viewpoint of QCD interactions, the processes pp W + jet an pp Z+ jet are quite similar to each other at p T,V M W,Z. Thus, the respective QCD uncertainties are expecte to be strongly correlate. However, ue to the presence of q qγ collinear singularities an the nee to suppress them with an appropriate photon-isolation prescription, higher-orer QCD contributions to γ+ jet prouction can behave in a significantly ifferent way as compare to the case of pp W/Z+ jet. In orer to reuce such ifferences, we aopt the ynamic photon isolation approach efine in Eq. (3. As iscusse in Section 3., this prescription reners the QCD ynamics of pp γ+ jet an pp Z/W + jet processes almost universal. As a result, QCD K-factors K (V N k (x an their uncertainties δ (i K (V (x epen only very N k weakly 9 on V at high p T, an in this situation the small process-epenent part of QCD K-factors can be use as an estimator of the egree of correlation across processes. To this en we consier the highest available term in the perturbative expansion, K (V (x = N k K(V N k (x/k(v (x, (37 N k an as estimate of unknown process correlation effects we take the ifference of the known QCD K-factors with respect to Z+ jet prouction, δ (3 K (V (x = N k K(V (x N k K(Z (x. (38 N k This process correlation uncertainty can be assesse using the central scale (28 throughout. Applying it to nominal preictions, i.e. replacing K (V N k K(V N k ± δ (3 K (V, amounts to oubling or removing K-factor N k ifferences between processes. The choice of Z+ jet prouction as reference process in Eq. (38 is arbitrary, but 9 For what concerns process correlations, it is crucial that (apart from the M V epenence all V + jet processes are evaluate using equivalent ynamical scales.