A Phase space parameterisation and integration

Higgs boson decay into four leptons at NLOPS electroweak accuracy


In view of precision studies of the Higgs sector at the Run II of the LHC, the improvement of the accuracy of the theoretical prediction is becoming a pressing issue. In this framework, we detail a calculation of the full Next-to-Leading Order (NLO) electroweak corrections to Higgs boson decay into four charged leptons, by considering the gold-plated channel , . We match the NLO corrections with a QED Parton Shower (PS), in order to simulate exclusive multiple photon emission and provide novel results at NLOPS electroweak accuracy. We compare our NLO predictions to those of the program Prophecy4f and present NLOPS phenomenological results relevant for Higgs physics studies, with particular attention to precision measurements of the Higgs boson mass, spin-parity assignment and tests of the Standard Model. Our calculation is implemented in a new code, Hto4l, which can be easily interfaced to any generator describing Higgs boson production. As an example, we provide illustrative results for Higgs production and decay in the process using POWHEG with NLOPS accuracy in the production mode.

a]Stefano Boselli, b]Carlo M. Carloni Calame, a]Guido Montagna, c]Oreste Nicrosini c]and Fulvio Piccinini \affiliation[a]Dipartimento di Fisica, Università di Pavia, and INFN, Sezione di Pavia, Via A. Bassi 6, 27100 Pavia, Italy \affiliation[b]Dipartimento di Fisica, Università di Pavia, Via A. Bassi 6, 27100 Pavia, Italy \affiliation[c]INFN, Sezione di Pavia, Via A. Bassi 6, 27100 Pavia, Italy \emailAddstefano.boselli@pv.infn.it \emailAddcarlo.carloni.calame@pv.infn.it \emailAddguido.montagna@pv.infn.it \emailAddoreste.nicrosini@pv.infn.it \emailAddfulvio.piccinini@pv.infn.it \keywordsHadronic colliders, Higgs physics, NLO computations

1 Introduction

With the announcement in 2012 of the discovery of a new particle in the search for the Standard Model (SM) Higgs boson by the ATLAS [Aad:2012tfa] and CMS [Chatrchyan:2012ufa] collaborations at the CERN LHC, particle physics entered a new era. The data collected at the centre-of-mass (c.m.) energies of 7 and 8 TeV have been analyzed by the two experiments in order to establish whether the newly discovered particle is actually the boson predicted in the SM as relic of the mechanism of electroweak symmetry breaking (EWSB) [PhysRevLett.13.321, Higgs1964132, PhysRevLett.13.508, PhysRevLett.13.585, PhysRev.145.1156, PhysRev.155.1554].

The mass of the observed particle has been precisely measured by studying the two cleanest decay channels given by the decays into a photon pair and into four charged leptons. The combination of the two channels and , which have excellent mass resolution and where excesses with large significance are observed [Aad:2014eva, Aad:2014eha, Aad:2014aba, Khachatryan:2014jba, Khachatryan:2014ira, Chatrchyan:2013mxa], presently provides a mass measurement of approximately 125 GeV for each experiment, with a relative uncertainty of better than 0.2% for the combined ATLAS-CMS measurement.

Concerning the main production mechanisms of the SM Higgs boson at hadron colliders, i.e. gluon-gluon fusion, vector boson fusion (VBF), associated production with a massive vector boson and associated production with top quarks, the studies performed at the LHC, based on the analysis of individual production signal strengths for various decay modes, have provided a clear observation of Higgs production through gluon fusion and an evidence for VBF production, with a significance above the level, and for associated production at about  [Khachatryan:2014jba, Aad:2013wqa].

Various tests of the couplings of the new particle to bosons and fermions have been carried out both by ATLAS and CMS collaborations. In particular, the measured ratio of the couplings of the Higgs particle to and bosons, which is an important probe of the EWSB mechanism as fixed by the custodial symmetry, is compatible with the SM expectation and, more generally, no significant deviation from the SM is observed from the coupling strength studies [Aad:2014eva, Khachatryan:2014jba, Aad:2013wqa]. Noticeably, evidence for the direct coupling of the Higgs boson to down-type fermions has been reported through the study of the challenging decay modes Higgs into bottom quarks and leptons [Aad:2015vsa, Chatrchyan:2014vua].

Last but not least, the spin and parity quantum numbers of the discovered particle have been assessed, by means of a systematic analysis of its production and decay processes. The data strongly favor the scalar nature of the observed particle, while rejecting other non-standard hypotheses () or possibility of mixtures at high confidence [Khachatryan:2014kca, Chatrchyan:2013mxa, Aad:2013xqa].

All these measurements marked the start of a new era of precision Higgs physics and were accompanied by an impressive theoretical effort summarized in three CERN reports by the LHC Higgs Cross Section Working Group [Dittmaier:2011ti, Dittmaier:2012vm, Heinemeyer:2013tqa]. These studies, as well as the related theoretical work, are in continuous progress and will continue during the Run II of the LHC at higher energies and luminosity.

In this paper, we focus on the Higgs boson decay into four charged leptons, i.e. , in order to provide novel precision predictions of interest for future studies of the Higgs sector at the LHC. This decay channel plays a particularly relevant rôle, as it provides the cleanest experimental signature, given by a peak in the four lepton mass spectrum on top of a flat and small background. Actually, the decay mode allows to derive a precise mass measurement in the different combinations of lepton final states, to assess the spin-parity quantum numbers using sensitive angular distributions and to perform precision tests of the SM at the level of differential cross sections [Aad:2014tca]. In the off-shell region, the data can also be used to put constraints on the total width of the Higgs boson [Aad:2015xua, Khachatryan:2014iha].

In the light of the above motivations, we compute the full set of next-to-leading order (NLO) electroweak corrections to , with . We match the NLO corrections to a QED Parton Shower (PS), in order to simulate multiple photon emission exclusively and provide final results at NLOPS electroweak accuracy. The calculation is available in an event generator, Hto4l1, which can be interfaced to any code describing Higgs boson production. The PS approach is based on the ideas first presented in Ref. [CarloniCalame:2000pz] for the simulation of the Bhabha scattering process at GeV-scale colliders and later applied to the study of single production in hadronic collisions [CarloniCalame:2003ux, CarloniCalame:2005vc]. The matching procedure is a generalization of the method developed in Refs. [Balossini:2006wc, Balossini:2008xr] for the precision calculation of processes in QED (as encoded in the program BabaYaga@NLO [Actis:2010gg, CarloniCalame:2003yt]) and also implemented in the event generator Horace for the calculation of single hadroproduction processes at NLOPS electroweak accuracy [CarloniCalame:2006zq, CarloniCalame:2007cd].

The NLO electroweak and QCD corrections to  fermions decay processes have been calculated in Refs. [Bredenstein:2006rh, Bredenstein:2006ha] and are available in the Monte Carlo (MC) program Prophecy4f [Bredenstein:2006nk, Bredenstein:2007ec], which is used in the context of Higgs studies at the LHC for the precision calculation of the branching ratios of the decays  fermions. In Prophecy4f higher-order photonic corrections are taken into account in terms of QED collinear Structure Functions. A preliminary study of the impact of the gauge-invariant NLO QED and PS corrections to the determination of the Higgs boson mass in the decay was performed in Refs. [Piccinini:2005iu, CarloniCalame:2006vr].

The article is organized as follows. In Section 2 we describe the details of our calculation, with particular emphasis on the method used for the matching of the NLO electroweak corrections with the QED PS. In Section 3 we present our phenomenological results: in Section 3.1 we show a sample of comparisons between our predictions and those of Prophecy4f as a benchmark of the NLO computation, in Section 3.2 we provide results for various observables at NLOPS EW accuracy, while in Section 3.3 we present the results for Higgs production and decay in the channel obtained in terms of POWHEG [Nason:2004rx, Frixione:2007vw] interfaced to Hto4l. In Section 4 we draw our conclusions.

2 Details of the calculation

2.1 Next-to-leading order (NLO) electroweak corrections

The NLO electroweak corrections to the Higgs boson decay into four charged leptons consist of QED and purely weak contributions. Since the decay is a neutral-current process, the two subsets are separately gauge invariant and can be computed separately as well.

The QED corrections are obtained by attaching a virtual or real photon to each charged lepton leg. They are expected a priori to provide the dominant contribution, as photons which are emitted collinear to a lepton give rise to large logarithmic corrections of the form , where is the lepton mass and some typical energy scale.

The QED virtual corrections comprise vertex and pentagon diagrams (in the on-shell renormalization scheme), while real photon corrections are induced by the bremsstrahlung process . The two contributions are separately infrared (IR) divergent but their sum is IR-finite. We treat the IR singularity according to the standard QED procedure of assigning a small fictitious mass to the photon in the computation of the virtual and real contributions. More precisely, the Higgs decay width associated to the bremsstrahlung correction is separated in two pieces and calculated as follows (in a shorthand notation)


where is the Higgs mass, is a soft-hard energy separator (), is the amplitude of the lowest-order (LO) process and is the matrix element of the radiative decay process , being the 4 leptons plus 1 photon phase space element. In Eq. (1) eikonal factor stands for the analytical expression of the real radiation correction in the soft limit . The integral in the first line can be done analytically (see e.g. ['tHooft:1978xw]) and the one in the second line is performed using standard MC integration with importance sampling.

The QED virtual counterpart is computed according to the following formula


where is the one-loop amplitude associated to the vertex and pentagon diagrams.

We perform the IR cancellation by taking the sum of Eq. (1) and Eq. (2) in the numerical limit . As a cross-check of the calculation, we tested that the inclusive NLO QED correction coincides with , which is correctly twice the inclusive final-state electromagnetic correction to the decay [Montagna:1993mf]. An important comment is in order here. The tree-level amplitude, as well as the amplitude for the real radiation process, contains poles in the phase space, corresponding to the points where the momenta of the pairs and of the system cross the zero of the inverse propagators: or . These poles are avoided considering that the boson is an unstable particle, i.e. its propagator contains the finite -width. This, however, would spoil the IR cancellation between real and virtual corrections of Eq. (1) and Eq. (2), respectively, unless in Eq. (2) the QED virtual corrections are calculated with unstable bosons. The scheme which we adopt for the introduction of the width in the boson propagator, without introducing gauge invariance violations, is the complex mass scheme [Denner:2005fg, Denner:2006ic], which also allows us to include weak loop corrections consistently.2

Concerning the basic features underlying the computation of the complete virtual corrections, we briefly describe the most important aspects in the following. Since we work in the ’t Hooft-Feynman gauge, all the particles present in the spectrum of the SM, including the Fadeev-Popov and Higgs-Kibble ghosts, are involved in the calculation. The corresponding Feynman diagrams include, in addition to two-point functions, rank-two tensor three-, four- and five-point functions. The related ultraviolet divergencies are regularized by means of dimensional regularization. The reduction of the tensor -point functions is carried out by means of the symbolic manipulation program FORM [Vermaseren:2000nd]. The necessary scalar form factors with complex masses are evaluated using Looptools v2.10 [Hahn:1998yk, Hahn:2000jm], which implements the evaluation of the reduction of tensor five-point integrals according to Refs. [Denner:2002ii, Denner:2010tr], as well as according to Passarino-Veltman reduction techniques [Passarino:1978jh]. The form factors are calculated with complex masses and real external squared momenta. This is sufficient for the implementation of the “simplified version of the complex renormalization”, as described in Refs. [Denner:2005fg, Denner:2006ic]. The complete expressions for the counterterms in the on-shell scheme and for the basic self-energy diagrams are taken from Ref. [Denner:1991kt]. Since the collinear singularities associated to the photon becoming collinear with one of the leptons are regulated by the finite lepton mass, the kinematics of the radiative process is calculated including exactly the contribution of lepton masses. In order to allow the cancellation of soft IR singularities, also the tree-level kinematics is calculated with complete lepton mass effects taken into account. In addition, this gives automatically the correct phase space integration boundaries for the diagrams of the virtual contribution where a virtual photon is connected to one external lepton pair. Although the kinematics is treated exactly, the non-IR virtual amplitudes are calculated in the approximation of neglecting finite fermion mass effects (with the exception of the quark Yukawa couplings, e.g. in the fermion-loop Higgs vertex corrections). These contributions are neglected in our calculation as they are irrelevant in view of a target theoretical accuracy of the order of % and their inclusion would make the numerical computation more time consuming.

In formulae, the Higgs width including one-loop weak corrections is obtained as


where is the one-loop amplitude associated to the full set of weak diagrams.

To check some relevant ingredients of our calculation of one-loop weak corrections, we compared our predictions for the Higgs decays at NLO electroweak accuracy with those of Ref. [Kniehl:1993ay], finding perfect agreement.

In conclusion, our predictions for the Higgs boson decay into four leptons at NLO EW accuracy are given by the sum of Eq. (1), Eq. (2) and Eq. (3), supplemented with the necessary renormalization conditions.

2.2 Matching NLO electroweak corrections to QED Parton Shower

In the present section, we sketch our scheme for the matching of the NLO EW corrections with a QED PS. We closely follow the approach already presented and successfully applied to QED processes at low energies and Drell-Yan production at hadron colliders [Balossini:2006wc, Balossini:2008xr, CarloniCalame:2006zq, CarloniCalame:2007cd].

On general grounds, the partial decay width corrected for the emission of an arbitrary number of photons in a PS framework can be written as follows:


where stands for the set of the final state lepton and photon momenta , , , , , (of order ) is the PS approximation to the squared amplitude for the decay , is the exact phase space for the decay and is the Sudakov form factor accounting for unresolved emission, i.e. soft (up to a cut-off energy ) and virtual corrections in the PS approximation. It is understood that the integral over the phase space has a lower limit for the photon energies set to , to ensure the cancellation of the IR divergencies.

The quantities and read explicitly


In Eq. (6), generates the soft/virtual collinear logarithms, including also interferences effects of radiation coming from different charged legs, and , the integral of the Altarelli-Parisi vertex for the branching , generates the infrared logarithms. It is explicitly given by:


In the definition of , the integral is performed over the angular variables of , and equals if is an anti-fermion or if it is a fermion.

The integral over the phase space as in Eq. (5) is performed after choosing a convenient set of independent variables and using multi-channel MC importance sampling techniques to improve the integration convergence and follow the peaking structure of the partial decay width of Eq. (4) to help events generation. The fully exclusive information on final state particles momenta is kept. Details of the implementation are given in appendix A.

Before discussing the inclusion of NLO corrections into Eq. (4), it is interesting to point out that the squared amplitudes with photon emissions are enhanced in regions of the phase space where the photons are soft and/or collinear or where the propagators are resonating. In this perspective, a good approximation to the exact matrix element can be written in the form3:


In the previous equation, is a shorthand for the coupling, label fermion and photon elicities, , and are the photon polarization vectors. is a -dimensional vector whose component is the index of the fermion to which the photon is attached and the sum over denotes all possible ways to share photons among the four fermions. Finally, is the sum of the momenta of the photons, for a given , attached to the electron current ( to the muon current).

Equation (8) is derived from the amplitude for the emission of photons in the soft limit but keeping the dependence on the photon momenta in the propagators. The sum over the elicities of the squared amplitudes of Eq. (8) gives an approximation of the exact squared matrix elements, coherently including also interferences among diagrams. The final step to obtain of Eq. (4) from Eq. (8) is to replace the photon energy spectrum with the Altarelli-Parisi distribution for a better treatment of hard collinear radiation.

Equation (4), with the building blocks described above, can then finally be improved to include exact NLO corrections according to our master formula:


The correction factors and carry the information of the exact NLO calculation: is the sum of the virtual corrections of Eq. (2) and Eq. (3) and soft real correction given by the first line of Eq. (1), is its PS approximation, i.e. the term without any real hard photon of Eq. (4), is the exact one-photon bremsstrahlung amplitude and is its PS approximation.

We want to remark that and are by construction free of collinear and/or infrared logarithms and that the expansion of Eq. (9) exactly coincides with the NLO calculation, without any double counting. Furthermore, Eq. (9) is still fully differential in the final state momenta and can be conveniently implemented in a MC event generator.

Finally, we remark that the NLO virtual and real corrections used in and are strictly defined only for or photon, while in Eq. (9) they are used also when there are additional photons: this requires a mapping of the -photons phase space to or photon phase space. The mapping is implemented in close analogy to the one described in appendix A.2 of Ref. [Balossini:2006wc], and here we do not discuss it in further detail.

3 Numerical results

In the present Section, we show and discuss the numerical results provided by our calculation, as obtained with the new tool Hto4l. First, we show some tuned comparisons with the predictions of the reference code Prophecy4f at the level of NLO electroweak corrections. Then, we present our best predictions for various observables at NLOPS electroweak accuracy, as well as for Higgs production and decay in the presence of NLO QCD and electroweak corrections matched to PS.

The results presented here are obtained using Prophecy4f v2.0.4 In both codes, we use the following set of input parameters

GeV  GeV
 GeV  GeV  GeV
 KeV  MeV  GeV
 MeV  GeV  GeV
 MeV  MeV  GeV
Table 1: Values of the input parameters used in the numerical calculations.

The and are the running-width PDG values which have to be converted to the fixed-width scheme adopted here through, for example, the relations of Eq. (7.2) of Ref. [Bredenstein:2006rh]. As we work in the scheme, for the electromagnetic coupling constant we use the expression


with , in the calculation of the LO width and NLO weak corrections, while we use for the coupling of the photon to the external charged particles.5 The top-quark width is set to the LO prediction in the SM, and a fixed width is employed in all the resonant propagators in the framework of the complex mass scheme.

3.1 NLO electroweak corrections: comparisons to Prophecy4f


[width=0.7]massscaneemm \includegraphics[width=0.7]massscanmmmm

Figure 1: Comparison between the NLO electroweak calculation of Prophecy4f and Hto4l for the decay widths (upper plot) and (lower plot), as a function of the Higgs mass in the range [125, 400] GeV. For each plot, upper panel: absolute predictions in KeV; lower panel: ratio between LO width and NLO EW corrected width.

A sample of the Prophecy4f vs. Hto4l comparisons at NLO electroweak accuracy is shown in Tab. 2 and in Figs. 1-3, in order to check the technical accuracy of our predictions in its different aspects sketched in Sect. 2.1. Generally speaking, we observe very good agreement between our predictions and the independent results of Prophecy4f.

In Fig. 1 we show the comparison for the NLO width in the leptonic decay channels and ,6 as a function of the Higgs mass in the range [125, 400] GeV, together with the relative contribution due to the NLO electroweak corrections where the effect of mass thresholds present in the loop computation is particularly visible. As can be seen, the two calculations perfectly agree. For the sake of clarity and completeness, we quote in Tab. 2 the predictions of the two codes for the decay channels and for three specific values of the Higgs mass: the level of agreement is within the statistical numerical uncertainty which is well below the 0.1% accuracy.

/Final State Prophecy4f Hto4l
125 GeV/ 0.24151(8) 0.24165(2)
140 GeV/ 1.2672(2) 1.2667(1)
200 GeV/ 825.9(1) 825.8(1)
125 GeV/ 0.13324(2) 0.13325(2)
140 GeV/ 0.6713(1) 0.6711(1)
200 GeV/ 413.02(7) 412.98(2)
Table 2: Comparison between the NLO electroweak predictions of Prophecy4f and our calculation (Hto4l) for the Higgs decay width (in KeV), for different values of the Higgs mass and final states. The numbers in parenthesis are the statistical uncertainty on the last digit due to MC integration.




Figure 2: Comparison between the NLO electroweak calculation of Prophecy4f and Hto4l for the invariant mass (in the Higgs rest frame), in the range [60, 100] GeV (upper plot) and in the range [85, 95] GeV (lower plot). Predictions for the decay at  GeV. Upper panels: absolute predictions for ; lower panels: relative effect of the NLO corrections.

In Fig. 2 a comparison between Prophecy4f and Hto4l is shown for the invariant mass (in the Higgs rest frame), in the range [60, 100] GeV (upper plot) and in the range [85, 95] GeV (lower plot). The results refer to the decay channel for  GeV. Also in this case, the agreement between the two codes is remarkable, in spite of the large effect due to the radiative corrections7. Actually, at and above the peak of the electron-pair invariant mass distribution the corrections are of the order of 30%, while for below they can reach 50%. The lowering of the peak and the raising of a tail can be mainly ascribed to the photon radiation off the leptons, as typical final-state radiation (FSR) effect observed around the peak of resonant processes [CarloniCalame:2006zq, CarloniCalame:2007cd, Barze:2012tt, Barze':2013yca].





Figure 3: Comparison between the NLO electroweak calculation of Prophecy4f and Hto4l for the angle distribution (in the Higgs rest frame) for the decay channels (upper plot) and (lower plot) at  GeV. Upper panels: absolute predictions in GeV/deg; lower panels: relative effect of the NLO corrections.

A further comparison is given in Fig. 3 for the distribution of the angle between the decay planes of the virtual bosons in the rest frame for the channels (upper plot) and (lower plot) for  GeV, which is the observable of main interest for spin-parity assignment. For the angle we use the definition


where and , , , are the three-momenta of the final-state leptons.

Again the predictions of the two codes nicely agree. The contribution of the NLO corrections is particularly visible at the edges of the distribution, where it can reach the 5% level for both the decay channels.

3.2 Predictions at NLOPS electroweak accuracy

Some illustrative results obtained according to a number of variants of the theoretical approach described in Sect. 2.2 are given in Figs. 4-6. In order to disentangle the impact of the different sources of correction, we consider the results obtained according to the following levels of accuracy:

  1. the pure PS approximation for the decay width as in Eq. (4), associated to multiple photon emission in the soft/collinear limit;

  2. the truncated approximation of Eq. (4), describing one photon radiation in the PS framework;

  3. the complete NLO electroweak calculation;

  4. the NLO QED calculation, given by the gauge-invariant subset of electromagnetic contributions within the full set of electroweak corrections;

  5. the NLO electroweak corrections matched to the QED PS, as in Eq. (9);

  6. the NLO QED corrections matched to the QED PS, i.e. the QED gauge-invariant realization of Eq. (9).

The comparison between approximations 1. and 2. is useful to quantify the higher-orders contribution due to photon emission beyond , while the difference between options 3. and 4. is a measure of pure weak loop corrections, the difference between approximations 2. and 3. is an estimate of non-logarithmic QED terms plus pure weak loop corrections. The comparison between approximations 3. and 5., as well as between 4. and 6., allows us to check that the NLOPS matching procedure correctly preserves the effect of QED exponentiation as given by the difference between options 1. and 2. Moreover, the results of 1. vs. those of 5. and of 3. vs. those of 5. provide an estimate of the accuracy of the predictions available in the literature for Higgs physics at the LHC, in particular of of the process-independent, widely used code PHOTOS [Barberio:1993qi], which describes multiple photon emission but does not include exact NLO electroweak corrections, and of Prophecy4f, that does not take into account the contribution of exclusive QED exponentiation.





Figure 4: Relative contribution of the QED/electroweak corrections to the (upper plot) and (lower plot) invariant mass in the Higgs rest frame, in the range [85, 95] GeV. Predictions for the decay at  GeV. The theoretical approximations corresponding to the different lines are explained in the text.

In Fig. 4 we show the relative contribution of the different theoretical approximations discussed above for the (upper plot) and (lower plot) invariant mass in the Higgs rest frame, in the range [85, 95] GeV. The results refer to the process for  GeV, according to a bare lepton definition. By inspection of Fig. 4 we can draw the following conclusions: the NLO corrections to the lepton invariant masses are quite large, since they amount to about 50% (30%) to the () invariant mass below the peak and about 30% (20%) at and above it. They are largely dominated by the enhanced leading logarithmic contributions of QED nature , as can be inferred from the comparison between the results of the pure PS algorithm and those of the NLO QED/electroweak calculations. From this comparison, one can also conclude that the non-logarithmic QED terms contribute at the some per cent level, both for the and invariant mass, whereas the pure weak loops have a much smaller effect, not exceeding the 1% level.

The large impact of NLO QED corrections, which significantly modify the shape of the invariant mass distribution, translates in a relevant contribution due to higher-order photonic corrections. Multiple photon emission is of the order of 10% for the final-state and at the level of some per cents for the case, as a consequence of the different magnitude of the lepton-photon collinear logarithm. It can also be noticed that QED exponentiation reduces the impact of NLO corrections and that the NLOPS matching correctly preserves the size of multiple photon emission.





Figure 5: The same as Fig. 4 for the angle distribution in the decay channels (upper plot) and (lower plot) at  GeV.

Quite different conclusions derive from the analysis of Fig. 5, which shows the relative corrections of the different theoretical recipes on the angle distribution for the and decays. For such an observable, the pure PS approximation significantly underestimates the contribution of NLO EW corrections for close to and , while it provides an overestimate around . Actually, it can be noticed that the angle distribution receives a non-negligible contribution from fixed-order non-logarithmic terms and that, more importantly, is particularly sensitive to pure weak corrections, which set the correct overall size and shape of the radiative corrections. On the other hand, the effect of QED exponentiation is moderate, varying between a few per mille to about 1%.





Figure 6: Relative contribution of the QED/electroweak corrections to the invariant mass (upper plot) and the angle distribution (lower plot) for recombined electrons and positrons. Predictions for the decay at  GeV in the Higgs rest frame. The theoretical approximations corresponding to the different lines are explained in the text.

For completeness, we show in Fig. 6 results for the invariant mass of the pair and the angle distribution (for the process ) under the more realistic experimental condition of calorimetric or recombined electrons and positrons. In this case, we replace the three-momentum of the with the effective momentum for each photon satisfying the condition , as typically done by LHC experiments, where is the lepton-photon separation angle in the transverse plane. As can be seen from Fig. 6 in comparison to Fig. 4 and Fig. 5, the contribution of the radiative corrections is largely reduced, as expected, when switching from bare to recombined electrons/positrons. For the invariant mass, the corrections are reduced by about a factor of three, almost independently of the considered theoretical approximation, and preserve their shape. However, non-negligible corrections still remain under the calorimetric condition, of about +15% in the left tail of the invariant mass and of the order of -10% around the peak of the distribution, when considering the most accurate matched predictions. In comparison to the case of bare electrons, the effect of QED exponentiation for dressed electrons reduces to about 1% in the tail and at the per mille level at and above the peak.

More interestingly, the QED relative corrections to the angle distribution are substantially modified both in size and shape by the recombination effects, whereas the full electroweak predictions receive a slight size reduction and a less pronounced shape modification. In particular, we checked through detailed numerical inspections that the especially visible difference in shape between the pure PS and the diagrammatic QED predictions is of virtual origin and has to be ascribed to the QED pentagons, which are exactly included in the Feynman diagram calculation and only (crudely) approximated in the soft/collinear limit in the PS calculation. To some extent, we expect that the rich angular correlations introduced by pentagon diagrams is only poorly reproduced by the PS approximation. All in all, the results shown in the lower plot of Fig. 6 reinforce the already noticed particularly relevant rôle played by loop contributions with complex topology, both of QED and weak nature, to obtain reliable predictions for the angle observable.

To summarize, the main conclusion of this Section is that both NLO electroweak and higher-order QED corrections, as well as their combination, are relevant for reliable simulations of the most important observables considered in precision studies of the Higgs sector at the LHC.

3.3 Interface to Powheg: results for production and decay

In order to facilitate phenomenological studies of Higgs boson production and decay in the presence of both QCD and electroweak contributions, we have implemented an interface which allows to use our code in association with any event generator describing Higgs production. In Figs. 7-9 we show a sample of illustrative results obtained by interfacing Hto4l with POWHEG [Frixione:2007vw] for the simulation of Higgs boson production in gluon-gluon fusion. We use the POWHEG version with NLOPS accuracy in QCD [Alioli:2008tz] from the POWHEG BOX framework [Alioli:2010xd] and we consider Higgs production in proton-proton collisions at a c.m. energy of 8 TeV8. The events generated by POWHEG are interfaced to Hto4l according to the following procedure:

  • generate unweighted events for the process in the Les Houches format using POWHEG, where is an on-shell Higgs boson and stands for the extra parton of the NLO QCD calculation;

  • the Les Houches file is read event by event by Hto4l and the particles momenta are stored in the generic common block structure introduced in Ref. [Boos:2001cv];

  • each event is decayed into the selected channel in the rest frame, using Hto4l. After boosting the decay products back to the laboratory frame, the events including production and decay are written in a file in the Les Houches format.

The Les Houches file can be finally passed to a shower event generator for QCD showering and hadronization. In our examples we use PYTHIA v6.4 [Sjostrand:2006za] as QCD PS. According to the above procedure, the process is treated in narrow width approximation, as it is the case for a 125 GeV Higgs boson, and factorized in on-shell Higgs production and decay.





Figure 7: Comparison between the results obtained with POWHEG + Hto4l (Born) + PYTHIA v6 (red dashed line) and POWHEG + Hto4l (NLOPS) + PYTHIA v6 (blue solid and black dash-dotted lines) for the transverse momentum (upper plot) and rapidity (lower plot) of the Higgs boson. In the lower panels the relative contribution of NLOPS electroweak corrections for bare (solid line) and recombined (dash-dotted line) leptons is shown.




Figure 8: The same as in Fig. 7 for the invariant mass of the subleading lepton pair (upper plot) and the cosine of the angle of the leading lepton pair in the four-lepton rest frame with respect to the beam axis (lower plot).

In our analysis we consider, for definiteness, the decay channel and the following observables: the transverse momentum and rapidity of the Higgs boson (Fig. 7), the invariant mass of the subleading lepton pairs and the magnitude of the cosine of the decay angle of the leading lepton pair in the four-lepton rest frame with respect to the beam axis (Fig. 8). The leading pair is defined as the lepton pair with invariant mass closest to the boson mass and its angle is obtained by summing the three-momenta of the two leptons. For the POWHEG calculation of Higgs production in gluon fusion, we use the PDF set MSTW2008nlo68cl [Martin:2009iq] with factorization/renormalization scale . The values of the other input parameters are the same as the ones given in Tab. 1. The results shown in the following refer to a sample of unweighted events and to the same selection cuts adopted in Ref. [Aad:2014tca] and correspond to bare (solid line) and recombined (black dash-dotted line) leptons. In the latter case, we recombine photons with both electrons and muons, in analogy to the selection criteria adopted in the experimental study of Ref. [Aad:2014tca], if the condition is satisfied.





Figure 9: Distribution of the transverse momentum of the hardest photon (upper plot) and the angular separation between the hardest photon and the closest lepton (lower plot) obtained using POWHEG + Hto4l (NLOPS) + PYTHIA v6. The minimum photon energy is  GeV.

In Fig. 7 and Fig. 8 we show the comparison between the predictions obtained using POWHEG interfaced to our code at LO and NLOPS electroweak accuracy. It can be noticed that the contribution due to NLOPS electroweak corrections is almost flat and of about % for , and when considering bare (recombined) leptons, while the invariant mass of the subleading lepton pairs receives a varying correction of size between % and % for bare (calorimetric) leptons, respectively.

In Fig. 9 we show the results for two observables which are fully exclusive over QED radiation and which can be easily treated in our approach. The results correspond to the process , with  GeV, for which we show the transverse momentum of the hardest photon and the angular separation between the hardest photon and the closest lepton, that exhibit the expected features of photon emission in radiative events.

4 Conclusions

In this work we have presented a precision calculation of the SM Higgs boson decay into four charged leptons, in view of improved measurements of the properties of the Higgs particle at the LHC Run II. Our approach is based on the computation of the full one-loop electroweak corrections supplemented with the contribution of multiple photon emission taken into account according to a fully exclusive QED PS algorithm. Our results, which have a NLOPS electroweak accuracy, are available in the form of a new event generator, Hto4l, that can be easily interfaced to any QCD program simulating Higgs production.

We have cross-checked our NLO electroweak corrections against the predictions of the reference code Prophecy4f and found perfect agreement. We have also shown that both NLO electroweak and higher-order QED corrections, as well as their interplay, are necessary for actually precise simulations of the variety of observables involved in Higgs physics at the LHC. This provides the main novel theoretical feature of our work, which goes beyond the presently available results limited to the fixed-order approximation or to a leading logarithmic QED modeling. The second relevant aspect is given by the possibility of interfacing Hto4l to any generator describing Higgs boson production, thus allowing simulations of Higgs production and decay in the presence of higher-order QCD and electroweak corrections. In this respect, we have shown some illustrative results obtained in terms of the combined usage of POWHEG and Hto4l.

Our results can find application in precision measurements of the Higgs boson mass, spin-parity determination and tests of the SM at the level of differential cross sections in the future run of the LHC. They can be generalized to other processes yielding four leptons in hadronic collisions, like e.g. or .


This work was supported in part by the Research Executive Agency (REA) of the European Union under the Grant Agreement number PITN-GA-2010- 264564 (LHCPhenoNet), and by the Italian Ministry of University and Research under the PRIN project 2010YJ2NYW. F.P. wishes to thank the CERN PH-TH Department for partial support and hospitality during several stages of the work. The authors thank the Galileo Galilei Institute for theoretical physics and INFN for partial support during the completion of this research.

Appendix A Phase space parameterisation and integration

The bodies phase space as in Eq. (5) is integrated according to standard multi-channel MC techniques, combined with importance sampling to reduce the variance of the integral and help event generation9. The first step is to generate a photon multiplicity and associate () photons to the electron (muon) current (), defining the channel of the multi-channel integration. The phase space is then conveniently split into two decaying objects to follow the propagators, namely


where () are the momenta of the virtual bosons.

We refrain from writing explicitly the simple decay phase spaces of Eq. (13) and we focus instead on the case where at least one photon is present. As discussed in appendix A.3 of Ref. [Balossini:2006wc], an efficient sampling of photons collinear to final state leptons is a non trivial task, because the directions of the leptons are known only after all the momenta are generated. In Ref. [Balossini:2006wc] we adopted a solution based on a properly chosen multi-channel strategy. Here we adopt a different and elegant solution, which consists in writing the phase space in the frame where the leptons are back-to-back, i.e. (see for example [Schonherr:2008av, JadachLect, Jadach:1999vf]).

Omitting overall numerical factors for brevity, the building block we are interested in is

where we defined , and contains the infinitesimal phase space element divided by the final state particle energies. It is usually understood that all the variables are expressed in the frame where is at rest, but we want to express them where is at rest. In order to do that, the previous equation can be further manipulated by inserting the following identities


which help to make explicit the Lorentz invariance of the phase space element.

With the help of Eq. (14) and appropriately rearranging the terms, we can write


In the cascade of identities (15) we used the result (see [JadachLect]) and we made use of Lorentz invariance. In the last identity it is understood that all the variables are expressed in the frame where is at rest and , , is the speed of particle and . The big advantage of the last equation is that the lepton momenta and lie on the same direction defined by and , hence all photons can be generated along this direction to sample the collinear singularities. Once all particle momenta are generated, they can be boosted back to the rest frame of the decaying Higgs boson.

One last remark concerns the integration limits of the phase space. As mentioned in Sect. 2.2, photon energies should be generated larger than the infrared cut-off in the Higgs frame, which is a non Lorentz invariant cut. Since the minimum photon energy can not be determined a priori in the frame where is at rest (because itself depends on the photons momenta), we decide to generate photon energies starting from to cover the whole phase space and then, once boosted back, cut the event if a photon enegy falls below . Finally, in order to flatten the infrared divergence, we choose to sample the photon energies according to the function


  1. The reference web page is http:/www.pv.infn.it/hepcomplex/hto4l.html
  2. Actually, the complex mass scheme is used in our calculation of the weak contributions due to the exchange of bosons and top quarks as well, where a complex top mass is introduced, in particular, to evaluate loop diagrams with internal top quarks when the Higgs mass is close to the threshold.
  3. For the sake of simplicity, we consider the decay , the generalization to or being straightforward.
  4. Available at http://omnibus.uni-freiburg.de/~sd565/programs/prophecy4f/prophecy4f.html
  5. This value is used for all the numerical results shown in the following, with the exception of the comparisons with Prophecy4f, where we use everywhere, to be consistent with the default choice of Prophecy4f.
  6. Analogous results are valid in the channel, which coincides for the integrated partial width with the final state (apart from negligible mass effects).
  7. For simplicity, in the present Section we provide results for bare electrons only, i.e. in the absence of lepton-photon recombination effects.
  8. However, as we are interested to study the relative impact of electroweak corrections dominated by contributions of the kind , the results shown in the following are in practice independent of the c.m. energy.
  9. Here we consider only the decay , the generalization to 4 identical leptons being straightforward.
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Comments 0
Request comment
The feedback must be of minumum 40 characters
Add comment
Loading ...

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description