The fully differential top decay distribution

We write down the four-dimensional fully differential decay distribution for the top quark decay t→Wb→ℓνb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t \rightarrow Wb \rightarrow \ell \nu b$$\end{document}. We discuss how its eight physical parameters can be measured, either with a global fit or with the use of selected one-dimensional distributions and asymmetries. We give expressions for the top decay amplitudes for a general tbW interaction, and show how the untangled measurement of the two components of the fraction of longitudinal W bosons – those with b quark helicities of 1 / 2 and -1/2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-1/2$$\end{document}, respectively – could improve the precision of a global fit to the tbW vertex.


Introduction
The detailed study of the properties of the top quark has a widespread interest as a probe of new physics beyond the Standard Model (SM) [1][2][3]. In hadron collisions top quarks are produced in pairs, mediated by the QCD interaction, and they are also singly produced via electroweak interactionssetting aside other sub-dominant production mechanisms in association with additional bosons γ , W , Z or H . In contrast, the decay of the top (anti-)quarks is almost completely dominated by the mode t → W b, with the subsequent decays of the W bosons into charged leptons, W → ν, = e, μ, τ , or into quarks, W → qq .
The top decay differential distribution in t → W b → ν is determined by four angles. The first two angles (θ, φ) are the spherical coordinates of the W boson momentum in the top quark rest frame, in some arbitrary reference system (x, y, z). The remaining two angles (θ * , φ * ) are the spherical coordinates of the charged lepton momentum in the W boson rest frame, in the reference system (x , y , z ) obtained by a "standard" boost from the previous one, such that thê z axis is in the direction (θ, φ) of the W boson momentum and theŷ axis is in the x y plane (see for example a e-mail: jaas@ugr.es Ref. [4] for a detailed discussion). One-dimensional distributions or asymmetries were already proposed long ago to measure some physical quantities involved in the top decay, such as the W helicity fractions, measured from the θ * distribution, or the top polarisation, measured from the θ and φ distributions [5]. Since then, W helicity fractions have been experimentally measured in a number of experiments, with the most precise measurements obtained by the ATLAS [6] and CMS [7] Collaborations, and the top polarisation in certain directions has also been measured in single top [8,9] and top pair [10,11] production. Other angular distributions, like the polar angle of the charged lepton in the top quark rest frame [12], the azimuthal angle [13], and the azimuthal angle in the laboratory frame [14], are also sensitive to top polarisation effects.
More generally, it has been shown that the full set of eight W boson spin observables can be measured from selected distributions in top quark decays [15], and preliminary measurements have been performed by the ATLAS Collaboration [9]. Moreover, the large top samples available at the LHC Run 1 have allowed one to perform more demanding measurements, such as the determination of the two-dimensional (θ * , φ * ) distribution by the ATLAS Collaboration in t-channel single top production [16]. With the increased statistics at Run 2, it is likely that measurements of the full four-angle distribution, until recently unconceivable, will be achieved. The aim of this work is to provide the framework for such measurements and point out their advantages for a global fit of the tbW vertex. We provide analytical expressions for the fully differential distribution, show how to extract the relevant information either directly or using suitable asymmetries, and relate the top decay amplitudes with W spin observables. Finally, we write the physical quantities involved in the top decay for a general tbW effective Lagrangian, showing that the untangled measurement of the two components of the fraction of longitudinal W bosons may improve the precision of the global fit.

The fully differential distribution
Using the helicity formalism of Jacob and Wick [17], the amplitude for the top decay t → W b → νb can be written, in the narrow width approximation, as with λ 1 , λ 2 , λ 3 and λ 4 the helicities of the W boson, b quark, charged lepton and neutrino, respectively, M the third spin component of the top quark, and = λ 1 − λ 2 , λ = λ 3 − λ 4 . The angular dependence is given by the well-known Wigner D functions [18] and a λ 1 λ 2 , b λ 3 λ 4 are constants. For the top decay, angular momentum conservation implies that there are only four nonzero reduced amplitudes (or just "amplitudes" for short) a . 1 The SM interaction of the W boson with charged leptons implies λ 3 = ±1/2, λ 4 = ∓1/2 for W ± decays, assuming massless charged leptons. Therefore, λ = 1 for top quarks, λ = −1 for top anti-quarks, and there is only one non-zero b constant in each case, which can be factored out and does not play any further role in the distributions. Let us introduce a top spin density matrix in terms of the top (anti-)quark polarisation in the three axes, P i = 2 S i , with S i the spin operators. The normalised differential decay width can be written, summing over b quark helicities, as determining the top decay distribution. One can define four untangled helicity fractions, where F + + and F − − equal the usual helicity fractions F + and F − , that is, the relative fractions of W bosons with helicity λ 1 = 1 or λ 1 = −1, respectively. We will thus drop the superscript for these in the rest of the paper. The fraction of W bosons with helicity The fractions F − + and F + − , which would correspond to the amplitudes a 1 − 1 2 and a −1 1 2 , vanish at leading order (LO) due to angular momentum conservation, and they are still extremely small at next-to-leading order [19]. The untangled W helicity fractions yield three independent parameters because their sum equals unity. The remaining two parameters can be taken as the phases of the only interference terms appearing in the sum (4), The property of the D functions and the compostion rule in terms of Clebsch-Gordan coefficients D guarantee that the four-dimensional distribution can be expanded as a finite combination of the set of functions we define as As it is shown in Appendix A, these functions are orthonormal. By writing the distribution as we find that the non-zero coefficients in the expansion are c 00 00 = 1 4π , Because of the orthonormality of the M functions, the different coefficients can be determined by projecting the dif- The Monte Carlo estimate of the integral is the average of the function (M j 1 j 2 m m ) * over a selected set of points, given by a probability density function g. Therefore, with real data, the estimate of the coefficients c j 1 j 2 m m is simply done by computing the average of (M j 1 j 2 m m ) * over the dataset [20], properly corrected for detector effects [21]. The statistical error on the coefficients, and the covariance matrix between different coefficients, can easily be estimated simply by taking averages of the product of M functions.

Physics parameters from asymmetries
The set of coefficients c j 1 j 2 m m allows one to extract the three top polarisation components and the five physical parameters in the decay of the top quark at the same time. Alternatively, one can measure them using asymmetries or one-dimensional distributions. Using the explicit expressions for the D functions, one can compactly write the fully differential distribution in terms of the four angles as Here, u L = (sin θ cos φ, sin θ sin φ, cos θ) is the unit vector in the direction of the W boson momentum in the top quark rest frame, and u T = (cos θ cos φ, cos θ sin φ, − sin θ), u N = (sin φ, − cos φ, 0) are two orthonormal vectors. As before, λ = 1 for top quarks and λ = −1 for anti-quarks. The explicit expression above suggests the observables one has to use to extract the desired quantities. First, as mentioned above, the (entangled) W helicity fractions are and they can be determined from the θ * distribution [5]. In order to measure the top polarisation, one can use the double forward-backward (FB) asymmetry defined in [22] and analogous ones for thex andŷ axes, since the factor between brackets is merely the sum of F + and F − . The two squared moduli of λ 1 = 0 amplitudes, whose sum appears in F 0 , can be disentangled by a forwardbackward-edge-central asymmetry, defined as Alternatively, we point out that the W boson spin analysing power α W or, equivalently, the FB asymmetry in cos θ , is also sensitive to the difference |a 0 1 2 | 2 − |a 0 − 1 2 | 2 . We point out that measuring the untangled helicity fractions F + 0 , F − 0 requires polarised top quarks, as is also seen from the expressions of c j 1 j 2 m m . The relative phases of the interfering amplitudes can be extracted from FB and double FB asymmetries, As is apparent from the above equations, the introduction of a FB asymmetry in cos θ * allows one to flip the sign of the latter terms between brackets, allowing to measure the real and imaginary parts of each product independently. Again, polarised top quarks are required to measure the quantities, as otherwise there is no privileged direction in the W boson rest frame other than theẑ axis. For completeness, we also give the relation between the eight W spin observables [15] and top decay amplitudes. They are with B 1 = B 2 = 0 due to angular momentum conservation.
Finally, let us stress that the global approach in Eqs. (10)-(12) and the use of selected observables (14)-(18) are formally equivalent. In particular, some of the asymmetries have a direct relation to coefficients in the expansion, 2π Im c 11 10 , However, the calculation of the correlation between observables, which is necessary to include all of them in a global fit, seems easier with the global approach. Whether one method or the other gives more precise results has to be determined with an analysis including all systematic uncertainties.

Physics parameters and the t bW interaction
The measurement of the top decay amplitudes a λ 1 λ 2 can be interpreted in terms of limits on anomalous tbW interactions.
The most general effective tbW interaction arising from the addition of dimension-six operators to the SM Lagrangian can be parameterised as [23] using standard notation, with g the electroweak coupling, M W the W boson mass and q ν its four-momentum. In the SM, V L equals the Cabibbo-Kobayashi-Maskawa matrix element V tb 1, and the rest of couplings V R , g L and g R vanish at the tree level. For this general vertex, expressions of the W boson spin density matrix have been obtained in Ref. [24]. Matching our general expressions for the top decay amplitudes, obtained by general angular momentum conservation arguments, with the explicit calculations there, we obtain for top quarks up to a global normalisation factor that is irrelevant. Here, m t is the top quark mass and | q| the modulus of the W boson three-momentum in the top quark rest frame. The form factors A 0,1 , B 0,1 , C 0,1 and D 0,1 depend on the couplings in (21) and are given in Appendix B for completeness. For top anti-quarks the decay amplitudes (denoted here by a bar) are related to the top quark ones by It is well known [25] that there is a cancellation between anomalous contributions to helicity fractions when M W V R m t g L , as well as when M W V L m t g R . This cancellation stems from the Gordon identities that one can write for onshell t and b quarks, Neglecting the b quark mass, the combinations of couplings with M W V R = m t g L or M W V L = m t g R are equivalent to an interaction of the type ( p t + p b ) μ P L ,R , which does not contribute for any element of the W boson spin density matrix except for λ 1 = λ 1 = 0, because the product of p b with the W polarisation vectors of helicities ±1 vanishes. 2 Therefore, the effect in the W spin observables is residual, and given by the change in the partial width to λ 1 = 0 states. The cancellation is apparent if we define new couplings η 1,2 , ζ 1,2 as unitary rotations of the ones in (21), The dependence of the helicity fractions on η 1 and ζ 1 , for V L = 1 and g R = 0, is depicted in Fig. 1, at LO. (Nextto-next-to-leading order calculations of the helicity fractions are available [26], but since the difference with LO is smaller than the experimental uncertainty, we use the latter for consistency.) The helicity fractions are rather insensitive to η 1 , or, in other words, a combination with gives a very small effect. The remaining W boson spin observables exhibit the same behaviour. This can also be seen analytically. Neglecting the b quark mass, the form factors in Appendix B read Dd , which were shown to be redundant in Ref. [23]. Therefore, the insensitivity to these combinations of couplings can be viewed as insensitivity to these effective operators.
We observe that, for ζ 1 = 0, the anomalous contributions from η 1 to the B 0,1 , C 0,1 and D 0,1 form factors vanish, in agreement with our previous argument, but for A 0 and A 1 , precisely the form factors involved in F + 0 and F − 0 , they do not. Therefore, the measurement of these two untangled helicity fractions, for example using A z,z FB,EC , can break the degeneracy. For illustration, we plot in Fig. 2 the dependence of A z,z FB,EC on η 1 and ζ 1 , for V L = 1 and g R = 0. There is a sharp difference with the helicity fractions, and for this asymmetry the variations with η 1 and ζ 1 are similar. The same information can of course be obtained, in the global analysis, by the measurement of c 10 00 and c 12 00 .

Summary
In this work we have set the framework to extract all the relevant physical quantities from the measurement of the fourdimensional top decay distribution, that is, the top polarisation in three orthogonal directions and the five parameters determining the top decay distributions. Such demanding measurements may be possible in the near future with the LHC Run 2 data. With the full distribution, or with a suitable forward-backward-edge-central asymmetry, it will be possible to measure the two λ 1 = 0 untangled W boson helicity fractions F + 0 and F − 0 . The precision of the global fit to the tbW vertex might be substantially improved if these untangled helicity fractions are accurately measured in upcoming analyses.
in the first D function is precisely the first index in the second D function. Explicitly, we have

Appendix B: Expressions for top decay form factors
We collect here the expressions for the eight dimensionless form factors involved in the decay t → W b for a general tbW vertex, from Ref. [24]. Defining x W = M W /m t , x b = m b /m t , with m b the b quark mass, they are