SFB/CPP1401
TTP14001
LPN14001
KATP0114
MCnet1401
Vectorboson pair production and electroweak corrections in HERWIG++
Stefan Gieseke, Tobias Kasprzik, Johann H. Kühn
[2ex] Karlsruhe Institute of Technology (KIT), Institut für Theoretische Physik (IThP)
D76131 Karlsruhe, Germany
[1ex] Karlsruhe Institute of Technology (KIT), Institut für Theoretische Teilchenphysik (TTP)
D76131 Karlsruhe, Germany
[2ex]
Abstract: The detailed study of vectorboson pair production processes at the LHC will lead to a better understanding of electroweak physics. As pointed out before, a consistent inclusion of higherorder 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 invariantmass and angular distributions at high energies, but may also significantly affect the cross section near threshold, as is the case e.g. for Zpairs. For this reason, we present an analysis of the nexttoleadingorder 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 fixedorder results, a simple and straightforward method is motivated to combine the electroweak corrections with stateoftheart Monte Carlo predictions, focusing on a meaningful combination of higherorder 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.
January 2014
1 Introduction
Vectorboson pair production processes play a central role in LHC phenomenology. These processes are not only of great importance with respect to background analyses in StandardModel (SM) Higgs production, they will also provide deeper insight into the structure of the electroweak interaction at highest energies. This is particularly true for the future highluminosity run of the LHC at a centerofmass (CM) energy of 13 TeV, which will allow for an unprecedented accuracy in the analysis of vectorboson interactions at the TeV scale. Consequently, theoretical prediction with high accuracy are needed for this important class of processes.
Given the experimental accuracy already achieved by LHC experiments in the 7 and 8TeV runs [1, 2, 3, 4, 5, 6], at least the nexttoleading order (NLO) QCD corrections are mandatory for a robust prediction of boson pair production processes (for selected references see Refs. [7, 8]). In addition, even first steps towards NNLO predictions of massive vectorboson pair production processes have been made [9, 10]. (As far as photonpair and Z production at the LHC are concerned, the full strong twoloop corrections are known even fully differentially [11, 12, 13].) In particular, approximate NNLO results for and WW production have been provided for hightransversemomentum observables [14, 15], as well as for WW production in the threshold limit [16].
However, at energies of a few hundred GeV, also electroweak (EW) corrections (and other related electroweak effects) are becoming more and more important, and a lot of activity has taken place also in this field. In particular, the interplay of EW corrections and anomalous couplings has been investigated in Ref. [17]. The corresponding EW corrections have been computed in Ref. [18] in the highenergy limit, including leptonic decays and offshell effects. Recently, also the full EW corrections to Wpair production, also taking into account mass effects, have been evaluated for the leptonic final state [19]. Leading twoloop effects at high transverse momenta were evaluated in Ref. [20] for Wpairs. A detailed analysis of onshell boson pair production () and production including EW corrections has been provided in Refs. [7, 8], consistently including all mass effects. Recently, a detailed review of NLO effects in pair production of massive bosons has been provided, emphasizing the importance of photoninduced contributions [21].
Expecting first results for the full NNLO QCD corrections to Wpair production in the near future [22], a natural next step would be the combination of EW and QCD predictions at accuracy on a consistent theory basis, as has been partially done for the Drell–Yan process [23] already. However, these multiscale twoloop calculations are beyond feasibility at present. Nevertheless, at least a pragmatic prescription of combining QCD predictions with EW corrections is desirable, aiming for a combination of EW precision with standard Monte Carlo (MC) tools.
In this work we extend the above results in two ways. In Sect. 2 a first study of EW corrections to massive vectorboson pair production processes is presented, including mass effects as well as leptonic decays to allow for a realistic event definition in the leptonic decay modes. We point out that our results are not restricted to the highenergy regime but are also valid at moderate energies of GeV, where already corrections of 5–10% may be observed. Applying these results, in a second step we propose a straightforward and simple strategy for the implementation of EW corrections in any stateoftheart MC generator (Sect. 4), relying on factors for unpolarized twobytwo 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.
2 Electroweak corrections
As stated in the introduction, the complete evaluation of the combined QCD and electroweak corrections to gaugeboson pair production of order 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 twobody collisions, described by effective (squared) partonic energies and energy transfers and . 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 oneloop 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 gaugeboson 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 factor which is given as a function of and 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 onshell 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 onshell WW, WZ and ZZ production at the LHC.
2.1 Onshell gaugeboson pair productions
Here, we summarize the combination of different electroweak effects in onshell pair production which have been discussed in detail in Refs. [7, 8, 21]. For the numerical implementation, we use the default setup defined in Refs. [7, 8].
In Wpair production, the invariantmass distribution (Fig. 1 (top)) receives wellknown logarithmically enhanced negative EW corrections () growing with energy. Positive contributions arise from the partonic subprocess () and the photonquark induced processes (), the latter evaluated applying the jet veto defined in Ref. [7]. As already pointed out in Ref. [8], the effect of massiveboson radiation () is moderate, however, strongly dependent on the event selection.
The above picture significantly changes if angular distributions of the Wpair 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 GeV. While the genuine EW corrections drastically reduce the differential cross section at high , corresponding to small rapidity gap, the photoninduced 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 photoninduced 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 Wpair cross section. A relative error of on the leadingorder (LO) cross section at 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 Wpair 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 Zpair production, however, the induced contributions are negligible [21], and the realradiation contributions are always below 10%. In total, particularly large negative corrections, reaching %, dramatically affect Zpair production at high invariant masses and transverse momenta (Fig. 3).
2.2 factors for the electroweak corrections
Let us start with the simplest reaction, inclusive Zboson pair production. The full set of electroweak corrections, including photon radiation, Zboson mass effects and virtual top quarks, has been discussed in Ref. [8], where the purely weak corrections were also evaluated for the 4lepton 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 Zboson 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 induced reaction which we denote as
(2.1) 
the others, and for the and induced reactions. All these factors depend on the partonic and only.
Similar considerations apply to diphoton 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 twobody approximation reproduces the full answer to better than 1 percent. The correction factors , and 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 Wpair 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 factor which depends on and only is nontrivial.
In order to return to the kinematics of a twototwo 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 and .) At the same time the kinematics of the resulting final states is now fully described by the partonic and . In this way we obtain again three correction factors , and which can be employed in the framework of the Monte Carlo generator, just as before. A similar approach is valid for WZ and WZ production, which have, of course identical correction factors, denoted . The corresponding endpoint contributions are also listed in the Appendix.
2.3 Gaugeboson polarization and fourlepton 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 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 region, separated according to longitudinal and transverse polarizations.
From these considerations it becomes clear that in case of Zpair production a single partonic 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 Wpair (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 transverselypolarized 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 factors to reproduce the full corresponding EW corrections with sufficient accuracy.
The corrections evaluated in Refs. [8, 7] and encoded in our factors were obtained for
onshell Z or W bosons. Any realistic simulation of fourfermion
production must, necessarily, include contributions from offshell
configurations. In the case of Z also digrams with offshell 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 , to suppress the admixture of virtual photons
and enhance closetomassshell gauge bosons. In this case the
neglection of virtual photons as implemented in HERWIG++ can be
justified
In our approach offshell effects and nonresonant 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 offshell effects in the EW corrections only amount to , as demonstrated in Ref. [19] for Wpair production.
From these considerations it can be expected that the dominant weak corrections to fourlepton production at the LHC are well described by process dependent factors which can be taken from the unpolarized results for the corresponding 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 gaugeboson 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 gaugeboson pair production as defined in Ref.[29]. Specifically, in the leptonic event selection we apply the following basic cuts
(2.2) 
for the chargedlepton transverse momenta and rapidities. If neutrinos are present in the final state, a minimal missing transverse momentum of
(2.3) 
is also required. The invariant mass of the lepton pair is restricted to
(2.4) 
to suppress nonresonant backgrounds. For the gauge boson widths we use the values
(2.5) 
and the weak coupling constant is defined in the 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 nonresonant and offshell 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 ComplexMass Scheme (CMS) [32] is at the permill level and hardly visible. All results presented here for the full LO cross sections therefore correspond to the naive fixedwidth implementation.
In addition to the full results we provide the results for Vboson pair production in the doublepole approximation (DPA) originally discussed in Ref. [18]. Here, the amplitudes for Vpair production and decays are evaluated onshell, but the Breit–Wigner shape of the resonance is included in the evaluation of the squared matrix elements to account for the dominant offshell effects. We apply the onshell projection procedure proposed in Ref. [33] to construct proper onshell momenta of the intermediate bosons from the fourparticle phasespace. Note that in addition to the physical cuts displayed above we impose a technical cut, , on the 4lepton invariant mass since the onshell projection suggested in Ref. [33] only gives sensible results above threshold.
As a third alternative, we work in the narrowwidth approximation (NWA), where the gauge bosons are strictly forced onshell from the beginning via the replacement
(2.6) 
for the resonant squared propagators.
In Figs. 5, 6 and 7 we present results for WW, WZ 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
(2.7a)  
(2.7b)  
(2.7c) 
Note that if intermediate Z bosons are present in the process, the contributions and all related interference contributions are taken into account in the full LO results, while those contributions are absent in DPA and NWA, respectively.
In the left panels of Figs. 5, 6 and 7 we present LO results for various typical differential distributions for processes (2.7) at LHC13, resepctively. Besides the full results, the respective approximate results in NWA and DPA are also shown, always taken in leading order.
The righthandside panels of the respective plots show the relative deviations of the NWA () and DPA () from the full results. For all pairproduction 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 Zpair production threshold where offshell 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 vectorboson 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 wellestablished FeynArts/FormCalc/LoopTools [34, 36, 35, 37, 38] setup already used in Refs.[8, 7], 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 finalstate photon radiation (FSR), which leads to significant distortions of the phasespace 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 nonfactorizable 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 nonfactorizable corrections largely cancel in sufficiently inclusive observables [42].
In addition to the full EW corrections (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 in the virtual+endpoint (V+E) approximation are presented in the same plots, where unpolarized onshell 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 factors in the V+E approximation to describe vectorboson pair production at the LHC at a sufficient accuracy.
3 Electroweak and QCD corrections – combined
As stated in the introduction, the complete treatment of combined QCD and electroweak corrections, involving e.g. twoloop terms of order is presently out of reach. Consequently neither an additive treatment of the corrections, , nor a multiplicative one, , 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 quarkantiquark 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 and dependent 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 reaction in the 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 noncollinear radiation, which leads to a diboson system of large transverse momentum. In this case only an approximate prescription can be formulated, since the original 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 and from the distorted kinematics for the evaluation of : The squared CM energy is calculated from the fourlepton final state via
(3.1) 
The momenta are boosted into the fourlepton CM frame (denoted by ). In this frame the unit directions of initialstate hadrons shall be denoted by
(3.2) 
The direction of the effective scattering axis in is now defined by
(3.3) 
and the effective scattering angle is, correspondingly, given by
(3.4) 
where the denotes the momentum direction of vector boson . The Mandelstam variable is then computed from assuming onshell 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 , 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 pair production is dominated by new topologies which are absent at lowest order in QCD. These topologies correspond to production with the radiation of an additional from the quark jet rather than being a correction to pair production, and spoil the perturbative series for the prediction of leptonic observables at high transverse momenta. They lead to huge QCD 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 observables, employing suitable restrictions on the transverse momentum of the 4lepton 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 ( vs ) 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 and 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.
4 Implementation into Herwig++
Our starting point of the implementation of the EW corrections in HERWIG++ is the electroweak 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 vectorboson pair production HERWIG++ delivers unweighted events, hence we can compute and use this directly as a reweighting factor for each event such that leads to a suppression of events with given while 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 are known.
Let us discuss the kinematics of our events in more detail. For EW corrections to the Born process the calculation of a factor is straightforward as we can directly access the kinematic setup of the hard process once this is generated. In this case 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 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 [44, 29]. (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
(4.1) 
Here, and denote the Born and radiative phase space variables, respectively, and is the minimum transverse momentum that is generated by the parton shower. only parametrizes the additional variables to specify a hard emission relative to the Born configuration. and are the Born and realemission matrix elements squared for the hard process under consideration. is the Born differential cross section while
(4.2) 
Here, is the (infrared finite) sum of virtual corrections and the divergent part of the real corrections, while is the (also finite) real correction matrix element squared with the divergent terms subtracted. Finally, the POWHEG Sudakov form factor is given by
(4.3) 
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 with weight . The first term gives the noemission probability, the second the contribution from a single emission. At the same time, upon expansion in , 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 . 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 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 nexttoleading 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 in the hard subprocess and hence we may easily compute the kinematic variables for every event. In effect, this assumes that the EW corrections are applied only on the level of Borntype 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 gaugeboson pairs are accompanied by additional hard quark or gluon jets, leaving the 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 gaugeboson 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 in the laboratory frame,
(4.4) 
Here, we consider all visible leptons from EW boson decays, i.e. runs to 2, 3 or 4 in case of WW, WZ or events. We additionally assume that the missing transverse momentum in the event solely stems from neutrinos from W decays. The lefthandside 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 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 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. 911. Let us compare LO and NLO predictions for the distribution in ZZ production (Fig. 9, left) which corresponds essentially to the transversemomentum distribution of the Zboson. Without cut the NLO distribution exceeds the LO distribution by a factor of two at GeV and more at larger transverse momenta. Introducing the cut (4.4) with removes most of this excess such that the difference between LO and NLO distributions is reduced to and remains relatively constant as function of . The ratio between the NLO and LO rapidity distribution, in contrast, is fairly constant with a ratio of NLO/LO , and is reduced to by the cut for . A similar behaviour is observed for WZ (Fig. 10) in the 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 by introducing the cut (4.4), and a similar behaviour is observed for WW production in the 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 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 . As an alternative, we reconstruct the kinematics from the vectorboson final state as outlined in the previous section. Here, no knowledge of the initialstate partons is needed for the computation of the kinematic invariants. is given by the invariant mass of the vectorboson pair and 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 due to the parton showering. We then find the four momenta of the outgoing vector bosons and compute 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 with . Then we compute the three different factors from these three sets of variables. Taking as a reference we compute the ratio 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 and . 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 factor.
In strong contrast, the EW factors after the kinematics reconstruction from the final state compared to the 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 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.
5 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 . 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 M 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 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 for small , reaching up to for close to 1 TeV or close to 500 GeV. The rapidity distributions, shown in the lower two plots of Fig. 13, receive the correction of typical for the low configuration.
A similar picture emerges in the case of WZ production which we show in Fig. 14 for 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 , and are again sizable for high transverse momenta of the Z bosons.
Finally, in Fig. 15 we consider selected observables for the case of Wpair 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.
6 Summary and Conclusions
We have computed the full NLO EW corrections to resonant vectorboson pair production at the LHC, taking into account leptonic decays and corresponding spin correlations. We propose a simple and straightforward method – relying on unpolarized factors – to implement our results in any stateoftheart 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 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 fourlepton 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 +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 boson pair production process is given by the general expression
(A.1) 
where the sum runs over all dipole contributions of charged external particles, i.e. , and the universal (i.e. processindependent) functions carry the respective endpoint singularities. The charge flow of the external particle is denoted by . Let us again state explicitly that the sum of the oneloop virtual corrections and the endpoint contribution, , is IR finite by construction and defined on a LO phase space solely parametrized by and . 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 and a massless IS spectator we find and
(A.2) 
with the auxiliary function
(A.3) 
Here, soft singularities are regularized by an infinitesimal photon mass , and a small quark mass 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 Wpair production, is given by Eq. (4.10) from Ref. [28], with . The endpoint contributions and for massless IS emitters (spectators) and massive FS spectators (emitters) can be found in Eq. (A.4) of the respective paper, where has been set to zero in the expression for . Here, the auxiliary momentum is given by .
Note that the finite terms (i.e. terms not proportional or ) 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 realradiation 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 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.
Acknowledgements
This work was supported by the DFG Sonderforschungsbereich/Transregio 9 “Computergestützte Theoretische Teilchenphysik”. SG acknowledges support from the EU Initial Training Network “MCnetITN” and the Helmholtz Alliance “Physics at the Terascale”.
Footnotes
 The HERWIG++ implementation of pair production [29] relies on the doublepole approximation, where only doublyresonant contributions are taken into account including offshell effects.
References
 G. Aad et al. [ATLAS Collaboration], JHEP 1303 (2013) 128 [arXiv:1211.6096 [hepex]].
 G. Aad et al. [ATLAS Collaboration], Phys. Rev. D 87 (2013) 112001 [arXiv:1210.2979 [hepex]].
 G. Aad et al. [ATLAS Collaboration], Eur. Phys. J. C 72 (2012) 2173 [arXiv:1208.1390 [hepex]].
 S. Chatrchyan et al. [CMS Collaboration], arXiv:1306.1126 [hepex].
 S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 721 (2013) 190 [arXiv:1301.4698 [hepex]].
 S. Chatrchyan et al. [CMS Collaboration], JHEP 1301 (2013) 063 [arXiv:1211.4890 [hepex]].
 A. Bierweiler, T. Kasprzik, H. Kühn and S. Uccirati, JHEP 1211 (2012) 093 [arXiv:1208.3147 [hepph]].
 A. Bierweiler, T. Kasprzik and J. H. KÃ¼hn, arXiv:1305.5402 [hepph].
 T. Gehrmann, L. Tancredi and E. Weihs, JHEP 1308 (2013) 070 [arXiv:1306.6344 [hepph]].
 G. Chachamis, M. Czakon and D. Eiras, JHEP 0812 (2008) 003 [arXiv:0802.4028 [hepph]].
 Z. Bern, A. De Freitas and L. J. Dixon, JHEP 0109 (2001) 037 [hepph/0109078].
 S. Catani, L. Cieri, D. de Florian, G. Ferrera and M. Grazzini, Phys. Rev. Lett. 108 (2012) 072001 [arXiv:1110.2375 [hepph]].
 M. Grazzini, S. Kallweit, D. Rathlev and A. Torre, arXiv:1309.7000 [hepph].
 F. Campanario and S. Sapeta, Phys. Lett. B 718 (2012) 100 [arXiv:1209.4595 [hepph]].
 F. Campanario, M. Rauch and S. Sapeta, arXiv:1309.7293 [hepph].
 S. Dawson, I. M. Lewis and M. Zeng, arXiv:1307.3249 [hepph].
 E. Accomando and A. Kaiser, Phys. Rev. D 73 (2006) 093006 [hepph/0511088].
 E. Accomando, A. Denner and A. Kaiser, Nucl. Phys. B 706 (2005) 325 [hepph/0409247].
 M. Billoni, S. Dittmaier, B. Jäger and C. Speckner, arXiv:1310.1564 [hepph].
 J. H. Kühn, F. Metzler, A. A. Penin and S. Uccirati, JHEP 1106 (2011) 143 [arXiv:1101.2563 [hepph]].
 J. Baglio, L. D. Ninh and M. M. Weber, arXiv:1307.4331 [hepph].
 G. Chachamis, arXiv:1312.4220 [hepph].
 R. Boughezal, Y. Li and F. Petriello, arXiv:1312.3972 [hepph]; W. B. Kilgore, C. Sturm, [arXiv:1107.4798 [hepph]]; A. Kotikov, J. H. Kühn, O. Veretin, Nucl. Phys. B788 (2008) 4762. [hepph/0703013 [HEPPH]]; Q. H. Cao and C. P. Yuan, Phys. Rev. Lett. 93 (2004) 042001 [hepph/0401026]; B. F. L. Ward, C. Glosser, S. Jadach and S. A. Yost, Int. J. Mod. Phys. A 20 (2005) 3735 [hepph/0411047]; B. F. L. Ward and S. A. Yost, Acta Phys. Polon. B 38 (2007) 2395 [arXiv:0704.0294 [hepph]]; A. Vicini et al., PoS RADCOR2007 (2007) 013; G. Balossini, G. Montagna, C. M. Carloni Calame, M. Moretti, O. Nicrosini, F. Piccinini, M. Treccani, A. Vicini, JHEP 1001 (2010) 013. [arXiv:0907.0276 [hepph]]; G. Balossini et al., Nuovo Cim. 123B (2008) 741.
 M. Bähr, S. Gieseke, M. A. Gigg, D. Grellscheid, K. Hamilton, O. LatundeDada, S. Plätzer and P. Richardson et al., Eur. Phys. J. C 58 (2008) 639 [arXiv:0803.0883 [hepph]].
 K. Arnold, L. d’Errico, S. Gieseke, D. Grellscheid, K. Hamilton, A. Papaefstathiou, S. Plätzer and P. Richardson et al., arXiv:1205.4902 [hepph].
 A. D. Martin, R. G. Roberts, W. J. Stirling, R. S. Thorne, Eur. Phys. J. C39 (2005) 155161. [hepph/0411040].
 R. D. Ball et al. [ The NNPDF Collaboration], arXiv:1308.0598 [hepph].
 S. Dittmaier, Nucl. Phys. B 565 (2000) 69 [hepph/9904440].
 K. Hamilton, JHEP 1101 (2011) 009 [arXiv:1009.5391 [hepph]].
 H. L. Lai, M. Guzzi, J. Huston, Z. Li, P. M. Nadolsky, J. Pumplin and C. P. Yuan, Phys. Rev. D 82 (2010) 074024 [arXiv:1007.2241 [hepph]].
 M. R. Whalley, D. Bourilkov and R. C. Group, in HERA and the LHC, eds. A. de Roeck and H. Jung (CERN2005014, Geneva, 2005), p. 575, hepph/0508110.
 A. Denner, S. Dittmaier, M. Roth and L. H. Wieders, Nucl. Phys. B 724 (2005) 247 [hepph/0505042]; A. Denner and S. Dittmaier, Nucl. Phys. Proc. Suppl. 160 (2006) 22 [hepph/0605312].
 A. Denner, S. Dittmaier, M. Roth and D. Wackeroth, Nucl. Phys. B 587 (2000) 67 [hepph/0006307].

J. Küblbeck, M. Böhm and A. Denner,
Comput. Phys. Commun. 60 (1990) 165;
H. Eck and J. Küblbeck, Guide to FeynArts 1.0, University of Würzburg, 1992.  T. Hahn, Comput. Phys. Commun. 140 (2001) 418 [hepph/0012260].
 T. Hahn and M. PérezVictoria, Comput. Phys. Commun. 118 (1999) 153 [hepph/9807565].
 T. Hahn, C. Schappacher, Comput. Phys. Commun. 143 (2002) 5468. [hepph/0105349].
 G. J. van Oldenborgh, J. A. M. Vermaseren, Z. Phys. C46 (1990) 425438.
 J. Alwall et al., JHEP 0709 (2007) 028. [arXiv:0706.2334 [hepph]].
 K. Hamilton and P. Richardson, JHEP 0607 (2006) 010 [hepph/0603034].
 D. R. Yennie, S. C. Frautschi and H. Suura, Annals Phys. 13 (1961) 379.
 P. Falgari, A. S. Papanastasiou and A. Signer, JHEP 1305 (2013) 156 [arXiv:1303.5299 [hepph]].
 M. Rubin, G. P. Salam and S. Sapeta, JHEP 1009 (2010) 084 [arXiv:1006.2144 [hepph]].
 S. Frixione, P. Nason and C. Oleari, JHEP 0711 (2007) 070 [arXiv:0709.2092 [hepph]].
 A. Buckley, J. Butterworth, S. Gieseke, D. Grellscheid, S. Höche, H. Hoeth, F. Krauss and L. Lönnblad et al., Phys. Rept. 504 (2011) 145 [arXiv:1101.2599 [hepph]].
 S. Gieseke, Prog. Part. Nucl. Phys. 72 (2013) 155.