QCD effects in electroweak W Zjj production at current and future hadron colliders

We present an update of an existing implementation of W Zjj production via vector-boson scattering in the framework of the POWHEG BOX program. In particular, previously unavailable semi-leptonic and fully hadronic decay modes of the intermediate vector bosons are provided, and operators of dimension six in an effective-field theory approach to account for physics beyond the Standard Model in the electroweak sector are included. For selected applications phenomenological results are provided to illustrate the capabilities of the new program. The impact of the considered dimension-six operators on experimentally accessible distributions is found to be small for current LHC energies, but enhanced in the kinematic reach of a potential future hadron collider with an energy of 100 TeV. The relevance of fully accounting for spin correlations and off-shell effects in the decay system is explored by a comparison with results obtained with the MadSpin tool that are based on an approximate treatment of the leptonic final state resulting from vector boson scattering processes. For selected semi-leptonic and hadronic decay modes we demonstrate the sensitivity of realistic signal selection procedures on QCD corrections and parton-shower effects.


Introduction
Vector boson scattering (VBS) processes are a particularly appealing class of reactions for exploring the electroweak (EW) sector of the Standard Model (SM) and possible extensions thereof.In the context of the SM, cross sections for the scattering of the longitudinal modes of an EW gauge boson are unitarised by Higgs-boson exchange contributions.The underlying cancellation mechanism is sensitive to both the Higgs-and the gauge-boson sector with deviations from the SM in particle content or properties immediately affecting the relevant cross sections.This makes VBS one of the most promising classes of processes for the discovery of physics beyond the SM in the EW sector.
In hadronic collisions the scattering of EW gauge bosons can be accessed in VBS processes which involve the scattering of hadronic constituents by EW gauge boson exchange.The experimental signature of such a reaction includes the gauge bosons' decay products and two distinctive jets resulting from the scattered partons.Because of the colour-singlet nature of the t-channel gauge-boson exchange these so-called tagging jets tend to be located in the far forward and backward regions of the detector with a large separation in rapidity.This feature of VBS reactions helps to identify the signal in the presence of QCD background processes with large production rates but rather different signatures.
VBS processes have received a lot of attention from the particle physics community.Tree-level predictions are available from the dedicated parton-level generator Phantom [1], and the multi-purpose tool Whizard [2].Next-to-leading-order (NLO) QCD corrections to VBS processes with massive gauge bosons have first been considered in refs.[3][4][5][6][7] and been implemented in the framework of the VBFNLO parton-level Monte-Carlo program [8].An alternative option is provided by the multi-purpose program MadGraph5 aMC@NLO [9].The matching of the NLO-QCD calculations with parton-shower (PS) programs according to the POWHEG formalism [10,11] has been considered in refs.[12][13][14][15] and made publicly available in the POWHEG BOX [16] repository.An independent implementation of VBS processes at NLO-QCD accuracy matched with PS in the context of the HERWIG7 Monte Carlo program [17] was presented in ref. [18].More recently, in addition to QCD corrections the NLO electroweak corrections have been considered in refs.[19][20][21][22][23][24][25][26].For a review of existing work at LO and NLO-QCD accuracy on the representative VBS W + W + jj channel and a comparison of the individual programs' features we refer the interested reader to ref. [27].A LO review of event generators for the VBS W Zjj channel is available from [28].
At the CERN Large Hadron Collider (LHC), the ATLAS collaboration reported the observation of EW W ± Zjj production at a center-of-mass energy of √ s = 13 TeV in final states with three identified leptons of electron or muon type in 2018 [29].The CMS collaboration observed the production of W ± Z pairs with leptonic decays in association with two jets at √ s = 13 TeV in 2020, and also provided constraints on anomalous quartic vector boson interactions [30].In refs.[31,32], semi-leptonic decay modes were considered.
In all of these experimental publications, however, the signal simulation was severely limited in accuracy: Refs.[29,30,32] resorted to a leading-order approximation of the signal process.In ref. [31] NLO-QCD corrections to the on-shell W Zjj production process were taken into account, but a factorised ansatz was used to simulate the subsequent decay of the gauge-boson system.This approximation works reasonably well in the resonance region, but fails to provide an accurate description in other regions of phase space, as will be illustrated below.
In this article we specifically consider the VBS-induced W ± Zjj process.Building on an existing implementation of the VBS W Zjj process [15] in the POWHEG BOX [16] that, however, was limited to fully leptonic decays of the gauge bosons, we provide an updated version of the program accounting for leptonic, semi-leptonic, and fully hadronic decay modes of the EW gauge bosons.For each mode, NLO-QCD corrections, off-shell effects, and spin correlations in the decay system are taken into account.Moreover, we provide an option to consider generic extensions of the SM.To illustrate the capabilities of the updated program, phenomenological results are presented for selected scenarios.
In detail, our work builds on previous developments for VBS-induced νe e − µ − µ + jj and ν e e + µ − µ + jj production in the context of the SM.In ref. [15] the NLO-QCD corrections for these reactions as calculated in [5] have been implemented in the POWHEG BOX.Here, we go beyond this existing implementation by adding semi-leptonic and hadronic decays of the gauge bosons, and considering an extension of the SM using the effective field theory (EFT) approach of ref. [33] that has already been used in the related ZZjj VBS process [14].Effects of dimension-six EFT operators in VBS have also been studied in refs.[34,35].In addition, we investigate the relevance of off-shell effects and spin correlations in the decays by an explicit comparison to approximate treatments of these effects.We would like to point out the existence of a complementary calculation [21] for EW W + Zjj calculation in the fully leptonic decay mode which focuses on perturbative corrections and provides results including the full NLO QCD and EW corrections.This calculation has, however, not been matched to a parton shower.
The article is structured as follows: In sec. 2 we describe the features of the updated implementation of EW W Zjj production in hadronic collisions in the POWHEG BOX.Using this program, we present some representative numerical results of EW W Zjj production at the LHC and a potential future circular hadron collider (FCC) operating at an energy of 100 TeV.We conclude in sec.4.

Details of the POWHEG BOX implementation
In order to provide a Monte-Carlo program for the simulation of EW W Zjj production in hadronic collisions at NLO+PS accuracy with the option for various leptonic, semileptonic, and hadronic final states in the context of the SM and a generic model supporting anomalous interactions in the gauge boson sector we provide appropriate extensions of the public POWHEG BOX implementation of Ref. [15].
We recall that this existing program resorts to tree-level and NLO-QCD matrix elements for the purely EW processes pp → ν e e + µ − µ + jj and pp → νe e − µ − µ + jj adapted from the VBFNLO parton-level Monte-Carlo generator [8].Even though in the following referred to as "EW W Zjj production" for the sake of brevity, it is implicitly understood that all resonant and non-resonant diagrams giving rise to a ν e e + µ − µ + jj or νe e − µ − µ + jj final state system are taken into account within the so-called VBS approximation.The VBS approximation only retains contributions from t-channel and u-channel diagrams, but not their interference, and disregards s-channel contributions.When selection cuts typical for an experimental VBS analysis are applied this approximation has been found [28] to reproduce the full result for the VBS cross section very well.For instance, for a representative setup at LO the authors of [28] report an agreement at the level of 0.6% between calculations based on full matrix elements and predictions of VBFNLO within the VBS approximation.We note, however, that the validity of this approximation deteriorates when more inclusive selection cuts are applied, see e.g.Ref. [27] for a comprehensive study of the VBS approximation.For this reason we only consider analysis setups with tight VBS cuts in sec.3. The Cabibbo-Kobayashi-Maskawa matrix is assumed to be diagonal, and contributions from external top or bottom quarks are not taken into account.For the updated POWHEG BOX implementation of the EW W Zjj production process presented here we resort to the same approximations.
As long as fully leptonic decays of the gauge bosons are considered, within the mentioned approximations the structure of the NLO-QCD corrections does not change if new interactions in the EW gauge-boson sector are taken into account.The original SM amplitudes for ν e e + µ − µ + jj and νe e − µ − µ + jj production are structured in a modular way with leptonic tensors for those building blocks of the relevant Feynman diagrams that only contain colour-neutral particles, and hadronic currents accounting for the scattering quarks and, in the real-emission contributions, gluons.The NLO-QCD corrections only affect the hadronic currents.Extensions of the SM in the EW sector thus merely require an appropriate replacement of the leptonic tensors.
We provide such an extension in the framework of the effective field theory (EFT) approach of Ref. [33] accounting for anomalous interactions in the EW gauge boson sector by an extension of the SM Lagrangian with operators of higher mass dimension, Here, d denotes the mass dimension of the operators O are the expansion coefficients.The sum over d includes contributions of all higher-dimensional operators starting from d = 6.The summation index i runs over all non-vanishing operators of a given mass dimension.The parameter Λ denotes the energy scale up to which the EFT is supposed to be valid.We assume Λ to be much larger than the EW scale and thus restrict ourselves to the contribution of operators up to dimension six.
Following Ref. [33] we consider three independent operators that conserve charge (C) and parity (P ), and two C and/or P violating operators, These operators are constructed from the Higgs doublet field Φ and the electroweak field strength tensors W a µν (a = 1, 2, 3) and B µν , with the U (1) and SU (2) gauge fields B µ and W a µ and their respective couplings g ′ and g.The σ a denote the Pauli matrices.The covariant derivative D µ is given by and the modified field strength tensors Ŵµν and Bµν are defined by while the modified dual field strength tensor is given by To simplify our notation we denote the coefficients of the operators of eqs.( 2)-( 6) that appear in the EFT expansion of eq. ( 1) up to dimension six as In the following, instead of a numbered index i we use the label of the corresponding operator to identify each operator coefficient.For instance, C W W W is the properly normalized coefficient of the O W W W operator.
In the actual calculation of scattering cross sections care has to be taken to ensure a consistent EFT expansion up to the desired order in 1/Λ 2 .Schematically, the operators of different mass dimension enter in the relevant matrix elements squared as Thus, if one truncates the EFT expansion at order 1/Λ 2 , in addition to the pure SM contribution one should only keep the interference term 2Re(M SM M ⋆ dim6 ), and disregard the quadratic term |M dim6 | 2 , which is part of the O(1/Λ 4 ) result.However, in the past in many applications this quadratic term was considered as part of the "dimension six" results.Below we therefore consider both options and refer to them as SM+lin and SM+quad.In addition to the ν e e + µ − µ + jj and νe e − µ − µ + jj final states provided in ref. [15], in this work we also implemented EW production of a νe e − ν µ νµ , a q q′ µ − µ + , a νe e − Q Q, a ν e e + Q Q, or a q q′ Q Q system, respectively, in association with two tagging jets (here q and Q refer to massless quarks of different types).Representative diagrams for the partonic channel uc → us du µ + µ − are shown in fig. 1.For each channel, not only resonant diagrams related to the leptonic or hadronic decay of a W or Z boson are taken into account, but also non-resonant diagrams resulting in the same final state.For simplicity we will refer to the previously listed processes (including off-resonant contributions) as fully leptonic, leptonic-invisible, semi-leptonic, and fully hadronic decay modes.In the case of semileptonic and fully hadronic decay modes we do not take QCD corrections to the decays into account, and we neglect QCD corrections connecting the W Zjj production with the decay part of the considered final state.The latter type of corrections are expected to be negligible.Corrections to the hadronic decays of the Z and W bosons are accounted for by the multi-purpose Monte-Carlo programs matched to our NLO-QCD calculation.
For each final state, at Born level singularities in the production cross section arise from diagrams with the t-channel exchange of a photon of very low virtuality Q 2 .Such contributions are entirely negligible after selection cuts on the tagging jets are applied and can thus be removed already at generation level by a cut, To further improve the numerical efficiency of the Monte-Carlo integration additionally a Born-suppression factor of the form can be employed.Here, the p T,i denote the transverse momenta of the final-state partons of the underlying Born configuration Φ, and Λ Φ is a technical parameter, by default set to 10 GeV.We remind the reader that contributions with a pair of same-type charged fermions (f ) in the final state cannot only stem from decays of the W Z system, but also from diagrams where a photon decays into a fermion pair.An additional type of singularity at Born level arises from diagrams where such a photon exhibits very low virtuality.For analyses that require a fermion pair with an invariant mass close to the mass of the Z boson, such singularities can easily be removed already at generation level by an invariant mass cut on the respective lepton pair.The requirement suffices to remove any potentially problematic contributions from photons of very low virtuality at generation level.

Phenomenological results
In the following we will provide some phenomenological results generated by our implementation in version 2 of the POWHEG BOX.For all of these results we use the following general settings: We set the Fermi constant to G µ = 1.1663787 • 10 −5 GeV −2 .For the masses and widths of the EW bosons we use: m H = 125.25 GeV, m W = 80.377 GeV, m Z = 91.1876GeV, Γ H = 0.0032 GeV, Γ W = 2.085 GeV, and Γ Z = 2.4952 GeV.The EW coupling, α em , is calculated therefrom via tree-level EW relations.
In addition to the POWHEG BOX we also use the tool-chain MadGraph5 aMC@NLO [9,36] including the program MadSpin [37] for comparison to the POWHEG BOX results.To simulate the parton shower we use PYTHIA8, version 8.245, with the Monash2013 tune [38].Hadronisation, MPI, and QED emissions are turned off in order to isolate the effect of the shower and matching.We note that in realistic simulations non-perturbative effects have to be considered [39].For reconstructing jets we resort to FastJet [40], version 3.3.4.For all results presented in this section jets are clustered via the anti-k T algorithm [41] with a radius-parameter of R = 0.4, and the NNPDF31 nlo as 0118 (ID 303400) set [42] of parton distribution functions (PDFs) is used as provided by version 6.3.0 of the LHAPDF library [43].

EFT results
In this subsection we explore the impact of SM extensions in the EFT framework introduced in sec.2. We individually set the coefficient of each EFT operator defined in eqs.( 2)-( 6) to the largest value compatible with the experimental limits of ref. [44] while setting the coefficients of all other EFT operators to zero.As it turned out that the impact of the O W W W operator is most pronounced, below we only display results obtained for non-vanishing values of the C W W W coefficient.For instance, results obtained with non-vanishing values of the C W operator coefficient are basically identical to the SM results and will thus not be further discussed here.
For the unitarisation of our EFT predictions we proceeded along the lines of ref. [45] where unitarity violations are avoided by using appropriate cuts on the invariant mass of the vector bosons produced in a VBS reaction.We calculated the limits beyond which unitarity violations are to be expected for the setups considered in this work using the tool calc-formfactor [46,47] that is available within the VBFNLO package [48,49].We found that unitarity violations would occur only beyond scales relevant for the results shown below.
Throughout this subsection we use a renormalisation scale, µ R = ξ R µ 0 , and factorisation scale, µ F = ξ F µ 0 , that is expected to optimally account for the region of high transverse momenta where the effects of the dimension-six operators are expected to have the largest impact (c.f.ref. [14]).The scale µ 0 is given by where the sum includes the transverse momenta p T,f of all n part final-state partons of a considered Born-type or real-emission configuration.In addition we define where p T,Z and p T,W are the transverse momenta of the muon pair and the positronneutrino pair of the fixed-order configuration, respectively.The factors ξ R and ξ F are varied between 0.5 and 2 with 7-point variation in our NLO+PS simulations and set to one for the fixed-order calculations.
For our representative numerical studies we use settings inspired by the ATLAS analysis described in ref. [50].In particular, we construct jets using the anti-k T algorithm with R = 0.4.In the following we will denote the jets with index j 1 for the hardest jet and j 2 , j 3 , ... for further jets ordered by their transverse momentum.The two hardest jets are identified as tagging jets and are required to have a transverse momentum, rapidity and invariant mass of Moreover, we only keep events where the tagging jets lie in opposite hemispheres, and have a large rapidity separation of To identify additional non-tagging jets we require them to be located in a rapidity range of No other cuts are applied on non-tagging jets unless specifically stated otherwise.
In our fixed-order results the final state contains two muons, one positron and one neutrino.Within the PYTHIA setup we consider, i.e. in the absence of QED radiation in the parton shower and without hadron decays, the events we simulate at NLO+PS level do not exhibit any additional leptons or neutrinos.For the charged leptons ℓ we demand p T,ℓ > 15 GeV, We do not apply any cuts on the neutrino.Furthermore, we require a clear separation of the tagging jets and the charged leptons, i.e. we require a separation in the rapidityazimuthal angle plane of For the muons which are stemming from the Z decay we additionally demand as well as that their reconstructed invariant mass, m Z inv , lies in a window around the physical Z-boson mass of 66 GeV < m Z inv < 116 GeV .
Finally, we also require all charged leptons to lie in the rapidity gap between the two tagging jets min(y tag j 1 , y tag j 2 ) < y ℓ < max(y tag j 1 , y tag j 2 ).
Let us now discuss results for the LHC with a center-of-mass energy of √ s = 13 TeV.In fig. 2 and fig. 3 we compare SM results at NLO+PS accuracy for selected distributions of the tagging jets and the leptons with those obtained in the EFT framework of sec.2. In particular, we set the operator coefficient C W W W to the maximal negative value compatible with current experimental limits, i.e.C W W W = −4.2TeV −2 , and consider separately the case where only the linear term of the EFT expansion sketched in eq. ( 15) is taken into account (dubbed SM+lin), and the case where additionally the quadratic term is retained (referred to as SM+quad).
For the SM+lin implementation we find only small differences to the SM results for all considered distributions.These effects are best visible for the azimuthal angle separation of the two tagging jets, Φ jj .In particular for Φ jj ≲ π/2 the shape is slightly different from the SM case.Larger differences to the SM case are found for the SM+quad implementation.However, we would like to remind the reader that the limits chosen for C W W W have been derived for the SM+lin implementation and that the SM+quad version is shown only for the purpose of comparison.The SM+lin and SM+quad predictions can barely be distinguished from the respective SM results for the transverse momentum and the rapidity distributions of the hardest tagging jet and the transverse momentum of the Z boson.The latter is reconstructed from the momenta of the muon pair closest in invariant mass to m Z .
In fig. 4 and fig. 5 we display the same observables as in fig. 2 and fig. 3 for a potential future circular collider (FCC) with proton-proton collisions at a center-of-mass energy of √ s = 100 TeV.For the FCC discussion we use the following setup inspired by ref. [51].We use the same renormalisation and factorisation scale as in eq. ( 19), but apply stronger cuts on the tagging jets.More precisely, we require Additionally, the tagging jets have to fulfill For the charged leptons we require as well as a separation of the tagging jets and the charged leptons in the rapidity-azimuthal angle plane of We keep the requirements of eqs.( 27)- (29).Differences between the SM and the SM+quad implementation are clearly enhanced at this energy for all distributions.However, the azimuthal angle separation of the tagging jets exhibits sensitivity to the effects of the SM+lin implementation.

Leptonic decays
The objective of this subsection is to explore the relevance of simulating the full leptonic final state of VBS-induced νe e − µ − µ + jj production as opposed to approximations where a W + Z boson pair is produced on-shell and combined with a simulation for the decays of these bosons into the desired leptonic final state.Such an approximation is  27)- (33).used, for instance, in the search for anomalous EW production of vector boson pairs in association with two jets by the CMS collaboration [32] and in the search for EW diboson production in association with a high-mass dijet system in semi-leptonic final states by the ATLAS collaboration [31].Since QCD corrections do not directly affect the leptonic decays, we conduct this discussion at LO.To that end we compare our POWHEG BOX implementation for νe e − µ − µ + jj production via VBS with two alternative ones: The first is using MadGraph5 aMC@NLO where a W + Zjj final state is produced onshell.The decays of the two bosons are afterwards simulated via MadSpin (this simulation is denoted by MG5+MadSpin below).The other implementation is using the LO version of MadGraph5 aMC@NLO and includes off-shell contributions and spin-correlations in the   27)- (33).lepton system (this implementation will be denoted by MG5-full in the following).This comparison serves as a consistency check for the correct usage of the two tools.It also shows the equivalence of the POWHEG BOX and the MG5-full implementation after the application of VBS cuts despite more approximations being used in the matrix elements entering the POWHEG BOX than the MG5-full implementation.
For the comparison presented in this subsection we use a fixed factorisation and renormalisation scale of As before, we consider proton collisions at the LHC with √ s = 13 TeV.For the selection of signal events we proceed along similar lines as in the previous subsection.We impose the cuts of eqs.( 21)-( 29) defined above.
In fig.6 we show several distributions related to the decay system of the VBS-induced ν e e + µ − µ + jj production process.The invariant masses of the W and Z systems, m W inv and m Z inv , are reconstructed from the momenta of the ν e e + and µ − µ + pairs, respectively.The transverse mass of the W Z system is defined by where p T,W and p T,Z are the transverse momenta of the reconstructed W and Z systems, and the transverse energies ẼT,i (i = W, Z) are given by We find that, for each considered distribution, the results of the MG5-full and the POWHEG BOX implementations are in very good agreement within their respective statistical uncertainties.In contrast, the results of MG5+MadSpin deviate from the implementations that retain full control on off-shell contributions and spin correlations in the leptonic decay system.While near the W and Z resonances MG5+MadSpin yields satisfactory results, further away from the peaks of the invariant mass distributions, and in particular in the tails of the transverse mass distribution, the on-shell approximation does no longer accurately reproduce the full results.Deviations can reach almost an order of magnitude for m T (W Z) ≲ 150 GeV and m T (W Z) ≳ 350 GeV.This should be kept in mind for the simulation of VBS processes when off-shell regions are of interest.

Semi-leptonic and hadronic decays
While in ref. [15] an implementation for VBS-induced W Z production with fully leptonic decays was developed and made available in the POWHEG BOX program package, semileptonic and fully hadronic final states were not considered before in that framework.We have closed this gap and are now able to simulate all possible decay modes of the W Z system in the VBS mode.
As an example for a semi-leptonic decay mode in this subsection we provide phenomenological results for the VBS process where a W + boson decays hadronically into an u d pair, and the Z boson into a muon pair.All off-shell diagrams giving rise to the same final state are taken into account, see fig. 1.
The decay quarks of the W boson give rise to jets.For the selection of the VBS signal in the presence of background processes it is important to distinguish these decay jets from the tagging jets of the production process.The representative numerical analysis below is designed to take that into account.
Throughout this subsection we use the dynamical scale of eq. ( 19).For the selection of events in the semi-leptonic mode at the LHC with an energy of √ s = 13 TeV we impose the following cuts: Charged leptons are required to fulfill the basic requirements p T,ℓ > 20 GeV, For the hardest lepton we additionally request Furthermore, the invariant mass reconstructed from the muon pair has to fulfill 83 GeV < m Z inv < 99 GeV.(39) Jets and the leptons have to be separated by We require all jets to have p T,j > 20 GeV, The jets in an event are then further classified.The decay jets are identified as those two jets with the invariant mass being closest to the mass of the W boson.These two jets are required to lie within a window around the W mass of 64 GeV < m W inv < 106 GeV , and to fulfill the transverse-momentum requirements where j 1 and j 2 denote the hardest and second hardest decay jets.
After the selection of the decay jets, the two hardest remaining jets are identified as the tagging jets.For the tagging jets we require Fig. 7 and fig.8 display some representative distributions related to the tagging jets and the decay jets at LO, NLO, and NLO+PS accuracy.As expected, tagging jets and decay jets exhibit entirely different properties.Both, transverse momentum and rapidity distributions look very different for the two types of jets.From the rapidity distributions one can understand that the decay jets are preferentially located a central rapidities, while the tagging jets peak in the forward and backward regions, analogous to the leptonic decay mode.We note that generally the curves for the NLO+PS results lie below the fixed-order predictions, indicating smaller event rates.As can be deduced from the invariant mass distribution in fig.8 (b) the shift of momenta beyond LO results in a considerable change of shape in m W inv .However this quantity is used as a selection criterion (c.f.eq. ( 42)).With m W inv shifted to values further away from m W , at NLO+PS level fewer events pass the selection criterion of eq. ( 42) resulting in a smaller value of the associated cross section.
Let us now consider a representative fully hadronic decay mode with the W + boson decaying into a u d pair and the Z boson into a pair of s quarks.The large number of jets emerging from this mode requires a dedicated analysis allowing for a proper classification of tagging and decay jets.
We require all jets to have p T,j > 10 GeV, To identify the jets associated with the decay of the W and Z boson we proceed in the following way: In the first step, from all pairs of jets we choose the one with its invariant mass closest to m W .These two jets are then considered to correspond to the decay of the W boson.In a second step, we proceed analogously for the Z boson (replacing m W with m Z ).In a third step, we check if a jet is contained in both the jet-pair associated with the W boson and the one associated with the Z boson.If this is the case we assign it to the boson with mass closer to the corresponding jet-pair.In a fourth step, we repeat the first or second step for the boson that was not chosen in the third step with the remaining jets not assigned to the other boson.The jets having thus been associated with the W and Z decays have to fulfill the transverse-momentum requirements of If an event does not pass these cuts, it is discarded.After the selection of the decay jets, the two hardest remaining jets are identified as the tagging jets.As in the semi-leptonic case, for the tagging jets we additionally require In fig. 9 and fig. 10 we show several distributions for the fully hadronic decay mode.We observe results that are qualitatively similar to the semi-leptonic case discussed above.However, the "smearing effect" in the NLO+PS results is even larger than in the semi-leptonic case, because now window cuts for both the Z and the W system have to be passed for an event to be accepted.

Conclusions and outlook
In this article we have presented new features for the implementation of VBS-induced W Zjj production in the framework of the POWHEG BOX V2.We are providing semileptonic and fully hadronic decays of the intermediate vector bosons that were missing in the previously existing implementation, and account for physics beyond the SM by the inclusion of dimension-six operators of a generic EFT expansion in the EW sector.
To illustrate the capabilities of the updated implementation we considered some selected applications.We explored the sensitivity of typical VBS observables to dimensionsix operators in an EFT framework, and found that, when the expansion is performed consistently, predictions with contributions from EFT operators compatible with current    experimental limits barely deviate from the SM case at LHC energies.Larger effects are found at higher energies which could be achieved, for instance, at a future FCC.
Using the leptonic decay mode as an example, we investigated the relevance of including off-shell contributions and spin correlations in the simulation.By comparing our results to those obtained with the MadGraph5 aMC@NLO+MadSpin tool that combines a calculation of VBS-induced W Z production with a simulation of the gauge-boson decays we could show that the on-shell approximation is appropriate when all final-state leptons stem from the resonant decay of a gauge boson, but deviates from the full result in regions away from the resonance.This limitation of the approximation should be kept in mind for ensuring its application is restricted to its region of validity.
Finally, we considered semi-leptonic and fully hadronic decay modes.We found that in these cases QCD corrections and PS effects can lead to a reshuffling of momenta such that they pass different selection cuts than the corresponding LO configurations.This kinematic effect results in a reduction of cross section beyond the LO, which becomes particularly pronounced once the NLO result is matched with a PS.
The new features of the VBS W Zjj code have been made available via the POWHEG BOX V2 repository, see https://powhegbox.mib.infn.it/.

Figure 2 :
Figure 2: NLO+PS predictions for pp → ν e e + µ − µ + jj at the LHC with √ s = 13 TeV within the cuts of eqs.(21)-(29) for the SM+lin (green) and the SM+quad case (orange) with C W W W = −4.2TeV −2 , and within the SM (purple).The upper panels show the transverse momentum of the hardest tagging jet (a), the rapidity of the hardest tagging jet (b), the reconstructed transverse momentum of the Z boson (c), and the azimuthal angle separation of the tagging jets (d).The respective lower panels show the ratios of the SM+lin and SM+quad predictions to the pure SM results.

Figure 3 :
Figure 3: NLO+PS predictions for pp → ν e e + µ − µ + jj at the LHC with √ s = 13 TeV within the cuts of eqs.(21)-(29) for the SM+lin (green) and the SM+quad case (orange) with C W W W = −4.2TeV −2 , and within the SM (purple).The upper panels show the transverse momentum of the 3rd jet (a) and the rapidity of the 3rd jet for a p T cut of p T,j 3 > 10 GeV (b).The respective lower panels show the ratios of the SM+lin and SM+quad predictions to the pure SM results.

Figure 4 :
Figure 4: Same as in fig.2, but for the FCC with an energy of √ s = 100 TeV and with the cuts of eqs.(27)-(33).

Figure 5 :
Figure 5: Same as in fig.3, but for the FCC with an energy of √ s = 100 TeV and with the cuts of eqs.(27)-(33).

Figure 6 :
Figure 6: LO predictions for pp → ν e e + µ − µ + jj at the LHC with √ s = 13 TeV within the cuts of eqs.(21)-(29) as obtained with the POWHEG-BOX (purple), MG5+MadSpin (green), and with MG5-full (orange).The upper panels show the reconstructed invariant mass of the W boson (a), the reconstructed invariant mass of the Z boson (b), the transverse mass of the W Z system (c), and the reconstructed transverse momentum of the Z boson (d).The respective lower panels show the ratios of the MG5+MadSpin and MG5-full predictions to the POWHEG-BOX results.

Figure 7 :
Figure 7: Predictions for VBS-induced W + Z production in the semi-leptonic decay mode at the LHC with √ s = 13 TeV within the cuts of eqs.(37)-(44) at LO (purple), NLO (green) and NLO+PS (orange).The upper panels show the transverse momentum of the hardest tagging jet (a), the rapidity of the hardest tagging jet (b), the transverse momentum of the hardest decay jet (c), the rapidity of the hardest decay jet (d).The respective lower panels show the ratios of the LO and NLO+PS predictions to the NLO results.

Figure 8 :
Figure 8: Predictions for VBS-induced W + Z production in the semi-leptonic decay mode at the LHC with √ s = 13 TeV within the cuts of eqs.(37)-(44) at LO (purple), NLO (green) and NLO+PS (orange).The upper panels show the transverse momentum of the hardest lepton (a) and the reconstructed invariant mass of the W -boson system (b).The respective lower panels show the ratios of the LO and NLO+PS predictions to the NLO results.

Figure 9 :
Figure 9: Predictions for VBS-induced W + Z production in the fully hadronic decay mode at the LHC with √ s = 13 TeV within the cuts of eqs.(45)-(48) at LO (purple), NLO (green) and NLO+PS (orange).The upper panels show the transverse momentum of the hardest tagging jet (a), the rapidity of the hardest tagging jet (b), the transverse momentum of the hardest decay jet from the W decay (c), the transverse momentum of the hardest decay jet from the Z decay (d).The respective lower panels show the ratios of the LO and NLO+PS predictions to the NLO results.

Figure 10 :
Figure 10: Predictions for VBS-induced W + Z production in the fully hadronic decay mode at the LHC with √ s = 13 TeV within the cuts of eqs.(45)-(48) at LO (purple), NLO (green) and NLO+PS (orange).The upper panels show the reconstructed mass of the W boson (a), the reconstructed mass of the Z boson (b).The respective lower panels show the ratios of the LO and NLO+PS predictions to the NLO results.