Techniques for the treatment of IR divergences in decay processes at NLO and application to the top-quark decay

We present the extension of two general algorithms for the treatment of infrared singularities arising in electroweak corrections to decay processes at next-to-leading order: the dipole subtraction formalism and the one-cutoff slicing method. The former is extended to the case of decay kinematics which has not been considered in the literature so far. The latter is generalised to production and decay processes with more than two charged particles, where new “surface” terms arise. Arbitrary patterns of massive and massless external particles are considered, including the treatment of infrared singularities in dimensional or mass regularisation. As an application of the two techniques we present the calculation of the next-to-leading order QCD and electroweak corrections to the top-quark decay width including all off-shell and decay effects of intermediate W\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm {W}} $$\end{document} bosons. The result, e.g., represents a building block of a future calculation of NLO electroweak effects to off-shell top-quark pair (W+W-bb¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm {W^+}} {\mathrm {W^-}} {\mathrm {b}} {\bar{\mathrm{b}}} $$\end{document}) production. Moreover, this calculation can serve as the first step towards an event generator for top-quark decays at next-to-leading order accuracy, which can be used to attach top-quark decays to complicated many-particle top-quark processes, such as for tt¯+H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm {t}} {\bar{\mathrm{t}}} +{\mathrm {H}} $$\end{document} or tt¯+jets\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm {t}} {\bar{\mathrm{t}}} +\text {jets}$$\end{document}.


Introduction
The calculation of radiative corrections to cross-section predictions at next-to-leading order (NLO) saw tremendous progress in the past decades leading to a high degree of automation of calculations. This advancement was also driven by the development of general techniques for the treatment of infrared (IR) singularities. Such divergences arise in the intermediate steps of the calculation involving massless a e-mail: ahuss@phys.ethz.ch states and cancel in the final result by virtue of the Kinoshita-Lee-Nauenberg (KLN) and factorisation theorems. The various methods that were proposed in the literature can be broadly classified into two types: the phase-space-slicing [1][2][3][4][5] and the subtraction formalisms [6][7][8][9][10][11][12][13][14][15][16][17][18][19][20]. Both approaches allow for the analytic cancellation of all IR singularities without jeopardising the flexibility of a fully differential numerical phase-space integration.
In the first part of the paper we consider the so-called dipole subtraction formalism, which is one of the most established subtraction methods. It was continuously extended since its original formulation in Ref. [12], which was restricted to the computation of NLO QCD radiative corrections involving only massless partons, to cover all possible cases that are relevant in scattering processes. In Ref. [16] the dipole formalism was developed for the calculation of NLO QED corrections where the emission of photons from both massless and massive fermions was considered. The extension of the original algorithm for NLO QCD calculations to the case with massive partons was worked out in Refs. [17,18]. The treatment of non-collinear-safe observables in QED corrections was presented in Ref. [19] where the formalism was further extended to cover the various fermion-photon splittings not considered in Ref. [16].
For the case of decay processes, subtraction terms were constructed in Ref. [21] for the calculation of NLO QCD corrections to the radiative decay of the top-quark. These subtraction terms were based on the results of Ref. [22] and constructed to fit into the dipole subtraction formalism. Although the discussion in Ref. [21] was restricted to the case of the top-quark decay, the subtraction terms given there should in principle be applicable to the decay of any strongly interacting particle, however, with the restriction that only massless partons appear in the final state.
So far, no general algorithm exists for the calculation of NLO QED corrections to decay processes using the dipole subtraction formalism. In this paper we fill this gap by explicitly constructing such a process-independent subtraction term considering both massless and massive final-state fermions and work out its extension to the treatment of non-collinearsafe observables. The further extension of the formalism to treat massive partons in NLO QCD corrections to decay processes is then fully straightforward.
The one-cutoff phase-space slicing (OCS) method, considered in the second part of the paper, was initially introduced in Ref. [1] for the calculation of NLO QCD corrections to the multi-jet production cross section in e + e − annihilation. As such, its application was limited to the case with only massless coloured partons in the final state. It was later extended to allow for incoming partons [2] and massive quarks [4]. The OCS method is distinguished by the feature that one single Lorentz-invariant cut is introduced to isolate the soft and/or collinear regions in the real-emission phase space. Another frequently used slicing method, which is not considered in detail in this paper, is the so-called twocutoff slicing method [5,23]. As a consequence of introducing two independent cuts on energies and angles, however, this approach is not manifestly Lorentz invariant.
In its original formulation, the OCS method is restricted to amplitudes that involve only two sources (particles) of IR singularities, a feature that is naturally supported by colourordered QCD amplitudes, but neither by sub-leading colour structures in QCD nor by IR singularities originating from U (1) gauge bosons, such as the photon in QED. Furthermore, no results employing mass regularisation are documented in the literature, although it is a common choice in the calculation of electroweak radiative corrections. In this paper we extend the OCS formalism in these two respects, i.e. presenting a formulation employing mass regularisation and dealing with an arbitrary number of soft or collinear singularities. To this end, we employ results from the aforementioned dipole subtraction formalism to obtain an approximation of the realemission matrix elements in the singular regions together with the associated phase-space parametrisations suitable for the subsequent analytic integration. Furthermore, we discuss the subtleties that arise in the case of processes that involve arbitrary many sources of IR singularities.
At its design luminosity in Run 2, the LHC will be a topquark factory, copiously producing top-quarks in pairs, as single particles, or in association with other particles. As the heaviest elementary particle known to date, the top-quark is a unique window to explore the mechanism of electroweak symmetry breaking and to test the Standard Model (SM) on the one hand, and to probe physics beyond the SM on the other. Apart from the top-quark mass, key observables in the top sector comprise the top-quark width and all accessible production cross sections, for which precise theoretical pre-dictions are required. In these predictions the top-quark width t plays a double role: as an observable in its own right and as an ingredient in cross-section predictions that take into account top-quark decays. In the latter case, t enters the Breit-Wigner propagator of the top-quark resonance, so that it should be known to the same perturbative order as the calculated cross section.
The corrections to the top-quark decay process t → bW + are known at NLO [24][25][26][27] and next-to-next-to-leading order (NNLO) [28][29][30][31][32] in QCD and up to NLO for the electroweak (EW) corrections [33][34][35][36][37][38]. However, all prior calculations either consider on-shell W bosons or only partially account for off-shell effects by performing an expansion about the intermediate W resonance. As an application of the techniques introduced in this work, we present the fully differential calculation of the QCD and EW corrections at NLO to the three-body decay of the top-quark, accounting for off-shell effects of the intermediate W boson within the complex-mass scheme [39][40][41] and retaining a finite bottomquark mass. We consider both the hadronic and the semileptonic decay channels and further allow for non-collinearsafe photon emission ("bare leptons") for the latter. The presented calculation of the NLO EW corrections to the decays t → bW + → b + ν /bqq , thus, provides an ingredient to the (not yet existing) calculation of NLO EW corrections to the tt production process pp → bbW + W − → 6 f including top-quark decays, similar to the respective NLO QCD corrections to the top-quark decay in the NLO QCD calculation to pp → bbW + W − → 6 f [42][43][44][45][46].
On the experimental side, the first direct experimental determination of the top-quark decay width was performed by the experiments at the Tevatron collider [47,48]. It, however, still suffered from large experimental uncertainties. Recently, the first measurement at the LHC was presented by the CMS experiment with the data collected at √ s = 8 TeV [49], significantly improving on the previous measurements, with a precision at the 10 % level. This paper is organised as follows: In Sect. 2 we briefly review the general structure of an NLO calculation to decay processes to set up the notation and conventions used throughout this work. Section 3 first outlines the dipole subtraction formalism in Sect. 3.1 and subsequently discusses its extension to deal with decay kinematics. Massive final-state fermions are treated in Sect. 3.2, and the special case of light fermions is covered in Sect. 3.3, where we also discuss the modifications for non-collinear-safe observables. The OCS method is presented in Sect. 4 with a general overview of the formalism in Sect. 4.1, explaining in particular the origin and treatment of so-called "surface terms" which arise in the case of processes with more than two sources of IR singularities. The various cases that arise from the different combinations of massless or massive particles in the initial or final states are covered in Sect. 4.2. The application of the two tech-niques to the top-quark decay at NLO is presented in Sect. 5 while the calculational details are summarised in Sect. 5.1. In Sects. 5.2 and 5.3 we present the numerical results for the corrections to the decay width and the differential distributions, respectively. The conclusions and an outlook are given in Sect. 6.

NLO corrections to decay processes -notation and conventions
The leading-order (LO) partial decay width of a particle with mass m a decaying into a specific final-state configuration comprising n particles is given by where M 0 is the lowest-order transition amplitude for the 1 → n decay process with representing the possible sum (average) over the final (initial) state degrees of freedom and S {n} is the symmetry factor for identical particles in the final state. Here we have introduced the abbreviation for the n-particle phase-space measure, where p i and m i denote the momentum and mass of the ith final-state particle, respectively, and p a the momentum of the decaying particle, with p μ a = (m a , 0) in its centre-of-mass frame. Furthermore, we have used the shorthand notation n ≡ {p 1 , . . . , p n } for the associated kinematics given by the set of momenta of the particles.
At NLO, the decay width additionally receives contributions from radiative corrections, where δM and M 1 denote the virtual and real-emission matrix elements, respectively. Here we assume that the virtual corrections already include the counterterm contributions, so that all ultraviolet singularities are properly cancelled after renormalisation.
The calculation of higher-order corrections involving the emission of massless particles, in general, leads to the occurrence of IR singularities in both the virtual and the real corrections which cancel in the sum of Eq. (2.4) by virtue of the KLN theorem. 1 Although the cancellation of IR singularities is well understood, in practise, the procedure is non-trivial, since the virtual and real corrections are defined on different phase spaces. Moreover, the divergences contained in the real corrections are implicitly hidden in the phase-space integration, whereas the singularities in the virtual corrections are made explicit using regulators when performing the integration over the loop momentum. To accomplish an analytic cancellation of all IR singularities therefore requires one to make the singularities of the real corrections explicit by means of regulators. On the other hand, it is desirable to retain the flexibility of a fully differential numerical approach for the phase-space integration. Both the slicing and the subtraction techniques provide a prescription to perform a fully differential calculation suitable for a numerical evaluation while having the full analytic control over the IR singularities and their cancellation. Both approaches exploit the factorisation properties of the real-emission amplitude in the various singular limits which we briefly review in the following for the case of photon radiation off fermions in order to set up our notation. For further details we refer to Ref. [16] which we follow closely. Throughout this work we regularise the IR singularities using mass regularisation where an infinitesimal photon mass m γ is introduced and the masses of the light fermions m f are retained in the divergent logarithms. A prescription for the conversion of the results obtained using mass regularisation to the corresponding expressions using dimensional regularisation is given in the "Appendix".
The singular behaviour of emission amplitudes in the softphoton limit (k → 0) can be described by the well-known eikonal approximation where the squared real-emission amplitude, summed over the photon polarisations λ γ , can be written as follows: where e is the elementary charge. Here the summation over f, f extends over all charged external particles (of any spin), and the reduced matrix element M 0 is evaluated with the set of momenta that is obtained by omitting the momentum k of the soft photon. We have further introduced the sign factors σ f = ±1 describing the charge flow, which are defined as σ f = +1 for incoming particles and outgoing anti-particles and σ f = −1 for incoming anti-particles and outgoing particles. Overall charge conservation can then be expressed as For the case of decay processes considered here, the mass m a of the decaying particle represents the maximal scale in the process and therefore will never be considered in the small-mass limit. As a consequence, only finalstate collinear singularities arise where the radiated photon becomes collinear to a light fermion in the final state. The asymptotic behaviour of the squared transition matrix in this limit takes the form where κ i denotes the helicity of particle i which might change to any other spin value τ κ i . If i is a spin-1 2 fermion, the spin can only potentially change sign, i.e. there is just a sum over τ = ±. In this case, the functions g (out) i,τ are given by with the dimensionless variable representing the momentum fraction carried away by the fermion after the collinear splitting f → f + γ . For the case of unpolarised fermions, Eq. (2.7) reduces to i,− . Note that Eq. (2.8a) and (2.8b) does not describe particles with a spin other than 1 2 . However, any subtraction formalism based on this asymptotics may still be used for particles with another spin (such as a W boson) as long those particles are heavy enough to avoid collinear singularities, because the asymptotic behaviour (2.5) in the soft-photon limit is spin independent.

The dipole subtraction formalism
As was pointed out in the introduction, so far, no general formalism exists yet in the literature for the calculation of NLO QED corrections to decay processes within the dipole subtraction formalism. We present the construction of such a general subtraction term in the following and begin in Sect. 3.1 with a brief overview of the formalism and of our approach for its generalisation to decay processes. The explicit form of the corresponding radiator function and its integrated counterpart is given in Sect. 3.2 for the case of a massive final-state fermion. Section 3.3 considers the case of light fermions and the extension to non-collinear-safe observables.
The described formalism was already successfully applied in the evaluation of the NLO EW and mixed NNLO QCD/EW corrections to Drell-Yan-like W-and Z-boson production in the resonance region presented in Refs. [50,51], where the corrections to the gauge-boson decay processes appear as building blocks of the full calculation.

Overview of the method
The general idea of the subtraction formalism can be summarised in the following formula: where M 1 denotes the transition amplitude of the bremsstrahlung process. The auxiliary term |M sub | 2 is specifically tailored to act as a local counterterm to the realemission corrections. The subtraction term |M sub | 2 pointwise mimics the real-emission contributions in all its singular limits, i.e. it satisfies the following properties: where p i denotes the four-momentum associated with a light fermion in the final state and λ γ indicates the summation over the photon polarisation. The explicit asymptotics in the respective limits are given in Eqs. (2.5) and (2.7). The property (3.2) ensures that the first term on the r.h.s. of Eq. (3.1) is free of IR singularities. As a result, the integration over for the decay of a particle a and a charged final-state particle i. The decay dipole merges the two configurations associated with photon radiation from the final (left) and initial (right) state into one subtraction term the full real-emission phase space d n+1 can be performed without regulators which makes it suitable for the evaluation using numerical methods. The singularities originally contained in the real-emission corrections are now fully encapsulated in the second term of Eq. (3.1) which is constructed in such a way to allow for a direct analytical integration of the singular contributions. To this end, the real-emission phase space is factorised into the singular subspace given by the photonic part [dk] and the remaining phase space integration d n associated with the non-radiative process. The symbol "⊗" indicates that this phase-space factorisation may contain additional convolutions or summations. Performing the analytical integration over [dk] results in an expression where the singularities are explicitly regularised in some appropriate regularisation scheme. These divergences can then be analytically combined with the respective singularities appearing in the virtual corrections, so that possible cancellations can be made explicit. The remaining integration d n is free of singularities and can be evaluated numerically.
The dipole subtraction formalism exploits the factorisation of λ γ |M 1 | 2 in the singular limits by partitioning the subtraction term |M sub | 2 into so-called dipoles M sub, f f 2 according to where f and f are called the emitter and spectator particles, respectively. The dipole terms are written as a product of the auxiliary functions g (sub) f f and the squared matrix element |M 0 | 2 , where the former represent universal functions that encapsulate the singularity structure and can be integrated analytically over the singular photon phase space once and for all. The associated squared matrix element |M 0 | 2 is evaluated on the reduced n-particle phase space n, f f that remains after isolating the singular one-particle subspace [dk] of the photon. These dipole mappings are required to reproduce the appropriate kinematic configuration in the singular regions and further respect exact momentum conservation as well as all mass-shell conditions.
A plain extension in the spirit of the dipole subtraction formalism for the decay process would construct two new subtraction terms: a dipole where the decaying particle a is the emitter, g (sub) ai , and second where a acts as the spectator g (sub) ia . However, since there are no collinear singularities associated with the initial-state particle a, we combine the two subtraction terms to a single term called d which is illustrated in Fig. 1 in terms of the two schematic diagrams that would be associated with g (sub) ia and g (sub) ai . Such a combined treatment reduces the number of independent dipoles that need to be computed and, moreover, allows part of the calculation to be taken over from the scattering case with a final-state emitter and spectator described in Ref. [16]. The asymptotic behaviour of the subtraction function d (sub) ia differs from those of g (sub) ia and g (sub) ai of Ref. [16] in the soft limit and explicitly reads In the limit of small fermion masses m i → 0 the asymptotic form in the collinear region is required to be with the functions g (out) i,± given in Eq. (2.8a) and (2.8b). The subtraction term |M sub | 2 in the dipole subtraction formalism is then given by the following expression in the case of decay processes, cf. Eq. (3.3): where the summation over the indices i, j extends over all charged particles in the final state, a denotes the decaying particle, and g (sub) i j are the known subtraction functions of the dipoles with a final-state emitter-spectator pair given in Ref. [16].

Subtraction function for massive final-state fermions
In the following we present the construction of the subtraction function d (sub) ia for the decay of a particle a of mass m a with a massive final-state fermion i of mass m i . The special case of light fermions can be obtained by a systematic expansion in m i and will be discussed in the next section together with the extension to non-collinear-safe observables.
In the construction of the phase-space mapping n,ia for the decay dipole it is convenient to keep the four-momentum of the decaying particle a fixed. This constraint can be maintained by the deformation of the momenta k l of all particles that are not involved in the respective dipole term via a Lorentz transformation. Such a transformation is necessary in order to compensate the recoil induced by the radiated photon so that four-momentum conservation and all mass-shell conditions are retained. To this end, we define the auxiliary momentum P ia of the "recoiling system", see Fig. 1, which absorbs the recoil of the radiated photon through the transformation where and λ(x, y, z) denotes the Källén function defined as We further introduce the quantitȳ (3.11) where m γ denotes the infinitesimal photon mass that is introduced as a regulator for the soft singularities for later convenience. The phase-space mapping of the decay dipole then explicitly reads where k l denotes the momenta of all remaining particles of the decay process, including all neutral particles. The Lorentz transformation μ ν can be expressed in terms of the momenta P μ ia and P μ ia as follows: This phase-space parametrisation is identical to the one used in the scattering case with a final-state emitter and a final-state spectator j of Ref. [16], where particle i plays the role of the emitter and the recoiling system the role of the spectator of mass m 2 j = P 2 ia = P 2 ia . In detail, the following substitutions have to be made in order to transfer the quantities of Sect. 4.1 of Ref. [16] to the kinematics considered here: (3.14) Introducing the dipole variables we construct the subtraction functions d (sub) ia,τ as follows: , (3.16) where we have defined the auxiliary functions Note the similarity of these subtraction terms to the final-final dipoles g (sub) i j of Ref. [16] up to modifications that establish the appropriate behaviour in the soft limit (3.4). In detail, the functions g i j,τ defined in Eq. (4.4) of Ref. [16] coincide with d (sub) ia,τ up to the first term in the second line of Eq. (3.16), which is new, and by an additional factor of (1 + 2m 2 i /P 2 ia ) in the eikonal term (first term in the first line). As a consequence, many parts of the calculations can be taken over from Ref. [16], considerably simplifying the analytic integration over the photonic part of the phase space. It is easy to verify that the radiator function d The analytic integration of the subtraction functions over the photonic part of the phase space gives rise to the integrated counterparts The integrated decay dipole can then be written as with α = e 2 /(4π) and (3.21) and the abbreviations

Light fermions and non-collinear-safe observables
The subtraction functions d (sub) ia,τ given in Eq. (3.16) are already defined in such a way to correctly reproduce the asymptotic form in the collinear limit as demanded in Eq. (3.5) for the special of case of small fermion masses m i → 0. The corresponding expressions for the integrated dipoles are obtained by expanding Eq. (3.21) for small m i and omitting mass-suppressed terms, where we made use of the auxiliary function which was introduced in Eq. (3.8) of Ref. [16] to describe the universal singular behaviour of the integrated dipole functions G (sub) . The corresponding result where the IR singularities are regularised using dimensional regularisation can be obtained from Eq. (3.23) by performing a transition between the regularisation schemes using the results of Ref. [52] which is briefly summarised in the "Appendix".
In non-collinear-safe observables, photons arbitrarily close to an outgoing charged particle are not treated inclusively, but it is assumed that such a configuration can be experimentally resolved. This is, for instance, the case for muons in the final state, which are observed in the muon chambers further out in the detector, whereas the photons are detected independently in the electromagnetic calorimeter. Such an observable definition is sensitive to the details of the collinear f → f + γ splitting. The subtraction formalism can be extended to cover this case by exploiting the information on the dipole variable z ia which corresponds to the momentum fraction z i of Eq. (2.9) carried away by the emitter fermion in the collinear limit. This is accomplished by treating the n-particle kinematics of the dipole phase space n,ia as an (n + 1)-particle event with momenta p i → z iapi and k → (1 − z ia )p i . For further details we refer to Ref. [19] where this extension was introduced. Owing to the explicit dependence of the observables on the dipole variable, the integration over z ia ≡ z cannot be performed analytically, but is left open for a numerical evaluation. Employing a plus prescription, the singular endpoint part can be isolated so that the extension to the non-collinear-safe case amounts to an additional term which vanishes for inclusive observables, where the plus prescription is defined by for any smooth test function g(x). The functionD (sub) ia is obtained after carrying out the analytical integration over the dipole variable y ia and explicitly reads (3.27) where y + (z) denotes the upper boundary of the y ia integration and is given by (3.28) We have further introduced the dimensionless variable (3.29) and the auxiliary functions (3.30)

The one-cutoff phase-space slicing method
The main idea of the phase-space slicing approach is the isolation of resolved hard emissions from the unresolved soft and/or collinear regions in the real-emission phase space upon introducing a technical cut parameter. The hardemission contributions are free of singularities, so that the corresponding phase-space integration can be performed numerically using Monte Carlo methods. The contributions from the unresolved regions encapsulate the IR singularities of the real corrections, and the integrations over the subspaces containing the singularities are performed analytically, exploiting the universal factorisation properties of amplitudes in the soft or collinear limits. This makes the IR singularities that were implicitly contained in the realemission phase space explicit by means of regulators and allows for an analytic cancellation of all IR singularities. Various approaches have been proposed for the introduction of the technical cut. The OCS method considered in this work imposes a single cut 2 p i p j > s on the invariant scalar products of the momenta p i , p j of two particles that can cause a soft singularity ( p i → 0 or p j → 0) or a collinear singularity ( p i p j → 0). The method, thus, isolates soft and collinear singularities simultaneously and results in a formalism that is manifestly Lorentz invariant. The technical cut parameter s should be chosen small enough to suppress all effects of order O( s) to a numerically irrelevant level. Despite the drawback of introducing an arbitrary resolution parameter, such as s, slicing methods offer the possibility to widely suppress negative weights and to avoid unbounded weights in Monte Carlo integrations, which considerably simplifies a later extension of the Monte Carlo integrator to a Monte Carlo event generator with unweighted events.
Although the OCS method has been extensively studied in NLO QCD calculations involving massless [1,2] or massive [4,23,53] cases, there is so far no exhaustive survey of results covering all relevant configurations of massive particles in production and decay processes. The OCS method is particularly suited to NLO QCD calculations in leadingcolour approximation where the corrections are decomposed into colour-ordered contributions that naturally involve only singularities associated with the two neighbouring hard radiators of the unresolved parton in the colour ordering. If more than two collinear singularities appear in one contribution, as for instance in the sub-leading colour contributions, the application of OCS becomes subtle. In this case, processspecific treatments are required, such as reconstructing the IR singularities from the leading-colour building blocks as described in Refs. [1,2].
In the following, we describe a generic method to apply OCS in the presence of arbitrarily many collinear singularities in phase space. Moreover, we are not aware of any application of the OCS method to QED radiation processes, where e.g. no "colour ordering" exists leading to amplitude contributions with only two collinear singularities. Furthermore, the results in the literature all employ dimensional regularisation with either vanishing masses or large mass parameters; in QED applications, however, the case of small, but finite masses frequently occurs as well (e.g. for electrons or muons emitting photons). It should be noted that the case of a small mass m cannot be obtained upon simply taking the limit m → 0 in the OCS method, since the hierarchy s m 2 is used in the integration over the unresolved regions for a radiating particle of mass m. In order to isolate the mass singularity regulated by a small mass m (which defines a mass regularisation), however, the inverted hierarchy m 2 s is required. In the following we describe OCS for QED radiation covering all relevant mass configurations (with large and small masses). The transfer to QCD calculations is straightforward and merely requires the rein-troduction of colour charges known from the massless QCD case. 2

Overview of the method
The photonic real-emission corrections are divided into the finite "hard" contribution I H and the singular "soft-collinear" contribution I SC as follows: (4.1) The two regions are separated with the resolution parameter s imposing conditions s f γ ≶ s on the Lorentz-invariant quantities s f γ defined as where p f and k are the momenta of a charged fermion f and the photon, respectively. The definition (4.2) holds irrespective whether the fermion is massless or massive. In order to formally partition the phase space into hard and soft-collinear regions, we introduce shorthands for the θfunctions For the hard region we require s f γ > s for all charged fermions f of the process, so that the corresponding contribution to the real correction reads The phase-space integration in I H can be performed numerically, producing a result that depends on s via terms proportional to ln s or ln 2 s for small s. The integral over the soft-collinear region is given by In this integral, at least one of the invariants s f γ is small, i.e. s f γ < s. This means that in the limit s → 0, non- vanishing contributions from the soft-collinear region to the cross section originate entirely from those terms in the matrix element that lead to soft or collinear singularities, because regular contributions tend to zero with s → 0. Thus, in the calculation of soft-collinear contributions to the cross section we can employ the universal factorisation form of the squared amplitude, as for instance described in Sect. 3. In fact it is convenient to directly use |M sub | 2 = f = f |M sub, f f | 2 as defined in dipole subtraction in terms of emitter-spectator pairs. We, thus, can evaluate I SC as follows: For each emitter-spectator pair f f we can write as Using θ f +θ f = 1, we can rewrite the factor (1 − ) to and decompose I SC into soft (S), collinear (C), and quasi-soft (QS) contributions as follows: ) The definition of the various sub-contributions is schematically illustrated in Fig. 2.
The integrals I S, f f and I C, f f are confined to the soft and collinear regions, respectively, and involve only the technical cuts on the invariants s f γ and s f γ of the emitter and spectator fermions. For this reason, the integrations in I S, f f and I C, f f over the singular subspaces can be performed analytically for small s, as outlined below. On the other hand, the quasi-soft integrals I QS, f f involve cuts on all charged particles, rendering a process-independent analytical treatment impossible; nevertheless their numerical evaluation is straightforward, as outlined at the end of this section. Note that the quasi-soft integrals I QS, f f do not occur if only two sources of collinear singularities exist in the bremsstrahlung phase space. It is the occurrence of these subtle terms that requires process-specific and non-generic manipulations of the formalism as described in Refs. [1,2] for the sub-leading colour contributions.
The analytical calculation of the singular phase-space integrals in I S, f f and I C, f f starts by using the factorised matrix element of Eq. (3.3) and the corresponding factorised (n+1)particle phase space, the construction of which is described in Ref. [16] in detail. It is convenient to directly calculate the sum of the two contributions, In the following we restrict ourselves to the case of unpolarised fermions with the radiator functions Note that the spin information of a massive decaying particle is still fully accessible in spite of this choice, because the introduction of a separate spin-flip function g (sub) f f ,− was motivated by the massless limit for an emitter fermion, which is discontinuous. For collinear singularities of initial-state fermions, the multiplication symbol "⊗" in Eq. (4.13) indicates that a convolution over the collision energy after photon emission is involved in the phase-space factorisation. In this case, the analytical integration over d [k] can only be performed partially, as described in Sects. 4.2.2 and 4.2.3 below. For final-state singularities, however, the squared matrix ele- so that the singular integration is entirely contained in the universal factor For initial-state singularities, this factor comprises an essential ingredient as well.
Results on G f f for all relevant mass configurations are provided in the following sections. Finally, we describe a simple way to numerically evaluate the quasi-soft integrals I QS, f f . Recall that the integral receives only support outside the soft-collinear region of the emitter-spectator pair f f on a phase-space volume of O( s), because the photon has to be collinear to at least one other charged fermion. Non-vanishing contributions to I QS, f f in the limit s → 0 can, thus, only result from an enhancement in the integrand lifting the s suppression in phase space. Power-counting in the photon momentum in the integrand reveals that this can only happen if the photon momentum becomes soft, i.e. the photon momentum has to be of order O( s) in order to lie outside the soft-collinear region of f f . The integral I QS, f f is therefore effectively a surface integral of the soft-collinear region of f f and can be evaluated upon taking into account only those parts in the integrand that produce soft singularities, i.e. the "eikonal terms" in g (4.16) After symmetrizing with respect to the interchange f ↔ f , we can identify the usual eikonal factor Since only the symmetric combination of g (sub) f f ( p f , p f , k) is relevant in the overall quasi-soft contribution (summed over all emitter-spectator combinations), I QS, f f can be evaluated via Note that this form is valid for any combination of incoming and outgoing emitters and spectators. In the form given in Eq. (4.18), each term I QS, f f involves the phase-space embedding n+1 → n, f f that is specific to each emitterspectator pair. Since, however, the integral receives contributions only from soft photon momenta k, we can simplify the evaluation even further. For the relevant contribution we can i p i γ k j p j Fig. 3 Generic diagram with a final-state emitter i and a final-state spectator j factorise the (n +1)-particle phase space into the corresponding n-particle phase space for all particles but the photon and the 1-particle photon phase space, leading to the alternative form which does not make use of any special phase-space mapping for f f pairs and is, thus, much simpler to evaluate than Eq. (4.18). In this approximation, some care has to be taken with respect to the upper bound of the photon energy k 0 , which has to be large enough to cover the full relevant soft region, but still small enough to avoid artefacts from terms of the order O( slnk 0 max ). In the following section we carry out the integrations over the soft-collinear regions in I sing, f f for the various emitterspectator pairs f f in the initial/final states.

Final-state emitter and final-state spectator
Following the notation and conventions of Refs. [16,19], in the case of a final-state emitter i and a final-state spectator j (illustrated in Fig. 3) the full integral over the soft-collinear regions is given by where with λ(a, b, c) defined in Eq. (3.10). The variables parametrise the photon phase space. Their boundaries are given by where Note that the infinitesimal photon mass m γ acts as regulator for the soft singularity appearing at y i j → 0. The integrand g (4.25) and the soft and collinear regions are characterised as follows: Depending on the masses of the final-state particles (massive or light emitter/spectator), we obtain different results for G i j (P 2 i j , s), which are listed below.
The full integral over the soft-collinear regions, as defined in Eq. (4.20), is sufficient to cover all kinematic configurations in which the photon is recombined with collinear fermions in the event selection, i.e. in the case of collinearsafe observables. If the photon and the emitter i become collinear, this photon recombination ensures that all icollinear configurations integrated over in G i j (P 2 i j , s) are treated inclusively, i.e. either the full contribution represented by G i j (P 2 i j , s) contributes or it is completely excluded by cuts. This even holds true for contributions in G i j (P 2 i j , s) where the photon is collinear to the spectator j, since the photon has to be soft in this case, i.e. s jγ = O( s). In non-collinear-safe observables, however, the photon can be separated from collinear fermions, which is for instance the case for muons detected in the muon chambers. In this case the differential information on the energy flow inside the i-collinear cone should be kept in the phase-space integral (4.20). This energy flow in the soft-collinear region is controlled by the variable z i j , where it is effectively given by In this situation, the integration over z i j in Eq. (4.20) should be included in the numerical phase-space integration, so that instead of the full integral G i j (P 2 i j , s) the following function of z = z i j appears: where the functions y 1,2 (z) result from inverting the parametrisation of the integral boundary defined in Eq. (4.23). The evaluation ofḠ i j (P 2 i j , s, z) essentially follows the analogous treatment described for dipole subtraction in the case of non-collinear observables in Ref. [19]. The procedure is greatly simplified upon introducing a plus distribution for the z-integration, which isolates the soft singularity in the endpoint function G i j (P 2 i j , s), so that [Ḡ i j (P 2 i j , s, z)] + can be evaluated for m γ = 0. The explicit evaluation ofḠ i j (P 2 i j , s, z) and of the full integral G i j (P 2 i j , s) requires a discrimination of the different mass patterns in the emitter-spectator pairs.
(i) Massive emitter and massive spectator For massive emitter and spectator fermions, we have the following hierarchy of parameters: The full integral G i j (P 2 i j , s) is tedious, but straightforward. The result involves a discrimination between the variables m 2 i , m 2 j , P 2 i j and is given by with the auxiliary functions and the variable To calculate the functionḠ i j (P 2 i j , s, z) for noncollinear-safe observables, we have to invert the functions z 1,2 (y) of Eq. (4.23) using the hierarchy (4.30). Making use of the fact that only y i j values of O( y) with y 1 are relevant, the functions y 1,2 (z) are given by  O( y). The integral G i j (P 2 i j , s), thus, only receives contributions from a tiny neighbourhood of the point (y i j , z i j ) = (0, 1), i.e. it does not involve an integration over the momentum flow in the collinear regions. In this case non-collinear-safe observable can be calculated using G i j (P 2 i j , s) without the need to introduceḠ i j (P 2 i j , s, z).

(ii) Massive emitter i and light spectator j
Due to the small mass of the spectator, the hierarchy between the parameters is The spectator mass can be set to zero, since no collinear singularity exists for the spectator. In fact the integral G i j (P 2 i j , s) can be obtained upon taking the limit m j → 0 from the result of the previous section, with the function L from Eq. (3.24). Similarly to the case of non-vanishing m j , G i j (P 2 i j , s) does not comprise an integral over the energy flow in the collinear regions, so that the quantityḠ i j (P 2 i j , s, z) is not needed either. Since the two limits m 2 i s and s m 2 i do not commute, the integral G i j (P 2 i j , s) does not result from the simple limit m i → 0 applied to the previous cases. An explicit integration yields the result It contains a term ln 2 s that results from a combination of the soft and the i-collinear singularities, which are both isolated by the cut parameter s in the integration over the hard photon emission which is based on matrix elements with m i = 0. In order to calculateḠ i j (P 2 i j , s, z) for non-collinear-safe observables, we first derive the functions y 1,2 (z) for the integration boundary of the y i j -integration in Eq. (4.28). The relevant range in y i j is characterised by O(m 2 i /P 2 i j ) y i j O( y), leading to the region defined by The functionḠ i j (P 2 i j , s, z) for the non-collinear-safe contribution is easily calculated tō

(iv) Light emitter i and light spectator j
In this case G i j (P 2 i j , s) can be calculated upon taking the massless limit m j → 0 of the previous case. The final result is withḠ i j (P 2 i j , s, z) as given in Eq. (4.40).

Final-state emitter and initial-state spectator, and vice versa
In the cases of a final-state emitter i and an initial-state spectator a, and of an initial-state emitter a and a final-state spectator i (depicted on the left-and on the right-hand sides of Fig. 4, respectively), again the variables and abbreviations of Ref. [16] are used. The integrals over the full soft-collinear regions are defined as where the variables . (4.44) are confined to the boundary determined by the variables , with x 0 being set by the least energy to trigger the process without photon emission. 3 Here we use with λ(a, b, c) defined in Eq. (3.10) and Note that we always haveP 2 ia < 0, but P 2 ia is not necessarily negative for massive particles i and a. The soft and collinear regions are characterised as follows: (4.48c) 3 For P 2 ia > 0 and 0 < P 2 ia < m a − m i , the lower limit has to obey , see Ref. [16].
By construction, the limits x 0,1 respect the condition x 0 < 1/(1 + x) < x 1 . Results for the integrals G ia (P 2 ia , s) and G ai (P 2 ia , s) for the various mass configurations of emitterspectator pairs are listed in the following sections.
Similar to the previous case of both emitter and spectator in the final state, the full integrals (4.42) and (4.43) are not sufficient to deal with non-collinear-safe observables if the emitter particle f = i or f = a is light, where m 2 f s. To keep track of the momentum flow in the collinear region of the emitter, the z ia or the x ia integration has to be done numerically in the respective cases. To this end, we definē using x = x ia and z = z ia as shorthands. In the case of a finalstate emitter, zp i is the fermion momentum after collinear photon emission from the final state; in the case of an initialstate emitter, x p a is the fermion momentum after collinear photon emission from the initial state, i.e. the momentum x p a enters the hard non-radiative process. The functions x 1,2 (z) appearing in Eq. (4.50) result from interchanging the order in the integrations over x ia and z ia and are derived below for the relevant cases. The soft singularities in the integrals again are extracted from the numerical integration upon introducing plus distributions in the integrals, (4.53)

(i) Massive final-state emitter i and massive initial-state spectator a, and vice versa
The case of massive initial-state particles is mostly interesting for decay processes, since the masses of colliding particles are typically much smaller than the scattering energies in high-energy collisions. If both final-state emitter i and initial-state spectator a (or vice versa) are massive, we have to respect the following hierarchy of parameters: No collinear singularities are encountered, and the effective integration region in G ia,ai (P 2 ia , s) is a neighbourhood of the point (x ia , z ia ) = (1, 1), so that the functions G ia (P 2 ia , s, z) and G ai (P 2 ia , s, x) are not required. The full integrals read for m i ≶ m a , (4.56) with the auxiliary functions (4.58) Here we partially used the results of Ref. [16], where we have identified (1 − x 0 ) with x in the limit x → 0, and (4.59) (ii) Light particle i and massive particle a If particle i is light and particle a is heavy, we have to consider the hierarchy so that we encounter a collinear singularity in the function G ia (P 2 ia , s), but not in G ai (P 2 ia , s). Therefore, we can obtain G ai (P 2 ia , s) from Eq. (4.56) upon taking the limit m i → 0, but to obtain G ia (P 2 ia , s) for small m i we have to perform a new integration. The explicit results are (4.62) In order to calculate the functionḠ ia (P 2 ia , s, z) for the noncollinear-safe contribution we have to perform the integration in Eq. (4.50) over the x ia range, which is bounded by

(4.63)
A simple integration yields (4.64) (iii) Light particle a and massive particle i If particle a is light and particle i is heavy, we have to consider the hierarchy so that we encounter a collinear singularity in the function G ai (P 2 ia , s), but not in G ia (P 2 ia , s). We can obtain G ia (P 2 ia , s) from Eq. (4.55) upon taking the limit m a → 0, but to obtain G ai (P 2 ia , s) for small m a we have to perform another integration. The explicit results are (4.67) In order to calculate the function G ai (P 2 ia , s, x) for the noncollinear-safe contribution we have to perform the integration in Eq. (4.51) over the z ia range, which is bounded by A simply integration yields (4.69)

(iv) Light particles a and i
Finally, we consider the case of both a and i being light, so that the mass hierarchy is Here we always have P 2 ia < 0. The full integrals G ia (P 2 ia , s) and G ai (P 2 ia , s) can be obtained from the results above where the respective emitter is light upon taking the limit of a massless spectator, i.e. upon taking m a → 0 in Eq. (4.61) and m i → 0 in Eq. (4.67). The results are The corresponding functionsḠ ia (P 2 ia , s, z) and G ai (P 2 ia , s, x) for non-collinear-safe observables are the same as given in Eqs. (4.64) and (4.69), respectively.

Initial-state emitter and spectator
Although not needed for the treatment of decays, for completeness here we collect the formulae with both emitter a and spectator b in the initial state with masses m a and m b . The corresponding structural diagram is shown in Fig. 5. Using again the kinematical variables introduced in Ref. [16], (4.74) the integrals over the full soft-collinear region is defined as with the integration boundary defined by where again x 0 is the lowest kinematically allowed value. The soft and collinear regions are characterised as follows: to keep track of the momentum flow in the collinear region of the emitter. Again we used x = x ab to shorten the notation.
The momentum x p a is the fermion momentum after collinear photon emission from the initial state, i.e. the momentum x p a enters the hard non-radiative process. The soft singularities in the integrals again are extracted from the numerical integration upon introducing plus distributions in the integrals,

(i) Massive initial-state emitter a and spectator b
If both initial-state emitter a and initial-state spectator b are massive, we have to respect the following hierarchy of parameters: No collinear singularities are encountered, and the effective integration region in G ab (s, s) is a neighbourhood of the point (x ab , y ab ) = (1, 0), so that the function G ab (s, s, x) is not required. The full integral reads G ab (s, s) =s √ λ ab ln(d 1 ) ln m 2 γ λ ab s 2 s with the auxiliary functions f ab,± (s) = ±s √ λ ab ln 2 s + 2m 2 a ± √ λ ab 2s Here we partially used the results of Ref. [16], where we have identified (1 − x 0 ) with 2 y in the limit y → 0, and There is no need to introduce the function G ab (s, s, x) if the emitter particle is massive.

(ii) Massive emitter a and light spectator b
The hierarchy between the parameters is Similarly to the case of non-vanishing m a and m b , G ab (s, s) does not comprise an integral over the energy flow in the collinear regions, so that the quantity G ab (s, s, z) is not needed either.

(iii) Light emitter a and massive spectator b
If the emitter a is light and the spectator b is heavy, we have to consider the hierarchy We encounter a collinear singularity in the function G ab (s, s), i.e. we have to perform another integration to obtain G ab (s, s) for small m a . The explicit result is (4.88) In order to calculate the function G ab (s, s, x) we have to perform the integration in Eq. (4.79) over the y ab range, which is bounded by

iv) Light emitter a and spectator b
Finally, we consider the case of both a and b being light, so that the mass hierarchy is The corresponding function G ab (s, s, x) is the same as given in (4.90).

Application to the top-quark decay
As an application of the techniques described in Sects. 3 and 4, in the following we present the calculation of the QCD and EW radiative corrections to the top-quark decay at NLO. First, we summarise our calculational setup in Sect. 5.1. Then, we discuss the different contributions to the top-quark decay width in Sect. 5.2. Finally, in Sect. 5.3 we present differential distributions in the top-quark rest frame for a set of kinematic observables.

Calculational setup and input parameters
Following Ref. [54], the set of input parameters used in our numerical evaluation is given by In order to describe the W-boson resonance, we employ the complex-mass scheme [39][40][41], a method that takes into account the effects of the instability and the off-shellness of the unstable particle in a gauge-invariant way, providing a consistent procedure at the one-loop level. The main idea of the complex-mass scheme is to consider the squared boson masses as complex quantities, defined through the gaugeinvariant pole mass M V and width V , and to introduce them directly at the level of the Lagrangian density. To apply the complex-mass scheme, the on-shell masses of the W and Z bosons given in Eq. (5.1) are converted to the corresponding pole masses according to [56][57][58] To preserve gauge invariance, the complex masses have to be introduced in a consistent manner. In particular, this implies a complex definition of the weak mixing angle, and a generalisation of the renormalisation procedure for unstable particles. Further details of the complex-mass scheme can be found in Refs. [40,41]. The dependence on the bottom-quark mass is fully accounted for without approximation. The value of m b given in Eq. (5.1) corresponds to the pole mass through the twoloop conversion formula given in Ref. [54]. Although all fermions with the exception of the top and bottom quarks are treated as massless in our calculation, a dependence on the mass of the light fermion arises when non-collinear-safe observables are considered (see Sect. 3.3). In the following, they are only considered for final-state muons where the experimental discrimination of a collinear muon-photon configuration is possible. Note that no event selection cuts are applied, so that the issue of collinear safety only arises in the event reconstruction. Thus, only differential distributions are concerned with this issue, while the integrated results, i.e. the (partial) decay widths, are collinear safe.
The EW coupling constant α is derived from the muon decay constant G μ via the relation This definition corresponds to the so-called G μ -scheme, which avoids large logarithms induced by the running of α from the Thomson limit (Q 2 = 0) to the EW scale (Q 2 ∼ M 2 W ) and further absorbs dominant corrections to the ρ parameter into the coupling definition. Using α G μ instead of α in LO predictions changes the charge-renormalisation constant by a term containing the quantity r [34,59] which quantifies the radiative corrections to muon decay (see, e.g., Ref. [60] for further details). For the strong coupling α s we employ the fixed value given in Eq. (5.1).

Details on the calculations
The implementation employing the slicing method of Sect. 4 was performed using matrix elements that were calculated using the Weyl-van-der-Waerden spinor formalism (as formulated in Ref. [61]) for the tree-level amplitudes and inhouse Mathematica routines for the virtual one-loop corrections. Optimised phase-space parametrisations, based on Refs. [39,62], were chosen for each of the pronounced structures in the integrand induced by collinear poles or Breit-Wigner resonances. The numerical integration was performed using the Vegas algorithm [63,64].
A second independent calculation was performed using the subtraction method described in Sect. 3. Here, the matrix elements were obtained with the combined packages FeynArts [65] and FormCalc [66,67]. An independent implementation of routines for the phase-space genera-tion was employed. As in the calculation based on the slicing approach, the numerical integration was performed with Vegas.
Both implementations resort to the Collier library [68], which is mainly based on the results of Refs. [69,70] for the numerical evaluation of loop integrals and supports both dimensional and mass regularisation for the treatment of IR singularities.

NLO corrections to the top-quark decay width
In Table 1 we collect the numerical results on the various contributions to the top-quark decay width, providing separate results on the semi-leptonic and hadronic decay channels.
The relative correction factors δ are defined with respect to the corresponding LO prediction LO t , where NLO QCD/EW t denotes the NLO prediction taking into account only QCD/EW corrections, while the full NLO QCD+EW result is denoted NLO t . The QCD corrections to the semi-leptonic decay of the top-quark are known in semi-analytic form [24] and are in full agreement with our numerical result. 4 The full NLO prediction 5 NLO t = 1.3693(2)GeV for the total width (see Table 1) can be compared to the most recent CMS measurement [49] of t = 1.36 ± 0.02 (stat) +0.14 −0.11 (syst) GeV, (5.7) assuming m t = 172.5 GeV, revealing good agreement between theory and experiment at the current accuracy level. A comparison between the results obtained with OCS and dipole subtraction is shown in Fig. 6.
As discussed in Sect. 4, the slicing approach introduces a cut parameter s in order to isolate the singular region of the real-emission phase space. The dependence on this cut parameter cancels between the individual parts of the calculation in a non-trivial manner when combined to the final result, constituting a powerful check on the correctness of the calculation. Although this implies a certain degree of arbitrariness in the choice of the specific value for the cut parameter, in practice, sensible values for s are constrained to a certain range. This interval of validity is limited by the approximation used in the analytic integration 4 A readjustment of the prefactor of the result given in Ref. [24] according to G μ → G μ · Re(s 2 w )/s 2 w is required in order to be consistent with the complex-mass scheme employed in our calculation. 5 The total width is calculated using lepton universality, i.e. the τ mass m τ is neglected, since those effects are of O(m 2 τ /M 2 W ).  on the upper, and by the stability of the numerical evaluation on the lower end. It is therefore necessary to perform a scan over the cut parameter to identify the region where both restrictions are simultaneously fulfilled and where the result can be trusted. Figure 6 illustrates such a variation of s covering three orders of magnitude, clearly depicting the on-set of the aforementioned breakdown of the calculation. Furthermore, we observe good agreement between the results obtained by the two techniques which overlap within their Monte Carlo errors where the s variation develops a plateau. The OCS results presented in the following are obtained for a specific value of s where a point in the middle of the plateau was chosen after performing such a scan. This procedure is repeated for each decay mode (semi-leptonic and hadronic) and each type of correction (QCD and EW). At this point, a comment on the relative performance of the two methods is in order: Although the OCSM has the advantage of avoiding negative weights to a large extent, the additional variation of the slicing cut that is required to locate the plateau constitutes an additional computational overhead. Furthermore, we require a higher degree of optimisation in the integrator to sample the regions of the real-emission phase space close to the technical cut to high accuracy, which is then cancelled against the large contribution that arises from the analytically integrated soft-collinear contributions. These are the well-known drawbacks of slicing approaches compared to local subtraction methods such as the dipole formalism which generally display better numerical convergence and stability. More specifically, in order to reach a comparable Monte Carlo error between the two implementations, the dipole subtraction formalism roughly requires an order of magnitude less statistics in the integration compared to the OCSM method.
Finally in the corrections to the partial decay widths, in order not to spoil the proper comparison to the results for the full topquark decays. Table 1 shows that the off-shell result and the NWA for the top decay widths differ by about 1.6 %, in good agreement with the naive expectation of O( W /M W ), which quantifies typical off-shell effects of the W boson in integrated quantities. The relative corrections are influenced by the off-shell effects only at the sub per-mille level. Off-shell effects of the W boson will, however, play a more important role in differential distributions.

Differential distributions in the top-quark rest frame
In the following, we present differential distributions for the decay of the top-quark in its rest frame. After describing the event reconstruction in Sect. 5.3.1, we present the numerical results for the semi-leptonic and hadronic decay channels in Sects. 5.3.2 and 5.3.3, respectively. The differential distributions and their corrections are shown as plots that are subdivided into three frames with the following conventions: The upper panel depicts the absolute distributions, while the bottom two frames show the relative EW (middle panel) and QCD (bottom panel) corrections.

Event reconstruction
It is not always possible to fully distinguish the particles of the process in the soft and collinear configurations. In our calculation, we try to be as close as possible to an experimental situation and distinguish the following scenarios: (i) Electroweak corrections to the semi-leptonic decay "Dressed leptons": The EW corrections contain IR singularities originating from soft photons and configurations where the photon becomes collinear to the final-state charged lepton. For sufficiently inclusive observables these singularities cancel in the final result. To this end, we apply the photon-recombination procedure described in Ref. [71] which attributes the momentum of a collinear photon to the charged lepton when close to it, i.e., when R( , γ ) = y( , γ ) 2 + φ( , γ ) 2 < 0.1, where y ( φ) is the rapidity (azimuthal angle) difference between two objects. Such a recombination procedure is mandatory for a realistic experimental description in the case of electrons. "Bare muons": In the case of muons, a collinear muon-photon pair can be experimentally resolved, so that no photon recombination is required. This results in non-collinear-safe observables, which are treated following the procedures described in Sects. 3.3 and 4.2.

(ii) QCD corrections to the semi-leptonic and hadronic decays
In contrast to the semi-leptonic case, where we assume that the bottom quark could be directly identified with a bottom jet, in the hadronic case we are confronted with the situation that experimentally indistinguishable partons are present in the final state. It is therefore necessary to define an IR-safe observable which groups soft and/or collinear partons into jets. In our calculation we adopt the inclusive generalised anti-k T algorithm [72,73] for two values of the clustering parameter, R = 0.6 and R = 1.0. The first value corresponds to the ATLAS default choice [74], the second one will be motivated in Sect. 5.3.3. We assume a bottom-tagging efficiency of 100 % and identify the jet containing the b quark as the b-tagged jet. In the following, we refer to the jet candidate with the highest transverse momentum 6 as the leading jet.
(iii) Electroweak corrections to the hadronic decay When computing the EW corrections to the hadronic decay, we do not distinguish the photon from the other QCD partons and apply the same jet algorithm as described in the previous paragraph.

Semi-leptonic decay of the top-quark
The numerical results for the NLO corrections to the semileptonic decay of the top-quark, t → b ν + , are summarised in Fig. 7.
They comprise differential distributions in the energy of the charged lepton (E ) and of the bottom quark (E b ), the invariant mass of the bottom quark and the charged lepton (M b ), and the cosine of the angle between the directions of the lepton and the bottom quark (cos θ b ).
Following the reconstruction described in Sect. 5.3.1, we distinguish the two cases of "dressed leptons" and "bare muons" for the EW corrections. For observables that are not sensitive to the collinear → γ splitting, i.e. which are collinear safe, both corrections must be nearly identical. This can be seen in the distributions in cos θ b and E b shown in the upper right and lower right plots in Fig. 7, respectively. In the former case in particular, the photon is indeed recombined only if collinear to the positron, i.e. it does hardly change the direction of the charged lepton. The EW corrections have a minimum in the proximity of the peak (cos θ b ∼ −0.8), are positive and smaller than 2 %. For comparison, the QCD corrections reach down to −12 %. The distribution in the b-jet energy is characterised by the fact that a 1 → 2 particle decay is effectively taking place when the W boson is on shell. In this case, the energy E b of the b-jet tends to the fixed value This explains the peak at E b ≈ 68 GeV. The asymmetric shape of the EW corrections is due to the typical radiative tail produced by photon radiation off a particle near a resonance energy. Most b quarks have an energy E b ≈ E OS b after the top-quark decay, so that photon radiation off b quarks, which is collinearly enhanced, reduces E b , leading to large positive corrections at energies below E OS b . This effect is further enhanced in the relative corrections, because the LO distribution is strongly suppressed there. The same effect occurs in the QCD correction as well, leading to a similar shape, but even more than an order of magnitude larger. Note that in this case corrections more than 100 % do not signal a breakdown of perturbation theory, they merely signal that the effective lowest order is provided by the radiative process with hard gluon (and partially photon) emission. A proper description of this tail, thus, would require the inclusion of multi-jet (and multi-photon) emission, as provided, e.g., by parton showers.
For the remaining distributions, i.e. in M b and E , we observe EW corrections that are larger in magnitude for "bare muons", stemming from the large mass-singular logarithm α ln m μ and reflecting the non-collinear-safe nature of the corresponding observable. The emitted photon in the → γ splitting removes some energy from the charged lepton. Contrary to the "bare muon" case, the "dressed lepton" effectively can reabsorb some photons in its redefinition when the photon-recombination procedure is applied. Hence, the EW correction to E , always monotonically decreasing, are steeper for the "bare muon" than for the "dressed lepton", as can be seen in the bottom left plot of Fig. 7 Fig. 7. Here, the EW corrections decrease from +6 to −4 % for "dressed leptons", while they drop from +8 to −10 % when "bare muons" are considered in the top-quark decay. The QCD corrections to this distribution are decreasing down to −45 % for large invariant masses.

Hadronic decay of the top-quark
Here, we present the results for the QCD and EW corrections to the top-quark hadronic decay width. In our setup, the corrections to the two processes t → b u N d and t → b cN s coincide and are given in Table 1 for the integrated decay width.
In the top-left plot of Fig. 8 we show the distribution in the invariant mass of the b-jet and the leading jet.
In the upper frame we notice that the QCD corrections evaluated for R = 0.6 turn the distribution negative at M b j ∼ m t . This is of course an unphysical effect: For a too low value of R, as in this case, the recombination between the real and the virtual parts is not sufficiently inclusive, so that the IR divergence leaves a trace as α s ln R corrections. This problem can be solved by enlarging the cone, e.g. to R = 1, as shown in the same plot. With a wider cone, there is indeed a larger probability for the particles to be merged. In the case where all the final-state particles (excluding the b quark) are combined into a single jet, the invariant mass M b j reproduces the top-quark mass value. This explains the appearance of a peak at M b j = m t . This obviously happens more frequently for a 1 → 3 kinematics (as the LO) rather than at NLO. At NLO, the extra gluon can populate M b j regions less probable at LO. This in turns explains the huge corrections in the lower frame of Fig. 8 (top left): The magnitude of the NLO correction, shown in the middle frame, gets very large due to the strong suppression of the LO in the normalisation. The EW correction is much smaller than the QCD correction, but is similar in shape. The curve for R = 0.6 does not suffer from the same issue as the QCD counterpart. The R = 1 case is nonetheless shown for comparison.
The distributions obtained choosing R = 0.6 or R = 1.0 are compared also for the other observables. The main difference is that the relative corrections induced by the choice of the larger radius are milder than those obtained when using R = 0.6, as expected from the inclusiveness argument. This effect can be clearly seen in the bottom-left and bottom-right plots in Fig. 8, where the NLO corrections to the energy of the leading jet and of the b-jet are presented, respectively. The energy of the leading jet receives positive QCD corrections due to gluon radiation, which turn negative above the energy peak at ∼65 GeV. As a consequence, the shape of the energy distribution at NLO is shifted to the left, more evidently for R = 0.6. In the low-energy tails the relative corrections are enhanced owing the normalisation to the small LO prediction in this region. The QCD and EW corrections to the distribution in the energy E b of the b-jet, displayed in Fig. 8 (bottom right), again show the distinctive radiative tails as already discussed for the semi-leptonic top decay in the previous section.
The different choice of the jet-algorithm parameter affects also the cosine of the angle between the b-jet and the leading jet, as displayed in the top-right plot in Fig. 8. Once again, a larger jet radius implies a larger probability to merge all final-state partons other than the b quark into a single jet. The latter will then be back-to-back to the b-jet since in this calculation the top-quark is considered at rest. This explains the peak appearing in the bin close to cos θ b j = −1 when R = 1 is employed. As in the previous case, the shape of the EW corrections match closely the QCD ones, however, are smaller in size by an order of magnitude. The strong rise of the QCD corrections toward larger values of cos θ bj again originates from the fact that gluon radiation populates the jet phase space more uniformly than at LO, where there is a strong tendency that the b quark and the hardest jet are back-to-back.

Conclusions
In this paper we have described extensions of two techniques for the treatment of infrared singularities in the computation of electroweak corrections at next-to-leading order: In Sect. 3 we have supplemented the construction of universal subtraction terms within the dipole subtraction formalism to cover decay kinematics. The one-cutoff slicing method described in Sect. 4 was derived on the basis of the dipole subtraction formalism and organises the singular contributions in terms of pairs of fermions. These results represent the first ones employing mass regularisation and further allow for an arbitrary number of charged particles, a case that has not been discussed in detail in the literature so far. When dealing with more than three charged particles in the process, the appearance of new "surface" terms has been emphasised. For both methods we have covered all possible cases of massless and massive fermions in the final state and further worked out the respective extensions for the treatment of non-collinear-safe observables. The extension of the described techniques to the calculation of NLO QCD corrections is straightforward.
Both methods have been applied to compute the NLO QCD and electroweak corrections to the top-quark decay width. In this computation, we have considered the b-quark as massive and evaluated the correction to the top-quark three-body decays fully accounting for the off-shellness of the intermediate W boson within the complex-mass scheme. Both the hadronic and the semi-leptonic decay channels are considered, further accounting for the case of non-collinearsafe observables in the latter. For the semi-leptonic decay, the NLO QCD corrections to the partial decay width amount to ∼−9.4 and ∼+1.3 % for the electroweak corrections. In the hadronic decay channel, the NLO QCD correction are more moderate with ∼−5.6 %, however, still much larger than the EW correction which amount to ∼+1.3 %. While off-shell effects of the W boson modify the integrated decay widths by about 1.5 %, the relative NLO QCD and EW corrections are modified only by few per mille by W off-shell effects. We have further presented differential distributions in the top-quark rest frame for different kinematic observables.
The final result for the total decay width of the top-quark at NLO QCD+EW is given by NLO t = 1.369 GeV, which corresponds to an overall correction with respect to the LO prediction of δ NLO = −5.5 %. This breaks down to an overall NLO QCD correction of ∼−6.8 %, and to a NLO EW correction of ∼+1.3 %, normalised to LO. As mentioned in the introduction, the total top decay width is known at NNLO QCD accuracy. The O α 2 s correction, evaluated with a finite b-quark mass and a finite W-boson width in a soft-collinear effective theory and normalised to the LO, amounts to a further ∼−2.1 % [29,32]. A comparison of our NLO QCD+EW prediction -with and without the additional NNLO correction of Refs. [29,32] -reveals good agreement with the most recent measurement by the CMS experiment However, the quoted experimental error is still one order of magnitude larger than the typical EW corrections of ∼+1.3 % addressed here. The experimental error is dominated by systematic uncertainties, such as the jet energy scale, where a significant improvement is difficult. The study of top-quark properties is an important part of the physics programme of future lepton colliders which aims for a precision at the order of 10−30 MeV in the experimental determination of the top width [75,76] and where the EW corrections are important. The codes developed in the context of this work can serve as a first step towards an unweighted Monte Carlo generator for the decay of the top-quark at NLO, which can easily be linked to any top-quark production process in order to simulate the decay at NLO accuracy. This is particularly important for processes with many particles in the final state when all intermediate particles are decayed. For instance, tt + V (V ) (V = γ, Z, W ± ), ttH, ttbb, and tttt, among others. Moreover, the top-quark decay width including both NLO QCD and EW corrections and W-boson off-shell effects comprises an ingredient in NLO predictions for the production of single top-quarks or tt pairs including top-quark decays and offshell effects, such as the process pp → bbW + W − → 6 f , for which NLO QCD corrections are known, but not yet the EW corrections.