Polarized Z bosons from the decay of a Higgs boson produced in association with two jets at the LHC

Investigating the polarization of weak bosons provides an important probe of the scalar and gauge sector of the Standard Model. This can be done in the Higgs decay to four leptons, whose Standard-Model leading-order amplitude enables to generate polarized observables from unpolarized ones via a fully-differential reweighting method. We study the Z-boson polarization from the decay of a Higgs boson produced in association with two jets, both in the gluon-fusion and in the vector-boson fusion channel. We also address the possibility of extending the results of this work to higher orders in perturbation theory.


Introduction
The discovery of the Higgs boson [1,2] has enabled a large number of tests of the Standard Model (SM) scalar and gauge sectors.
In the SM, electroweak bosons (W and Z) acquire their mass and, as a consequence, an additional longitudinal polarization, through their coupling to the Higgs field. Therefore the interplay between the Higgs boson and weak-boson polarizations represents a crucial probe for the SM electroweak symmetry breaking mechanism as well as an ideal framework for searches of possible modifications due to new-physics effects.
The gold-plated channel for these searches is given by vector-boson scattering (VBS), as in the high-energy regime the unitarity-violating behaviour of the scattering among longitudinal bosons is regularized by the inclusion of Higgs exchange contributions [3][4][5]. Studying the Higgs coupling to polarized weak bosons can provide relevant insight also at the Higgs-boson resonance itself [6], in particular with the aim of discriminating between the SM and models with modified Higgs sectors. Since vector-boson polarizations are best accessed through the angular distribution of their decay products, the cleanest channel is the Higgs decay to four charged leptons.
At the Large Hadron Collider (LHC), the Higgs boson is mostly produced in the gluon-fusion (GGF) and in the vector-boson-fusion (VBF) channels. Experimental analyses can be either inclusive or exclusive in the production mode of the Higgs boson, depending on the specific target of the investigation. For the purposes of this work, we only consider the production of a Higgs boson in association with two jets. In this channel the GGF and VBF contributions are comparable. The SM prediction for an on-shell Higgs boson produced in GGF in association with two jets is known at leading order (LO) in QCD, including also dependence on the top-quark mass [42][43][44]. Predictions in the large-m t limit are known up to NLO QCD [45]. The full mass dependence has also been studied in the presence of high-energy jets [46]. The predictions for VBF are known at NLO QCD+EW accuracy [47], and up to N 3 LO QCD [48] in the structure-function approximation [49]. It has been shown [42,43] that for an on-shell Higgs produced in GGF in association with two jets, the large-m t approximation gives a very good description of the loop-induced SM process with full m t dependence, of order O(α 4 s α), provided that the transverse-momenta of jets and the Higgs-boson mass are smaller than the topquark mass. The logarithmic structure of Higgs production in association with up to two jets and the kinematic configurations where finite top-mass effects become relevant are also known [50,51]. The large-m t approximation works well also in the presence of a two-jet system with large invariant mass [42,43,52], which is the typical VBF phase-space region. The interference between VBF and GGF signals has been proved to be negligible [53].
The phenomenology of weak-boson polarizations has been widely investigated both in experimental analyses and in theoretical studies.
From the theory side, the W-boson polarization at the LHC has been studied in ref. [67] in the absence of lepton cuts. Realistic selection cuts has been introduced in ref. [68] both in V + jets and in other multi-boson production processes. The interference between amplitudes for different polarizations has been investigated in ref. [69].
In ref. [70] a simple and natural method to define cross sections corresponding to vector bosons of definite polarization has been proposed. This method has been applied to study the polarization of W and Z bosons in vectorboson scattering at LO [70][71][72] and in di-boson production at NLO QCD [73,74] with purely leptonic final states. Polarizations in W ± Z production, including NLO QCD and EW corrections have been analyzed in refs. [75,76]. Recently MadGraph has introduced the possibility of generating polarized amplitudes [77]. Ref. [78] has suggested that a study of vector-boson polarizations in gluoninduced ZZ production could be sensitive to the Ztt coupling.
A study of the vector-boson polarization effects at the Higgs resonance in VBF has been performed in the W + W − channel, using effective Higgs couplings to longitudinal and transverse bosons and including a general dimension-six EFT interpretation [6]. The possibility of measuring the coupling of the Higgs to polarized bosons is studied in ref. [64].
In this work we perform a phenomenological study of polarized electroweak bosons from the decay of a SM Higgs boson produced in association with two jets at the LHC. In ref. [79] a first assessment of polarized bosons from Higgs decays has been performed in the LO production of a Higgs in gluon fusion, that is, for a Higgs with vanishing transverse momentum. In this paper we consider the case in which the Higgs boson is produced in association with two jets, both in VBF and in GGF. Although the focus is on the Higgs signals, a number of comments are made concerning the impact of the QCD background and of the pure electroweak contributions to VBF.
The phenomenological analysis is limited to the fourlepton Higgs decay, with two pairs of opposite-charge leptons and different flavours. Despite a very small crosssection, due to the small branching ratio, the considered decay channel has two advantages. First, it enables the complete reconstruction of the Higgs-boson kinematics, thanks to the absence of neutrinos in the final state. Second, it allows for an unambiguous determination of the kinematics of each Z-boson, thanks to the different lepton flavours.
VBF and GGF Higgs production with two jets has already been widely investigated [80][81][82][83], mostly with the purpose of finding kinematic regimes and observables that discriminate between the two signals. In this work, however, we do not aim to separate the Higgs signals, instead, we address the possibility to perform a polarization study of vector bosons in Higgs decay, rather independently of the Higgs-production mechanism. Given the very small fraction of Higgs bosons decaying into four charged leptons, being able to sum over different production channels is crucial to enhance the experimental sensitivity to the polarization structure.
As a last comment, we choose to work in the large-m t limit for the GGF signal. Given the purposes of this work, it represents a satisfactory approximation.
This paper is organised as follows. In Sect. 2 we show the details of the matrix-element reweighting method that we use to generate polarized cross-sections and distributions. The fiducial setup we employ for numerical simulations is given in Sect. 4. In Sect. 5 we validate the reweighting technique comparing its results with those obtained directly simulating polarized events [70]. The polarized results for VBF and GGF signals are shown and discussed in Sect. 6. In Sect. 7 we draw our conclusions.

Matrix-element reweighting
In Fig. 1 we show the general structure of Higgs-mediated amplitudes contributing to a reaction producing the Higgs-decay products (four charged-leptons) and an arbitrary number of jets.
We notice that these amplitudes (which we call signal contributions) are only a sub-set of all amplitudes for the production of four charged leptons in association with jets. Therefore, they are not gauge invariant per sè. Nonetheless, we have checked numerically that with a sharp but realistic cut on the invariant mass of the four leptons, e.g. |M 4 − M h | < 5 GeV, the s-channel Higgs contributions are largely dominant and sufficient to describe the fullmatrix-element calculation up to few-permille effects.
We observe that for this work it is important that the Higgs-decay products can be unambiguously identified and distinguished from the additional jets produced in association with the Higgs boson in order to minimize the impact of irreducible backgrounds and reconstruction effects.
Provided we define vector-boson polarizations in the Higgs rest frame, we can easily parametrize the Higgsdecay amplitude as a sum of polarized and interference terms [79]: where the first three terms in the sum correspond to a definite polarization state λ for both bosons, longitudinal (λ = L), right-handed (λ = +) and left-handed (λ = −), while the fourth and fifth terms are contributions from longitudinal-transverse interference, and the last one comes from left-right interference. Note that, since polarization vectors are defined in the Higgs-boson rest frame there are no mixed-state contributions, i.e. A λλ = 0 for λ = λ . The analytic expression for each term of the sum of Eq. 1 reads In Eqs. 2-7 we have defined the P (Q i ) propagator factor as depending on the Higgs-to-gauge-boson coupling (g hZZ ), the Z pole mass and width, and the invariant mass of the i-th Z boson (Q i ). We have also introduced the factor where Q is the Higgs-boson invariant mass. The c L , c R parameters represent the SM left-and right-chirality coupling of the Z boson to massless leptons. The variables θ * 1(2) , φ * 1(2) are the positively-charged-lepton decay angles computed in the corresponding Z-boson rest frame, with respect to the boson flight direction in the Higgs rest frame, and ∆φ * 12 = φ * 1 − φ * 2 . Instead of considering the left-handed and righthanded contributions separately, we combine them in a single transverse (T) contribution which includes also the left-right interference term, The definition of the transverse mode as a coherent sum of the left and right modes minimizes the interference effects, which now come only from longitudinal-transverse terms in Eqs. 5-6. Moreover, it is known [79] that in the Higgs decay to four charged leptons the longitudinal-transverse interference is only a few percent of the total result, much smaller than the interference between the right and left amplitudes included in Eq. 10. Introducing the following longitudinal or transverse weights, it is possible to compute any polarized distribution, with arbitrary cuts on the lepton kinematics, multiplying the weight of each unpolarized event by the factor in Eq. 11. In this way it becomes unnecessary to generate separately the individual polarized contributions. All the subtleties related to sampling events with negative weights are avoided. We stress that this procedure is accurate only because of the special (and relatively simple) Higgs-decay analytic structure. Applying the same approach to general multiboson processes would require building weights that depend on non-factorized amplitudes, which would be much more time-consuming and equivalent to directly generate events with polarized amplitudes.
This method is designed for the LO Higgs decay into four leptons. However, it could be extended to NLO, provided that it is possible to disentangle polarized contributions to the Z boson propagators. The NLO QCD corrections would only affect the production mechanism, therefore this entire formalism can be extended with no modifications to this perturbative order. However, we do not include QCD radiative corrections to the specific production processes that are considered in this paper, as this would not affect the main results of this work. The extension of this reweighting method to NLO EW corrections is more involved, in particular for the virtual contributions. These corrections are of the order of 2% [24]. The description of weak boson polarizations is not trivial if EW corrections are included, since a standard narrow-width or double-pole approximation [84] is not viable for vector bosons from a 125-GeV Higgs-boson decay. Below the 2M V threshold the NLO EW corrections can be computed via an improved-Born approximation [85], which is only valid for M VV ≈ 2M V − nΓ V , with n = 2 ÷ 3. However, a rigorous description of weak boson polarizations in Higgs decay at NLO EW is possible in the single-pole approximation, namely projecting on mass shell only one of the two vector bosons. This is discussed in Sect. 3.

Single-pole approximation
In this section we briefly address the possibility of studying Z boson polarizations in Higgs decay including NLO EW effects. The presence of non-factorizable corrections both in the real and the virtual contributions makes the separation of polarizations of off-shell Z bosons not well defined, as several diagram topologies do not feature two s-channel Z-boson propagators. Such diagrams cannot be simply dropped as the result would not be gauge invariant. The pole approximation [86,87] can help in addressing this issue. Since on the Higgs resonance one of the two Z bosons is off its mass shell, it is only possible to project on-shell one of the two bosons. This method consists in: selecting only diagrams which feature at least a resonant Z boson decaying into two leptons, projecting the momenta of the s-channel Z-boson decay products such that their sum gives the momentum of an on-shell Z momentum, computing the numerator of the amplitude with the on-shell-projected kinematics, while evaluating the denominator with the original kinematics, to retain the Breit-Wigner modulation of an off-shell Z boson.
The method is accurate in the vicinity of the resonance pole mass [85][86][87], and is characterized by an intrinsic un- . This procedure requires a careful treatment of non-factorizable QED corrections [85,88,89], because of spurious infrared singularities that arise when setting the Z-boson momentum on its mass shell. Furthermore, the presence of photon contributions in the virtual corrections makes it essential to constrain the decay leptons to have an invariant mass that is as close as possible to the Z pole mass. This requirement has the advantage of making the Z-resonant contributions more dominant over the non-resonant ones, giving a more accurate pole-approximated description of the process. This induces the loss of signal events, due to the tighter constraint on the two leptons decaying from the almost-on-shell Z boson. We have calculated that with realistic kinematic selections (see LEP setup of Sect. 4), the additional constraint |M e + e − − M Z | < 5 GeV reduces the VBF total cross-section by 70%. Most of the effect comes from the fact that only half of the events contain an almost-on-shell e + e − pair. However, one should include also the contribution from almost-on-shell µ + µ − pairs. The choice of the specific on-shell projection is not unique [70,71,84,85], and introduces an artificial modification of the momenta, which depends on the physical quantities that one chooses to conserve. In Higgs decay, the choice of such a projection is even more delicate than in other processes, as changing the momenta of one Z boson could induce a shift in the total momentum of the Higgs boson itself, resulting in a bad phase-space sampling in the region of the Higgs resonance and in large discrepancies w.r.t. the full calculation. Therefore, an essential requirement for the on-shell projection is to preserve the total momentum of the four leptons. A viable projection for the process h → Z 1 (→ e + e − ) Z * 2 (→ µ + µ − ) preserves: the spatial direction of e ± in the rest frame of the e + e − system, the spatial direction of µ ± in the rest frame of the µ + µ − system, the spatial direction of the e + e − system in the rest frame of the four-lepton system, the invariant mass of the µ + µ − system, the four-momentum of the four-lepton system. This choice preserves the decay angles of the leptons in the corresponding Z-boson rest frame, minimizing the effect of the pole approximation on polarization-sensitive variables. This specific on-shell projection only works if the invariant mass of the Higgs boson is larger than the sum M Z + M µ + µ − , giving a decrease in the total crosssection, which is of the same order of magnitude as the intrinsic uncertainty of the pole approximation.
This approach, despite some technical details which must be properly taken care of, is expected to give a reliable description of the Higgs-boson decay into four leptons at NLO EW, in the case where one of the two lepton pairs is close to the Z-boson mass shell. The single-pole approximation allows to select in a gauge invariant way only resonant contributions and therefore to reliably separate the polarizations of a single boson at the amplitude level. Since in the SM the Higgs couples to two weak bosons with like-wise polarization mode, selecting the polarization mode of the on-shell Z boson intrinsically gives important information about the polarization of the off-shell boson, which can be then studied by means of the usual angular observables of its decay leptons.
We conclude that the single-pole approximation represents a viable procedure to extend the LO polarization studies on the Higgs resonance to higher orders, in particular for NLO EW corrections. We leave this for future investigations.

Setup
We now proceed to the results for the parton-level process at the LHC with 13 TeV centre-of-mass energy. We have computed VBF and GGF signals at LO, the former simulated in the SM, the latter in the large-m t approximation, with the same numerical setup. A unit CKM matrix is assumed in both processes. We use NNPDF3.0 parton distribution functions (PDF) [90] computed at LO with α s (M Z ) = 0.118, via the LHAPDF interface [91]. The complex-mass scheme [92,93] is understood for the treatment of electroweak-boson masses and couplings in the SM. The pole masses and widths of weak bosons and of the Higgs are set to the following values: The electroweak coupling α is computed in the G µ scheme, with the Fermi constant set to G µ = 1.16637·10 −5 GeV −2 .
For both GGF and VBF, we work in the five-flavour scheme, i.e. including partonic processes with external massless b-quarks. These contributions account for less than 2% of the total cross-section. We use a dynamical factorization and renormalization scale, that is a typical choice for GGF production [42,43]. We have verified that the typical scale choice used for on-shell Higgs production [94] gives results for VBF which differ by a few percent (both at the integrated and differential level) from those obtained with the dynamical choice in Eq. 14. We stress that the scale choice does not affect the main results of our work, as the reweighting procedure detailed in Sect. 2 is independent of the Higgs-production mechanism and its corresponding central-scale choice. We have employed two different setups. The first setup (label INC) is inclusive in lepton kinematics, and is used for validation purposes: -jets with minimum tranverse momentum p T,j > 25 GeV and maximum pseudo-rapidity |η j | < 4.5; -a two-jet system with rather large invariant mass, M jj > 300 GeV, large pseudo-rapidity separation, |∆η jj | > 2.5, and such that η j1 · η j2 < 0; -a four-lepton invariant mass close to the Higgs pole mass, |M 4 − M h | < 5 GeV; -two pairs of opposite-charge leptons of different flavours, with M + − > 10 GeV, for = e, µ.
The requirements on the leading-jet pair are somewhat milder than those used for on-shell Higgs production in VBF [94] and slightly stronger than those used in ZZ scattering studies [95]. The second setup (label LEP) includes, in addition to the ones of the INC setup, the following cuts on the lepton kinematics: -for all charged leptons, p T, ± > 20 GeV, |η ± | < 2.5.
Note that the effect of the lepton p T and η cuts is dramatic, as they decrease the signal cross-section roughly by a factor of 10, with respect to the INC setup. These cuts are the same used in the study of polarized Z-boson scattering performed in ref. [71]. In recent experimental results on Higgs decay to four charged leptons [20,96] the fiducial lepton cuts are slightly looser than those used in this work, therefore we expect a milder effect on polarized results in a realistic experimental analysis than in the present phenomenological setup.
For the simulation of the unpolarized VBF signal we have used the Phantom Monte Carlo [97], which enables the selection of diagrams with an s-channel Higgs exchange from the complete set of tree-level, electroweak diagrams in the SM [98]. The same process has been simulated also with MadGraph [99], with an agreement to better than 0.5% both in the total cross-section and in all analyzed differential distributions. The GGF signal has been simulated with MadGraph (version 2.7.3 [77]), in the large-m t approximation [29], using the spin-correlated decay chain for the s-channel Higgs boson. We observe that the off-shell-ness of the Higgs boson is preserved both in MadGraph and in Phantom. The numerical integration is carried out by means of a Breit-Wigner phase-space mapping restricted to the region |M 4 − M h | < 5 GeV, which is also the kinematic region determined by the selection cuts.

Validation of the method
In order to validate the matrix-element-reweighting method (MERM), we first compare the polarized cross sections computed with the MERM with those extracted from Phantom in VBF. In the first case, we reweight each unpolarized event as described in Sect. 2. In the second case, we use the approach of polarized amplitudes that has been already applied to VBS and di-boson production [70][71][72][73][74]. Differently from previous studies in VBS   Table 1.
In the absence of lepton cuts, polarized cross sections can also be extracted from unpolarized angular distributions of the charged leptons, via appropriate projections onto Legendre polynomials [70]. The corresponding results are presented in the third column of the INC section of Table 1. They agree with those computed by other means at the sub-percent level. In the LEP setup the cross sections computed with the MERM and those obtained from polarized amplitudes are in excellent agreement. This check ensures that separating polarizations of off-shell Z bosons from Higgs decay is well defined.
The reweighting method works very well also at the differential level, as can be appreciated in Fig. 2, where we consider four kinematic variables, in the LEP setup.
The polarized distributions are perfectly reproduced by the reweighting both for the decay angle θ * e + [ Fig. 2(a)] that appears in Eqs. 2-7, and for the azimuthal difference ∆φ e + µ + computed in the Higgs rest frame [ Fig. 2(b)] with the z-axis defined by the Higgs direction of flight in the laboratory. Notice that this angular separation does not coincide with ∆φ * e + µ + in Eqs. 2-7, where the z-axis is defined by the Z boson direction of flight. The invariantmass variables considered in Figs. 2(c)-2(d) are also perfectly described by the reweighting procedure. We stress that the impressive agreement is motivated by the fact that the reweighting factors of Eq. 11 are fully-differential. Therefore not only the variables appearing in the amplitude parametrization of Eq. 1, but any variable depending on the lepton momenta is well modelled, independently of the selection cuts that are applied.
We have successfully validated the method also for other observables, giving us confidence that the matrix- element reweighting procedure furnishes a useful tool to generate polarized observables from unpolarized ones.

Results
In this section we only consider the LEP setup defined in Sect. 4. Before starting the discussion on the VBF and GGF signals, it is worth commenting on the impact of some irreducible EW and QCD backgrounds. The evaluation of such backgrounds to the Higgs signal is in fact an important step in any polarization study. We have computed with Phantom the full process which, at LO order [O(α 6 )], receives contributions from several diagram topologies in which no s-channel Higgs propagator is present. However, the tight but realistic cut on the four-lepton invariant mass (|M 4 − M h | < 5 GeV), introduced in sec4, suppresses these non-signal contributions. The contribution from Higgs-strahlung is suppressed by the large jet-pair invariant-mass cut. The impact of all non-signal contributions on the total crosssection is below 0.5%. The fiducial cross-section of the QCD background, computed at LO [O(α 2 s α 4 )] is 2.6% of the VBF signal. The signal cross-sections in GGF and VBF are shown in Table 2. The two LO cross-sections are of the same size. Given the very low statistics of the process it is essential to sum over them, provided that they show the same polarization structure. The sizable scale uncertainties that characterize the GGF channel [6.988(5) +37.7% −25.4% · 10 −6 pb] make it important, in the lights of a realistic analysis, to include higher-order QCD corrections. The LO scale uncertainties in VBF production are much smaller [8.566(5) +3.9% −3.8% · 10 −6 pb] and the NLO QCD corrections are expected to reduce them further. As already pointed out, the MERM applies, in the case of a Higgs-boson decay to four leptons, in exactly the same way whether radiative QCD corrections are included or not. Table 2 shows that the polarization fractions (polarized cross-sections over the unpolarized one) amount to 65% for the longitudinal mode, 35% for the transverse one. These fractions are almost the same for the two signals, showing that in the SM the polarization content only depends on the decay of the Higgs, and not on its production mechanism. The small differences in the fractions can be traced back to the slightly different kinematics for the two signals, which implies different effects of the selection cuts. The interferences are very small, in spite of quite tight lepton cuts.
The results of Table 2 are not directly related to the SM Higgs coupling to vector bosons (g hVV ), which is independent of the polarization mode of the weak bosons. Therefore, a modification of g hVV due to new-physics effects would result in a enhanced (or diminished) crosssection at the unpolarized level, but one would not expect a strong modification of the relative weight of the longitudinal and the transverse mode. However, as shown in ref. [6], if a new-physics model allows for different Higgs couplings to longitudinal and transverse vector bosons, the effect on the polarization fractions would be relevant. In addition, since the GGF cross-section depends on g 2 hZZ , while the VBF one depends either on g 4 hZZ or on g 2 hZZ g 2 hWW , polarization-dependent Higgs couplings to weak bosons could imply sizeable differences in the polarization structure of the GGF and the VBF signals. Therefore, a good understanding of both signals within the SM is the first step towards the search of beyond-SM effects via the polarization of weak-bosons.
The integrated results are not enough to fully characterize the polarization structure in the Higgs-boson decay. The study of differential distributions is essential to identify LHC observables that are capable of discriminating among polarization states. Furthermore, the interference terms could give noticeable shape distortions, in spite of an integrated contribution which is close to zero.
In Figs. 3-5 we compare the distributions for different polarization states, separating the GGF and VBF production mechanisms. We present the results in terms of differential cross-sections (upper panels) and normalized shapes (lower panels).
In Fig. 3(a) we consider the cosine of the positron decay angle θ * e + in the corresponding Z-boson rest frame, that is directly related to the polarization mode of the weak boson (see Eq. 1). This angular variable can be directly reconstructed at the LHC, thanks to the final state with four charged leptons. The VBF and GGF distributions have the same shape, both in the LL and in the TT component. The interferences play a negligible role for this distribution, as at the integrated level. The symmetric LL shape features a maximum at cos θ * e + = 0 and a minimum in the (anti)collinear regime, similar to the corresponding distribution in the INC setup (dσ/d cos θ * e + ∝ f LL (1−cos 2 θ * e + )). The TT distributions has constant convexity which is very close to zero, but with opposite sign w.r.t. the LL one. The noticeable difference in shape between the LL and the TT distributions makes this variable well suited for polarization discrimination.
In Fig. 3(b) we consider the azimuthal difference between the two tagging jets, computed in the laboratory frame. The GGF distribution shape is sizeably different from the VBF one, as expected from the different kinematics of the forward-backward jets in the two signals. In fact, this variable has been used to discriminate the GGF signal from its largest EW and QCD backgrounds in the h → WW * decay channel [80]. The LL and TT distributions in a given signal do not show relevant dif-  ferences, as the kinematics of the decay leptons is just mildly affected by the different kinematics of the production part of the amplitudes. This azimuthal difference can be useful for polarization measurements only in combination with other variables that discriminate among polarization modes, and if it is needed to separate different signal processes.
In Fig. 4 we consider two azimuthal-difference variables that concern the charged-lepton kinematics.
The difference between the two azimuthal decay angles of the positively-charged leptons ∆φ * e + µ + is considered in Fig. 4(a). This variable appears in the squaredamplitude parametrization of Eq. 1, and gives a modulation to the interferences terms of Eqs. 5-7. In the INC setup the LL distribution is flat, while the TT one it is dominated by a cos(2∆φ * e + µ + ) functional dependence, since the interference between left and right modes is included in the cross-section [79]. In the LEP setup, the behaviour of the TT distribution is quite similar to the one in the INC setup, while the LL shape is no more flat. The non-symmetric character of the unpolarized distributions is due to the longitudinal-transverse interferences that are negative for ∆φ * e + µ + < π/2, positive otherwise. This effect (at most of 2%) reflects, even in the presence of lepton cuts, the cos(∆φ * e + µ + ) functional dependence of the interference terms of Eqs. 5-6. The GGF and VBF signals show exactly the same behaviour in all polarized contributions and interference terms. The ∆φ * e + µ + variable can be easily reconstructed at the LHC, with the considered final state, in the same fashion as cos θ * e + .
In Fig. 4(b) we consider the azimuthal difference ∆φ e + µ + between the positron and the antimuon computed in the Higgs-boson rest frame. This angular variable is related to ∆φ * e + µ + , with the difference that ∆φ e + µ + does not depend on angles computed in the Z-boson rest frame, but only on the kinematics of the four leptons in the Higgsboson rest frame. The shapes of the polarized distributions is markedly different from those for ∆φ * e + µ + . The interferences are very small and negative for ∆φ e + µ + < 3π/4, positive and slightly larger in size for ∆φ e + µ + > 3π/4. Also this variable enables a clear discrimination between the LL and TT modes. The LL shape is monotonically increasing from 0 to π, while the TT one has an absolute minimum at π/2, an absolute maximum at π and another maximum at around 0. The impressive similarity of the SM polarized shapes for the GGF and VBF signals confirms that, if the focus is put on angular variables describing the Higgs-boson decay, it is safe to sum over production mechanisms.
As pointed out in ref. [79], not only angular variables but also invariant-mass observables are suitable for polarization discrimination. In Fig. 5 we consider the invariant mass of the e + e − and e + µ + pairs.
In the first case ( Fig. 5(a)), the reconstructed Z-boson mass is peaked below 40 GeV and at the pole mass, as at the Higgs resonance one Z is typically on-shell while the other is far off-shell. Interferences are almost negligible for this observable. The difference between the LL and the TT normalized shapes concerns only the off-shell region below 40 GeV: the longitudinal curve peaks between 20 and 25 GeV, while the transverse one has a narrower peak around 30 GeV. The GGF and VBF curves for a definite polarization state are almost identical, apart from mild differences in the TT curve around its two peaks.
The invariant mass of the two positively-charged leptons is shown in Fig. 5(b). The GGF and VBF signals behave in the same way, as for most of the other analyzed observables, giving distribution shapes that are almost independent of the production mode. The unpolarized and polarized distributions feature a maximum around M Z /2, but the LL shape is narrower than the TT one, making this variable quite suitable for polarization discrimination. The polarization fractions (both for GGF and for VBF) show that the transverse mode gives a larger cross-sections than the longitudinal in the soft region (M e + µ + < 20 GeV), while the LL mode is larger in the rest of the spectrum. The interference pattern is slightly more evident than for M e + e − . Its size is at most 2% of the total at moderate masses (M e + µ + ≈ 70 GeV). Comparing Fig. 5(b) with Fig. 9 in ref. [79] one notices that different selection cuts can produce significant distortions in the observed distributions.
As a last comment, we note that in the SM not only the polarized shapes but also the relative fraction of longitudinal and transverse modes are almost independent of the production channel. This results in unpolarized distributions shapes (gray curves in bottom panels of Figs. [3][4][5] that are also independent of the production channel.

Conclusions
In this paper we have studied the polarization of Z bosons decaying from a Higgs boson produced in association with two jets, in a vector-boson-fusion kinematic regime.
We have proved that, thanks to the simple analytic structure of the Higgs-decay Standard-Model amplitude, it is possible to avoid generating separately polarized event samples by simply reweighting unpolarized events with fully-differential weights.
We have considered the two main channels that give contribution to Higgs+2j production, namely gluon-gluon fusion and vector-boson fusion. The polarized signals show the same behaviour in the two channels, both at the level of polarization fractions, and at the level of the shape of polarized distributions. This allows to sum over production channels. This holds for a SM measurement, and for modelling the background for the search of beyond-the-Standard-Model effects.
The possibility of extending this work to higher-orders in perturbation theory is also addressed. The usage of a single-pole approximation enables the description of the Higgs decay to polarized bosons in the presence of EW radiative corrections.
The results presented in this paper provide general techniques to study polarizations of vector bosons from the decay of a scalar Higgs boson produced in any channel at hadron colliders. The extension of such methods to an effective-field-theory framework allows for a modelindependent assessment of the polarization structure at the Higgs resonance in the presence of new-physics effects.