Higgs Pseudo Observables and Radiative Corrections

We show how leading radiative corrections can be implemented in the general description of $h\to 4\ell$ decays by means of Pseudo Observables (PO). With the inclusion of such corrections, the PO description of $h\to 4\ell$ decays can be matched to next-to-leading-order electroweak calculations both within and beyond the Standard Model (SM). In particular, we demonstrate that with the inclusion of such corrections the complete next-to-leading-order Standard Model prediction for the $h\to 2e2\mu$ dilepton mass spectrum is recovered within 1% accuracy. The impact of radiative corrections for non-standard PO is also briefly discussed.


Introduction
The decays of the Higgs particle, h(125), can be characterised by a set of Pseudo Observables (PO) that describes, in great generality, possible deviations from the Standard Model (SM) in the limit of heavy New Physics (NP) [1].
The Higgs PO are defined from a momentum expansion of the on-shell electroweak Higgs decay amplitudes. More precisely, the PO relevant to h → 4 are defined by the momentum expansion around the physical poles (due to the exchange of SM electroweak gauge bosons) of the following three-point correlation function where J µ (x) are generic leptonic currents. This expansion encodes in full generality the short-distance contributions to the decay amplitudes in extensions of the SM with no new light states [1]. However, in order to compare this amplitude decomposition with data, also the long-distance contributions due to soft and collinear photon emission (i.e. the leading QED radiative corrections) must be taken into account.
Soft and collinear photon emission represents a universal correction factor [2,3] that can be implemented, by means of appropriate convolution functions (or, equivalently, showering algorithms such as those adopted in PHOTOS [4], PYTH-IA [5], or SHERPA [6]) irrespective of the specific shortdistance structure of the amplitude. 1 In this paper we illustrate how this works, in practice, in the h → 4 case.
We focus our analysis to the h → 2e2µ case, that is particularly interesting for illustrative purposes: the effect of radiative corrections can be implemented by simple analytic formulae, allowing a transparent comparison with numerical methods. As we will show, the inclusion of the universal QED corrections is necessary and sufficient to reach an accurate theoretical description of the Higgs decay spectrum, that recovers the best up-to-date SM predictions in absence of NP.

QED corrections for the h → 4 dilepton spectrum
In this section we describe how leading QED radiative corrections affect the dilepton spectrum of h → 4 decays assuming a generic PO decomposition of the amplitude. As anticipated, we focus our discussion to the case of two lepton pairs with different flavor (h → 2e2µ) and, more precisely, on the double differential lepton-pair invariant-mass distribution The emission of soft and collinear photons leads to infrared (IR) divergences in the h → 4 spectrum. The full structure of such divergences is rather complicated. However, as we have checked by means of an explicit calculation at O(α), such divergences can be factorized in F(m ee , m µ µ ) and can be analyzed separately for each dilepton system. This happens because each fermion current in Eq. (1) carries an overall neutral electric charge.
Working in the limit of massless leptons, we need to introduce two independent IR regulators for soft and collinear divergences. We choose them to be: i) the minimal fraction of invariant mass lost by the dilepton invariant-mass system; ii) the minimal invariant mass of a single lepton plus (collinear) photon (m * ).
We then define the radiator ω(x, x * ), that represents the probability density function (PDF) that a dilepton system retains a fraction √ x of its original invariant mass after bremsstrahlung for a given x * ≡ 2m 2 * /m 2 0 , where m 2 0 is the initial dilepton invariant mass (pre bremsstrahlung). By construction, the kinematical range of x is Keeping only the leading terms for (1 − x) 1 and x * 1, the radiator is The first term, ω 1 , describes the real emission of a photon such that the lepton pair retains a fraction √ x of its invariant mass; the θ -function implements the corresponding IR cutoff. The second term, ω 2 , describes the events in which the soft radiation is below the IR cutoff, as well as the effect of virtual corrections.
We have determined the structure of ω 1 by means of an explicit O(α) calculation of the real emission, while ω 2 has been determined by the condition 1 0 dxω(x, x * ) = 1. The latter condition implies a redefinition of O(α/π), not enhanced by large logs, of the PO characterizing the non-radiative amplitude.
Denoting by m 01 and m 02 the invariant masses of the two dilepton systems before bremsstrahlung, defining fur- where F 0 (m 01 , m 02 ) denotes the non-radiative (tree-level) spectrum [1]. Starting from Eq. (6) we can extract the double differential spectrum after radiative corrections. To this purpose, we first trade x 1,2 for m 1,2 , obtaining: From Eq. (7) we then explicitly extract the double differential decay width by integrating over all the possible physical m 01,02 combinations, determined by the conditions m 01 + m 02 ≤ m h and m 01,02 ≥ m 1,2 / √ x max . In this way we finally obtain: We stress that the result in Eq. (8) includes both real and virtual QED corrections. The latter have been indirectly determined by the normalization condition for ω(x, x * ), that is the same condition applied in showering algorithms [4]. As anticipated, this implies a O(α/π) redefinition of the PO compared to their tree-level values (both within and beyond the SM). In the context of next-to-leading order (NLO) effective field theory (EFT) calculations [8,9], this procedure provides a well-defined condition for the matching between the full EFT calculation of the amplitude and the PO decomposition.

Comparison with full NLO electroweak corrections
In this section we present a comparison of the SM predictions for the h → 2e2µ dilepton invariant mass spectrum obtained using full NLO electroweak corrections [10], and the PO decomposition "dressed" with leading QED corrections, as described above.
The complete SM NLO electroweak corrections to h → 4 have been computed in [10], and the results have been implemented in the Monte Carlo event generator Prophecy4f [11]. We have used Prophecy4f version 2.0 to generate 200 millions weighted events for the recombination mass parameter m * = 1 GeV. We have used the default Prophecy4f SM inputs except for setting the Higgs boson mass to 125 GeV. Prophecy4f adopts the dipole subtraction formalism [12] for  the treatment of soft and collinear divergences, and the socalled "photon-recombination" is applied. In particular, if the invariant mass of a lepton and a photon is smaller than m * , the photon momentum is added to the lepton momentum [10] . As a result, m * coincides with the collinear cut-off introduced in the previous section. In Fig. 1 (left) we show the decay distribution as a function of the dilepton invariant mass normalized to the total decay width for h → 2e2µ in the SM (upper plot) and the ratio between NLO and leading-order (LO) predictions (lower plot). Shown in solid black is our improved prediction obtained by convoluting the leading order distribution, shown in dashed black, with the radiator function as described in the previous section. The PO have been fixed to their SM treelevel reference values (κ ZZ = 1, ε i = 0 [1]). The Prophecy4f predictions within MC uncertainty are shown with red and blue bands for µ + µ − and e + e − invariant mass spectra, respectively.
We list here a series of conclusions that can be derived from this numerical comparison.
-The spectrum obtained with the PO decomposition of the amplitude, "dressed" with leading QED corrections, provides an excellent approximation (within 1% accuracy) to the spectrum obtained with full NLO EW corrections. 2 -The effect of the leading QED corrections can be large, exceeding 10% in specific regions of the phase space. It therefore must be included, in view of a precise datatheory comparison, also when fitting beyond-the-SM parameters. -The PO "dressed" spectrum is obtained setting ε i = 0 (i.e. to their LO SM values). The good agreement with the complete NLO calculation confirms that the O(α/π) redefinition of the ε i is a small effect, with no observable consequences for the h → 2e2µ dilepton invariant mass spectrum.

Implications for New Physics
As shown in Fig. 1 (left), radiative corrections can be sizable and must be included also when going beyond the SM.
Having demonstrated the validity of our QED improved predictions to describe such effects, we are in position to apply the method in the presence of an arbitrary new physics contribution to h → 2e2µ decay as parameterised by generic PO [1]. As an illustrative example, we consider the impact of the leading QED corrections for non-standard values of κ ZZ , ε Ze L , ε Ze R , ε Zµ L , and ε Zµ R .
To draw some general conclusions we analyse three benchmark points, chosen such that the deviations of the total h → 2e2µ decay rate from the SM prediction are always small, 3 but the impact on the spectrum are quite different. The results of the inclusion of QED corrections are shown in Fig. 1 (right). As in the left panels, we plot the dilepton invariant mass distribution normalized to the total rate (upper plot) and the ratio between NLO and LO (lower plot).
The definition of the benchmarks, and the consequences following from the analysis of radiative corrections, are listed below.
Here the deviation from the SM point in the Higgs PO parameter space is small: this benchmark point is compatible with naive power counting in the linear EFT expansion. As a consequence, small deformations in the spectrum are obtained (upper panel) and the relative QED corrections are SM-like (lower panel). In this regime, the leading QED corrections can be directly extracted from the SM result (via an appropriate NLO/LO re-weighting).
Here the deviation from the SM point is sizable, beyond the naive power counting within a generic EFT (both linear and non-linear). However, the PO configuration is such that the deviations from the SM in the spectrum are small. This implies that the relative impact of QED corrections is still SM-like.
In this example we observe a sizable distortion of the dilepton shape (upper panel). As a consequence, the relative impact of the QED corrections is quite different from the SM case (a description of radiative corrections by 3 The dependence of the total rate on the PO can be found in Ref. [13]. NLO/LO re-weighting of the SM result would not provide a good approximation).

Conclusions
The dominant electroweak corrections to h → 4 decays are due to the universal soft and collinear photon emission. As shown in Fig. 1, these can lead to distortions of the dilepton invariant spectrum of O(10%) is specific regions of the phase space. These effects are of the same order as the expected modifications from the SM under the assumption of underlining linear EFT [13]. It is then mandatory to properly incorporate these corrections in a consistent way both within and beyond the SM.
As we have shown in this paper, this can be achieved in general terms within the framework of the Higgs PO. In particular we have shown that: (i) the QED corrected predictions for the h → 2e2µ dilepton invariant mass spectra, with PO fixed to their SM LO values, are in agreement with the full NLO electroweak SM predictions within 1% accuracy; (ii) the QED corrections in the presence of NP can be sizable and significantly different from the SM case.