Higgs boson decay into four leptons at NLOPS electroweak accuracy
Abstract
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 NexttoLeading Order (NLO) electroweak corrections to Higgs boson decay into four charged leptons, by considering the goldplated 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, spinparity 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 centreofmass (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 ATLASCMS measurement.
Concerning the main production mechanisms of the SM Higgs boson at hadron colliders, i.e. gluongluon 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 downtype 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 nonstandard 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 spinparity 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 offshell 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
nexttoleading
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, Hto4l
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 higherorder photonic corrections are taken into account in terms of QED collinear Structure Functions. A preliminary study of the impact of the gaugeinvariant 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 Nexttoleading 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 neutralcurrent 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 onshell renormalization scheme), while real photon corrections are induced by the bremsstrahlung process . The two contributions are separately infrared (IR) divergent but their sum is IRfinite. 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)
(1)  
where is the Higgs mass, is a softhard energy separator (), is the amplitude of the lowestorder (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
(2) 
where is the oneloop 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 crosscheck of the
calculation,
we tested that the
inclusive NLO QED correction coincides with ,
which is correctly twice
the inclusive finalstate electromagnetic correction to the
decay [Montagna:1993mf].
An important comment is in order here. The treelevel 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.
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 HooftFeynman gauge, all the particles present in the spectrum of the SM, including the FadeevPopov and HiggsKibble ghosts, are involved in the calculation. The corresponding Feynman diagrams include, in addition to twopoint functions, ranktwo tensor three, four and fivepoint 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 fivepoint integrals according to Refs. [Denner:2002ii, Denner:2010tr], as well as according to PassarinoVeltman 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 onshell scheme and for the basic selfenergy 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 treelevel 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 nonIR 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 fermionloop 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 oneloop weak corrections is obtained as
(3) 
where is the oneloop amplitude associated to the full set of weak diagrams.
To check some relevant ingredients of our calculation of oneloop weak corrections, we compared our predictions for the Higgs decays at NLO electroweak accuracy with those of Ref. [Kniehl:1993ay], finding perfect agreement.
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 DrellYan 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:
(4) 
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 cutoff 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
(5) 
(6) 
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 AltarelliParisi vertex for the branching , generates the infrared logarithms. It is explicitly given by:
(7) 
In the definition of , the integral is performed over the angular variables of , and equals if is an antifermion 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 multichannel 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 form
(8) 
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 AltarelliParisi 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:
(9) 
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 onephoton 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.
GeV  GeV  

GeV  GeV  GeV 
KeV  MeV  GeV 
MeV  GeV  GeV 
MeV  MeV  GeV 
The and are the runningwidth PDG values which have to be converted to the fixedwidth 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
(10) 
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.
3.1 NLO electroweak corrections: comparisons to Prophecy4f
A sample of the Prophecy4f vs. Hto4l comparisons at NLO electroweak accuracy is shown in Tab. 2 and in Figs. 13, 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 ,
/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) 
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 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 spinparity assignment. For the angle we use the definition
(11)  
(12) 
where and , , , are the threemomenta of the finalstate 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. 46. In order to disentangle the impact of the different sources of correction, we consider the results obtained according to the following levels of accuracy:

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

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

the complete NLO electroweak calculation;

the NLO QED calculation, given by the gaugeinvariant subset of electromagnetic contributions within the full set of electroweak corrections;

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

the NLO QED corrections matched to the QED PS, i.e. the QED gaugeinvariant realization of Eq. (9).
The comparison between approximations 1. and 2. is useful to quantify the higherorders 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 nonlogarithmic 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 processindependent, 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.
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 nonlogarithmic 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 higherorder photonic corrections. Multiple photon emission is of the order of 10% for the finalstate and at the level of some per cents for the case, as a consequence of the different magnitude of the leptonphoton 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.
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 nonnegligible contribution from fixedorder nonlogarithmic 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%.
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 threemomentum of the with the effective momentum for each photon satisfying the condition , as typically done by LHC experiments, where is the leptonphoton 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, nonnegligible 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 higherorder 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. 79
we show a sample of illustrative results obtained
by interfacing Hto4l with POWHEG [Frixione:2007vw] for the simulation of Higgs boson
production in gluongluon 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 protonproton collisions at a c.m. energy of 8 TeV

generate unweighted events for the process in the Les Houches format using POWHEG, where is an onshell 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 onshell Higgs production and decay.
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 fourlepton 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 threemomenta 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 dashdotted 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.
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 oneloop 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 crosschecked 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 higherorder 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 fixedorder 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 higherorder 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, spinparity 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 PITNGA2010 264564 (LHCPhenoNet), and by the Italian Ministry of University and Research under the PRIN project 2010YJ2NYW. F.P. wishes to thank the CERN PHTH 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 multichannel MC techniques, combined with importance sampling
to reduce the variance of the integral and help event
generation
(13) 
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 multichannel strategy. Here we adopt a different and elegant solution, which consists in writing the phase space in the frame where the leptons are backtoback, 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
(14) 
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
(15) 
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 cutoff 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
Footnotes
 The reference web page is http:/www.pv.infn.it/hepcomplex/hto4l.html
 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.
 For the sake of simplicity, we consider the decay , the generalization to or being straightforward.
 Available at http://omnibus.unifreiburg.de/~sd565/programs/prophecy4f/prophecy4f.html
 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.
 Analogous results are valid in the channel, which coincides for the integrated partial width with the final state (apart from negligible mass effects).
 For simplicity, in the present Section we provide results for bare electrons only, i.e. in the absence of leptonphoton recombination effects.
 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.
 Here we consider only the decay , the generalization to 4 identical leptons being straightforward.