QCD and electroweak corrections to WZ scattering at the LHC

We present the first computation of the full next-to-leading-order QCD and electroweak corrections to the WZ scattering process at the LHC. All off-shell, gauge-boson-decay, and interference effects are taken into account for the process $\mathrm{p} \mathrm{p} \to \mu^+\mu^-\mathrm{e}^+\nu_\mathrm{e} \mathrm{j} \mathrm{j} + X$ at the orders $\mathcal{O}{\left( \alpha_\mathrm{s} \alpha^6 \right)}$ and $\mathcal{O}{\left( \alpha^7 \right)}$. The electroweak corrections feature the typical Sudakov behaviour towards high energy and amount to $-16\%$ relative to the electroweak contribution to the integrated cross section. Moreover, the corrections induce significant shape distortions in differential distributions. The next-to-leading-order analysis of the quark- and gluon-induced channels is supplemented by a leading-order study of all possible contributions to the full $4\ell+2\mbox{jets}$ production cross section in a realistic fiducial phase-space volume.


Introduction
The accumulation of experimental data during Run II of the Large Hadron Collider (LHC) allows to measure some rare Standard Model (SM) processes for the first time. Vector-boson scattering (VBS) processes constitute a prime example of processes that have not been measured before Run II. While the scattering of like-sign W-boson pairs, the golden VBS channel, has been measured first [1][2][3][4][5], the WZ channel comes in second [6,7]. It features a lower cross section than W ± W ± scattering, but has only one neutrino in the final state, allowing thus for better reconstruction and a better study of its properties.
As experimental errors (both statistical and systematic) will shrink in the next few years, precise theoretical predictions should be carefully prepared. In particular, higher-order corrections of both QCD and electroweak (EW) type should be incorporated. The inclusion of next-to-leading-order (NLO) QCD corrections has become standard for LHC analyses, but not yet the inclusion of EW corrections, which are known to increase at high energies owing to Sudakov logarithms. For the class of VBS processes, EW corrections are expected to be particularly large [8]. This expectation was confirmed in the first complete NLO QCD+EW calculation presented in Ref. [9] for like-sign WW scattering where it turned out that the genuine EW corrections of order O α 7 are even the largest NLO contribution.
In this article, we present results for the first calculation of the full NLO QCD+EW corrections to the WZ scattering process at the LHC with the final state µ + µ − e + ν e jj. An analysis of the LO contributions to the WZjj production mode was presented in Ref. [10], where also different Monte Carlo programs were compared. In our NLO analysis, we include the whole set of contributing diagrams in the relevant orders, instead of only VBS configurations. The QCD and especially the EW corrections are rather involved, as the process features seven charged external particles. This is the first time that EW corrections are computed for a process involving so many charged particles. The leptonic final state with a single net charge gives rise to a larger number of partonic channels as compared to like-sign WW VBS, which complicates the calculation further. We also note that our calculation of QCD corrections is based on the full set of NLO diagrams including all interferences without approximation, i.e. we do not employ the so-called VBS approximation used in previous QCD calculations [11,12], which neglects colour exchange between the two incoming protons. While for the current experimental precision such approximations are most likely sufficient, in the future they might actually be inadequate, because they can fail at the level of 10% in differential distributions, as shown in Ref. [13] for like-sign WW scattering.
In addition to the contributions to the NLO cross section of orders O α s α 6 and O α 7 , we also provide predictions for all LO processes relevant for the µ + µ − e + ν e jj final state. These include the orders O α 6 (EW contribution), O α s α 5 (interference), and O α 2 s α 4 (QCD contribution). Contributions including photons in the initial state or external bottom quarks are discussed separately.
All these results are presented in the form of cross sections and differential distributions for realistic experimental cuts. Specifically, the event selection chosen is the so-called loose fiducial region presented by the CMS collaboration in Ref. [7]. It has the advantage to be simple enough to be implemented easily in a Monte Carlo program. Such an experimental effort is particularly welcome by theorists as it allows for a direct use of state-of-the-art theoretical predictions in experimental analyses.
Finally, we would like to mention that all results have been produced by two independent Monte Carlo programs, matrix elements providers, and loop libraries: One is the Monte Carlo program BONSAY with matrix elements from OpenLoops [14,15] and loop integrals evaluated with the DD mode of the Collier [16,17] library. The other Monte Carlo is MoCaNLO with matrix elements from Recola [16,17] and loop integrals evaluated with the COLI mode of the Collier library. The two independent calculations ensure a exhaustive validation of all results presented in this paper.
This article is organised as follows: In Section 2, the various contributions to the NLO cross section of the considered process are described. In Section 3 the details on the implementations of the computation as well as the checks performed to validate the results are presented. Section 4 is devoted to the description and the analysis of the results. Finally, Section 5 contains a summary of the article as well as concluding remarks.

Leading-order contributions
As for all quark-quark-initiated processes characterised by four leptons and two jets in the final state, two types of amplitudes occur for the quark-induced processes qq → µ + µ − e + ν e qq, where q generically stands for a quark or antiquark: These are diagrams of order O g 6 (some sample diagrams are shown in Figures 1a-1f) and diagrams of order O g 2 s g 4 (an example is depicted in Figure 1g), with g s and g generically denoting the strong and electroweak gauge couplings, respectively. Besides VBS the former diagrams involve also the production of three vector bosons as well as singly-resonant and non-resonant diagrams. Consequently, three different orders contribute to the LO cross section: O α 6 (EW contribution), O α s α 5 (interference), and O α 2 s α 4 (QCD contribution). At the order O α 6 , where all couplings are of EW origin, there are in addition contributions from γγ → µ + µ − e + ν e qq (Figure 1h provides an example). For the quark contributions, one can further distinguish the cases in which an external quark is a bottom quark or a light one. In our predictions, we show separately the O α 6 contributions with an external bottom quark, as the corresponding partonic channels can develop a top-quark resonance (see Figure 1f). The photon-induced contributions are also shown separately.
The order O α s α 5 contributions are obtained by interfering amplitudes of O g 6 and O g 2 s g 4 in the channels qq → µ + µ − e + ν e qq. Further contributions in this order result from squares of amplitudes of order O g s g 5 of the channels gγ → µ + µ − e + ν e qq and qγ → µ + µ − e + ν e qg (see Figures 1i and 1j for examples). Those contributions are shown together with the order O(α 6 ) photon-induced contributions.
Finally, the O α 2 s α 4 contributions result from channels with either four external quarks or two external quarks and two gluons (see Figures 1k and 1l for examples). The contribution with two gluons in the initial state is particularly large due to the large gluonic parton-distribution function (PDF) at the LHC. The order O(α 2 s α 4 ) contributions

Virtual corrections
We compute the NLO corrections of orders O(α 7 ) and O(α s α 6 ) for the process pp → µ + µ − e + ν e jj + X.   Figure 2e with gluon exchange within one quark line contribute. In both types of NLO corrections, partonic channels with initial-state photons are not taken into account, since their contribution is already strongly suppressed at LO. Channels with external bottom quarks are excluded as well. Those could only significantly contribute via singly-resonant top quarks, which corresponds to a different experimental signature. In total, 40 partonic channels must be taken into account at each coupling order with up to ∼ 83,000 1-loop Feynman diagrams contributing per channel. Tensor integrals appear up to 8-point functions with tensor ranks of up to 4.
In the VBS approximation, as employed in previous QCD calculations, only QCD corrections of the type Figure 2e with gluon exchange within one quark line are taken into account. With ∼ 1000 Feynman diagrams per partonic channel and up to 5-point functions this approximation reduces the complexity drastically in comparison to the calculation presented in this article.

Real corrections
At the order O α 7 , there are two types of real EW corrections: One is due to photon radiation, which results from radiating a photon from one of the charged particles of the LO processes of order O α 6 . The other type comprises photon-induced channels, which we do not take into account at NLO, as already the corresponding LO contribution turned out to be very small. 1 Therefore, only real photon radiation from the qq → µ + µ − e + ν e qq channels, i.e. the process qq → µ + µ − e + ν e qqγ, is considered. Some relevant Feynman diagrams are shown in Figures 3a and 3b. The related infrared (IR) divergences are subtracted using QED dipole subtraction [18,19]. At the order O α s α 6 a mixture of two types of real radiation contributes, because this NLO contribution comprises both QCD corrections to the order O α 6 and EW corrections to the order O α s α 5 . The EW corrections are obtained by attaching a photon to each LO diagram of order O g 2 s g 4 in all possible ways (see Figure 3c) and interfering the resulting diagrams with all photon emission diagrams of O g 7 . The QCD corrections are obtained by attaching a gluon to each LO diagram of order O g 6 in all possible ways (a sample diagram is given in Figure 3d) resulting in the process qq → µ + µ − e + ν e qqg, and squaring the corresponding amplitude. Of course, there are also real QCD corrections of the same order with the gluon crossed into the initial state, gq/qg → µ + µ − e + ν e qqq. The phase-space integration for the real corrections of O α s α 6 leads to both QCD and EW IR divergences  in the limits of soft and/or collinear gluon or photon emission, or via forward branchings of QCD partons in the initial state. Figure 3a shows a q → qγ splitting of QED type, Figure 3d displays a typical q → qg splitting of QCD type.
For real radiation of order O α s α 6 further subtleties arise. Some diagrams with external real gluons involve singularities associated with soft/collinear photons. One example is given in Figure 4a which has an initial-state collinear singularity and requires both QED and QCD dipoles to subtract all IR divergences. Another subtle case arises from final states involving a qq pair. This pair can result from a QED splitting γ * → qq where the off-shell photon has a very low virtuality (see Figures 4b and 4c). In the singular limit where this virtuality goes to zero, a collinear singularity develops with a universal singular structure factorising from the hard matrix elements of the underlying process with a real photon instead of the qq pair. Note, however, that the physical final state is still a jet, or at least some hadronic activity, emerging from the photon initiating the splitting. Technically, the collinear singularity can, e.g., be separated via dipole subtraction as described in Ref. [19], i.e. an auxiliary subtraction function is subtracted from the original integrand, rendering the resulting contribution integrable over the singular region. The formerly subtracted contribution is added back after integration over the singular region with the help of some regularisation, either by switching from four to D space-time dimensions or by employing small quark masses. Either way, the resulting singular contribution is not yet described in a physically meaningful way, since the splitting contains non-perturbative contributions. In the case of low-virtuality γ * → qq splittings, this contribution can be obtained from a dispersion integral for the R ratio of the cross sections for e + e − → hadrons/µ + µ − , as will be further detailed in Ref. [20]. As described there, this contribution can be tied to the quantity ∆α had , which is derived from experimental data. In our calculation we follow this approach, i.e. we separate the singularity via dipole subtraction [19] and add the nonperturbative fragmentation-like "photon-to-jets conversion part" from the collinear region based on ∆α had . Conceptually, it is quite important to properly treat this non-perturbative contribution, but in the present case the overall contribution matters only at the level of 10 −4 relative to the EW LO cross section.
Note that the previous discussion also applies in principle to the real corrections of order O(α 7 ) with an additional photon. However, the extra collinear singularities coming from the matrix elements similar to the ones in Fig. 4, but with external photons instead of gluons, are cut off due to our process definition. In particular, we treat a final-state photon in the real matrix elements always as a photon and never as a jet (or photon-jet). This choice implies that phase-space points with such collinear singularities have either zero or one jet. For final-state singularities, the two collinear quarks are clustered into a single jet while for an initial-state singularity, the collinear quark is along the beam pipe making it undetectable. Therefore our requirement of having at least two jets (see Section 4.1) cuts away such singularities, rendering the real corrections finite, so that no additional terms of type γ → qq are required.

Implementations
In order to ensure the correctness of the results, two independent Monte Carlo programs have been developed based on two entirely different sets of matrix elements constructed by independent matrix-element providers. One calculation is based on the combination BONSAY+OpenLoops, the other on MoCaNLO+Recola.
The program BONSAY is a general-purpose Monte Carlo integrator which can be used to calculate arbitrary NLO EW, QCD, and mixed corrections with matrix elements from an external provider. It has already been used before in Ref. [13] to calculate QCD corrections of like-sign WW scattering. It employs many different phase-space mappings that are combined via multi-channel techniques [21], similar to the Lusifer Monte Carlo program [22], but allows to run the integration in parallel on clusters using MPI [23].
MoCaNLO is also a generic Monte Carlo program, designed to compute arbitrary cross sections in the SM at NLO QCD and/or EW accuracy. The efficient integration is ensured by using phase-space mappings similar to the ones of Refs. [22,24,25]. It has already been used to compute NLO QCD and EW corrections for several high-multiplicity processes [26][27][28][29], including the like-sign W-boson VBS process [8,9,13]. Moreover, it has also been tested against other independent codes for the computation of EW corrections to di-boson production in Ref. [10].
In both Monte Carlo programs, IR divergences in the real radiation are handled with the dipole-subtraction method for QCD [30] and its extension to QED [18,19]. Although the same algorithms are used, the two implementations are completely independent. The library LHAPDF [31] provides PDFs in both codes.
Both OpenLoops [14,15] and Recola [32,33] use the Collier library [16,17] to obtain numerically stable results for the one-loop scalar [34][35][36][37] and tensor integrals [38][39][40]. In order to ensure independence, the two different modes of Collier have been used: the DD mode in BONSAY+OpenLoops and the COLI mode in MoCaNLO+Recola. The intermediate W/Z-boson resonances are treated in the complex-mass scheme [24,41,42] to ensure gauge independence of all LO and NLO amplitudes. In OpenLoops, we set the option use_cms=2 to switch to the same conventions for the complex-mass scheme as used by Recola. This allows us to verify the agreement of the results of the two matrix-element providers for individual phase-space points.
The numerical results presented in this article are obtained from BONSAY, 2 which agrees with the other implementation within integration errors, which are typically of the size of a per mille with respect to the LO prediction.

Validation
The first and strongest validation is that the final results (at the level of cross sections and for each bin of the differential distributions) of the two calculations agree within statistical errors. This constitutes a very solid check as the two Monte Carlo programs as well as the matrix-element providers are different and independent. This ensures the correct implementation of the event selection, input parameters as well as the subtraction on the one hand. On the other hand, it also ensures the validity of the matrix elements used.
The α dipole parameter [43] allows one to restrict the phase space to the singular regions, where α dipole = 1 corresponds to the full phase space (within the acceptance defined by selection cuts) without additional restrictions. Varying α dipole allows then for a robust check of the subtraction procedure. Representative contributions have been checked between BONSAY+OpenLoops and MoCaNLO+Recola for α dipole = 1 for both orders O α s α 6 and O α 7 . The final results have been obtained with α dipole = 1 for BONSAY+OpenLoops and α dipole = 10 −2 for MoCaNLO+Recola, showing agreement at the level of a per mille. This constitutes a strong check on the subtraction procedure used.
In addition, point-wise comparisons of (squared) matrix-element contributions have been carried out for the virtual corrections. At the order O α 7 , which comprises the numerically most delicate loop amplitudes, for 1000 phase-space points chosen in the fiducial volume described above, more than 99% of the points show at least 6 digits of agreement. In total, the level of agreement spans from about 6 to 12 digits.
Finally, 1000 points have been generated to check the real QCD corrections. This ensures the correct implementation of the event selection for both the real radiation and the dipoles. In that way, the correct implementation of the dynamical scale is ensured as well.

Input parameters and event selection
The results presented are for the LHC operating at a centre-of-mass (CM) energy of 13 TeV. We use the NLO NNPDF 3.1 QED set [44,45] with the photon PDF determined by the LUXqed method [46,47] and α s (M Z ) = 0.118 (LHAPDF ID 324900), employing the fixed N F = 5 flavour scheme throughout. We use the same PDFs for LO and NLO predictions. Both QCD and QED singularities from collinear initial-state splittings are factorised into redefined PDFs using the MS factorisation scheme.
The central renormalisation and factorisation scales, µ ren and µ fact , are set to the geometric average of the transverse momenta of the jets, The choice of this scale is motivated by the results of Ref. [48] on like-sign WW scattering, where it was shown that this choice reduces the difference between LO and NLO QCD predictions at large transverse momenta significantly. In the following, we perform a 7-point scale variation of the renormalisation and factorisation scales, i.e. apart from the "diagonal" variations µ ren = µ fact = µ 0 , µ 0 /2, 2µ 0 we set each of the two scales to µ 0 /2, 2µ 0 while keeping the other scale fixed.
Regarding the electromagnetic coupling, the G µ scheme (see, e.g. Refs. [49,50]) is used, i.e. the coupling is obtained from the Fermi constant G µ as The masses and widths of the massive particles read [51]  The bottom quark is considered massless and is neglected in the initial state by default. The width of the top quark is set to zero as it is never resonant, except for the O α 6 contributions with external bottom quarks, which we consider separately; there we set the top-quark width to Γ LO t = 1.449582 GeV [52]. The Higgs-boson mass is taken according to the recommendation of the Higgs cross section working group [53] with its corresponding width. The pole masses and widths entering the calculation are determined from the measured on-shell (OS) values [54] for the W and Z bosons according to The set of acceptance cuts is taken from the recent CMS measurement, more precisely, the ones of the loose fiducial region defined in Ref. [7]. Experimentally, the final state of the process is required to have three charged leptons and at least two jets. QCD partons are clustered into jets using the anti-k T algorithm [55] with jet-resolution parameter R = 0.4. Similarly, photons from real radiation are recombined with the final-state quarks into jets or with the charged leptons into dressed leptons, in both cases via the anti-k T algorithm and a resolution parameter R = 0.4.
In MoCaNLO only partons with rapidity |y| < 5 are considered for recombination, while particles with larger |y| are assumed to be lost in the beam pipe. In BONSAY all partons are considered for recombination, regardless of their rapidities. This difference turns out to be numerically irrelevant in our set up.
The pseudo-rapidity η and the transverse momentum p T of a particle are defined as where |p| is the absolute value of the three-momentum p of the particle, p z the component of its momentum along the beam axis, and p x , p y the components perpendicular to the beam axis. The charged leptons are required to pass the acceptance cuts In addition, an invariant-mass cut on the decay products of the Z boson is applied: A recombined QCD parton system is called a jet if it obeys the jet-identification criteria where the last condition requires a minimal distance between a jet and each of the charged leptons. The identified jets are then ordered according to the magnitude of their transverse momenta p T,j,i , where p T,j,1 denotes the largest p T,j value in the event and p T,j,2 the second largest. The distance ∆R ij between two particles i and j in the pseudo-rapidity-azimuthalangle plane reads ∆R ij = (∆φ ij ) 2 + (∆η ij ) 2 , (4.9) with ∆φ ij = min(|φ i − φ j |, 2π − |φ i − φ j |) being the azimuthal-angle difference and ∆η ij = η i − η j the rapidity difference. On the invariant-mass and rapidity separation of the leading and sub-leading jets, i.e. on the two jets with largest transverse momenta, the following VBS cuts are applied: M jj > 500 GeV, |∆y jj | > 2.5. (4.10)

Cross sections
We start our discussion of numerical results by reporting LO cross sections in the fiducial region. In Table 1 the cross sections at the orders O α 6 , O α s α 5 , and O α 2 s α 4 are shown for the central scale. In contrast to the like-sign WW channel, where the EW contribution largely dominates over the QCD one, here the EW contribution is smaller than the QCD contribution by about a factor four. The LO interference contribution of O α s α 5 only amounts to 0.5 % and is, thus, phenomenologically unimportant.
Taking into account scale variation as defined after (4.1), the LO cross section for the quark-induced EW contribution (often referred to as signal in experimental analysis) is:   α 4 ), and involve bottom quarks in the initial and/or final state. All contributions are given in fb as well as relative to the LO EW cross section of order O α 6 (in percent). The digits in parentheses indicate the integration errors.
Note that this order does not involve any strong coupling, which explains the relatively low scale dependence. We do not show the scale dependence of the LO contribution at the orders O α s α 5 and O(α 2 s α 4 ), since the corresponding NLO contributions balancing their scale dependence are not part of this calculation.
In addition to these quark-induced EW contributions, we have also computed all LO contributions featuring a photon in the initial state. This includes contributions with initial states gγ, qγ as well as γγ at orders O α 6 or O α s α 5 . As can be seen from Table 2, these contributions are phenomenologically negligible. In addition, the LO contributions at order O α 6 involving bottom quarks either in the initial state or in the final state are also reported. While the contributions with two bottom quarks in the initial state are negligible due to their PDF suppression, the contributions with one light quark and one bottom quark in the initial state are rather large. The latter are usually referred to as tZ + jet contributions in experimental analyses (see Figure 1f). These contributions are enhanced due to resonant top-quark contributions. In the final state they have one b-jet and one light jet and can therefore be suppressed in experimental analyses using b-jet tagging techniques. Note that these contributions also contain VBS contributions (for instance diagram Figure 1a with Table 3  On the other hand, the EW corrections of order O α 7 amount to −16 % and represent the dominant NLO contribution. This is in line with the findings of Ref. [8] for like-sign WW scattering and supports the expectation that large EW corrections are an intrinsic feature of VBS at the LHC. Following Ref. [56], one can derive a leading logarithmic approximation for the EW corrections to the process pp → µ + µ − e + ν e jj + X based on the logarithmic corrections to the sub-process WZ → WZ. Taking the mixing of photon and Z boson into account and using M WZ→Wγ ≈ − sw cw M WZ→WZ , one arrives at the approximation already given in Ref. [8] for WW → WW. This approximation holds in fact for all scattering processes of EW bosons owing to the fact that these scattering processes result from the same SU(2) w coupling. The approximation reads dσ LL = dσ LO (1 + δ EW,LL ) , (4.12) where  showing a significant reduction of scale uncertainty. This is mainly due to the O(α s ) PDF redefinition included in the O(α s α 6 ) NLO correction that cancels the factorisation scale dependence of the LO O(α 6 ) contribution. As shown in Table 3, the full NLO correction is about −17.8 % with respect to the LO of order O α 6 . Finally, for completeness, we also provide cross sections at NLO for the LHC running at a centre-of-mass energy of 14 TeV in Table 4. While the LO cross section increases by 17.2% with respect to 13 TeV, the relative NLO corrections are rather stable. These numbers can be important for future operation of the LHC at high luminosity [57] and serve as benchmarks.

Differential distributions
In this section, LO predictions and NLO corrections for several differential distributions are discussed. We start with a few LO predictions in Figure 5. The upper panels show the absolute predictions of order O α 6 (EW), O α s α 5 (interference), and O α 2 s α 4 (QCD). In the lower panels, the relative contributions are displayed with respect to the sum of the three contributions. Note that the contributions featuring external bottom quarks or initial-state photons are not included here. The first two distributions are the invariant mass and pseudo-rapidity difference of the two tagging jets in Figures 5a and 5b. These observables are often used to separate EW and QCD contributions in experimental analysis. This is perfectly justified by the fact that at higher invariant mass or larger pseudo-rapidity, the EW contribution is becoming dominant. The effect of the event selection for |∆η j 1 j 2 | > 2.5 and |M j 1 j 2 | > 500 GeV is clearly visible. In Figure 5c, the transverse momentum of the second hardest jet is shown. Around 500 GeV, both the EW and QCD contributions become of the same size, as the QCD contribution is falling much more steeply than the EW contribution. Interestingly, this behaviour is not visible in other transverse-momentum and invariant-mass distributions, like the transverse-momentum distribution of the leading jet, where the EW contributions are always larger than the QCD contributions. The comparably steep fall of the distribution in the transverse momentum of the subleading jet is due to the fact that QCD contributions are dominated by contributions with at least one jet with small transverse momentum. Finally, Figure 5d displays the distribution in the rapidity-azimuthal-angle distance between the leading jets which also shows a good discriminating power, as already noticed in Ref. [10]. Note that in all distributions, the interference contribution is very much suppressed reflecting its overall small cross section.
The following figures show our results on NLO differential distributions. In the upper panels, the LO contribution of order O α 6 is shown along with the NLO predictions including orders O α 7 or/and O α s α 6 . For simplicity, these are often denoted by EW and QCD corrections, respectively, in the following. We stress again that while the order O α 7 comprises genuine EW corrections, the order O α s α 6 contains both QCD and EW corrections. We adopt this assignment in order to facilitate the reading. The lower panels show the relative contributions where the bands in the plots reflect the variation of the numerator with the (renormalisation and/or factorisation) scale µ while keeping the scales in the denominator fixed to µ 0 . In Figure 6, several distributions in transverse momenta are presented. We start with those for the hardest and second hardest jet in Figures 6a and 6b, respectively. For both distributions, the EW corrections display the typical Sudakov behaviour in the high-energy region, where they become large in size and negative. The QCD corrections are positive for low transverse momentum of the leading jet, but steadily decrease towards high transverse momentum, becoming negative above 150 GeV and stabilising in the range 200−300 GeV. This behaviour is due to the reshuffling of jets from higher to lower transverse momentum by the real QCD corrections and has already been observed in like-sign WW scattering [9]. For the transverse momentum of the second leading jet, the QCD corrections turn again positive towards high transverse momentum. The enhanced corrections for small transverse momentum of the leading jet are due to the phase-space suppression of the LO when all jet transverse momenta are required to be small. This causes corrections above 20% for small transverse momenta of the hardest jet, while the corrections almost vanish for small transverse momenta of the second hardest jet. For the distributions in the missing transverse momentum (Figure 6c), which is identified with the neutrino transverse momentum p T,νe , and in the transverse momentum of the positron p T,e + (Figure 6d), the EW corrections follow the same trend as for the jet transverse momenta. The QCD corrections are almost independent of p T,miss and p T,e + until about 400 GeV. Figure 7 displays some angular distributions. For the distribution in the azimuthalangle difference of the two tagging jets (Figure 7a) QCD and EW corrections follow a similar trend and inherit mostly the corrections to the total cross section. This does not hold for the other angular distributions, where the corrections show differences in shape. For the rapidity-azimuthal-angle distance between the two jets (Figure 7b), the QCD corrections reach a minimum around ∆R j 1 j 2 = 4, while the EW corrections tend to increase slightly towards increasing ∆R j 1 j 2 , resulting in an increase of the combined NLO prediction. The QCD corrections are generally flat and small for distributions in leptonic angular variables, resulting in combined predictions very close to the EW ones. The distribution in the rapidity-azimuthal-angle distance of the muon and anti-muon (Figure 7c) displays increasing EW corrections with increasing ∆R µ − µ + , varying from −45% at ∆R µ − µ + → 0 to about −10% at ∆R µ − µ + = 4. The distribution in the cosine of the angle between the muon and anti-muon (Figure 7d) receives only a mild shape distortion towards cos θ µ − µ + → 1 from the EW corrections.   Figure 8 shows pseudo-rapidity distributions. The first two concern the hardest ( Figure 8a) and second hardest jet (Figure 8b). Both QCD and EW corrections are rather similar in shape and differ mainly by some offset. For the hardest jet, the corrections peak in the peripheral region, while for the second hardest jet they increase in the central region as well. The distribution in the pseudo-rapidity of the anti-muon (Figure 8c) receives flat corrections over almost the whole range with only a slight increase in the peripheral region. Finally, both QCD and EW corrections increase with growing pseudo-rapidity difference of the two leading jets.
Finally, we show the distributions in the invariant mass of the two jets ( Figure 9a) and in the transverse mass of the WZ system (Figure 9b)  the two jets displays a similar behaviour than the one in the transverse momentum of the jets: There is a steady increase in size of the negative EW corrections towards more and more negative values with increasing invariant masses due to EW high-energy logarithms. The QCD corrections are positive at 500 GeV and decrease slowly towards higher invariant masses. This is particularly interesting, as this observable is used to define fiducial regions in measurements. The transverse mass M T,W + Z , of the WZ system is defined as where is running over the four leptons (including the neutrino). In the region below

Conclusion
The process pp → µ + µ − e + ν e jj + X is of great interest at the LHC, because its EW contribution of order O α 6 to the cross section contains vector-boson scattering (VBS) as a subprocess. In this article we have reported on a calculation of NLO corrections  of order O α s α 6 and O α 7 to the EW process. This is the first time that the EW corrections of O α 7 are computed for such a final state. While the QCD corrections of order O α s α 6 have already been computed in the VBS approximation [11], for the first time their full computation (including interference contributions of EW type) is performed. The combination of these two NLO contributions constitutes the complete NLO prediction for the EW component of the pp → µ + µ − e + ν e jj + X process.
The EW corrections turn out to be relatively large, in accordance with similar observations made already for like-sign W scattering in Ref. [8]. This confirms the expectation that large EW corrections are indeed an intrinsic feature of VBS at the LHC. The corrections reach −16% in the chosen fiducial region and are driven by Sudakov logarithms that are large and negative and grow in size in the high-energy limit. The large EW corrections for the integrated fiducial cross section can be well reproduced by a simple logarithmic approximation. The corrections of order O α s α 6 (sometimes simply called QCD corrections) are small and negative. They amount to about −2% at the level of the integrated cross section and do not exceed 20% in differential distributions.
To complete the picture for the VBS process, we have also reported on several LO contributions, including some suppressed channels. The interference contribution of order O α s α 5 is below 1%, while the QCD contribution of order O α 2 s α 4 is larger than the EW process by a factor 4. This shows how challenging the measurement of the EW process of order O α 6 is, highlighting therefore the need for precise predictions in this context. Moreover, we have computed all the LO contributions with a photon in the initial state. These turn out to be rather small and can safely be neglected in future analyses. Finally, we have computed the contributions involving bottom quarks (either both in the initial state or one in the final and one in the initial state) which turn out to be non-negligible.
These are enhanced by contributions with a singly-resonant top quark, but can, in principle, be experimentally suppressed using b-tagging techniques.
In addition to the phenomenological relevance of the presented calculation, it is worth stressing that this also constitutes a non-trivial extension of previous computations. The complications arising here are manifold: (i) This is the first time that a process is computed at NLO EW accuracy with seven external charged particles. This increases the complexity of the real contribution as well as of the virtual corrections. (ii) The computation of the real contribution entails additional complications. QCD real radiation contains singular contributions arising from soft/collinear photons and gluons. This requires an advanced automation of computation of the real corrections. (iii) The number of partonic channels is very large. With respect to the like-sign W channel, the number of partonic channels increases by more than a factor three. This means that the computing time is considerably increased, and efficient book-keeping and parallelisation become decisive.
To come as close as possible to the situation realised in experimental analyses, our results are given in terms of integrated cross sections and differential distributions in the so-called loose fiducial region presented in Ref. [7] by the CMS collaboration. The event selections for this phase space are simplified with respect to the ones used for the actual measurement. Such efforts are particularly welcome to theorists as they allow a direct use of state-of-the-art theoretical progress in experimental analyses. Therefore the predictions provided in the present article should be particularly useful for the VBS program at the LHC.