BSM WW production with a jet veto

We consider the impact on WW production of the unique dimension-six operator coupling gluons to the Higgs field. In order to study this process, we have to appropriately model the effect of a veto on additional jets. This requires the resummation of large logarithms of the ratio of the maximum jet transverse momentum and the invariant mass of the W boson pair. We have performed such resummation at the appropriate accuracy for the Standard Model (SM) background and for a signal beyond the SM (BSM), and devised a simple method to interface jet-veto resummations with fixed-order event generators. This resulted in the fast numerical code MCFM-RE, the Resummation Edition of the fixed-order code MCFM. We compared our resummed predictions with parton-shower event generators and assessed the size of effects, such as limited detector acceptances, hadronisation and the underlying event, that were not included in our resummation. We have then used the code to compare the sensitivity of WW and ZZ production at the HL-LHC to the considered higher-dimension operator. We have found that WW provides complementary sensitivity with respect to ZZ, provided one is able to control theory uncertainties at the percent-level. Our method is general and can be applied to the production of any colour singlet, both within and beyond the SM.


Introduction
Di-boson production at the Large Hadron Collider constitutes a promising window into physics beyond the SM. This is particularly true for di-boson pairs with high invariant mass, which have been already probed by a number of recent experimental analyses [1][2][3][4][5][6][7][8][9][10][11]. On the one hand, their production through gluon fusion receives contributions from an off-shell Higgs boson [12][13][14]. In particular, the interference of the contribution of an off-shell Higgs boson and di-boson continuum background makes it possible to access the Higgs width in a model-independent way [15]. On the other hand, contact interactions arising from higher-dimensional effective field theory operators [16][17][18][19][20] could give rise to spectacular effects in the tails of di-boson differential distributions, due to the fact that their contribution increases with energy. Technically, in the SM, di-boson production via gluon fusion is a loop-induced process. At low di-boson invariant masses, top quarks in the loops behave as very heavy particles, thus giving rise to effective contact interactions. At high invariant masses, the two bosons probe virtualities that are much larger than the masses of the top quarks running in the loops, hence suppressing their contribution and enhancing the effect of BSM contact interactions. Such a feature has been already used to constrain the coefficient of a number of higher-dimensional operators, see e.g. [21] for a recent study.
In this article we restrict ourselves to considering the unique dimension-six operator coupling gluons to the Higgs boson, given by [17] L ⊃ c gg Λ 2 G a µν G a,µν φ † φ, (1.1) with G a µν the gluon field strength and φ the Higgs field. This operator can be used to represent contributions to SM Higgs production from particles with mass of order Λ m H . This operator has previously been considered in high-invariant-mass ZZ production with a fully leptonic final state in [22,23]. However, the leptonic final state for W W has larger cross section and so W W could give complementary or better sensitivity than leptonic final states for ZZ. However, in W W production, a tight jet veto is employed by experiments to suppress background from top-pair production. Such a veto "forbids" the radiation of jets from the initial-state partons, with the effect of suppressing not only the background, but also the operator-mediated signal. In the present case, the signal occurs through gluon fusion, whereas W W production is mainly driven by quark-antiquark annihilation. Since gluons radiate more than quarks, one expects the suppression due to a jet veto to be stronger for the signal than for the background. It is therefore important to address the general question of how BSM searches with W W production compare to ZZ in the presence of a jet veto. 1 The aim of this paper is to quantify in a simple way how the significance of such a BSM signal is affected by the presence of a jet veto. The same procedure can be applied to any BSM scenario that modifies the production rate of a colour singlet, for instance dimension-8 operators [31]. A similar study [32] investigates the impact of a jet veto in the determination of the Higgs width using interference. To be more specific, we veto all jets that have a transverse momentum (with respect to the beam axis) above p t,veto . First we observe that, at the level of the matrix element squared, a generic BSM signal mediated by a single higher-dimensional operator consists of an interference piece and a quadratic piece: The last piece is of higher order 1/Λ 4 . Therefore, if the interference piece is not suppressed or vanishing for some reason, then, to a first approximation, we can neglect it relative to the 1/Λ 2 interference piece. 2 The presence of a jet veto induces large logarithms of the ratio of p t,veto and the invariant mass of the W W pair M W W . Such logarithms arise at all orders in QCD, and originate from vetoing soft-collinear parton emissions. Considering just the leading logarithms, and neglecting the quadratic piece |M BSM | 2 , the deviation of a BSM signal that proceeds from gluon fusion from the SM prediction is approximately given by where C A = 3, α s is the strong coupling, and L gg (M W W ) is the gluon-gluon luminosity corresponding to a partonic centre-of-mass energy equal to M W W . The effect of the jet veto is an exponential (Sudakov) suppression with respect to a naive Bornlevel estimate. Note also that, for fixed p t,veto , such a suppression becomes more and more important, the higher the invariant mass of the W W pair. This is precisely where the contribution of the BSM operator in eq. (1.1) has the most impact on the signal. For the SM background, dominated by quark-antiquark annihilation, we have instead a contribution proportional to For α s = 0.1, M W W = 1 TeV, p t,veto = 20 GeV, the above factor is about 0.2. Therefore, despite the gain in the number of events one has in W W production with respect to ZZ, the significance of the signal might be reduced due to jet-veto effects. This is why it is crucial to have an estimate of jet-veto effects that is as accurate as possible.
The first question we address is what accuracy we can aim for in the description of a BSM signal and a QCD background involving the production of a colour singlet. In the absence of large jet-veto corrections, a generic BSM signal can be predicted at Born-level, or leading order (LO), in QCD, whereas any QCD background is nowadays known at least at next-to-leading order (NLO). In the presence of a jet veto, the production of a system of invariant mass M is affected by logarithms of the ratio p t,veto /M , which make fixed-order predictions unreliable. After the all-order starting point is to observe that, in eq. (1.6), the factor multiplying leading logarithms is in fact a new perturbative series, whose coefficients are functions of α s L. As stated previously, NLL corrections have the same structure as Born-level contributions, while NNLL corrections closely resemble NLO contributions. Therefore, NLL resummation could just be obtained by an event-by-event reweighting of a Born-level generator by keeping only the functions g 1 and G 2 in eq. (1.6). This is enough to estimate jet-veto effects to the BSM production of a colour singlet. Including NNLL corrections, needed for a precise estimate of the corresponding SM background, is also possible in a general way. In fact, resummation effects originate from soft and/or collinear emissions in such a way that NNLL corrections share the same phase space with Born-level contributions, but are of relative order α s . In all NLO calculations there is always a contribution that lives in the same phase space as the Born, and is of relative order α s . This is the subtraction term that cancels the infrared singularities of virtual corrections. Therefore, to implement NNLL effects, we can just modify the appropriate subtraction term in the NLO event generator. Having done this, all other NNLL effects factorise, and can be accounted for by an event-by-event reweighting, so as to reproduce eq. (1.6). The whole procedure requires generating Born-level events only, and hence is much faster than a full NLO calculation.
In the following two sections we give a detailed description of this procedure for the specific case of BSM effects induced by the operator in eq. (1.1). In section 2, we study the effect of such an operator on W W production with a jet veto. As discussed above, this operator induces a modification of the cross section of W W production through gluon fusion. We denote the (differential) cross section for gluon fusion, potentially including an additional BSM contribution, with dσ gg . The main result of this section is a recipe to compute cross sections for W W production with a jet veto at NLL accuracy, fully exclusive in the decay products of the W bosons. In section 3 we compute the cross section for the dominant contribution to the SM background, which is W W production via quark-antiquark annihilation, again in the presence of a jet veto. We denote the cross-section for this process with dσ qq , and compute exclusive cross sections in the decay products of the W bosons, while resumming ln(M W W /p t,veto ) at NNLL accuracy. The main result of this section is a general recipe to modify a NLO event generator for the production of any colour singlet so that it produces resummed cross-section with a jet veto at NNLL accuracy. In section 4 we present some numerical results for a simplified model derived from the Lagrangian in eq. (1.1), corresponding to a realistic experimental setup. We compare our resummed predictions with parton-shower event generators, and assess the size of effects, such as limited detector acceptances, hadronisation and the underlying event, that are not included in our resummation. In section 5 we perform some basic sensitivity studies to investigate the exclusion potential of the HL-LHC for the parameters of the simplified model of section 4. Finally, section 6 presents our conclusions.

Gluon fusion (including BSM effects)
Let us first consider W W production via gluon fusion, possibly with a modification of the amplitude induced by the BSM operator in eq. (1.1). For simplicity, we consider here the decays W W → e + ν e µ −ν µ and W W → e −ν e µ + ν µ . As explained in the introduction, if we impose that all jets have a transverse momentum below a threshold value p t,veto , the distribution in M 2 W W , differential in the phase space of the leptons, is affected by the presence of large logarithms ln(M W W /p t,veto ), that have to be resummed to all orders to obtain sensible theoretical predictions. Specifically, we consider jets obtained by applying the anti-k t algorithm [40] with a given radius R. At NLL accuracy, the best we can achieve for gluon fusion, the aforementioned observable is given by [35,41] where L = ln(M W W /p t,veto ), α s = α s (M W W ), and explicit expressions for the functions g 1 (α s L) and g 2 (α s L) can be found, for instance, in ref. [35]. In particular, they are the same for any colour singlet that is produced via gluon fusion (e.g. Higgs production). Note that, at NLL accuracy, the resummed distribution in eq. (2.1) does not depend on the radius R of the jets [41]. The phase space of the leptons is given by 2) with p = (E , p ) is the four-momentum of lepton = e, µ, ν e , ν µ , and p i = x i P i , i = 1, 2 are the momenta of the incoming partons, carrying each a fraction x i of the incoming proton momentum P i .
Last, we have a process dependent "luminosity" factor L (0) gg , given by 4 The two main ingredients entering L (0) gg are: • the SM amplitude M (gg) SM for the production of a W W pair (and its decay products) through gluon fusion, which can be supplemented with an additional contribution M (gg) BSM accounting for BSM effects; • the gluon density in the proton f g (x, µ F ) at the factorisation scale µ F = p t,veto .
This value of µ F reflects the fact that the factorisation scale is the highest scale up to which the considered observable is inclusive with respect to multiple collinear emissions from the initial-state partons. Since all collinear emissions with a transverse momentum above p t,veto are vetoed, the factorisation scale has to be p t,veto (see e.g. [33] for a formal derivation).
By comparing eq. (2.1) to eq. (1.6), we obtain the function G 2 (α s L) resumming all NLL contributions: So far, with the exception of ref. [29], such resummations have been obtained by devising process-dependent codes that produce numerical results for L For instance, the program JetVHeto [38] returns NNLL resummations integrated over the full phase space of the decay products of a Higgs or a Z boson. However, the luminosity in eq. (2.3) can be obtained by running any Born-level event generator. In fact, any such program will compute a Born-level cross-section in W W production via gluon fusion (possibly with BSM contributions) starting from the formula: where µ F here is the default factorisation scale in the considered Born-level generator. Therefore, to obtain the differential distribution in eq. (2.1), it is enough to set that factorisation scale µ F to p t,veto , and multiply the weight of each phase-space point by exp[Lg 1 (α s L) + g 2 (α s L)]. Note that, if the programs returns event files with information on M W W for each event, or if one produces histograms binned in M W W , the reweighting can be performed without any need to touch the Born-level generator code.

Quark-antiquark annihilation (SM only)
Since SM background processes are typically known at least to NLO, in the presence of a jet veto, the SM cross-section for W W production can be computed at NNLL accuracy. The corresponding NNLL resummed expression is given by ) × e Lg 1 (αsL)+g 2 (αsL)+ αs π g 3 (αsL) , (3.1) where again L = ln(M W W /p t,veto ), α s = α s (M W W ), and dΦ leptons is the lepton phase space defined in eq. (2.2). The functions g 1 , g 2 and g 3 are reported in [35], and are the same as for Drell-Yan production. The dependence on the jet radius R appears for the first time at NNLL accuracy in the functions F clust (R), F correl (R), whose explicit expressions can be found in [41]. At NNLL accuracy we have two process-dependent "luminosities" L The only difference with respect to L ij , which is different from zero only if i, j is a quark-antiquark pair with the same flavour.
At NNLL accuracy we need to add the luminosity L qq , which is of relative order α s with respect to L (0) qq , and is given by Here new ingredients appear: • one-loop virtual corrections to W W production. They are included in the term H (1) , the coefficient of α s (M W W ); • coefficient constants arising from real collinear radiation. They are included in the terms C ik (z), whose explicit expressions can be found in ref. [35], and are the same as for Drell-Yan production. They multiply α s (p t,veto ), which reflects the fact that the characteristic scale of collinear radiation in jet-veto cross sections is p t,veto .
With reference to eq. (1.6), the function G 2 resumming NLL contributions is whereas the function G 3 resumming NNLL contributions is As explained in the previous section, the function L (0) qq can be obtained from an appropriate Born-level program. The function L (1) qq instead represents a correction to L (0) qq of relative order α s , that cannot be obtained from a LO calculation. A viable possibility to perform NNLL resummation would be to modify eq. (2.3) so that it includes the convolutions over the variable z in eq. (3.3), and implement the modification in a Born-level generator. This is the approach taken in ref. [29], and in some way underlying the current implementation of the JetVHeto program [38].
Here we want to present an alternative procedure. First, let us consider how the NLO W W cross section is calculated in a NLO event generator: The first term in the sum is the LO SM cross section dσ qq,r /dΦ leptons dM 2 W W , represents NLO corrections coming from the emission of an extra parton. They include the counterterms needed to ensure their finiteness in four space-time dimensions. The second term, dσ (1) qq,v+ct /dΦ leptons dM 2 W W , gives NLO corrections arising from the sum of virtual corrections, and the counterterms integrated over the full extra-parton phase space. This contribution lives in the same phase-space as the Born contribution, and is of relative order α s . It has the form (3.8) In the above equation, µ R , µ F are the renormalisation and factorisation scales used by the NLO generator,H (1) represents virtual corrections to qq → W W , andC (1) ik (z) the integrated counterterms. The explicit expressions ofH (1) andC (1) ik (z) depend on their actual implementation in the NLO generator, in particular on the employed subtraction scheme. However, the form of eq. (3.8) is the same as that of the NNLL luminosity L and by evaluating the parton distribution functions at the factorisation scale µ F = p t,veto . Finally, in order to obtain the resummed distribution in eq. (3.1), we need to reweight each phase space point by This rescaling can also be performed when constructing histograms, as long as one has access to M W W for each bin, or for each event in an event record.
We have implemented this procedure in the code MCFM-RE [43], a suitable modification of the NLO program MCFM [44]. The actual implementation is richer than what has been discussed so far, because it allows a user to change the default renormalisation and factorisation scales, and contains additional features. Since these details are not relevant for a general discussion, we have omitted them here. The interested reader is referred to appendix A for the actual formulae we implement, and to appendix B for a short manual of the code.
In the following two sections, we use this implementation to produce numerical results and sensitivity studies for an explicit BSM model.

Numerical results
Let us discuss first our results for W W production via qq annihilation. We consider W pairs produced at the LHC with √ s = 13 TeV, specifically W + W − → e + ν e µ −ν µ and W + W − → µ + ν µ e −ν e , and select the final state according to a simplified version of the experimental cuts of ref. [3], reported in table 1. Jets are reconstructed according to the anti-k t algorithm [40] with a jet radius R = 0.4. In table 1 we encounter the newly introduced observable E / T,Rel , which is defined as follows [45]: , and φ e , φ µ and φ MET the azimuthal angle of the electron, the muon and the missing transverse energy respectively. The simplified cuts in table 1 slightly differ from the ATLAS ones. First, ATLAS vetoes only jets with |η| < 4.5. This does not cause problems for our resummed calculation because, according to the arguments in ref. [46], it just limits its validity to the Fiducial selection requirement Cut value Table 1. Definition of the W W → eµ fiducial phase space, where p T , η are the transverse momentum and rapidity of either an electron or a muon, M eµ is the invariant mass of the electron-muon pair, E / T is the missing transverse energy, and E / T,Rel is defined in eq. (4.1).
range ln(M W W /p t,veto ) < 4.5, which is within the region we consider. However, ATLAS employs an additional cut on the jets, vetoing also jets with p T > 25 GeV and |η| < 2.5. If we compute dσ/dM W W with the cuts in table 1, we miss a (formally NNLL) contribution of order exp[−C(α s /π)∆η ln(30 GeV/25 GeV)], with C = C F or C = C A according to whether we have quarks or gluons in the initial state and ∆η the size of rapidity region in which the jet veto cuts differ, in this case ∆η = 5. Last, the definition of E / T,Rel used to define the cuts in table 1 considers only leptons, whereas ATLAS considers all reconstructed particles, including jets. This leads to small NNLL corrections that depend on the area in the η-φ plane occupied by the rejected jets. We will investigate the impact of these effects using parton-shower event generators. We also omit b quark-initiated contributions to pp → W W . At LO, the bb scattering subprocess contributes only 1% to the cross section. The gb and gb subprocesses, which enter at NLO QCD increase the NLO cross section by a factor 1.5. This large increase is due to graphs like gb → W − (t → W + b). Such graphs feature a resonant top quark propagator, which effects an enhancement of O(m t /Γ t ) = O(10 2 ), which compensates the O(1%) suppression due to the b PDF, and altogether an O(1) contribution is obtained. This contribution is commonly attributed to W t production and decay (at LO QCD) [47], and hence has to be omitted in the NLO QCD corrections to W W production, which we consider here.
Given these considerations, we produce both NLO, NNLL resummed, and matched NLO+NNLL (with the matching procedure explained in appendix A.3) predictions for the differential distribution dσ/dM W W using PDF4LHC15 parton distribution functions (PDFs) at NLO [48], accessed through LHAPDF6 [49], corresponding to α s (M Z ) = 0.118, and we set both renormalisation and factorisation scales at M W W /2, as customary in Higgs precision studies [50]. Fig. 1 shows the differential cross section in the invariant mass M W W of the W W pair. We first note that both NLO and NNLL+NLO are both smaller than the LO, as expected due to the presence of a jet veto, with the suppression with respect to LO increasing with M W W . This implies that, in this situation, a naive Born-level calculation fails to capture this effect and that, in the absence of a resummation, one should use at least a NLO prediction. NNLL+NLO gives a mild extra suppression with respect to NLO, revealing that logarithms are not particularly large in the considered kinematical region. However, we note that the difference between pure NNLL resummed and matched NNLL+NLO (the so-called "remainder"), which contains the part of the NLO which is not enhanced by logarithms, is basically negligible. This means that the resummation alone is very close to the best prediction we have at this order. This is remarkable in view of the fact that to obtain NNLL predictions we need to perform a calculation with Born-level kinematics. On the contrary, the computational cost of the NLO calculation is larger due to the presence of an extra emission, without any significant gain in accuracy compared to the NNLL prediction.
To complete our discussion of the qq channel, we compare our predictions to those obtained from SCET via the program aMC@NLO-SCET of ref. [29]. The comparison is shown in Fig. 2. Our results contains theoretical uncertainties evaluated both with the most recent jet-veto efficiency (JVE) method [51] at the relevant accuracy (the wider, lighter band), and pure scale variations (the tighter, darker band). The details of both prescriptions can be found in appendix A. The SCET prediction  corresponds to the default scale choices, and is at the boundary of scale variation uncertainties and well within JVE uncertainties. We remark that we do not expect perfect agreement, because, although both methods share the same formal accuracy, they differ in the treatment of subleading effects.
A last comment on uncertainties is in order here. Within MCFM, we do not have access to NNLO calculations for di-boson production, so we cannot match our resummed predictions to NNLO. As a result of this, the JVE method may be overly conservative, due to the largish (∼ 1.5) K-factor of the W W inclusive total cross-section, which propagates in the evaluation of the uncertainty according to the JVE method. If we could match to NNLO, the JVE uncertainty would be reduced and, as happens for Higgs production [35], would probably get closer to plain scale uncertainties.
In order to have a specific example of a BSM theory that implements the effective operator of eq. (1.1), we consider the following modification of the SM Lagrangian [52]: with t, h, G a µν the top field, the SM Higgs field, and the gluon field strength re-spectively. The SM corresponds to (κ t , κ g ) = (1, 0), and in this section we will only explore BSM scenarios such that κ t + κ g = 1, which ensures that the Higgs total cross section stays unchanged (modulo quark-mass effects, which give a correction of a few percent [51]). Such modifications of the SM Lagrangian only affect the gluon-fusion contribution to di-boson production. Their effect has been investigated before for the case of ZZ production [23], where one does not need to impose a jet veto to suppress unwanted background. Here we wish to study how the presence of a jet veto, required for studies of W W production, affects the relative size of a BSM contribution with respect to the SM background. We consider the three benchmark scenarios studied in ref. [23], i.e.
First, in Fig. 3 we compare the loop-induced gluon fusion contribution to the M W W distribution at LO, which is what is given by default by any automated Born-level event generator, with the NLL analytic resummation, which gives the best modelling of jet-veto effects at the currently available accuracy. Our best qq prediction is also shown for comparison. We see that, if we include resummation effects, the cross   section for each benchmark point is reduced by almost an order of magnitude in the tail of the distribution, where BSM effects start to become important. We then investigate more quantitatively how this impacts the deviations we might observe with respect to the SM, by plotting the quantity is set to zero, and dσ qq follows from eq. (3.1). Fig. 4 (left) shows δ(M W W ) for the benchmark point (0.7, 0.3). We first note the growth of this quantity with energy, as expected from the effective nature of the ggH vertex. Fortunately, the growth persists after including jet-veto effects through NLL resummation, however the deviation from the SM reduces from the 1% that one would obtain using fixed-order calculations (see fig. 3) to fractions of a percent. The same quantity shown in the right panel of fig. 4 for the benchmark point (0.0, 1.0) displays qualitatively the same behaviour, although the deviation is a factor ten bigger. We see that, in the presence of jet-veto restrictions such as the one in ATLAS cuts [3], one is bound to use a theoretical tool that resums large logarithms. This could be either resummed predictions, or simulations with parton-shower event generators.  The variable δ(M W W ) is of theoretical interest only, because we do not have access to the momenta of the neutrinos. To have experimentally accessible observables, we consider differential distributions in M T 1 [53], M T 2 [54] and M T 3 [53], three measurable variables that are strongly correlated with M W W In the above equations, p T,eµ = p T,e + p T,µ , and M 2 eµ = (p e + p µ ) 2 . The vector p / T is the missing transverse momentum, defined as minus the vector sum of all detectable particles. Note that, if no jets are present, as at Born-level and in NNLL resummed predictions, p / T = − p T,eµ . Last, ∆φ eµ,miss is the azimuthal angle between p T,eµ and p / T . The corresponding results for δ are shown in Fig. 5. We note that M T 2 gives rise (κ t ,κ g )=(0.7,0.3) (κ t ,κ g )=(0.0,1.0) Figure 5. The relative difference between BSM and Standard Model W W production, differential in M T 1 (left) and M T 2 (right).
to considerably larger deviations with respect to M T 1 . This is because low values of M T 2 are correlated to larger values of M W W , so M T 2 effectively probes the M W W distribution in the high-mass tail, where BSM effects are appreciable. However, this also means that the differential cross section in M T 2 is much smaller than that in M T 1 , as can be seen from Fig. 6. Therefore, the discriminatory power of M T 2 is only of use if we have a very large number of events. We have also studied the variable M T 3 defined again in ref. [54] and first devised in ref. [55]. The distribution in this variable looks very similar to that of 2M T 1 , so the same discussion as for M T 1 applies here.
We now compare our results to parton-level predictions from parton-shower event generators, using existing tunes. In particular, for qq we consider POWHEG [56][57][58][59] matched to the AZNLO [34] tune of PYTHIA v8.230 [60], and aMC@NLO [61][62][63][64][65] matched to PYTHIA, this time with the default parameters. To investigate the dependence on the shower algorithm, we also consider the parton shower HERWIG v7.1.0 [66,67] matched as POWHEG+HERWIG, and aMC@NLO+HERWIG, both with the default parameters. For POWHEG+PYTHIA, we use the PDF set by the AZNLO tune, i.e. CT10 [68] for POWHEG and CTEQ6L1 [69] for the parton shower. For consistency, we use CT10 everywhere for POWHEG+HERWIG. For POWHEG+HERWIG, we also performed runs with default shower PDFs, and noted no significance difference in the resulting distributions. For all the aMC@NLO runs we use PDF4LHC15 PDFs, both for the generation of the hard configurations and the shower.
The comparison of resummation with event generators is shown in Fig. 7 for the SM (for qq → W W and gg → W W separately), and in Fig. 8 for the two BSM scenarios considered above. Resummed predictions include an estimate of theory uncertainties at the appropriate accuracy, as explained in appendix A.3. Note that, due to the missing NLO total cross-section for the incoming gg channel, JVE and scale uncertainties for gg → W W are of comparable size, with the JVE ones slightly larger. We first observe that, both for qq-and for gg-initiated W W production, all event generators agree with the resummation within its uncertainties. For qq, where we can match parton-shower predictions to NLO, POWHEG+PYTHIA shows a remarkable agreement with the resummation, but other event generators give comparable results. We note that predictions obtained with aMC@NLO show a slightly different trend with M W W . In particular aMC@NLO+PYTHIA is slightly above our central prediction at low M W W , and a bit lower at high M W W , whereas aMC@NLO+HERWIG shows the same trend but is everywhere lower than our predictions.
In the gg case, both for the SM and the considered BSM scenarios, we can only compare to unmatched parton-showers results, as no NLO calculation is available. We observe that PYTHIA is in better agreement with our predictions at large values of M W W , whereas HERWIG's predictions have the same shape as ours, but are systematically lower by about 10%. Overall, there is agreement between our predictions and parton showers within uncertainty bands, so the latter can be reliably used for this process. We remark that parton-shower predictions not only have lower formal accuracy, but are also much more expensive computationally. Hence it might be lengthy to assess with those tools if a range of BSM parameters leads to sizeable deviations from the SM, whereas with our numerical implementation such analyses could be performed at the cost of an unshowered Born-level calculation.
Last, we wish to investigate the impact of actual ATLAS cuts on the jets with respect to the simplified jet-veto cuts needed for the validity of our resummations,         as well as the impact of non-perturbative corrections due to hadronisation and underlying event (UE). In fig. 9 we investigate the effect of different cuts on the jets on dσ/dM W W , using parton shower event generators at parton level, in particular we use POWHEG+PYTHIA for qq and plain PYTHIA for gg. We observe that the rapidity cut |y| < 4.5 has essentially no effect. On the contrary the ATLAS cuts give a sizeable but constant extra suppression. This is reasonable given that the jet veto cut imposed by ATLAS is more stringent, since jet veto in the central region |y| < 2.5, is taken to be p t,veto = 25 GeV. In the case of gg, the suppression is larger with respect to qq due to the larger colour factor of the initial-state gluons with respect to the quarks. In fig. 10 we investigate the effect of hadronisation and underlying event on dσ/dM W W , using parton shower event generators. Again we make use of POWHEG+PYTHIA for qq and plain PYTHIA for gg. We observe that hadronisation corrections are essentially negligible, which is expected since they scale like inverse powers of the hard scale, in this case M W W . Corrections arising from the underlying event are a few percent, smaller than the typical theoretical uncertainties of our predictions.
To summarise, the effect with the greatest impact is the different jet veto procedure employed by ATLAS. This could be modelled better, either by making use of an effective p t,veto , or by performing the appropriate resummation. Both are beyond the scope of the present work.

Sensitivity studies
In this section, we compare the sensitivity of W W and ZZ production at HL-LHC ( √ s = 14 TeV, with 3 ab −1 of integrated luminosity) to the BSM operator considered in eq. (1.1). Here we consider only the decay ZZ → e + e − µ + µ − . First we present the best predictions that could be obtained with the theoretical tools considered here, for a given choice of observables for the two processes. For W W we choose M T 1 in eq. (4.5a), and our best prediction is NNLL for qq → W W and NLL for gg → W W . For ZZ we consider M ZZ , and our best prediction is NLO for qq → ZZ and LO for gg → ZZ. Note that the accuracy of the predictions for qq annihilation for both W W and ZZ production can be improved to include the most recent NNLO calculations of refs. [70,71]. For gluon fusion, full NLO corrections have yet to be calculated, although approximate results are available [72][73][74][75][76][77]. While the inclusion of NNLO corrections to ZZ is straightforward, and can be obtained by running the code MATRIX [78][79][80], the use of NNLO corrections to W W requires matching of fixed-order predictions to the NNLL resummation. Although this can be achieved by interfacing the NNLL resummation to MATRIX, it is technically more involved than the simple procedure described in section 3. Therefore, we leave matching to NNLO to future work. The differential distributions in M T 1 and M ZZ are shown in figure 11. We observe that, in the qq channel, the cross section dσ/dM T 1 with a jet-veto is comparable to the cross section dσ/dM ZZ where no jet veto is applied. We  Figure 11. Our best predictions for the differential distributions dσ/dM T 1 for W W production with the experimental cuts in table 1 (left) and dσ/dM ZZ for ZZ production with the cuts in ref. [81] (right) for qq and gg processes.
note that, even with a jet veto, the qq background is much larger in the W W case. Therefore, we naively expect W W to perform slightly worse than ZZ for exclusion of BSM effects.
To be more quantitative, we generate exclusion plots for a range of values of the parameters κ t and κ g entering the Lagrangian of eq. (4.2). To do this we ask ourselves how likely it is that predictions corresponding to different values of (κ t , κ g ) are compatible with data that agree with the SM. Quantitatively, given a value of (κ t , κ g ), we compute n i (κ t , κ g ), the expected number of events in bin i of the distribution in a suitable leptonic observable. Specifically, we choose M T 1 for W W production and M ZZ for ZZ. Given a set of data points {n i } i=1,...,N , and a given value of (κ t , κ g ), we define and from that we construct our test statistic where (κ t ,κ g ) are the values of (κ t , κ g ) that minimise χ 2 (κ t , κ g ). This test statistic is a good approximation to the usual log-likelihood ratio for counting experiments [82] in the limit of a large number of events, and in the assumption that there are no correlations between bins. Assuming n i (κ t , κ g ) is the expected number of events, in the denominator of eq. (5.1) we can approximate n i n i (κ t , κ g ). Therefore, ∆χ 2 (κ t , κ g ) is asymptotically distributed according to a chi-squared distribution with two degrees of freedom (see e.g. [83]), which we denote by f (∆χ 2 (κ t , κ g ) | κ t , κ g ).
We now consider data {n i } i=1,...,N generated in such a way that the expected number of events in each bin is the "central" SM prediction, corresponding to µ R = µ F = Q = M W W /2 for W W and µ R = µ F = M ZZ /2 for ZZ, which we denote with n i (1, 0). This constitutes our "background-only" hypothesis. We now set exclusion limits in the (κ t , κ g ) plane using the median significance [82,84], assuming those data, with which one reject the hypothesis corresponding to each value of (κ t , κ g ) (our "signal" hypothesis). More precisely, for each value of (κ t , κ g ), we construct the distribution in ∆χ 2 (κ t , κ g ) under the assumption of the background-only hypothesis, which we denote by f (∆χ 2 (κ t , κ g ) | 1, 0). We then compute the median of that distribution, which we denote with ∆χ 2 med (κ t , κ g ). The p-value for each (κ t , κ g ) is given by and we exclude at the 95% confidence level all (κ t , κ g ) such that p(κ t , κ g ) < 0.05. In practice, we have binned the variables M T 1 and M ZZ in such a way that, when computing ∆χ 2 med (κ t , κ g ), in the denominator of eq. (5.1) we can always approximate n i withn qq i , the number of events obtained using central scales and the qq subprocess only.
We first consider the case in which the expected number of events for the signal hypothesis corresponds ton i (κ t , κ g ). We have examined two cases, both corresponding to di-boson invariant masses above the Higgs mass, so as to ensure to have complementary information with respect to Higgs cross sections. In one case, we have considered only two bins, a low-mass bin (200 GeV < M T 1 , M ZZ < 400 GeV) and a high-mass bin containing the rest of the distributions. The low-mass bin is more sensitive to κ t , and the high-mass bin to κ g . The corresponding exclusion regions in the (κ t , κ g ) plane are bounded by the dashed contours in Fig. 12. We see that W W is complementary to ZZ for low values of κ t , whereas the sensitivity to κ g of ZZ is larger. This can be understood from figure 11. Note that, despite the fact that the W W cross-section is larger, the presence of the jet veto kills a good fraction of the gg signal, with the net effect that its cross-section decreases with increasing M T 1 . In the ZZ case, where there is no suppression due to a jet veto, the contact interaction driven by κ g is fully effective, and makes the gg signal flatter with respect to the qq background, thus giving a larger sensitivity to κ g . We gain sensitivity by considering a greater number of bins. For instance, we have considered 60 bins equally spaced from 200 GeV to 1400 GeV, and an extra bin containing the distribution with larger values of M T 1 or M ZZ . The corresponding exclusion contours are the solid lines in Fig. 12. For reference, we also plot the line κ t + κ g = 1, and three points corresponding to the SM, and the scenarios BSM 1 and BSM 2 considered in the previous section. We also draw bands corresponding to 95% confidence-level bounds on κ t + κ g and κ t obtained from ref. [85]. These give more stringent constraints than our observables, which have nevertheless complementary sensitivity, since the anal- ysis of ref. [85] probes regions of di-boson invariant masses that we do not consider here. Also, having full control of theoretical predictions for both the signal and the background, our procedure is suitable for optimisation of both the observables and the binning procedure, and is open to improvements of the theoretical predictions. The exclusion contours we have obtained so far do not take into account theoretical uncertainties. Including theoretical uncertainties, the true theory value n i (κ t , κ g ) will differ from its central predictionn i (κ t , κ g ) by some theoretical error δ i , taken to lie in some interval ∆ i . In every bin, n i (κ t , κ g ) will be the sum of a contribution n (qq) i arising from quark-antiquark annihilation, and a contribution n  , and considering the fact that these predictions correspond to completely uncorrelated processes, we take the theoretical uncertainty on n i (κ t , κ g ) to be given by (5.4) Therefore, the χ 2 corresponding to a given value of (κ t , κ g , δ ≡ {δ 1 , δ 2 , . . . }) is given by In order to estimate the impact of theoretical uncertainties on our sensitivity contours, we adopt the approach of ref. [83], and add to χ 2 exp a Gaussian "theory term", with a width ∆ i (κ t , κ g )/2, as follows: The test statistic corresponding to (κ t , κ g ) is then obtained by profiling with respect to δ, i.e. computing For χ 2 exp and χ 2 th as in (5.5) and (5.6) this gives In other words, for a Gaussian theory term our treatment is equivalent to the common prescription to combine theoretical and experimental errors in quadrature. 5 With our choice of bins, we can approximate ∆ i (κ t , κ g ) ∆ (qq) i . Before presenting sensitivity contours including theory uncertainties, it is worth comparing the impact of statistical and theoretical uncertainties. In the case of W W production, theory uncertainties differ according to whether we use the efficiency method described in appendix A.3, or we just perform 9-point scale variations in the resummed cross section. In the former case, as can be seen from Fig. 7, relative theory uncertainties are of order 40%, whereas in the latter they are of order 10%, with a mild dependence on M W W . In both cases then ∆ (qq) i roughly scales like n i . Therefore, by looking at the denominator of eq. (5.8), we see that in the bins with larger n i , theory uncertainties will dominate over statistical uncertainties (∼ √ n i ), and hence these bins have very little power to constrain (κ t , κ g ). In the case of ZZ, theory uncertainties are smaller, around 5%, so all bins retain their constraining power. This is illustrated in Fig. 13. All contours have been obtained with 61 bins, as explained above. The outer contour (dotted) corresponds to W W production with theory uncertainties estimated with the JVE method. As explained in sec. 4, the method is 5 In fact, (5.6) itself can similarly be obtained as follows: (i) introduce separate δ  probably overly conservative, and the corresponding contour cannot compete with the constraints from ZZ production. Note in particular that large theoretical uncertainties affect mostly the bins with lowest values of M T 1 , which are the most sensitive to κ t . This explains why the JVE contour is so wide compared to the others. The solid contours correspond to uncertainties obtained with the appropriate scale variations, both for W W and for ZZ. Based on previous works on Higgs production with a jet-veto [35,41,51], we believe that scale variations for W W give a realistic estimate of the best theoretical uncertainties that could be obtained with a matching to NNLO with the JVE method. We see that, taking into account theory uncertainties at the currently achievable accuracy, W W does not have complementary constraining power with respect to ZZ. However, the dashed curves, corresponding to all predictions fixed at their central value without theory uncertainties, show that W W might compete with ZZ. We have therefore determined the necessary accuracy on W W production such that one obtains a comparable sensitivity with ZZ. First, we have observed that, in the case of ZZ, adding the NNLO contribution to qq does not improve the overall theory accuracy, due to missing higher orders in the gg channel.
So we assume that the uncertainties on ZZ production will remain the NLO ones, i.e. around 5%. The solid contour for W W in Fig. 14  theoretical uncertainty of 3% in every bin, which is approximately the one you need for W W to be competitive with current ZZ predictions. Based again on previous work on Higgs production [51], such an uncertainty could be reached by matching NNLL resummation to a future NNLO calculation for W W plus one jet, and maybe even further decreased after an N 3 LL resummation. We note that improving ZZ predictions hardly offers any stronger constraint. However, improved predictions for the gg channel, both for W W and ZZ, might move the central prediction, and may open up further space for constraints. We conclude this section with a comment on the actual implementation of the calculation of χ 2 (κ t , κ g ). If we consider the numerator of χ 2 (κ t , κ g ) in eqs. (5.1) and (5.8), we see that it involves n (gg) i (κ t , κ g ). This quantity is a second-order polynomial in κ t and κ g , arising from the square of the matrix element the remaining contributions, giving rise to the so-called "continuum" background. The fact that we have full control over M (gg) allows us to compute the coefficient of each power of κ t and κ g separately, and once and for all. This is crucial for an accurate calculation of χ 2 (κ t , κ g ), because a naive implementation of this quantity might involve cancellations between large numbers, whose control requires Monte Carlo samples with large statistics.

Conclusions
We have studied the impact of a veto on additional jets on setting limits on the coupling of a dimension-6 operator affecting W W production. In the presence of such a veto, large logarithms of the ratio of the maximum allowed jet transverse momentum p t,veto and the invariant mass of the W W pair M W W have to be resummed at all orders in QCD. These logarithmically enhanced contributions give rise to the so-called Sudakov suppression of cross sections with respect to naive Born-level predictions. The dimension-6 operator we considered affects W W production via gluon fusion, but does not affect W W production via quark-antiquark annihilation, which stays unchanged with respect to the SM. At Born level, the effect of this operator amounts to a growth of the cross section at large values of M W W . Unfortunately, the suppression due to the jet-veto gets larger with increasing M W W . Also, such suppression affects gluon fusion more than quark-antiquark annihilation due to the fact that gluons radiate roughly twice as much as quarks, so vetoing radiation off gluons cuts a larger portion of cross sections. Therefore, enhancement due to a contact interaction and Sudakov suppression are in competition.
To investigate quantitatively the impact of a jet veto on W W production, we have devised a new method to interface the most accurate resummed predictions for the gg and qq channels to fixed-order QCD event generators. This procedure provides events that are fully differential in the decay products of the W W pair, so that suitable acceptance cuts can be applied. The method involves minimal modifications of the ingredients already present in fixed-order event generators, and can be applied to the production of any colour singlet. In particular, we have implemented the procedure in the fixed-order program MCFM, which resulted in the code we called MCFM-RE, a Resummation Edition of MCFM.
Our program MCFM-RE has been used to produce differential cross sections for W W production with a simplified version of the ATLAS acceptance cuts, both in the SM, and including BSM effects induced by the aforementioned dimension-6 contact interaction. The main message is that, with the value of p t,veto used in current analyses, Sudakov suppression effects dominate over the enhancement produced by a contact interaction, so that deviations from the Standard Model are in general quite small for reasonable values of the strength of the contact interaction.
We have compared our results with those obtained from a number of partonshower event generators, and we have found very good agreement. We have used parton-shower event generators to estimate effects that cannot be not be taken into account by our analytical calculation, and found that they have a small impact, well within our theory uncertainties. We emphasise that our predictions have the computational cost of a Born-level event generator, and provide full analytical control of theoretical uncertainties. Our predictions are also in agreement, within uncertainties, with those obtained by interfacing a SCET calculation with the same formal accuracy with aMC@NLO.
We produced projections for the sensitivity to the considered BSM effects for HL-LHC, and compared with what could be obtained using ZZ production, which is not affected by the presence of a jet veto. We have found that W W has complementary sensitivity, provided it is possible to reduce theory uncertainties below 3%. This could be achieved by both matching current resummed predictions with a future calculation of W W plus one jet at NNLO, and improving the resummation to achieve N 3 LL accuracy. We hope this work encourages further theoretical work in both directions. We remark that the main advantage of using MCFM-RE for such studies compared to parton-shower event generators is that we have access to amplitudes, so we can compute separately all terms contributing to square matrix elements, in particular interference terms which can be computed separately with an arbitrary numerical accuracy.
We have found that, with the current acceptance cuts, the observables we have considered are not yet competitive with Higgs total cross sections, although they do provide additional information. However, our code does provide an accurate and fast tool to explore different choices of cuts and observables, so could be used for further studies in this direction. Also, it makes it possible to implement other models of new physics affecting the production of a colour singlet.
Last, our code is the only implementation of the jet-veto resummation of ref. [51] that is fully exclusive in the decay products of a colour singlet, so it can be used for precision determinations of Standard Model parameters, notably those characterising the Higgs boson.

A Collection of relevant formulae
In this appendix we report the explicit expressions that we have implemented in MCFM to achieve NLL and NNLL resummation of the cross section for the production of a colour singlet with a jet-veto. This discussion is of a technical nature, and we assume that the reader is familiar with the details of the jet-veto resummations performed in refs. [35,41,51].
In general, we consider the production of a colour singlet of invariant mass M , for instance a Higgs, a Z boson, or a pair of W bosons. At Born-level, this proceeds via either qq annihilation or gluon fusion. We then compute the cross section dσ i.s. /(dM 2 dΦ n ), with i.s. = qq, gg, fully differential in the phase space of the decay products of the colour singlet. Given their momenta q 1 , q 2 , . . . , q n , and incoming momenta p 1 and p 2 , the phase space dΦ n is defined as with E i and q i the energy and three-momentum of particle q i . Any prediction for dσ i.s. /(dM 2 dΦ n ) depends on the renormalisation scale µ R at which we evaluate the strong coupling α s , as well as the factorisation scale µ F at which we evaluate the PDFs. Both scales are typically set at values of order M . Furthermore, in the presence of a jet-veto, dσ (i.s.) /(dM 2 dΦ n ) is affected by large logarithms L ≡ ln(M/p t,veto ), with p t,veto the maximum allowed transverse momentum of the observed jets. When resumming such logarithms at all orders, our predictions become functions ofL, defined as The quantityL is such that for large p t,veto ,L → 0, which implements the fact that, in this regime, there are no large logarithms to be resummed. Also, at small p t,veto ,L ln(Q/p t,veto ), so in fact we resum logarithms of the ratio of p t,veto and the so-called resummation scale Q. The three scales µ R , µ F , Q are handles that we will use to estimate theoretical uncertainties, as explained in app. A.3. The power p determines how fast the resummation switches off at large p t,veto . We choose p = 5, as in refs. [35,41,51].

A.1 NLL resummation
At NLL accuracy, the distribution dσ i.s. /(dM 2 dΦ n ) is given by Explicit expressions for the functions g 1 , g 2 can be found in the supplemental material of ref. [35]. The NLL "luminosity" L (A.4) In the above expression, M (i.s.) ij is the Born-level amplitude for the production of the colour singlet via annihilation of the two partons i and j, and f i,j is the density of parton i, j in the proton.
Given any Born-level event generator, the recipe to implement the NLL resummation of eq. (A.3) is straightforward: 1. change the factorisation scale µ F provided by the generator to µ F e −L ; 2. multiply the weight of every event by a factor exp L g 1 (α sL ) + g 2 (α sL ) .
Note that, if p t,veto is fixed, and we do not integrate over different values of M 2 , both operations can be performed without touching the Born-level generator code. In fact, many programs allow a change in the factorisation scale by a constant factor. Also, the rescaling of the weight can be performed by the analysis routines that produce histograms for physical distributions. In our implementation, since we do want to integrate over M 2 , we have implemented the change in factorisation scale inside the MCFM code.
Another advantage we have in using MCFM is that it gives us access to the matrix elements in a form that is human readable. This is particularly useful in case one wishes to separate contributions from different parts of the matrix element, for instance a possible BSM contribution from that of the SM background. We consider here the case of W W production via gluon fusion, but the argument applies to other processes as well. There, the Born-level matrix element has the form M (gg) = M where we have used the notation with i = t, g, c, and with ij = tg, tc, gc. Using these luminosities we can interpret L (0) gg as a polynomial in the various κ i , and compute each coefficient separately. All one has to do then is to reweight each phase-space point using the Sudakov exponent exp L g 1 (α sL ) + g 2 (α sL ) . In doing so, we have used the fact that the Sudakov exponent depends only on the colour and kinematics of the incoming partons, and therefore is the same for every single contribution to the luminosity.

A.2 NNLL resummation
At NNLL accuracy, the cross section dσ i.s. /(dM 2 dΦ n ) with a jet veto is given by where the function g 3 can be found in ref. [35]. The functions F clust (R), F correl (R) depend on the jet radius R. Their expressions can be found in ref. [41]. As for the NLL resummation, α s = α s (µ R ). The remaining new ingredient for NNLL resummation is the luminosity L with H (1) the finite part of one-loop virtual corrections to the process in question, e.g. W W production through qq annihilation. The coefficients C ij depend on whether incoming partons i and j are quarks/antiquarks (q) or gluons (g), and are given by: As explained in the previous section, the NLL luminosity L i.s. of relative order α s . Therefore, its implementation requires at least a NLO generator. Any NLO event generator includes the calculation of virtual corrections, as well as integrated counterterms. This contribution, which we denote by dσ (1) i.s.,v+ct /(dΦ n dM 2 ), has the same form as the luminosity L i.s. and C (1) ij (z). Its expression in general depends on the way each process is implemented in the NLO event generator. For instance, the implementation of W W production in the NLO program MCFM follows from the general coding of the production of a colour singlet, whose details can be found in ref. [86]. Schematically, After direct inspection of the MCFM code, we realised that the term H MCFM,i.s. does not contain just the finite part of the virtual corrections H (1) , but also the terms −(π 2 /12) δ(1 − z) in the coefficients C i.s. through MCFM, we had to perform the following changes to the MCFM code: (A.14) 2. modify the integrated counterterms as follows 3. change the factorisation scale in all PDFs from µ F to µ F e −L .

A.3 Matching to fixed order and theoretical uncertainties
Our MCFM implementation includes the matching of resummed predictions with NLO calculations. In particular, we have implemented the relevant contributions to the two multiplicative matching schemes introduced in refs. [35,51]. At NLO, the total cross section σ NLO for the production of a colour singlet, satisfying a set of kinematical cuts for its decay products, is given by with σ (0) its Born-level contribution, and σ (1) a correction of relative order α s . Similarly, at NLO, the corresponding cross section with a jet-veto Σ NLO (p t,veto ) is given by Σ NLO (p t,veto ) = σ (0) + Σ (1) (p t,veto ) . which implies Σ NLO (p t,veto ) = σ NLO +Σ (1) (p t,veto ). We also denote by Σ N k LL (p t,veto ) the resummed jet-veto cross section at N k LL accuracy, again satisfying the chosen set of kinematical cuts for the decay products of the considered colour singlet. At this order, it has the following expansion in powers of α s : As in refs. [35,51], the matching is performed at the level of the jet-veto efficiency (p t,veto ), the fraction of events that survives the jet veto. This quantity is matched to exact NLO, as follows: N k LL (p t,veto ) σ 0 (1 + δL N k LL (p t,veto )) .
(A.21b) compile and run the MCFM code, in all its operation modes. If not, the interested reader should consult the MCFM manual [44].

B.1 Overview
MCFM-RE (an acronym for Resummation Edition) is a modification of MCFM-8.0 to include the resummation of jet-veto effects in colour-singlet processes up to NNLL+LL R accuracy. The modifications are modular, as most of the resummation effects are included through an interface to the code JetVHeto [38], suitably modified to become a library linkable to MCFM. Although a small number of modifications require us to directly change the MCFM code, these do not interfere with its usual modes of operation. The program is available at [43]. Included in the package are a README file and an example input card.
To run MCFM-RE, one must simply provide a suitably modified MCFM input card. We list here the new parameters we have added or changes made to existing parameters, described with the same conventions and terminology as the MCFM manual.
• file version number. This should match the version number that is printed when mcfm is executed.
{blank line} [Flags to specify the mode in which MCFM is run] • part ll. Jet-veto resummation at LL accuracy, i.e. each event produced by MCFM is reweighted with exp[Lg 1 (α sL )].
nnll. Jet-veto resummation at NNLL accuracy, with or without the inclusion of small jet radius resummation (LL R ), see eq. (A.9). -nllexp1. Expansion of the NLL resummation at order α s (for matching).
-nnllexp1. Expansion of the NNLL resummation at order α s (for matching).
ptj. The default mode of the resummation, resum logarithms of the jetveto.
-ptj+small-r. Available for NNLL resummations only. Include the effect of resumming the jet radius at leading logarithmic accuracy.
• Qscale. This parameter may be used to adjust the value of the resummation scale Q introduced in eq. (A.2). It behaves in the same way as the MCFM parameters scale and facscale do, i.e. if dynamicscale is .false., Q is set to Qscale, otheerwise Q = Qscale × µ 0 , with µ 0 the dynamic scale specified by the parameter dynamicscale.
• Rscale. This parameter may be used to adjust the value of the jet-radius resummation scale.
• ptjveto. The value of the jet-veto cut p t,veto in units of GeV.
{blank line} [Coupling rescaling in the kappa formalism] • kappa t. The parameter κ t of the Lagrangian in eq. (4.2), a.k.a the anomalous top Yukawa coupling.
• kappa g. The parameter κ g of the Lagrangian in eq. (4.2).
• interference only. Flag to control whether to compute just the interference terms, e.g. the coefficient of κ t κ g arising from squaring the amplitude in eq. (5.9). All other coefficients can be determined by setting a single κ i , i = t, g, b to zero.
Normally, MCFM identifies whether a process is qq-or gg-initiated, and running MCFM-RE in resummation mode does not lead to any problems. However, in cases like process 61, in fact W W production, MCFM includes in the NLO correction to a qq-initiated process formally higher-order gg-initiated contribution. As a consequence, not specifying the colour of the initial state leads to an ambiguity that is impossible to resolve. To avoid such problems, we have decided that, when running MCFM-RE in any resummation mode for ambiguous processes, the user must impose that a process is either qq-or gg-initiated, by making use of the MCFM flags omitgg and ggonly. Failure of doing so will result in MCFM-RE stopping and returning an error message.

B.2 Details of MCFM implementation
We modify MCFM version 8.0 to include the resummation of jet-veto effects. To this end there are two pieces that we must include, the computation of the luminosities L i.s. , L i.s. , and the Sudakov form factor combined with the functions F clust , F correl . The computation of the luminosities requires structural changes to MCFM whereas we are able to include the Sudakov form factor through an interface in src/User/usercode.f90.

end function end interface
The user should not normally make changes to this function. The reweighting is applied to all histograms, including the default MCFM ones, as wt and wt2 are intent(inout), so our reweighting is applied globally. The cost of doing the reweighting here is that the cross section returned by the main MCFM program is wrong, or rather it includes only the contribution of the luminosities and not the Sudakov exponent. To that end we include the extra histogram xsec, a single-bin histogram to record the correct total cross section for runs with the jet-veto.
To include the luminosities we have to modify the factorisation scales of the PDFs. Instead of adding lots of switches to the default MCFM integration routines, we create our own special routines resmNLL.f (based on lowint.f) and resmNNLL.f (based on virtint.f), which we include in the src/Procdep directory along with the other default integration routines. The changes made in resmNLL.f are modest with respect to lowint.f, schematically function resmNLL(r,wgt) use rad_tools, only: Ltilde implicit none include 'types.f' real(dp):: resmNLLint ! resummation include 'jetvheto.f' real(dp) :: facscaleLtilde real(dp) :: L_tilde_arr(1) At the beginning of each event we determineL, and the modified facscale which we call facscaleLtilde. We then use this scale in the computation of the PDFs. The simplicity here is that at NLL accuracy all we need to do is change the factorisation scale and reweight, so these changes are very modest.
To perform the same calculation at NNLL is much more involved, since there are three separate actions that must be performed to compute the luminosity. First, we need to cast the virtual matrix element into the correct form for the resummation. We do this with a utility function in the file src/Procdep/virtfin.f, which performs the replacement detailed in eq. (A.14). This is carried out by the subroutine subroutine virtfin(p,msq,msqv) real(dp) :: p(mxpart, 4) real(dp) :: msq(-nf:nf,-nf:nf), msqv(-nf:nf,-nf:nf) end subroutine virtfin where one must provide the array of momenta p(mxpart,4), the tree level matrix element squared msq(-nf:nf) and the matrix element of the virtual corrections msqv(-nf:nf) (using the conventions of MCFM).
The second contribution to the luminosities comes from the convolution of the coefficient functions. To include this coefficient function we modify the integrated dipole functions located inside src/Need/dipoles.f, adding switches to choose between the different types of "dipoles" that we have added as well as the default MCFM subtraction dipole.