# Higgs pseudo observables and radiative corrections

• Marzia Bordone
• Gino Isidori
• David Marzocca
• Andrea Pattori
Open Access
Letter

## Abstract

We show how leading radiative corrections can be implemented in the general description of $$h\rightarrow 4\ell$$ decays by means of pseudo observables (PO). With the inclusion of such corrections, the PO description of $$h\rightarrow 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 SM prediction for the $$h\rightarrow 2e2\mu$$ dilepton mass spectrum is recovered within $$1\,\%$$ accuracy. The impact of radiative corrections for non-standard PO is also briefly discussed.

## Keywords

Invariant Mass Radiative Correction Effective Field Theory Standard Model Prediction Dilepton Invariant Mass
These keywords were added by machine and not by the authors. This process is experimental and the keywords may be updated as the learning algorithm improves.

## 1 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\rightarrow 4\ell$$ 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
\begin{aligned} \langle 0 | {\mathscr {T}} \{ J_{\ell }^\mu (x), J_{\ell ^\prime }^\nu (y), h(0) \} | 0 \rangle , \end{aligned}
(1)
where $$J_{\ell }^\mu (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], PYTHIA [5], or SHERPA [6]) irrespective of the specific short-distance structure of the amplitude.1 In this paper we illustrate how this works, in practice, in the $$h\rightarrow 4\ell$$ case.

We focus our analysis to the $$h\rightarrow 2e 2\mu$$ 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.

## 2 QED corrections for the $$h\rightarrow 4\ell$$ dilepton spectrum

In this section we describe how leading QED radiative corrections affect the dilepton spectrum of $$h\rightarrow 4 \ell$$ 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\rightarrow 2e 2\mu$$) and, more precisely, on the double differential lepton-pair invariant-mass distribution
\begin{aligned} F(m_{ee},m_{\mu \mu })=\frac{\text {d}^2\Gamma (h\rightarrow 2e 2\mu ) }{\text {d}m_{ee} \text {d}m_{\mu \mu } }. \end{aligned}
(2)
The emission of soft and collinear photons leads to infrared (IR) divergences in the $$h\rightarrow 4 \ell$$ spectrum. The full structure of such divergences is rather complicated. However, as we have checked by means of an explicit calculation at $$O(\alpha )$$, such divergences can be factorized in $$F(m_{ee},m_{\mu \mu })$$ 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 $$\omega (x, x_*)$$, that represents the probability density function that a dilepton system retains a fraction $$\sqrt{x}$$ of its original invariant mass after bremsstrahlung for a given $$x_* \equiv 2 m_*^2/ m_0^2$$, where $$m_0^2$$ is the initial dilepton invariant mass (pre bremsstrahlung). By construction, the kinematical range of x is
\begin{aligned} 0 < x < x_\mathrm{max} = 1 - x_* . \end{aligned}
(3)
Keeping only the leading terms for $$(1-x) \ll 1$$ and $$x_* \ll 1$$, the radiator is
\begin{aligned} \omega (x,x_*) = \omega _1(x,x_*) \theta (1 - x_* - x) + \omega _2(x,x_*) \delta (1-x),\nonumber \\ \end{aligned}
(4)
where
\begin{aligned} \omega _1(x,x_*)= & {} - \frac{\alpha }{\pi } \left( 1 + x - \frac{2}{1-x}\right) \log \left( \frac{2(1-x) - x_*}{x_*} \right) , \nonumber \\ \omega _2(x,x_*)= & {} 1 + \frac{\alpha }{2\pi } \left[ \frac{\pi ^2}{3} - \frac{7}{2} - 3 \log \left( \frac{x_*}{2}\right) - 2 \log \left( \frac{x_*}{2}\right) ^2 \right] .\nonumber \\ \end{aligned}
(5)
The first term, $$\omega _1$$, describes the real emission of a photon such that the lepton pair retains a fraction $$\sqrt{x}$$ of its invariant mass; the $$\theta$$-function implements the corresponding IR cutoff. The second term, $$\omega _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 $$\omega _1$$ by means of an explicit $$O(\alpha )$$ calculation of the real emission, while $$\omega _2$$ has been determined by the condition $$\int _0^1 \text {d}x \omega (x,x_*) =1$$. The latter condition implies a redefinition of $$O(\alpha /\pi )$$, 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 further $$x_{i}=(m_{i} / m_{0i})^2$$, it is easy to show that
\begin{aligned} \frac{\text {d}^4 \Gamma }{\text {d}m_{01} \text {d}m_{02} \text {d}x_1 \text {d}x_2} = F_0(m_{01},m_{02}) \omega (x_1, x_{1*}) \omega (x_2, x_{2*}), \end{aligned}
(6)
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:
\begin{aligned} \frac{\text {d}^4 \Gamma }{\text {d}m_{01} \text {d}m_{02} \text {d}m_1 \text {d}m_2}= & {} \frac{4 m_1 m_2}{m_{01}^2 m_{02}^2} F_0(m_{01},m_{02}) \nonumber \\&\times \omega \left( \frac{m_1^2}{m_{01}^2}, \frac{2m_*^2}{m_{01}^2}\right) \omega \left( \frac{m_2^2}{m_{02}^2}, \frac{2m_*^2}{m_{02}^2}\right) .\nonumber \\ \end{aligned}
(7)
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}\le m_h$$ and $$m_{01,02} \ge m_{1,2}/\sqrt{x_\mathrm{max}}$$. In this way we finally obtain:
\begin{aligned} F(m_1,m_2)= & {} \int _{\sqrt{m_1^2+2m_*^2}}^{m_h} \text {d}m_{01} \int _{\sqrt{m_2^2+2m_*^2}}^{m_h - m_{01}} \text {d}m_{02} ~\frac{4 m_1 m_2}{m_{01}^2 m_{02}^2} \nonumber \\&\times F_0(m_{01},m_{02}) \omega \left( \frac{m_1^2}{m_{01}^2}, x_{1*}\right) \omega \left( \frac{m_2^2}{m_{02}^2}, x_{2*}\right) . \nonumber \\ \end{aligned}
(8)
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 $$\omega (x,x_*)$$, that is the same condition applied in showering algorithms [4]. As anticipated, this implies a $$O(\alpha /\pi )$$ 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.

## 3 Comparison with full NLO electroweak corrections

In this section we present a comparison of the SM predictions for the $$h\rightarrow 2e 2\mu$$ 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\rightarrow 4\ell$$ 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 so-called “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\rightarrow 2e2\mu$$ 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 tree-level reference values ($$\kappa _{ZZ}=1$$, $$\epsilon _i=0$$ [1]). The Prophecy4f predictions within MC uncertainty are shown with red and blue bands for $$\mu ^+\mu ^-$$ 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 data-theory comparison, also when fitting beyond-the-SM parameters.

• The PO “dressed” spectrum is obtained setting $$\epsilon _i=0$$ (i.e. to their LO SM values). The good agreement with the complete NLO calculation confirms that the $$O(\alpha /\pi )$$ redefinition of the $$\epsilon _i$$ is a small effect, with no observable consequences for the $$h\rightarrow 2e2\mu$$ dilepton invariant mass spectrum.

## 4 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\rightarrow 2e 2\mu$$ 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 $$\kappa _{ZZ}$$, $$\epsilon _{Ze_L}$$, $$\epsilon _{Ze_R}$$, $$\epsilon _{Z\mu _L}$$, and $$\epsilon _{Z\mu _R}$$.

To draw some general conclusions we analyse three benchmark points, chosen such that the deviations of the total $$h\rightarrow 2e 2\mu$$ 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.
• Benchmark I [$$\kappa _{ZZ}=1.3$$, $$\epsilon _{Ze_L}=\epsilon _{Z\mu _L}=-0.05$$, $$\epsilon _{Ze_R}=\epsilon _{Z\mu _R}=0.05$$ (dot–dashed blue)]. 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).

• Benchmark II [$$\kappa _{ZZ}=0$$, $$\epsilon _{Ze_L}=\epsilon _{Z\mu _L}=0.26$$, $$\epsilon _{Ze_R}=\epsilon _{Z\mu _R}=0$$ (dotted red)]. 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.

• Benchmark III [$$\kappa _{ZZ}=0.3$$, $$\epsilon _{Ze_L}=-0.45$$ and $$\epsilon _{Z\mu _L}=\epsilon _{Ze_R}=\epsilon _{Z\mu _R}=0$$ (solid green)]. 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 NLO/LO re-weighting of the SM result would not provide a good approximation).

## 5 Conclusions

The dominant electroweak corrections to $$h\rightarrow 4\ell$$ 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\rightarrow 2e2\mu$$ 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.

## Footnotes

1. 1.

For a discussion about the implementation of universal QED corrections in a general EFT context see also Ref. [7].

2. 2.

The $$\sim$$2 % deviations at the border of the phase space are expected due the breakdown of the approximation $$m_{\ell \ell } \gg m_*$$ employed in the analytic evaluation of the radiation function.

3. 3.

The dependence of the total rate on the PO can be found in Ref. [13].

## Notes

### Acknowledgments

This research was supported in part by the Swiss National Science Foundation (SNF) under Contract 200021-159720.

## References

1. 1.
M. Gonzalez-Alonso, A. Greljo, G. Isidori, D. Marzocca, Eur. Phys. J. C 75(3), 128 (2015). arXiv:1412.6038
2. 2.
D.R. Yennie, S.C. Frautschi, H. Suura, Ann. Phys. 13, 379 (1961)
3. 3.
S. Weinberg, Phys. Rev. 140, B516 (1965)
4. 4.
N. Davidson, T. Przedzinski, Z. Was. arXiv:1011.0937
5. 5.
T. Sjostrand, S. Mrenna, P.Z. Skands, Comput. Phys. Commun. 178, 852 (2008). arXiv:0710.3820
6. 6.
T. Gleisberg et al., JHEP 0902, 007 (2009). arXiv:0811.4622
7. 7.
G. Isidori, Eur. Phys. J. C 53, 567 (2008). arXiv:0709.2439
8. 8.
C. Hartmann, M. Trott. arXiv:1505.02646
9. 9.
M. Ghezzi, R. Gomez-Ambrosio, G. Passarino, S. Uccirati. arXiv:1505.03706
10. 10.
A. Bredenstein, A. Denner, S. Dittmaier, M.M. Weber, Phys. Rev. D 74, 013004 (2006). arXiv:hep-ph/0604011
11. 11.
12. 12.
S. Dittmaier, Nucl. Phys. B 565, 69 (2000). arXiv:hep-ph/9904440
13. 13.
M. Gonzalez-Alonso, A. Greljo, G. Isidori, D. Marzocca, Eur. Phys. J. C 75(7), 341 (2015). arXiv:1504.04018

Funded by SCOAP3.

## Authors and Affiliations

• Marzia Bordone
• 1
• 1
• 3
• Gino Isidori
• 1
• 2
Email author
• David Marzocca
• 1
• Andrea Pattori
• 1
1. 1.Physik-InstitutUniversität ZürichZurichSwitzerland
2. 2.INFN, Laboratori Nazionali di FrascatiFrascatiItaly
3. 3.Faculty of ScienceUniversity of SarajevoSarajevoBosnia and Herzegovina
{"extraction_info": {"found_math": true, "script_math_tex": 0, "script_math_asciimath": 0, "math_annotations": 0, "math_alttext": 0, "mathml": 0, "mathjax_tag": 0, "mathjax_inline_tex": 0, "mathjax_display_tex": 2, "mathjax_asciimath": 0, "img_math": 0, "codecogs_latex": 0, "wp_latex": 0, "mimetex.cgi": 0, "/images/math/codecogs": 0, "mathtex.cgi": 0, "katex": 0, "math-container": 0, "wp-katex-eq": 0, "align": 0, "equation": 0, "x-ck12": 0, "texerror": 0, "math_score": 0.985479474067688, "perplexity": 1825.4434292039298}, "config": {"markdown_headings": true, "markdown_code": true, "boilerplate_config": {"ratio_threshold": 0.18, "absolute_threshold": 10, "end_threshold": 15, "enable": true}, "remove_buttons": true, "remove_image_figures": true, "remove_link_clusters": true, "table_config": {"min_rows": 2, "min_cols": 3, "format": "plain"}, "remove_chinese": true, "remove_edit_buttons": true, "extract_latex": true}, "warc_path": "s3://commoncrawl/crawl-data/CC-MAIN-2019-09/segments/1550247484928.52/warc/CC-MAIN-20190218094308-20190218120308-00199.warc.gz"}