Vector-boson pair production and electroweak corrections in HERWIG++

The detailed study of vector-boson pair production processes at the LHC will lead to a better understanding of electroweak physics. As pointed out before, a consistent inclusion of higher-order electroweak effects in the analysis of corresponding experimental data may be crucial to properly predict the relevant phenomenological features of these important reactions. Those contributions lead to dramatic distortions of invariant-mass and angular distributions at high energies, but may also significantly affect the cross section near threshold, as is the case e.g. for Z-pairs. For this reason, we present an analysis of the next-to-leading-order electroweak corrections to WW, WZ and ZZ production at the LHC, taking into account mass effects as well as leptonic decays. Hence, our predictions are valid in the whole kinematic reach of the LHC and, moreover, respect the spin correlations of the leptonic decay products at NLO accuracy. Starting from these fixed-order results, a simple and straight-forward method is motivated to combine the electroweak corrections with state-of-the-art Monte Carlo predictions, focusing on a meaningful combination of higher-order electroweak and QCD effects. To illustrate our approach, the electroweak corrections are implemented in the HERWIG++ generator, and their phenomenological effects within a QCD environment are studied explicitly.


Introduction
generator (Sect. 4), relying on K-factors for unpolarized two-by-two scattering and the assumption of factorization of EW and QCD effects, as will be motivated in Sect. 3. Our method will be tested against an alternative implementation method especially tailored for the HERWIG++ [24,25] MC generator. In Sect. 5 selected numerical results obtained within the HERWIG++ setup will be discussed.

Electroweak corrections
As stated in the introduction, the complete evaluation of the combined QCD and electroweak corrections to gauge-boson pair production of order αα s is presently out of reach. As a first step we, therefore, consider a factorized ansatz, where the kinematics of events with additional QCD radiation is mapped to effective two-body collisions, described by effective (squared) partonic energies and energy transfersŝ andt. We assume that the bulk of the QCD corrections arises from events with soft or collinear emission, where the factorized ansatz is expected to work well. The EW corrections are taken directly from the result of the one-loop calculation, evaluated at the same kinematical point. This approach is expected to fail for events with large transverse momenta of the gauge boson pair recoiling against a jet with large transverse momentum. These events, however, are of lesser relevance for the study of gauge-boson dynamics and can be eliminated by suitable cuts, as discussed in Sect. 4. The motivation of this approach will be discussed in Sect. 3, and details of the implementation are given in Sect. 4. In the present section we concentrate on the electroweak corrections and motivate that indeed the bulk of the electroweak corrections can be collected in a K-factor which is given as a function ofŝ andt only (Sect. 2.1). Photonic corrections (which evidently lead to a more complicated kinematic situation) can be split off such that the corresponding modifications of the electroweak corrections are small. This aspect will be investigated in Sect. 2

.2.
A second simplification is introduced by applying a correction factor which does not depend on the helicities of the gauge bosons. In Sect. 2.3 we argue that this approximation still preserves the proper angular distributions and correlations of the Z and W decay products and investigate the phenomenological implications of this approximation in detail. Finally, the corrections as derived for on-shell gauge boson production are applied for the cases where W or Z are slightly off mass shell (In the case of Z bosons we also include the amplitude where fermion pairs are produced through the virtual photon). The quality of this approximation at lowest order will also be discussed in Sect. 2.3. First, however, let us briefly recall the most important phenomenological features of various EW effects in on-shell WW, WZ and ZZ production at the LHC. In W-pair production, the invariant-mass distribution ( Fig. 1 (top)) receives wellknown logarithmically enhanced negative EW corrections (δ EW ) growing with energy. Positive contributions arise from the partonic subprocess γγ → WW (δ γγ ) and the photonquark induced processes (δ veto γq ), the latter evaluated applying the jet veto defined in Ref. [7]. As already pointed out in Ref. [8], the effect of massive-boson radiation (δ 3V ) is moderate, however, strongly dependent on the event selection.
The above picture significantly changes if angular distributions of the W-pair are studied at high invariant masses. This can be seen in Fig. 1 (bottom) where distributions of the rapidity gap of the two Ws are shown for M WW > 1000 GeV. While the genuine EW corrections drastically reduce the differential cross section at high p T,W , corresponding to small rapidity gap, the photon-induced contributions significantly increase the rate at small scattering angles, corresponding to large rapidity gap. As a result, a dramatic distortion of the angular distribution is visible which might easily be misinterpreted as signal of anomalous couplings. We point out that the photon-induced corrections presented above (which are obtained using the MRST2004qed PDF set [26] for the photon density) suffer from a large systematic error stemming from our ignorance of the photon content of the proton. This becomes obvious looking at Fig. 25 of Ref. [27], where the NNPDF2.3QED [27] set has been used to estimate the error on the γγ-induced W-pair cross section. A relative error of ±50% on the leading-order (LO) cross section at M WW = 1000 GeV can be deduced, solely induced by the photon PDF error. This picture indicates that a significant improvement in the determination of the photon PDFs is mandatory to reliably predict the W-pair production cross section at high energies.
Turning to WZ production, the situation is qualitatively similar to WW production, though here the γγ process is absent and the genuine EW corrections are smaller (Fig. 2). In Z-pair production, however, the γq-induced contributions are negligible [21], and the real-radiation contributions are always below 10%. In total, particularly large negative corrections, reaching −40%, dramatically affect Z-pair production at high invariant masses and transverse momenta (Fig. 3).

K-factors for the electroweak corrections
Let us start with the simplest reaction, inclusive Z-boson pair production. The full set of electroweak corrections, including photon radiation, Z-boson mass effects and virtual top quarks, has been discussed in Ref. [8], where the purely weak corrections were also evaluated for the 4-lepton final state in the pole approximation. The real and virtual QED corrections which can be considered as gauge invariant subset were included in the onshell analysis. However, their contribution is relatively small, in general below 1%. This is demonstrated in Fig. 4 for four characteristic distributions: the distributions in Z-boson rapidity and transverse momentum, in the invariant mass of the Z pair and, furthermore, in the difference between the rapidities of the two Z bosons. Note that in fact three correction functions are required, one for the uū induced reaction which we denote as  Figure 4: Various differential distributions for on-shell Z-pair and γ-pair production at LHC14. Left: LO predictions; right: relative weak corrections (δ V V weak ) and the full set of electroweak corrections (δ V V EW ), including QED contributions. All results are obtained in the default setup defined in Ref. [8].
the others, K ZZ dd and K ZZ bb for the dd and bb induced reactions. All these factors depend on the partonicŝ andt only.
Similar considerations apply to di-photon production, and results analogous to ZZ production are also shown in Fig. 4. Again the relative contribution from purely photonic corrections is small and the two-body approximation reproduces the full answer to better than 1 percent. The correction factors K γγ uū , K γγ dd andK γγ bb are defined in an obvious way. The results presented in Fig. 4 were obtained using the default setup as defined in Ref. [8].
The situation is more involved in the case of W-pair and WZ production where weak and electromagnetic corrections are intimately intertwined. As is well known, virtual photonic corrections are required to arrive at an ultraviolet finite answer. The resulting infrared divergencies are then canceled by real radiation, the remaining collinear singularities are finally absorbed into a redefinition of the PDFs. Since one is now dealing partially with a three body final state, the direct use of a K factor which depends onŝ andt only is nontrivial.
In order to return to the kinematics of a two-to-two body reaction the corrections from real radiation are now replaced by just subtracting the endpoint singularities as defined originally in Ref. [28] and described in the Appendix. As will be shown explicitly in Sect. 2.3, the physical predictions remain practically unaffected by this simplification. (The agreement at the level of around one percent is demonstrated explicitly in Figs. 6 and 7 through the comparison between δ full EW and δ V+E EW .) At the same time the kinematics of the resulting final states is now fully described by the partonicŝ andt. In this way we obtain again three correction factors K WW uū , K WW dd and K WW bb which can be employed in the framework of the Monte Carlo generator, just as before. A similar approach is valid for W + Z and W − Z production, which have, of course identical correction factors, denoted K WZ . The corresponding endpoint contributions are also listed in the Appendix.

Gauge-boson polarization and four-lepton production
Our corrections are presented for unpolarized gauge bosons, which in the case of W and Z are observed through their decay products. The distributions of the decay products, however, are affected by the boson polarization. Hence an additional complication could result from the fact that the polarization pattern of the gauge bosons is modified by the radiative corrections. In principle one would have to employ K factors for the full set of helicity amplitudes. However, as demonstrated in Ref. [8] for ZZ production, in practice a fairly simple pattern emerges. Let us first consider the case of Z pairs: For small transverse momenta the electroweak corrections are small (about -4%) and of similar magnitude for all four combinations of transverse and longitudinal polarizations. For large transverse momenta one single configuration dominates completely and corrections for the subdominant combinations are irrelevant. This feature has been demonstrated in Table 7 of Ref. [8], where the cross sections and the corrections are displayed in the low-, intermediate-and large-p T region, separated according to longitudinal and transverse polarizations.
From these considerations it becomes clear that in case of Z-pair production a single partonic K factor is sufficient for large as well as for small transverse momenta, and polarizations of the gauge bosons and, correspondingly, the correlations between the decay products as predicted in Born approximation are maintained even after inclusion of the electroweak corrections. The situation is slightly different for W-pair (and also WZ) production. Here, roughly 10% of the LO cross section are given by longitudinally polarized Ws even at high transverse momenta, and the corresponding relative EW corrections are substantially different compared to the transversely-polarized case (see Fig. 7 of Ref. [20]). However, the numerical effects from the longitudinal polarization on the radiative corrections are small, and it is still sufficient -as will be shown later -to apply unpolarized K-factors to reproduce the full corresponding EW corrections with sufficient accuracy.
The corrections evaluated in Refs. [7,8] and encoded in our K factors were obtained for on-shell Z or W bosons. Any realistic simulation of four-fermion production must, necessarily, include contributions from off-shell configurations. In the case of Z also digrams with off-shell Z replaced by virtual photons would be required for a description away from the Z peak. However, for the experimental analysis of gauge boson production the invariant mass of the decay products (lepton pairs or jets) must be restricted to an interval around the nominal mass, say |M ll − M Z | < 25 GeV, to suppress the admixture of virtual photons and enhance close-to-mass-shell gauge bosons. In this case the neglection of virtual photons as implemented in HERWIG++ can be justified 1 .
In our approach off-shell effects and non-resonant contributions are consistently accounted for in the LO predictions, while at NLO the NWA is applied to compute the relative corrections. This approximation is well motivated, since off-shell effects in the EW corrections only amount to ∼ 0.5%, as demonstrated in Ref. [19] for W-pair production.
From these considerations it can be expected that the dominant weak corrections to four-lepton production at the LHC are well described by process dependent K-factors which can be taken from the unpolarized results for the corresponding 2 → 2 production process. The validity of the approximations discussed above will now be studied in detail.
For this purpose we give a concise presentation of the computation of EW corrections to massive gauge-boson pair production at the LHC, consistently taking into account leptonic decays and, related to this, spin correlations. In particular, we discuss the validity of the various approximations discussed in the previous subsection. To allow for a sensible comparison with the HERWIG++ results presented in Sect. 5, we generally stick to the default HERWIG++ setup for gauge-boson pair production as defined in Ref. [29]. Specifically, in the leptonic event selection we apply the following basic cuts for the charged-lepton transverse momenta and rapidities. If neutrinos are present in the final state, a minimal missing transverse momentum of is also required. The invariant mass of the lepton pair is restricted to  Figure 6: Various differential distributions for e + ν e µ − µ + production at LHC13. Left: The full LO prediction as well as NWA and DPA are shown; right: full relative EW corrections δ full EW evaluated in the NWA including spin correlations; EW corrections evaluated with unpolarized 2 → 2 K-factors in the V+E approximation (δ V+E EW ); relative deviations of NWA and DPA w.r.t. the full LO are also shown. and the weak coupling constant is defined in the G µ scheme to systematically absorb universal corrections related to the running of α to the weak scale in the LO predictions. All remaining SM input parameters are directly carried over from Ref. [8].
For the evaluation of hadronic cross sections we use the CT10NLO PDF set [30] in the LHAPDF framework [31], and the CKM dependence in the WZ production channels is included at leading order, while in the computation of EW corrections the CKM matrix is set to unity.
At leading order we present full results, including non-resonant and off-shell effects, as well as results in two different approximations. As far as the full LO cross sections are concerned, we have checked that the difference between a naive fixed width implementation and results obtained in the Complex-Mass Scheme (CMS) [32] is at the per-mill level and hardly visible. All results presented here for the full LO cross sections therefore correspond to the naive fixed-width implementation.
In addition to the full results we provide the results for V-boson pair production in the double-pole approximation (DPA) originally discussed in Ref. [18]. Here, the amplitudes for V-pair production and decays are evaluated on-shell, but the Breit-Wigner shape of the resonance is included in the evaluation of the squared matrix elements to account for the dominant off-shell effects. We apply the on-shell projection procedure proposed in Ref. [33] to construct proper on-shell momenta of the intermediate bosons from the four-particle phase-space. Note that in addition to the physical cuts displayed above we impose a technical cut, m 4l > M V 1 + M V 2 , on the 4-lepton invariant mass since the on-shell projection suggested in Ref. [33] only gives sensible results above threshold.
As a third alternative, we work in the narrow-width approximation (NWA), where the gauge bosons are strictly forced on-shell from the beginning via the replacement for the resonant squared propagators. In Figs. 5, 6 and 7 we present results for WW, W + Z and ZZ production at the LHC. Since in our approach the EW corrections to the partonic subprocesses are insensitive to the leptonic decay mode, in each case we concentrate on one specific decay channel, namely The right-hand-side panels of the respective plots show the relative deviations of the NWA (∆ NWA LO ) and DPA (∆ DPA LO ) from the full results. For all pair-production channels one observes that the NWA overshoots the full LO predictions at the level of 5%, while the DPA results give a good approximation valid at the 2-3% level. From the upper right plot of Fig. 5 it becomes obvious that both approximate predictions become crude near the Z-pair production threshold where off-shell effects apparently become more important.
Let us now turn to the EW corrections. As stated before, we compute the full EW corrections to the respective polarized vector-boson pair production processes, and include the spin correlations in the decays to leptons, which are treated at leading order. The actual computation is carried out in the well-established FeynArts/FormCalc/Loop-Tools [34][35][36][37][38] setup already used in Refs. [7,8], and Madgraph [39] was useful for internal checks. We do not take into account the EW corrections to the leptonic decay processes for various reasons. On the one hand, the bulk of the EW corrections to Z and W decay are given by final-state photon radiation (FSR), which leads to significant distortions of the phase-space distributions of the leptonic decay products. These FSR corrections, however, are included in HERWIG++ [40] in the YFS framework [41] and will therefore not be considered here. On the other hand, electroweak corrections to the inclusive boson decay widths are implicitly included in the experimentally determined values for the branching ratios used in the HERWIG++ framework. Additionally, we strictly stick to the NWA for the computation of the EW corrections. In this simplified approach, the corrections completely factorize into corrections either to the production or the decay process, and no non-factorizable corrections, which connect production and decay, have to be considered. Those contributions have to be taken into account using the DPA as demonstrated in Ref. [18]. However, it is known that the non-factorizable corrections largely cancel in sufficiently inclusive observables [42].
In addition to the full EW corrections δ full EW (which contain proper spin correlations and, in case of WW and WZ production, also the full set of QED corrections to the respective production process) approximate results δ V+E EW in the virtual+endpoint (V+E) approximation are presented in the same plots, where unpolarized on-shell K-factors have been used to obtain the relative corrections, as detailed in Sect. 2.2. One observes that the approximate ansatz gives an almost perfect approximation for the full result, in general better than 1%. The best agreement is observed for ZZ production, while for WW production a slight discrepancy is visible. This can be understood recalling that photon radiation and related QED corrections are largest for WW production, as has been demonstrated in Ref. [8], while they remain small in case of Z pairs. In general, the agreement between the full result and the V+E approximation is even better than expected. As a conclusion one finds that it is justified to use unpolarized K-factors in the V+E approximation to describe vector-boson pair production at the LHC at a sufficient accuracy.

Electroweak and QCD corrections -combined
As stated in the introduction, the complete treatment of combined QCD and electroweak corrections, involving e.g. two-loop terms of order α s α weak is presently out of reach. Consequently neither an additive treatment of the corrections, (1 + δ QCD + δ weak ), nor a multiplicative one, (1 + δ QCD )(1 + δ weak ), will lead to the correct result. However, the bulk of the QCD corrections for the processes under consideration is related to soft or collinear gluon radiation from the incoming quark-antiquark system. Soft radiation leaves the initial state practically unchanged, the quarks remain close to their mass shell and the weak corrections can be taken from theŝ andt dependent K factor.
A similar line of reasoning is applicable to hard collinear radiation which leads to final states with vanishing or small transverse momentum of the diboson system. Let us assume that a collinear gluon is radiated from the incoming (anti-)quark, which subsequently initiates the boson pair production. Also in this case the (anti-)quark stays close to its mass shell. The scattering angle of the q 1q2 → V 1 V 2 reaction in the V 1 V 2 rest frame can again be directly identified with the scattering angle of the partonic reaction. A similar line of reasoning applies to reactions with quarks or antiquarks originating from gluon splitting.
The situation becomes more involved in the case of hard non-collinear radiation, which leads to a diboson system of large transverse momentum. In this case only an approximate prescription can be formulated, since the original 2 → 2 kinematics is distorted. However, this approximate prescription must coincide with the previous one in the limit of vanishing transverse momentum. To be specific we advocate the following strategy to compute the effective partonic Mandelstam variablesŝ andt from the distorted kinematics for the evaluation of K(ŝ ,t ): The squared CM energy is calculated from the four-lepton final state viaŝ = m 2 4l . (3.1) The momenta are boosted into the four-lepton CM frame (denoted by Σ * ). In this frame the unit directions of initial-state hadrons shall be denoted by The direction of the effective scattering axis in Σ * is now defined bŷ and the effective scattering angle is, correspondingly, given by where the v * i denotes the momentum direction of vector boson V i . The Mandelstam variablet is then computed from θ * assuming on-shell kinematics.
Diboson events with large transverse momenta necessarily require the presence of at least one hard quark or gluon jet, and electroweak corrections would have to be evaluated separately for this class of processes. As long, as they can be treated as a small admixture to the diboson sample, suppressed by an additional factor α s , the distortion of the weak corrections should not lead to a significant error for the inclusive sample. As long as one is interested specifically in the analysis of the diboson process, a cut on the transverse momentum of the dibosons system will eliminate the pollution with events of a very different nature. Let us discuss this important issue in some more detail. As pointed out by several groups [7,14], at large transverse momenta V -pair production is dominated by new topologies which are absent at lowest order in QCD. These topologies correspond to V + jet production with the radiation of an additional V from the quark jet rather than being a correction to V -pair production, and spoil the perturbative series for the prediction of leptonic observables at high transverse momenta. They lead to huge QCD K-factors together with large residual scale uncertainties. To improve the corresponding theory predictions the authors of Ref. [14] have provided approximate NNLO QCD predictions for these particular observables in WZ production, applying the LoopSim method [43]. They observe pronounced shifts of the predictions going from NLO to NNLO and at the same time a significant reduction of the theory uncertainties. In this work, however, we follow another strategy to stabilize the problematic high-p T observables, employing suitable restrictions on the transverse momentum of the 4-lepton system, as will be detailed in Sect. 5.
In the present section a fairly general prescription has been suggested how to implement weak corrections into a generic Monte Carlo program which includes QCD corrections and hadronisation already. Indeed this prescription allows for an a posteriori application of EW corrections to any QCD Monte Carlo, at least as long as the partonic origin (uū vs dd) of the final state remains under control. In the next section a slightly different approach will be described which is specifically tailored to the program HERWIG++. In this case a partonicŝ andt is introduced in the program from the very beginning and, in the limit of diboson events with small transverse momenta, coincides with the prescription described above. In the following section this correspondence will be investigated in more detail.

Implementation into HERWIG++
Our starting point of the implementation of the EW corrections in HERWIG++ is the electroweak K-factor. Following the strategy of multiplicative QCD and EW corrections we reweight the events that have been generated in HERWIG++ as the hard processes. For vector-boson pair production HERWIG++ delivers unweighted events, hence we can compute K(ŝ,t) and use this directly as a reweighting factor for each event such that K(ŝ,t) < 1 leads to a suppression of events with given (ŝ,t) while K(ŝ,t) > 1 leads to an enhancement. The actual implementation of a reweighting factor for a given hard event is straightforward in HERWIG++ and ThePEG once the variables (ŝ,t) are known.
Let us discuss the kinematics of our events in more detail. For EW corrections to the Born process the calculation of a K factor is straightforward as we can directly access the kinematic setup of the hard process once this is generated. In this case (ŝ,t) can be computed uniquely. As soon as we want to apply the EW corrections to an event that is generated from an NLO QCD matrix element this is no longer the case. In this case we face the complication that an event with real radiation is not described by 2 → 2 kinematics anymore.
Fortunately, the situation in the HERWIG++ case is largely simplified by the fact that the NLO QCD corrections are applied in the POWHEG scheme [29,44]. (See [45,46] for more details.) The POWHEG scheme provides a consistent matching of NLO corrections and parton shower emissions without double counting. In this scheme, the differential cross section for the hard process can be written as Here, Φ B and Φ R denote the Born and radiative phase space variables, respectively, and k ⊥,min is the minimum transverse momentum that is generated by the parton shower. Φ R only parametrizes the additional variables to specify a hard emission relative to the Born configuration. B(Φ B ) and R(Φ R ) are the Born and real-emission matrix elements squared for the hard process under consideration. B(Φ B ) is the Born differential cross section whilē Here, V (Φ B ) is the (infrared finite) sum of virtual corrections and the divergent part of the real corrections, while R S (Φ B , Φ R ) is the (also finite) real correction matrix element squared with the divergent terms subtracted. Finally, the POWHEG Sudakov form factor is given by which, opposed to the Sudakov form factor in a common parton shower, contains the full real emission matrix element squared. The two terms in Eq. (4.1) are constructed to resemble the result of a single parton shower emission applied to a configuration Φ B with weightB(Φ B ). The first term gives the no-emission probability, the second the contribution from a single emission. At the same time, upon expansion in α s , Eq. (4.1) reproduces the differential cross section at NLO QCD.
In the HERWIG++ implementation we find exactly this prescription. In a first step we generate an event with kinematical configuration Φ B . Technically, this is already the hard process. Only in a second step, already as part of the parton shower algorithm, the potential hard emission with relative kinematics Φ R is generated according to the Sudakov form factor Eq. (4.3). Once this is done, the default parton shower is used as at this point, all terms will become formally higher than next-to-leading order and the default parton shower is computationally much simpler than the Sudakov form factor (4.3).
As the hard emission from POWHEG is technically already part of the parton shower, we have direct access to the Born type variables Φ B in the hard subprocess and hence we may easily compute the kinematic variables (ŝ,t) for every event. In effect, this assumes that the EW corrections are applied only on the level of Born-type kinematics and are not strongly influenced by the hard emission. In fact, we later also focus our analysis to regions where the transverse momenta generated by hard gluon emissions are not too large by applying a suitable veto. This veto will suppress events where gauge-boson pairs are accompanied by additional hard quark or gluon jets, leaving the qq events largely unaffected. This also enforces the kinematics to be reasonably close to a Born configuration in order to justify our approach.
Events with strong QCD activity, e.g. jets with large transverse momentum give rise to large QCD corrections. In order to suppress these enhanced corrections [43] we apply an additional cut on the final state in our analysis. Focussing on the leptonic final state we have to make sure that the gauge-boson pairs or its decay products, the four leptons, are not produced with too much transverse momentum from the recoil against a system of strongly interacting particles. We achieve this by applying a condition on the leptons' transverse momenta i,T in the laboratory frame, Here, we consider all visible leptons i from EW boson decays, i.e. i runs to 2, 3 or 4 in case of WW, WZ or ZZ events. We additionally assume that the missing transverse momentum in the event solely stems from neutrinos from W decays. The left-hand-side of Eq. (4.4) is small whenever the system has a small transverse momentum, as the leptons from the EW system recoil against each other. After some experimenting, we find that ρ = 0.3 gives a good selection of events while leaving enough events for our analysis at the same time.
In Fig. 8 we show histograms of the ratio ρ of vector and scalar sums of lepton transverse momenta that we finally apply the cut on. We show runs with and without NLO QCD corrections. We find that the chosen value ρ = 0.3 is a sensible choice for all three combinations of boson pairs. A good number of events is still available for the analysis while at the same time we find that the events with large distortions from hard radiation are vetoed with our selection. These are peaked at large values of the ratio in all three cases. In a full experimental analysis the value of ρ might be subject to optimization for the individual cases of vector boson pairs.
In addition we show actual leptonic observables with and without application of the lepton veto in Figs. 9-11. Let us compare LO and NLO predictions for the p T,ll distribution in ZZ production (Fig. 9, left) which corresponds essentially to the transverse-momentum distribution of the Z-boson. Without cut the NLO distribution exceeds the LO distribution by a factor of two at p T,ll = 300 GeV and more at larger transverse momenta. Introducing the cut (4.4) with ρ = 0.3 removes most of this excess such that the difference between LO and NLO distributions is reduced to O(20%) and remains relatively constant as function of p T,ll . The ratio between the NLO and LO rapidity distribution, in contrast, is fairly constant with a ratio of NLO/LO ∼ 1.2, and is reduced to ∼ 1 by the cut for ρ = 0.3. A similar behaviour is observed for WZ (Fig. 10) in the e + ν e µ + µ − mode with a giant correction factor of 4 for transverse momenta of the Z boson of 300 GeV and a correction between 1.4 and 1.6 for the rapidity distribution. As before, one obtains a significant reduction of the QCD correction down to O(20%) by introducing the cut (4.4), and a similar behaviour is observed for W + W − production in the e −ν e µ + ν µ channel (Fig. 11). In total, after applying the cut, LO and NLO results are close over the whole range in transverse momentum and rapidity. In order to assess the sensitivity to the choice of partonic kinematic variables we have studied the impact of K(ŝ,t) with two additional choices for the kinematical variables for ZZ production, as in this case we have full access to the kinematics of the final state.
Let us call the variables from the direct access to the kinematics within HERWIG++ before the parton shower (ŝ hard ,t hard ). As an alternative, we reconstruct the kinematics from the vector-boson final state as outlined in the previous section. Here, no knowledge of the initial-state partons is needed for the computation of the kinematic invariants.ŝ rec is given by the invariant mass of the vector-boson pair andt rec is reconstructed from the scattering angle in the centre of mass frame as proposed in Eqs. (3.1)- (3.4). For illustration, we also consider a choice of variables that is very likely to be wrong. We take the initial state partons after the termination of the parton shower. This parton pair will have a much larger invariant massŝ PS due to the parton showering. We then find the four momenta of the outgoing vector bosons and computet PS with respect to this initial state. This last choice will demonstrate us the sensitivity to the parton shower emissions. For each event produced by HERWIG++ we construct three sets of kinematic variables (ŝ i ,t i ) with i ∈ {hard, rec, PS}. Then we compute the three different K factors K i from these three sets of variables. Taking K hard as a reference we compute the ratio K i /K hard for each event. If the reconstruction would be insensitive to the choice of kinematics this ratio ought to be unity all the time.
In Fig. 12 we show a histogram of this ratio for the two different ratios K rec /K hard and K PS /K hard . As expected, the results from the reconstruction after the application of the parton shower are quite different from unity and show a very broad peak. This is reasonable and shows us that a wrong choice of kinematic variables gives a sensitive difference in the computation of the K factor.
In strong contrast, the EW K factors after the kinematics reconstruction from the final state compared to the K factor from the variables of the hard process inside the event generator are very close and their ratio shows a very strong peak around unity, showing us that the two prescriptions lead to nearly identical results.
In the same figure 12 we compare the results from runs with and without NLO QCD corrections applied. In both cases we get very similar results for the EW K factor. The additional hard gluon from the real emission graph distorts the kinematics that was reconstructed in the leading order case only slightly. As expected the QCD corrections give a slightly bigger difference between the two reconstruction schemes.
We conclude that our computation of the EW correction is very robust against small variations of the kinematics as long as the variables are sensibly chosen and a sensible veto is applied. The reconstruction of kinematics from the final state alone, as proposed in the previous section, is a viable choice. One should, however, be careful, as e.g. naively ignoring the parton shower in the reconstruction of kinematics will lead to significantly different results.

HERWIG++ results
Having established the role of the veto on the leptonic final state we present results of differential distributions in leptonic observables. We apply the same selection of final states as described above in Sect. 2.3 and Sect. 4. Now, the leptons are, however, selected from a full hadronic final state with an additional isolation cut of R = 0.2. Furthermore, the final state is of course modified by parton showers and hadronization as well as additional soft and collinear photon radiation. We left the underlying event switched off as we expect only a small effect for the observables presented here. In all cases we have generated 10M unvetoed events. The leptonic veto has only been applied at the analysis level.
In Fig. 13 we show a number of observables for the final state of ZZ production. In all cases we show four lines: the leading order result (LO), results with only electroweak (NLO EW) or QCD corrections (NLO QCD) applied and, finally, the result with the combined EW and QCD corrections in the multiplicative scheme as outlined above. In both cases, with and without NLO QCD corrections, the additional EW correction is as sizable as in the partonic case, i.e. of the order of −5 % for small m 4l , reaching up to −20 % for m 4l close to 1 TeV or p T,ll close to 500 GeV. The rapidity distributions, shown in the lower two plots of Fig. 13, receive the correction of −4% typical for the low-ŝ configuration.
A similar picture emerges in the case of WZ production which we show in Fig. 14 for √ s = 13 TeV. Here, the EW corrections are smaller than in the ZZ case. In every observable we find that the EW corrections act quite similar on the final states with and without NLO QCD corrections. While the QCD corrections are moderate, at the level of 20 %, the EW corrections vary between zero and −20 %, and are again sizable for high transverse momenta of the Z bosons.
Finally, in Fig. 15 we consider selected observables for the case of W-pair production at 13 TeV. Here, the EW corrections are slightly larger than in the WZ case but the overall picture remains the same. The QCD corrections are quite large but tamed by our veto on the leptonic final state, typically of the order of 20 %. The EW corrections are typically of the order of 5 % but again sizable in the case of large lepton transverse momentum, where they can completely compensate the QCD corrections.

Summary and Conclusions
We have computed the full NLO EW corrections to resonant vector-boson pair production at the LHC, taking into account leptonic decays and corresponding spin correlations. We propose a simple and straight-forward method -relying on unpolarized 2 → 2 K-factors -to implement our results in any state-of-the-art MC generator. The EW corrections are combined with precise QCD predictions in a multiplicative approach. We emphasize that our method also allows for an a posteriori implementation of EW K-factors into MC samples that have already been generated. To demonstrate the practicability of our approach, we have included our corrections in the HERWIG++ MC generator and presented various distributions for four-lepton production at the LHC obtained in the HERWIG++ setup, including EW corrections, NLO QCD corrections matched to parton showers, as well as hadronization effects. In the future our method could also be applied to other process classes, such as V +jet production at the LHC, allowing for phenomenological studies that combine EW precision on the one hand and an adequate treatment of dominating QCD effects on the other.

A Endpoint contributions from dipoles
According to Ref. [28] the endpoint contribution compensating soft and collinear QED singularities from the virtual corrections of the unpolarized V -boson pair production process qq → V 1 V 2 is given by the general expression where the sum runs over all dipole contributions (IJ) of charged external particles, i.e. I, J = q,q , W ± , and the universal (i.e. process-independent) functions G (sub) IJ (P IJ , m I , m J ) carry the respective endpoint singularities. The charge flow of the external particle J is denoted by σ J . Let us again state explicitly that the sum of the one-loop virtual corrections and the endpoint contribution, σ V qq →V 1 V 2 + σ E qq →V 1 V 2 γ , is IR finite by construction and defined on a LO phase space solely parametrized byŝ andt. In our computation the quark masses are neglected whenever possible and only introduced to regularize collinear singularities, whereas a finite W mass is kept throughout the calculation. Therefore, the dipole formulas for massive FS particles and massless IS particles have to be applied here.
In case of a massless IS emitter a and a massless IS spectator b we find P ab = √ŝ and with the auxiliary function L(P 2 , m 2 ) = ln m 2 P 2 ln λ 2 P 2 + ln Here, soft singularities are regularized by an infinitesimal photon mass λ, and a small quark mass m a is kept to regularize collinear singularities, as mentioned before. For the case of a massive FS emitter and a massive FS spectator, which contributes to W-pair production, G (sub) ij (P 2 ij , m i , m j ) is given by Eq. (4.10) from Ref. [28], with P 2 ij =ŝ. The endpoint contributions G (sub) ai (P 2 ia , m a , m i ) and G (sub) ia (P 2 ia , m i ) for massless IS emitters (spectators) and massive FS spectators (emitters) can be found in Eq. (A.4) of the respective paper, where m a has been set to zero in the expression for G (sub) ia (P 2 ia , m i ). Here, the auxiliary momentum is given by Note that the finite terms (i.e. terms not proportional ln(λ) or ln(m q )) of the endpoint contributions in Eq. (A.1) do not have any physical meaning, since they are tailored to cancel certain subtraction contributions from the real-radiation processes, which are neglected anyway in the V+E approximation. Nevertheless, our approach is still justified since in the V+E approximation applied in this paper we do not claim to control the QED part of the EW corrections, and, more importantly, we have explicitly shown in Sect. 2.3 that this approximation works remarkably well in V -pair production processes. However, we clearly point out that this might be different for different process classes, and the validity of the V+E approximation has to be carefully checked if applied to other processes.    Figure 15: Results from WW events at √ s = 13 TeV. We show observables that can be reconstructed from the visible leptonic final state after we applied the veto Eq. (4.4).