NLO QCD+EW predictions for \boldsymbol{HV} and \boldsymbol{HV}+jet production including parton-shower effects

# NLO QCD+EW predictions for HV and HV+jet production including parton-shower effects

F. Granata,    J. M. Lindert,    C. Oleari,
###### Abstract

We present the first NLO QCD+EW predictions for Higgs boson production in association with a or pair plus zero or one jets at the LHC. Fixed-order NLO QCD+EW calculations are combined with a QCD+QED parton shower using the recently developed resonance-aware method in the POWHEG framework. Moreover, applying the improved MiNLO technique to and production at NLO QCD+EW, we obtain predictions that are NLO accurate for observables with both zero or one resolved jet. This approach permits also to capture higher-order effects associated with the interplay of EW corrections and QCD radiation. The behavior of EW corrections is studied for various kinematic distributions, relevant for experimental analyses of Higgsstrahlung processes at the 13 TeV LHC. Exact NLO EW corrections are complemented with approximate analytic formulae that account for the leading and next-to-leading Sudakov logarithms in the high-energy regime. In the tails of transverse-momentum distributions, relevant for analyses in the boosted Higgs regime, the Sudakov approximation works well, and NLO EW effects can largely exceed the ten percent level. Our predictions are based on the POWHEG BOX RES+OpenLoops framework in combination with the Pythia 8.1 parton shower.

###### Keywords:
\preprint

IPPP/17/46

ZU-TH 14/17

Università di Milano-Bicocca and INFN, Sez. di Milano-Bicocca, Piazza della Scienza 3, 20126 Milano, Italy Institute for Particle Physics Phenomenology, Durham University, South Rd, Durham DH1 3LE, UKUniversità di Milano-Bicocca and INFN, Sez. di Milano-Bicocca, Piazza della Scienza 3, 20126 Milano, Italy Physik-Institut, Universität Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland

## 1 Introduction

The discovery of the Higgs boson [1, 2] has opened the door to the direct experimental investigation of the Higgs and Yukawa sectors of the Standard Model. While present measurements of Higgs boson properties and interactions are consistent with the Standard Model [3], the full set of data collected during Run II and in subsequent runs of the LHC will provide more and more stringent tests of the mechanism of electroweak symmetry breaking.

In this context, the associated production of a Higgs and a vector boson, with and , plays a prominent role. In spite of the fact that the total cross sections for these so-called Higgsstrahlung processes are subleading as compared to Higgs boson production via gluon fusion and vector-boson fusion, the possibility to reconstruct the full final state and the clean signatures that result from leptonically decaying vector bosons offer unique opportunities of testing Higgs boson interactions with vector bosons and heavy quarks (see Refs. [4, 5, 6] and references therein). The associated production makes it possible to disentangle Higgs boson couplings to and bosons from one another and to measure them in a broad kinematic range. In addition, the presence of the associated vector boson allows for an efficient suppression of QCD backgrounds. In particular, is the most favorable channel for measurements of the branching ratio, and thus for determinations of the bottom Yukawa coupling. In production with decay, the boosted region, with Higgs boson transverse momentum above 200 GeV, plays a particularly important role, both in order to achieve an improved control of the QCD backgrounds [4] and for the sensitivity to possible anomalies in the couplings. Higgsstrahlung processes permit also to probe invisible Higgs boson decays, both through direct measurements of with invisible Higgs decays and through indirect bounds based on measurements of the branching ratio.

The accuracy of present and future measurements of production, at the level of both fiducial cross sections and differential distributions, calls for increasingly accurate theoretical predictions. The inclusion of higher-order QCD corrections is crucial, both for total rates and for a precise description of the QCD radiation that accompanies the production of the system. The role of QCD corrections can be particularly important in the boosted regime or in the presence of cuts and for observables that are sensitive to QCD radiation.

In general, in order to account for experimental cuts and observables, higher-order QCD and EW predictions should be available for arbitrary differential distributions, and experimental analyses require particle-level Monte Carlo generators where state-of-the-art theoretical calculations are matched to parton showers. Finally, when QCD and EW higher-order effects are both sizable, also their combination needs to be addressed.

Theoretical calculations for the associated-production processes are widely available in the literature. Among the numerous studies on production at next-to-leading order (NLO) QCD we quote here Refs. [7, 8, 9]. Predictions for inclusive and production at next-to-next-to-leading order (NNLO) QCD have first been obtained in Refs. [10, 11] and are implemented in the VH@NNLO program [12]. Besides contributions of Drell–Yan (DY) type, where the Higgs boson results from an -channel subtopology, Higgsstrahlung at NNLO QCD involves also extra contributions where the Higgs boson couples to heavy-quark loops. Such non-DY contributions arise via squared one-loop amplitudes in the channel [13] and through the interference of one-loop and tree amplitudes in the and crossing-related channels. Studies of possible anomalous coupling in the channel can be found in Ref. [14, 15]. Heavy-quark loop contributions to the channel are known up to in the limit of the mass of the bottom quark going to zero, and the mass of the top quark going to infinity [16]. Their impact, especially in the boosted regime, can be quite significant [17].

Fully differential NNLO calculations for production with off-shell vector-boson decays were first presented in Refs. [18, 19, 20], including all DY contributions plus heavy-quark-loop contributions to . More recently, a NNLO QCD calculation that includes also the small heavy-quark loop contributions in the channel and in the crossing-related and channels became available [21] and also production with NNLO QCD corrections both in the production and in the decay part of the process [22].

Analytic resummations have been discussed in Refs. [23, 24, 25, 26], while leading-logarithmic resummation can be routinely achieved through the matching of NLO QCD calculations to parton showers (PS). The first NLO+PS generators in the MC@NLO [27] and POWHEG frameworks [28, 29, 30] have been presented in Refs. [31] and [32], respectively. More recently, new generators that provide an NLO accurate description of and + jet production became available. The first generator of this kind was presented in Ref. [33] based on the MiNLO method [34, 35], while a simulation of  + 0 and 1 jet, based on the MEPS@NLO multijet merging technique [36, 37], was presented in Ref. [38]. Concerning fermion loops, the POWHEG BOX generator of Ref. [33] can account for all NLO contributions of DY and non-DY type to  + jet and also for the finite loop-induced contributions, with the possibility of studying anomalous couplings in the “kappa framework”. A more general study, which uses an effective field theory approach and introduces generic six-dimensional operators, can be found in Ref. [39].

The heavy-quark loop-mediated production was first studied in Ref. [40]. More recently, the Sherpa generator of Ref. [38] has included also NNLO-type squared quark-loop contributions in the , , and plus crossing-related channels. Lately, a NNLO+PS generator for  [41] that combines the NNLO QCD calculation of Ref. [18] with the parton shower using the method of Refs. [35, 42] was presented.

Electroweak corrections to , including off-shell - and -boson decays, are known at NLO [43, 44] and are implemented in the parton-level Monte Carlo program HAWK [45]. These corrections are at the level of 5% for inclusive quantities, but in the high-energy regime they can reach various tens of percent due to the presence of Sudakov logarithms [46, 47, 48, 49, 50, 51, 52, 53]. For this reason, especially in boosted searches, the inclusion of EW corrections is mandatory. An interesting aspect of these corrections in production is that they induce also a dependence on the Higgs sector, and in particular on the trilinear coupling . Thus, precise measurements of Higgsstrahlung processes can be exploited for setting limits on  [54, 55, 56, 57]. To date, none of the existing NLO+PS generators implement EW corrections.

In this paper, for the first time, we present NLO QCD and NLO EW calculations for the production of a Higgs boson in conjunction with a or leptonic pair, plus zero or one jet, at the LHC. While, for convenience, the above-mentioned processes will often be denoted as / production (with and ) in the rest of the paper, all the results we are going to present always correspond to the complete decayed final-state processes, with spin effects, off-shell and non-resonant contributions taken into account.

At NLO QCD we include the full set of contributions to and contributions to . Although terms of non-DY type are implemented in our codes, we have not included them in our simulations. In addition, we do not include NNLO-like loop-induced contributions to plus 0 and 1 jet production.

Besides showing fixed-order NLO QCD+EW predictions at parton level for typical observables, we also present full NLO+PS simulations for and production. To this end, we have implemented our NLO calculations for and production into four separate codes (, , and ) in the POWHEG BOX framework. In this way, we have consistently combined the radiation emitted at NLO QCD+EW level with a QCD+QED parton shower. In this context, photon radiation from the charged leptons can lead to severe unphysical distortions of the - and -boson line shapes, if not properly treated. This problem was first pointed out in the context of NLO QCD+PS simulations of off-shell top-quark production and decay, and was solved in the context of the POWHEG BOX framework by means of the so-called resonance-aware method [58]. The first application of this method and its variants, in the context of electroweak corrections, has appeared in Refs. [59, 60]. In this paper, we exploit the flexibility of the resonance-aware method to perform a fully consistent NLO QCD+EW matching in the presence of non-trivial EW resonances. To this end, our NLO calculations and generators are implemented in the new version of the POWHEG BOX framework, known as POWHEG BOX RES. In this recent version, the hardest radiation generated by POWHEG preserves the resonance virtualities present at the underlying-Born level. At the same time, the resonance information can be passed on to the parton shower, which in turn preserves the virtualities of intermediate resonances of the hard process in subsequent emissions.

Similarly to what was done in Ref. [33] for production at NLO QCD, we have applied the improved MiNLO [34, 35] approach to production in order to get a sample of events that has simultaneously NLO QCD accuracy for plus 0 and 1 jet. In the MiNLO framework, also the NLO EW corrections to and production have been consistently combined in the same inclusive sample. This can be regarded as an approximate treatment of corrections in observables that are very sensitive to QCD radiation and receive, at the same time, large EW corrections. Moreover, although we do not present a rigorous proof, based on considerations related to unitarity and factorization of soft and collinear QCD radiation, we will argue that our MiNLO predictions should preserve full NLO QCD+EW accuracy in the phase space with zero or one resolved jets. As we will see, this conclusion is supported by our numerical results.

While our NLO EW results are exact (apart from the treatment of photon-initiated contributions), we also present approximate NLO EW predictions in the so-called Sudakov limit, where all kinematic invariants are well above the electroweak scale. Specifically, based on the general results of Refs. [49, 61], we provide explicit analytic expressions for all logarithmic EW corrections to + 0 and 1 jet in next-to-leading-logarithmic (NLL) approximation. Based on the observed accuracy of the NLL Sudakov formulas, this approximation can be exploited both in order to speed up the evaluation of EW corrections at NLO and in order to predict the dominant EW effects beyond NLO.

All needed matrix elements for + 0 and 1 jet at NLO EW have been generated using the OpenLoops program [62, 63], which supports the automated generation of NLO QCD+EW scattering amplitudes for Standard Model processes [64, 65, 66]. The implementation in the POWHEG BOX RES framework was achieved exploiting the generic interface developed in Ref. [67]. For what concerns NLO QCD corrections, on the one hand we implemented in-house analytic expressions for the virtual corrections. On the other hand, following the approach of Ref. [33], for real-emission contributions we used MadGraph4 [68] matrix elements, via the interface described in Ref. [69].

The paper is organized as follows. In Sec. 2 we introduce the various ingredients of and production at NLO QCD+EW. In particular, in Sec. 2.2 we present a schematic proof of the NLO QCD+EW accuracy of MiNLO predictions for inclusive observables. Further technical aspects of the calculation as well as input parameters and cuts are specified in Sec. 3. Fixed-order NLO QCD+EW predictions are discussed in Sec. 4, while in Secs. 5 and 6 we present NLO+PS QCD+EW results for production and MiNLO QCD+EW results for production, respectively. The predictions of the NLO+PS and MiNLO+PS generators are compared in Sec. 7. Our main findings are summarized in Sec. 8. In the appendices we document the validation of EW corrections in production against HAWK (App. A), detailed NLO EW formulas in the Sudakov approximation (App. B), a reweighting approach that we employ in order to speed up the evaluation of EW corrections (App. C), and technical aspects of the interface between the POWHEG BOX RES and Pythia 8.1 (App. D).

## 2 NLO QCD and EW corrections to HV and HVj production

In this section we describe the QCD and EW NLO corrections to the production of a Higgs boson in association with a or leptonic pair plus zero or one additional jets. For convenience, these Higgsstrahlung processes will be denoted as associated and production, with or . However, all results presented in this paper correspond to the complete processes

 pp→HW+(j)→Hℓ+νℓ(j), pp→HW−(j)→Hℓ−¯νℓ(j), (1) pp→HZ(j)→Hℓ+ℓ−(j),

including all spin-correlation and off-shell effects. The combination of / and / Higgsstrahlung will be denoted as / production. In our calculations, we have considered only one leptonic generation, and all leptons are treated as massless.

### 2.1 NLO QCD+EW matrix elements

In this section we describe the various tree and one-loop amplitudes that have been assembled to form a NLO QCD+EW Monte Carlo program based on the POWHEG BOX framework [30].

Associated production proceeds through quark–antiquark annihilation at leading order, which corresponds to . In production, where the leading order corresponds to , additional (anti)quark–gluon initiated processes contribute. All NLO QCD corrections to production have been computed analytically, since they simply affect the vertex, and the calculation of the real and virtual corrections is trivial. In production, the virtual NLO QCD corrections have been computed analytically [70]. The color- and spin-correlated Born amplitudes and the real contributions at have been computed using the automated interface [69] between the POWHEG BOX and MadGraph4 [68]. The real contributions involve tree diagrams with either an additional gluon or an external gluon replaced with a -pair. Example diagrams are shown in Fig. 1 (a, b).

The virtual EW corrections to and production comprise loop amplitudes up to pentagon and hexagon configurations, respectively. Example diagrams for production are shown in Fig. 2. All the internal resonances have been treated in the complex-mass scheme [71, 72] throughout.

As pointed out in the Introduction and illustrated in Fig. 2 (c), the virtual NLO EW amplitudes induce a dependence on the Higgs trilinear coupling . This dependence arises both from the bare virtual amplitudes and from the Higgs boson self-energy entering the Higgs boson wave-function renormalization. In view of the possibility of exploiting precision measurements of Higgsstrahlung processes for an indirect determination of , we allow to be set independently of the Higgs boson mass.111The corresponding parameter can directly be set in the POWHEG BOX RES input file.

The real NLO EW corrections to and production comprise QED radiation off all charged particles, i.e. they have an additional photon in the final state, as illustrated in Fig. 1 (c). Photon-induced real radiation contributions, where the photon is crossed to the initial state, are, on the other hand, not considered here, as they are suppressed by the small photon density in the proton. These corrections for production have been computed for the first time in Ref. [44] and are included in the HAWK [45] Monte Carlo generator. Interestingly, they reach several percent for inclusive production, but remain at the level when leptonic selection cuts are applied, and are negligible for production [73]. For production, photon-induced contributions enter already at Born level: however, they are of and thus formally subleading with respect to the leading order. Still, the NLO QCD corrections to these photon-induced processes are of and thus formally of the same order as the NLO EW corrections to the quark–antiquark and (anti)quark–gluon initiated channels in production. Also not considered here are mixed QCD-EW bremsstrahlung contributions to production at . These tree-level contributions are finite and can easily be investigated separately. Similar contributions in the NLO EW corrections to +jet production are known to yield relevant contributions only in jet observables at very large transverse momentum [74, 65]. Finally, also virtual QCD corrections to production contribute formally at and are thus of the same perturbative order as the NLO EW corrections to production. However, if a photon isolation is applied, as is done in this paper (see Sec. 3.4), production can be considered as a separate process and thus excluded from the definition of production.

All the electroweak real and virtual corrections have been computed using a recent interface of the POWHEG BOX RES to OpenLoops [67].

In this study we combine NLO QCD and EW corrections in an additive way, i.e. corresponding perturbative contributions are simply added. At fixed order, an improved description can easily be obtained via a factorized ansatz, where differential NLO QCD cross sections are multiplied with relative EW correction factors. Such a multiplicative combination can be motivated from the factorization of soft QCD radiation and EW Sudakov logarithms, which can be tested comparing relative NLO EW corrections for and production.

### 2.2 MiNLO approach for HVj production at NLO QCD+EW

In order to obtain an optimal description of QCD radiation, both in the hard and soft regime, all NLO QCD+EW calculations for have been performed using the “Multiscale improved NLO” (MiNLO[34] method. This approach effectively resums logarithmic singularities of soft and collinear type to NLL accuracy, thereby ensuring a finite cross section in all regions of phase space, even when the extra jet becomes unresolved. In the MiNLO approach, NLL resummation is achieved by means of a CKKW scale setting [75, 76] for the strong coupling factors associated with each QCD vertex, together with an appropriate factorization-scale choice and NLL QCD Sudakov form factors. These are applied to all internal and external lines corresponding to the underlying-Born skeleton of each event. In addition, improving the MiNLO resummation as described in Ref. [35], we have obtained a fully inclusive description of production with NLO QCD accuracy in all phase space regions. In other words, besides providing kinematic distributions that are NLO accurate and also finite when the hardest jet goes unresolved, the improved MiNLO predictions for are NLO accurate also for distributions in inclusive variables such as the rapidity or the transverse momentum of the pair.

All NLO QCD+EW predictions for presented in this paper, both at fixed order and including matching to the parton shower, are based on the MiNLO approach, which is applied to all contributions of NLO QCD and NLO EW type. Technically, the MiNLO Sudakov form factors and scale choices are implemented at the level of underlying-Born events that correspond to so-called terms in the POWHEG jargon.222Real-emission events of NLO QCD and NLO EW type are related to underlying-Born events of type via FKS mappings [77]. Note that the MiNLO procedure resums only logarithms associated with soft and collinear QCD singularities that result form the presence of QCD radiation at Born level, while QED radiation is not present at Born level. Thus there is no need to introduce NLO EW effects in the MiNLO Sudakov form factors. This implies that, in contrast to the case of NLO QCD, the NLO EW corrections to do not need to be matched to the MiNLO form factors. In practice, for what concerns the EW corrections, the MiNLO procedure is applied in a way that is equivalent to Born level.

For observables where QCD radiation is integrated out, the MiNLO improved NLO EW contributions assume the form

 dσMiNLOEWHVjdΦHV=∫dΦj¯BEWHVj(ΦHV,Φj)Δ(kT(Φj)), (2)

where and denote the factorized phase spaces of the system and the jet, respectively. The term includes corrections333Since Born contributions are part of the usual QCD term, in the term we only include corrections. of virtual and real type, and the latter are integrated over the corresponding emission phase space. The MiNLO approach is implemented through an implicitly understood CKKW scale choice for the term in , and through the NLL Sudakov from factor in Eq. (2). For later convenience, together with the Sudakov form factor, we introduce a corresponding emission kernel that is formally related to via

 (3)

In the following, based on the factorization properties of soft and collinear QCD radiation, encoded in the kernel , and using the unitarity relation

 ∫dΦjK(Φj)Δ(kT(Φj))=1, (4)

we will argue that the inclusive MiNLO predictions of Eq. (2) are not only NLO QCD accurate, but also NLO EW accurate. More precisely, we will prove (in a schematic way) that

 dσMiNLOEWHVjdΦHV=dσNLOEWHVdΦHV+O(αEMαS), (5)

where

 dσNLOEWHVdΦHV=¯BEWHV(ΦHV). (6)

We first demonstrate the Born-level version of Eq. (5), which corresponds to

 dσMiLOHVjdΦHV=dσLOHVdΦHV+O(αS), (7)

where MiLO denotes the Born (or LO) version of the MiNLO approach. The above identity can be written as

 ∫dΦjBHVj(ΦHV,Φj)Δ(kT(Φj))=BHV(ΦHV)+O(αS), (8)

where and are the Born counterparts of the and terms in Eqs. (2) and (3). The meaning of Eqs. (7) and (8) is that the MiNLO approach at Born level guarantees LO accuracy for observables that are inclusive with respect to the extra jet. In order to demonstrate this property, we split the Born term into an IR divergent and a finite part,

 BHVj(ΦHV,Φj)=BHV(ΦHV)K(Φj)+BfinHVj(ΦHV,Φj). (9)

Here the singularities associated with QCD radiation in the soft and collinear limits are factorized444In this schematic derivation we assume a simple factorization of multiplicative type, while the factorization of initial-state collinear singularities takes the form of a convolution. into the Born term times the NLL kernel , while the remainder is free from singularities. Thus, upon integration over the jet phase space, the remainder yields only suppressed contributions with respect to , while using the unitarity relation (4) it is easy to show that the singular term in Eq. (9) leads to Eq. (8).

Thanks to the fact that applying the MiNLO approach to NLO EW contributions is largely equivalent to applying MiNLO at Born level, the NLO EW accuracy property (5) can be proven along the same lines as for the LO accuracy property (7). As sole additional ingredient, the NLO EW proof requires certain factorization properties of soft and collinear QCD radiation. More precisely, the factorization properties of Eq. (9) must hold also in the presence of EW corrections, i.e.

 ¯BEWHVj(ΦHV,Φj)=¯BEWHV(ΦHV)K(Φj)+¯BEW,finHVj(ΦHV,Φj). (10)

Here the remainder should be free from QCD singularities, so that it yields only -suppressed contributions relative to , when the extra jet is integrated out. Based on this natural assumption, in full analogy with the LO case, we easily arrive at

 ∫dΦj¯BEWHVj(ΦHV,Φj)Δ(kT(Φj)) = ¯BEWHV(ΦHV)∫dΦjK(Φj)Δ(kT(Φj))+O(αEMαS) (11) = ¯BEWHV(ΦHV)+O(αEMαS),

which is equivalent to the hypothesis (5).

In summary, based on unitarity and factorization properties of QCD radiation, we expect that the improved MiNLO procedure applied to NLO QCD+EW matrix elements for should preserve its full QCD+EW accuracy when the jet is integrated out. As we will see, this conclusion is well supported by our numerical findings in Secs. 47. Nevertheless, due to the schematic nature of the presented derivations and related assumptions, the above conclusions should be regarded as an educated guess that deserves further investigation.

### 2.3 Sudakov approximation at NLO EW

In the Sudakov high-energy regime, where all kinematic invariants are of the same order and much larger than the electroweak scale, the NLO EW corrections are dominated by soft and collinear logarithms of Sudakov type. Based on the general results of Refs. [49, 61] we have derived analytic expressions for the NLO EW corrections to and production in NLL approximation. Details and scope of this approximation are discussed in App. B.

The Sudakov approximation at NLO provides us with qualitative and quantitative insights into the origin of the dominant NLO EW effects. Moreover, it can be easily extended to the two-loop level [51, 52], thereby opening the door to approximate NNLO EW predictions based on the combination of exact NLO EW corrections with Sudakov logarithms at two loops. From the practical point of view, the Sudakov approximation at NLO permits to obtain the bulk of the EW virtual corrections at much higher computational speed as compared to an exact NLO EW calculation.

In Sec. 5, we will assess the quality of the Sudakov approximation555 As explained in more detail in App. B, the Sudakov approximation is applied only to the virtual part of EW corrections, while real QED radiation is always treated exactly. through a detailed comparison against exact NLO EW corrections. Finally, in App. C, we show how the NLL EW approximation can be used in order to speed up the Monte Carlo integration, while keeping full NLO EW accuracy in the final predictions.

## 3 Technical aspects and setup of the simulations

### 3.1 The Powheg Box Res framework at NLO QCD+EW

The QCD+EW NLO calculations for and production have been matched to parton showers using the POWHEG method. To this end, we used the recently-released version of the POWHEG BOX framework, called POWHEG BOX RES. The major novelty of this new version is the resonance-aware approach [58], which guarantees a consistent treatment of intermediate resonances at NLO+PS level. This is achieved by generating the hardest radiation in a way that preserves the virtuality of resonances present at the underlying-Born level. At the same time, the resonance information can be passed on to the parton shower, which in turn preserves the virtuality of intermediate resonances of the hard process in subsequent emissions. This method was introduced in order to address the combination of NLO QCD corrections with parton showers in the presence of top-quark resonances. However, since it is based only on general properties of resonances and infrared singularities, the resonance-aware approach is applicable also to the combination of EW corrections with QED parton showers. In fact, this method has already been applied in the context of electroweak corrections in Refs. [59, 60].

In the POWHEG BOX RES jargon [58], a radiated parton (or photon) can be associated to one or more “resonances” present in the process, or to the “production” part, if it cannot be associated to a particular resonance. The POWHEG BOX RES framework automatically finds all the possible so-called “resonance histories” for a given partonic process. For the processes at hand, considering QED radiation, only two resonance histories are detected: a production history, where the photon can be emitted by any quark (both in the initial and in the final state), and a vector-boson decay history, where the photon is radiated off a final-state charged lepton, and the virtuality of the intermediate vector boson needs to be preserved. Soft photons that are radiated from a resonance are attributed either to the production subprocess or to the decay, consistently with the virtualities of the quasi-resonant propagators “before” and “after” the photon emission.

The treatment of QED radiation was first introduced in the POWHEG BOX for the calculation of the EW corrections to Drell–Yan processes [78, 79, 59, 60]. In this context, leptons were considered as massive particles, and QED subtraction in the POWHEG BOX was implemented accordingly. In the study at hand, leptons are treated as massless, and we have implemented the treatment of photon radiation off massless charged particles (both leptons and quarks). To this end, we have adapted the QCD soft and virtual counterterms already present in the POWHEG BOX to the QED case. Moreover, we have computed a new upper-bounding function for the generation of photon radiation with the highest-bid method, as described in Ref. [29].

By default, in the POWHEG BOX RES framework, only the hardest radiation out of all singular regions is kept, before passing the event to shower Monte Carlo programs like Pythia or Herwig. In this way, for each event, at most one of the decaying resonances (or the production part of the process) includes an NLO-accurate radiation. Moreover, in case of combined QCD and EW corrections, QED emission occurs in competition with the QCD one. The POWHEG BOX RES uses the highest-bid method to decide what kind of radiation (QED or QCD, initial- or final-state) is generated. Due to the larger center-of-mass energy available in the production stage, initial-state radiation is enhanced with respect to final-state radiation, and since the QCD coupling is larger than the QED one, initial-state quarks tend to radiate gluons rather than photons. Thus, QED emission from the decay of a resonance would hardly be kept at the Les Houches event (LHE) level, and the QED radiation would mainly be generated by the shower Monte Carlo program.

The resonance-aware formalism implemented in the POWHEG BOX RES framework offers the opportunity to further improve the POWHEG radiation formula. With this improvement, first introduced in Ref. [80], radiation from each singular region is generated and, instead of keeping only the hardest overall one, the hardest from each resonance is stored. As a result, the LHE file contains a radiated particle for each decaying resonance, plus possibly one emission from the production stage. In this way NLO+LL accuracy is ensured for radiation off each resonance. The subsequent shower from each resonance generated by the Monte Carlo shower program has to be softer than each corresponding POWHEG radiation.666This multiple-radiation mode can be activated by setting the flag allrad to 1 in the input file. All NLO+PS results presented in this paper are based on this multiple-radiation scheme.

As a final remark, we note that in the POWHEG BOX RES framework both the and processes can be computed at NLO or NLO+PS level with only QCD corrections, with only EW corrections, or with combined NLO QCD+EW corrections.777The flag qed_qcd controls this behavior in the input file. The values it can assume are: 0, to compute only QCD corrections, 1, to compute only EW corrections or 2, for both.

### 3.2 OpenLoops tree and one-loop amplitudes

All needed amplitudes at NLO EW have been generated with OpenLoops [62, 63] and implemented in the POWHEG BOX RES framework through the general interface introduced in Ref. [67]. Thanks to the recursive numerical approach of Ref. [62] combined with the COLLIER tensor reduction library [81], or with CutTools [82], the OpenLoops program permits to achieve high CPU performance and a high degree of numerical stability. The amplitudes employed for the EW corrections in this paper are based on the recently achieved automation of EW corrections in OpenLoops [64, 65, 66].

Within OpenLoops, ultraviolet and infrared divergences are dimensionally regularized in dimensions. However, all ingredients of the numerical recursion are handled in four space-time dimensions. The missing -dimensional contributions, called rational terms, are universal and can be restored from process-independent effective counterterms [83, 84, 85]. The implementation of the corresponding Feynman rules for the complete EW Standard Model in OpenLoops is largely based on Refs. [86, 87, 88, 89]. Relevant contributions for and production have been validated against independent algebraic results in dimensions. UV divergences at NLO EW are renormalized in the on-shell scheme [90] extended to complex masses [71].

### 3.3 Input parameters, scales choices and other aspects of the setup

In our (+jet) simulations at NLO QCD+EW, we have set the gauge-boson masses and widths to the following values [91]

 MZ=91.1876 GeV, MW=80.385 GeV, ΓZ=2.4955 GeV, ΓW=2.0897 GeV.

The latter are obtained from state-of-the-art theoretical calculations. Assigning a finite width to the Higgs boson in the final state would invalidate EW Ward identities: we then consider the Higgs boson as on shell with and set its mass to  GeV. The top-quark mass and width are set respectively to  GeV and  GeV. All other quarks and leptons are treated as massless. In the EW corrections, the top-quark contribution enters only at loop level, the dependence of our results on is thus completely negligible.

For the treatment of unstable particles we employ the complex-mass scheme [71, 72], where finite-width effects are absorbed into complex-valued renormalized masses

 μ2k=M2k−iΓkMkfor k=W,Z,t. (13)

The electroweak couplings are derived from the gauge-boson masses and the Fermi constant, , and the electromagnetic coupling is set accordingly to

 αEM=∣∣ ∣∣√2s2wμ2WGμπ∣∣ ∣∣, (14)

where and the squared sine of the weak mixing angle

 s2w=1−c2w=1−μ2Wμ2Z, (15)

are complex-valued.888By default we use the scheme throughout. However, in the POWHEG BOX RES framework, there is the option to evaluate the virtual EW corrections using computed in the scheme, and use the Thomson value in the evaluation of the contribution due to photon radiation.

The absolute values of the CKM matrix elements are set to

 |VCKM|=d  s  buct⎛⎜⎝0.974280.22530.003470.22520.973450.04100.008620.04030.999152⎞⎟⎠. (16)

Our default set of parton-distribution functions (PDF) is the NNPDF2.3_as_0119_qed set [92], that includes QED contributions to the parton evolution and a photon density.999It corresponds to the PDF set 244800, in the LHAPDF6 [93] numbering scheme. The value of the strong coupling constant corresponding to this PDF set is .

Finally, in production, the renormalization and factorization scales are set equal to the invariant mass of the pair at the underlying-Born level,

 μR=μF=MHV,M2HV=(pH+pℓ1+p¯ℓ2)2, (17)

where and are the final-state leptons, while in the improved MiNLO [34, 35] procedure is applied, and the scales are set accordingly.

Predictions at NLO+PS generated with the POWHEG method are combined with the Pythia 8.1 QCD+QED parton shower using the “Monash 2013” tune [94]. Effects due to hadronization, multi-particle interactions and underlying events are not considered in this paper.

### 3.4 Physics objects and cuts in NLO+PS simulations

In the following we specify the definition of physics objects and cuts that are applied in the phenomenological NLO+PS studies presented in Secs. 57.

All leptonic observables are computed in terms of dressed leptons, which are constructed by recombining the collinear photon radiation emitted within a cone (in the plane) of radius from charged leptons, and the recombined photons are treated as unresolved particles. Observables that depend on the reconstructed vector bosons are defined by combining the momenta of the dressed charged leptons and the neutrino associated with their decay. The latter is taken at Monte Carlo truth level.

Jets are constructed with FastJet using the anti- algorithm [95, 96] with . The jet algorithm is applied in a democratic way to QCD partons and non-recombined photons, with the exception of photons that fulfill the isolation criterion of Ref. [97] with a cone of radius and a maximal hadronic energy fraction . The hardest of such isolated photons is excluded from the jet algorithm and is treated as resolved photon.

The following standard Higgsstrahlung cuts are applied. For every dressed charged lepton we require

 pℓT≥25 GeV,|yℓ|≤2.5. (18)

In / production, we also impose

 ⧸ET≥25 GeV, (19)

where is the transverse momentum of the neutrino that results from the -boson decay at Monte Carlo truth level. In / production, the invariant mass of the dressed-lepton pair is required to satisfy

 60 GeV≤Mℓ+ℓ−≤140 GeV. (20)

Besides these inclusive selection cuts, we also present more exclusive results in the boosted regime. In this case, we impose the following additional cuts on the transverse momentum of the Higgs and vector bosons

 pHT≥200 GeV,pVT≥190 GeV. (21)

Such a selection of events with a boosted Higgs boson improves the signal-over-background ratio in the decay channel.

## 4 Results for HV and HVj production at fixed NLO QCD+EW

In this section we present fixed-order NLO QCD+EW predictions for and at 13 TeV. For production the improved MiNLO approach [34, 35] is applied. Higgs boson production in association with and bosons is discussed in Secs. 4.1 and 4.2, respectively. Predictions based on exact NLO EW calculations (apart from photon-initiated contributions that have been neglected) are compared against the Sudakov NLL approximation (see App. B), which includes virtual EW logarithms supplemented by an exact treatment of QED radiation.

The fixed-order results presented in this section are not subject to the cuts and definitions of Sec. 3.4. No acceptance cut is applied, and differential observables are defined in terms of the momenta of the Higgs and vector bosons. The latter are defined in terms of the momenta of their leptonic decay products at the level of underlying-Born events, i.e. before the emission of NLO radiation. Photons and QCD partons are clustered in a fully democratic way using the anti- algorithm with . Effectively this procedure corresponds to an inclusive treatment of QED radiation.

Besides total cross sections, we consider various differential distributions, focusing on regions of high invariant masses and transverse momenta, where EW corrections are enhanced by Sudakov logarithms. Such phase-space regions play an important role for experimental analyses of production in the boosted regime.

An in-depth validation of our fixed-order NLO EW results for production against the ones implemented in the public Monte Carlo program HAWK [44, 45], is presented in App. A.

### 4.1 HW and HWj production

In this section we focus on NLO results for and . In Tab. 1 we report inclusive NLO cross sections. In the case of production, the improved MiNLO approach yields finite cross sections without imposing any minimum transverse momentum on the hardest jet. For comparison, we report also  MiNLO cross sections for the case where a minimum of 20 GeV is required for the hardest jet. In this case, the MiNLO Sudakov form factor plays hardly any role, since it damps the cross section only at of the order of a few GeV, i.e. far below the imposed cut. Thus, at fixed order, the MiNLO procedure only affects the choice of scales, as described in Sec. 3.3. The EW corrections lower the inclusive NLO QCD cross section by roughly for production and for production, while they amount to only when a resolved jet with  GeV is required in the calculation.

Inclusive cross sections in the NLO QCD+NLL EW approximation differ by several percent from the exact NLO QCD+EW results. This is expected, since the NLL approximation is only valid in the high-energy regime.

In the following we investigate the impact of EW corrections and the validity of the NLL approximation in differential distributions for and production. Results for production (not shown) are very similar.

In Fig. 3 we plot the invariant mass of the reconstructed pair, both for and production. The three curves represent predictions at NLO QCD, NLO QCD+EW and in NLO QCD+NLL EW approximation. While EW corrections have a moderate impact on the total cross sections, they affect the tail of the distribution in a substantial way. At large we observe the typical Sudakov behavior, with increasingly large negative EW corrections that reach the level of () for () production at 2 TeV. The Sudakov NLL approximation captures the bulk of these large EW corrections as expected. In the tail it agrees at the percent level with the exact result for both processes, while for moderate invariant masses it overestimates EW correction effects by up to 5%.

In Fig. 4 we investigate the transverse momentum of the Higgs boson. Also in this case EW corrections become negative and large in the tail, exceeding in the TeV region. For both processes the Sudakov approximation agrees at the percent level with exact NLO EW results for  GeV.

The EW corrections have a sizable impact also on the missing transverse momentum distribution, shown in Fig. 5. Size and shape of these corrections are very similar to the ones observed for the Higgs boson distribution.

In Fig. 6 we present predictions for the distribution in the of the leading jet. At low (left plot) the MiNLO Sudakov form factor damps soft and collinear singularities at zero transverse momentum yielding finite cross sections below the Sudakov peak, which is located around 3 GeV. Concerning EW effects, the NLL approximation converges to the exact NLO results already for values of around 200 GeV. In the region of moderate transverse momentum, NLO EW corrections are nearly constant, and in the limit of vanishing jet- they converge towards an EW -factor that is very close to the one of the NLO QCD+EW calculation for the inclusive cross section (see Tab. 1). This observation is consistent with the theoretical considerations presented in Sec. 2.2, namely with the fact that EW corrections are insensitive to soft and collinear QCD radiation, and that MiNLO predictions for production preserve NLO QCD+EW accuracy when the extra jet is integrated out. In fact, in the inclusive distributions of Figs. 35, we observe that the EW corrections obtained from and calculations are very similar, with small differences that can be attributed to NNLO effects.

Finally, in Fig. 7 we see that the EW corrections affect the rapidity of the leading jet in a rather uniform way over the whole phase space. We have observed a similar behavior of EW corrections in several angular distributions.

### 4.2 HZ and HZj production

In line with the discussion of and production, we present in this section fixed-order results for and production. In Tab. 2 we collect the inclusive cross sections at NLO QCD, NLO QCD+EW and NLO QCD+NLL EW. The EW corrections decrease the total NLO QCD cross section for production by about 4%, and by about 1% for inclusive production. In the presence of a jet threshold of  GeV, the EW corrections are positive and amount to about .

In Figs. 8 and 9 we show distributions of the invariant mass of the reconstructed pair and of the transverse momentum of the Higgs boson. Similarly as for production, the inclusion of EW corrections is essential in the tails of these distributions, where the NLL Sudakov approximation agrees well with the exact NLO EW predictions.

## 5 Results for HV production at NLO+PS QCD+EW

In this section we present NLO QCD+EW predictions for production completed by the Pythia 8.1 QCD+QED parton shower using the “Monash 2013” tune [94]. All predictions are subject to the cuts and physics object definitions specified in Sec. 3.4, and NLO EW corrections are treated exactly throughout, except for photon-initiated processes, that have been neglected. The NLL Sudakov approximation is only used in order to speed up the Monte Carlo integration, as detailed in App. C.

In Sec. 5.1 we compare predictions at fixed-order NLO QCD+EW against corresponding predictions at the level of Les Houches events, which include only the hardest emission generated in the POWHEG BOX RES framework, and at NLO+PS level, where the full QCD+QED parton shower is applied. The effect of EW corrections is studied in Sec. 5.2 in the case of fully showered NLO+PS simulations.

By default, at NLO+PS level, the full QCD+QED parton shower is applied, both for NLO QCD+EW and for pure NLO QCD simulations. Occasionally, we also present NLO QCD simulations with a pure QCD shower, where QED radiation is switched off. Such predictions are labeled “QCD (no QED shower)”.

The consistent combination of the NLO radiation to the parton shower requires the vetoing of shower emissions that are harder than the radiation generated in the POWHEG BOX RES framework. Since no standard interface is available in a multi-radiation scheme, we have implemented a dedicated veto procedure on the Pythia 8.1 showered events, as described in App. D. This veto procedure is applied in case of NLO QCD+EW simulations. Instead, in case of NLO QCD simulations combined with the Pythia 8.1 QCD+QED shower, only QCD radiation is restricted by the POWHEG BOX RES hardest scale, while arbitrarily hard QED radiation can be generated by the shower.

We have verified that inclusive cross sections at NLO+PS QCD and NLO+PS QCD+EW agree within statistical uncertainties with the corresponding fixed-order results reported in Tabs. 1 and 2. Thus, in the following we will focus on differential distributions.

### 5.1 From fixed NLO QCD+EW to NLO+PS QCD+EW

In this section we compare NLO QCD+EW predictions at fixed order with NLO+PS ones at LHE level and completed with the Pythia 8.1 shower. Since the various Higgsstrahlung processes behave in a very similar way, we will focus on production.

In Fig. 10 we plot the rapidity of the reconstructed pair, which is NLO accurate, and its transverse momentum, which is only LO accurate. Due to the inclusiveness of the rapidity of the pair, we find, as expected, very good agreement, within the integration errors, among the three predictions. The fixed-order curve for the transverse momentum displays the typical divergent behavior at low . At LHE level, instead, the divergence is tamed by the Sudakov form factor. The effect of the parton shower is modest in the tail of this distribution, while at low it slightly shifts the position of the Sudakov peak.

In Fig. 11 we plot the pseudorapidity and the transverse momentum of the Higgs boson: thanks to the inclusiveness of this variable, we find again very good agreement among the three predictions.

### 5.2 Impact of the EW corrections in NLO+PS events

In this section we investigate EW correction effects at the level of fully showered NLO+PS predictions.

In Fig. 12 we show the rapidity (left) and the transverse momentum (right) of the charged dressed lepton in production. In the rapidity distribution, the impact of NLO EW effects is constant and amounts to about . The shape of the distribution, instead, changes drastically due to EW Sudakov logarithms in the high- region, where differences with respect to the pure QCD predictions reach around 1 TeV.

In Fig. 13 we plot the transverse mass of the reconstructed boson

 MWT=√2pℓT⧸ET(1−cosΔϕ), (22)

where is the azimuthal angle between the charged lepton and the missing transverse momentum. Similarly, as for the lepton rapidity, the EW corrections do not change the shape, but lower the differential cross section by about 7% with respect to the pure QCD corrections. If no QED shower is activated when Pythia 8.1 showers QCD-corrected events, the curve that is obtained is very similar to the QCD one, i.e. the impact of the QED shower is small for this distribution and no radiative tail can be observed.

In Fig. 14 we show the rapidity and the transverse momentum of the Higgs boson in the boosted regime, as defined by the cuts of Eq. (21). The EW corrections have a constant negative impact around 10% on the rapidity distribution, and reach around 1 TeV. Similar conclusions can be drawn for the rapidity and transverse momentum of the boson.

We conclude this section by presenting kinematic distributions for production in Figs. 1517.

In Fig. 15 we show the distribution in the rapidity and the transverse momentum of the dressed electron. The EW corrections give a constant contribution of about  in the plotted rapidity range, while in the high-energy tail of the distribution the EW corrections decrease the differential cross section by roughly 30% due to Sudakov logarithms.

In Fig. 16 we plot the invariant mass of the reconstructed leptonic pair in the region around the resonance. In spite of the fact that the shape of the resonance it known to receive very large radiative corrections (see e.g. Refs. [98, 79]), NLO EW effects turn out to be almost constant and as small as % when we compare showered NLO+PS predictions at NLO QCD+EW versus NLO QCD. This is due to the fact that the bulk of the radiation is correctly described by the QED shower in Pythia 8.1. The importance of radiation becomes evident when we switch off the QED shower (“no QED shower”) in the NLO QCD simulation. This results in a radiative tail with distortions of up to in the region below the peak.

In production, the momentum of the vector boson can be fully reconstructed. Thus, in Fig. 17 we display the rapidity and transverse-momentum distributions of the boson in the boosted regime, as defined by the cuts of Eq. (21). These results are very similar to the ones obtained for the Higgs boson in production in Fig. 14. While EW corrections have a constant impact of about  on the rapidity distribution, the tail of the distribution is dominated by large negative EW Sudakov logarithms, and we observe differences with respect to the pure QCD result of the order of  for  TeV.

## 6 Results for HVj production at NLO QCD+EW with MiNLO+Ps

In this section we study at NLO+PS accuracy in the MiNLO approach, denoted in the following as MiNLO+PS. Similarly as in the previous section, in Sec. 6.1 we first compare NLO QCD+EW predictions obtained with MiNLO at fixed order against corresponding results at the LHE level or including also the full QCD+QED parton shower. The effect of EW corrections is studied in Sec. 6.2 in the case of fully showered MiNLO+PS simulations. The cuts and physics object definitions of Sec. 3.4 are applied throughout, and we do not impose any cut that requires the presence of jets.

### 6.1 From fixed-order MiNLO to MiNLO+PS at NLO QCD+EW

In Fig. 18 we analyze the rapidity and the transverse momentum of the reconstructed system. As a result of the MiNLO prescription, the rapidity distribution, as well as any other inclusive observable, is finite. For the rapidity distribution we observe that the three predictions are very close to each other. At variance with Fig. 10, here fixed-order predictions for the distribution are finite at small , since soft and collinear divergences are suppressed by the MiNLO Sudakov form factors. Moreover, the NLO accuracy in the spectrum of the system leads to an improved agreement between fixed-order and NLO+PS results.

In Fig. 19 we show the pseudorapidity and the transverse-momentum distributions of the Higgs boson, finding again very good agreement among the three predictions.

We refrain from presenting results for and production since they behave qualitatively very similar as the results shown here for production.

### 6.2 Impact of EW corrections at MiNLO+PS level

The impact of EW corrections at the level of fully showered MiNLO+PS predictions for production is illustrated in Figs. 2022.

For the distributions in the rapidity and transverse momentum of the Higgs boson in the boosted regime (Fig. 20) we find that the EW corrections induced by the boosted cut, , are nearly independent of and around , while they grow up to and beyond when enters the TeV regime. These results closely agree with the corresponding ones shown in Fig. 14 for the NLO+PS simulation of inclusive production. Consistently with the fixed-order findings discussed in Sec. 4, also this observation supports the theoretical considerations of Sec. 2.2, where we have argued that MiNLO improved predictions for production should preserve NLO QCD+EW accuracy when the extra jet is integrated out. Also other inclusive observables, such as the distribution in the missing transverse momentum shown in Fig. 21, confirm this observation.

The EW corrections to the leading-jet distribution, shown in Fig. 22, do not feature the standard Sudakov behavior. In this distribution, EW effects remain rather small, at the level of , in the entire plotted range, i.e. from very low jet- up to 400 GeV. This is not surprising, since a similar “non-Sudakov” behavior for inclusive jet spectra was already observed in Ref. [74, 65] for the case of  jets production. Another important feature of Fig. 22 is that EW corrections are nearly constant in the region where the jet approaches zero. Again, this confirms the considerations made in Sec. 2.2 regarding the factorization of EW corrections in the presence of soft or collinear QCD radiation, and the NLO QCD+EW accuracy of inclusive MiNLO simulations. To be more precise, in the left panel of Fig. 22 we see that EW corrections effects are nearly constant at small with the exception of the first bin. This effect can be attributed to photonic contributions to the jet transverse momentum, and to the fact that the Sudakov peak associated with the damping of QCD radiation is located well above the one associated with the damping of QED radiation. This mismatch tends to enhance the relative importance of QED radiation in the region between the QCD and QED Sudakov peaks. In any case, this effect cancels upon integration over the soft region of the jet spectrum. Thus, it should not spoil the expected NLO QCD+EW accuracy of inclusive MiNLO predictions.

We conclude this section by discussing the impact of NLO EW effects in production, illustrated in Figs. 23 and 24. The distribution in the -boson in the boosted regime (Fig. 23) features the typical Sudakov EW behavior, with negative EW corrections that exceed in the tail. In the leading-jet distribution (Fig. 24) we observe relatively small and rather constant EW corrections. Both distributions behave similarly as the corresponding distributions for production.

We refrain from showing further plots for or production, since EW correction effects are quite similar to the ones already discussed.

## 7 Comparison between the HV and HVj generators

In this section, we discuss and compare NLO+PS predictions for against MiNLO+PS predictions for , both at NLO QCD+EW accuracy. A similar comparison at NLO QCD accuracy was presented in Ref. [33]. Since the various Higgsstrahlung processes behave in a very similar way, we will focus on the associated production. The comparison between the and generators is motivated by the fact that the improved MiNLO prescription [35] applied to production provides NLO accuracy also for inclusive quantities, i.e. for observables where the associated jet is not resolved. While this is well known at NLO QCD level, in Sec. 2.2 we have argued that also NLO EW accuracy should be preserved when the jet is integrated out.

We also study the dependence of our results on scale variations. To this end we apply standard seven-point variations obtained by multiplying the central value of the renormalization and factorization scales and , defined in Eq. (17) for production, by the factors and , respectively, chosen among the seven pairs

 (KR,KF)=(12,12),(12,1),(1,12),(1,1),(2,1),(1,2),(2,2). (23)

Scale-variation bands in the following plots are based on the envelope of the seven-point variations. In production, where the scale setting is based on the improved MiNLO prescription, the scaling factors of Eq. (23) are applied to the coupling constants at each interaction vertex and to the scale entering the Sudakov form factor.

For the fully inclusive NLO+PS and MiNLO+PS cross sections at NLO QCD+EW we find

 σNLO+PSHW−=55.29+0.80−0.74 fb,σMiNLO+PSHW−j=55.25+1.25−2.57 fb, (24) σNLO+PSHZ=24.41+0.27−0.38 fb,σMiNLO+PSHZj=24.9+0.6−1.1 fb,

where uncertainties correspond to scale variations. These results are well consistent, within statistics, to the corresponding ones reported in Tabs. 1 and 2 at fixed-order NLO. Moreover, it turns out that, in the presence of EW corrections, cross sections obtained from and simulations agree at the one-percent level, confirming again the expectation of inclusive NLO QCD+EW accuracy for MiNLO improved simulations.

Scale variations are in general larger in with respect to production. This is due to the fact that, in standard POWHEG BOX simulations, the scale associated to the emission of the hardest jet is kept fixed at the corresponding transverse momentum, while scale variations are applied only at the level of the so-called term, where QCD and QED radiation are integrated out. For this reason, scale variations in MiNLO+PS simulations provide a more realistic estimate of scale uncertainties associated with QCD radiation.

Figures 2527 display differential distributions subject to the cuts of Sec. 3.4. Red bands correspond to scale variations for and production. We do not show the statistical uncertainties associated to the integration procedure on these bands, since they are much smaller than the width of the bands. When plotting instead the blue curves for the distributions computed at the central scales, we display the statistical uncertainties of the integration procedure as an error bar. The plots on the left-hand side show the uncertainty band for the process, while the ones on the right-hand-side show the uncertainty band for production.

In Fig. 25 we display the rapidity distribution of the system. Since this inclusive quantity is predicted at NLO QCD accuracy by both simulations, we find very good numerical agreement between the two curves at NLO QCD+EW level. The uncertainty band is larger in the case. This is due to the fact that for production there is no renormalization-scale dependence at LO, while in such dependence is already present at leading order.

In Figs. 26 and 27 we compare the transverse momentum of the pair in two different ranges. Here we observe significant differences due to the fact that this distribution is only computed at leading order in the simulation, while it is NLO accurate in the case. Since we included also EW corrections, in our plots these differences are slightly more pronounced than in the pure QCD implementation of Ref. [33]. The fact that such differences emerge in the region below the QCD Sudakov peak (Fig. 26) is consistent with the observation of enhanced EW effects in that region (Fig. 22) as discussed in Sec. 6.2. We also note that the uncertainty band for the generator is smaller than the one. This is due again to the fact that, at Born level, production does not depend upon , while production does, and this dependence amplifies the scale-variation band.

## 8 Summary and conclusions

In this paper we have presented the first NLO QCD+EW calculations for and production, with , at the LHC. Specifically, we have considered complete Higgsstrahlung processes corresponding to Higgs boson production in association with off-shell or leptonic pairs plus zero or one jet. In addition to fixed-order predictions we have presented realistic simulations obtained by combining NLO QCD+EW calculations with a QCD+QED parton shower. This was achieved by means of the POWHEG BOX RES generator, a recent extension of the POWHEG BOX V2 framework, that allows for consistent NLO+PS simulations in the presence of resonances. In the case of production, using the improved MiNLO approach, we have extended the applicability of NLO QCD+EW predictions to the full phase space, including regions where the hardest jet is unresolved. This is the first application of the MiNLO and POWHEG BOX RES approaches in combination with NLO EW corrections.

We have studied several kinematic distributions for and production in proton-proton collisions at 13 TeV, and we have discussed predictions at fixed-order NLO, at the level of POWHEG BOX RES Les Houches events, and at NLO+PS level using Pythia 8.1. Particular care has been taken in combining the QCD+QED shower of Pythia 8.1 with the POWHEG BOX-generated events, since no standard interface is available, at present, to deal with multiple NLO emissions that can arise at production and decay level in resonant processes.

Electroweak corrections typically lower NLO+PS QCD predictions by 5 to 10% at the level of integrated cross sections and in angular distributions. We have observed quantitatively similar and rather constant EW corrections also in the jet- spectrum, as well as in the reconstructed -mass and transverse -mass in the vicinity of the corresponding resonances. In contrast, due to Sudakov logarithms, EW corrections can be much more sizable in the tails of transverse-momentum and invariant-mass distributions. For example, in the Higgs and vector-boson distributions, EW corrections reach up to around 1 TeV. In this respect, the and Higgsstrahlung processes behave similarly, i.e. the emission of a jet does not have a sizable impact on EW corrections.

We have studied theoretical uncertainties associated with standard factor-two variations of the renormalization and factorization scales. In the context of the POWHEG formalism, scale variations are performed only at the level of the underlying-Born cross section, while the scale of the strong coupling constant associated with NLO radiation is kept fixed at the corresponding transverse momentum. Thus the resulting scale-variation bands are typically smaller as compared to the ones obtained in fixed-order NLO calculations. In the total cross sections for and production we have found scale uncertainties around 1-2% and 2-4%, respectively, while scale variations in kinematic distributions are typically at the 10% level.

Thanks to the improved MiNLO prescription, simulations based on NLO QCD+EW matrix elements for production can be applied to inclusive observables and compared against more conventional simulations based on NLO QCD+EW matrix elements for production. At NLO QCD, the observed agreement between and predictions confirms that, as is well known, the improved MiNLO approach guarantees NLO QCD accuracy also when the extra jet is integrated out. A similarly good level of agreement was found also at NLO QCD+EW level in a variety of observables. In this regard, based on unitarity and factorization properties of soft and collinear QCD radiation, we have sketched a proof of the fact that the improved MiNLO approach, applied to QCD jet radiation computed with NLO QCD+EW matrix elements, should provide NLO QCD+EW accuracy in the full phase space.

All relevant matrix elements at NLO EW have been computed using a recent interface of the POWHEG BOX RES framework with the OpenLoops matrix-element generator. The other QCD amplitudes have been computed in part analytically and in part using the standard interface to MadGraph4. We have also presented simple analytic expression that approximate the virtual EW amplitudes in the Sudakov regime at next-to-leading-logarithmic accuracy. This approximation captures the bulk of EW corrections and reproduces exact NLO EW results with reasonable accuracy. Moreover it can be exploited in the combination of the reweighting approach that permits to speed up NLO QCD+EW simulations while providing full NLO EW accuracy in the final results.

The POWHEG BOX RES code together with the generators that we have implemented for and production can be downloaded following the instructions at the POWHEG BOX web page: http://powhegbox.mib.infn.it

\acknowledgments

J.L. wishes to thank S. Dittmaier and S. Kallweit for many detailed answers on the implementation of EW corrections to Higgsstrahlung processes in HAWK. This research was supported in part by the Swiss National Science Foundation (SNF) under contract PP00P2-128552 and by the Research Executive Agency (REA) of the European Union under the Grant Agreements PITN–GA–2010–264564 (LHCPhenoNet), and PITN–GA–2012–316704 (HiggsTools).

## Appendix A Validation of the fixed-order NLO EW corrections in HV production

In this section we compare our fixed-order NLO EW predictions for and production with predictions obtained with the Monte Carlo program HAWK [45].

### Setup for the comparison

In order to make a comparison between the results generated by HAWK and our results, we switched off photon-initiated contributions in HAWK, since these contributions are currently not included in the POWHEG BOX RES generators. Similarly, -initiated contributions have been discarded in the POWHEG BOX RES, since this sub-process is not included in HAWK. The CKM matrix elements have been set to

 ∣∣VCKMud∣∣=∣∣VCKMcs∣∣=0.974,∣∣VCKMus∣∣=∣∣VCKMcd∣∣=√1−∣∣VCKMud∣∣2, (25)

omitting mixing with third-generation quarks. The renormalization and factorization scales are set to the default values used in HAWK, i.e. to the sum of the Higgs and the vector boson masses

 μR=μF=MV+MH,V=W,Z. (26)

All other input parameters are chosen in accordance with Sec. 3.3.

Photons are recombined with collinear charged leptons if , where is the angular separation variable in the plane. If more than one charged lepton is present in the final state, the eventual recombination is performed with the lepton having the smallest value of . After photon recombination, we apply the following cuts on the charged dressed leptons

 pℓT>20GeV,|yℓ|<2.5, (27)

while for production we also require a missing transverse momentum of

 ⧸ET>25GeV. (28)

### Results

In Figs. 28 and 29 we compare NLO EW predictions obtained with POWHEG BOX RES (solid line) and HAWK (dashed line) for selected observables in and production.

Figure 28 displays the Higgs boson transverse-momentum and pseudorapidity distributions. Within statistical uncertainties the two predictions fully overlap.

As a further example, in Fig. 29 we plot the transverse momentum of the neutrino, i.e. the missing transverse momentum. Again, we observe perfect agreement between the fixed-order NLO POWHEG BOX RES and HAWK predictions, and a similar level of agreement was found in all considered observables.

As examples for the validation of production, in Fig. 30 we present the transverse momentum and the rapidity of the Higgs boson, and in Fig. 31 the rapidity of the produced electron. Again we find a perfect overlap between the POWHEG BOX RES and HAWK predictions, and the same level of agreement was found for all kinematic distributions that we have examined.

## Appendix B The virtual EW Sudakov approximation

The calculation of EW virtual corrections is typically more complex than in the case of QCD. This is due to the nontrivial gauge-boson mass spectrum, the presence of Yukawa and scalar interactions, the fact that EW corrections enter also in leptonic vector-boson decays, as well as subtleties related to the treatment of unstable particles. For these reasons, the numerical evaluation of NLO EW virtual corrections can be time consuming. Motivated by this practical issue, in this appendix we present compact analytic formulas that provide a decent approximation of the bulk of NLO EW effects, based on the Sudakov approximation. Besides speeding up the numerics, this approximation provides also valuable insights into the origin of the bulk of the EW corrections.

The largest EW corrections originate in the Sudakov regime, where all kinematic invariants are of the same order and much larger than the electroweak scale. In this high-energy regime, the EW corrections are dominated by Sudakov logarithms [99, 100, 46, 101, 48, 49] of the form

 L(s)=αEM4πlog2sM2V,l(s)=αEM4πlogsM2V, (29)

i.e. by leading (LL) and next-to-leading logarithms (NLL) involving the ratio of the partonic center-of-mass energy to the electroweak-boson masses, . Sudakov EW logarithms originate from virtual gauge bosons that couple to one or two on-shell external particles in the soft and/or collinear limits.

General factorization formulas for LLs and NLLs that apply to any Standard Model process at one loop have been derived in Refs. [49, 61, 102]. For a generic -particle scattering processes with all particles and momenta incoming101010In the following we adopt the notation of Refs. [49, 61, 102].

 φ1(p1)φ2(p2)…φn(pn)→0, (30)

high-energy EW logarithms in one-loop matrix elements assume the general factorized form

 δMφ1…φn({λi},p1…pn) = ∑λiδλi∂Mφ1…φn0∂λi({λi},p1…pn) (31) +

Here the first term is related to the running of the dimensionless coupling parameters in the Born amplitude, while the second term consists of process-independent correction factors that contain all LL and NLL terms and multiply Born matrix elements for the process at hand. Note that the correction factors are matrices in space. In general they act on one or two external particles, requiring the evaluation of -transformed matrix elements .

The logarithmic EW corrections of Eq. (31) can be schematically split into five contributions

 δM=(δ