Research Article  Open Access
Angular and CPViolation Analyses of Decays at Hadron Collider Experiments
Abstract
The branching fractions ratio between muon and tau lepton decay modes has shown intriguing discrepancies between the Standard Model prediction and measurements performed at BaBar, Belle, and LHCb experiments, a possible sign of beyond the Standard Model physics. Theoretical studies show how observables related to the differential decay distribution can be used to further constrain New Physics contributions, but their experimental measurements are lacking to date. This article proposes the measurement of angular and violating observables at hadron collider experiments, by exploiting approximate reconstruction algorithms using information from detectable finalstate particles only. The resolution on the phase space variables is studied using decays simulated in a forward detector geometry like LHCb. A method to correct the observable values for the reconstruction inaccuracies based on detector simulation is successfully tested on simulated data and the decrease in precision with respect to a perfect reconstruction is evaluated. The longitudinal polarization fraction and the violating observable can be measured losing a factor of 2 and 5 in precision, respectively. The extraction of angular distributions from the template fit selecting decays and associated systematic uncertainties are also discussed.
1. Introduction
Semileptonic decays, in which stands for one of the three charged leptons, have shown intriguing discrepancies between the Standard Model predicted ratio of branching fractions between muon and tau lepton decay modes [1], indicated as , and the measured values at BaBar [2], Belle [3–5], and LHCb [6, 7] experiments. This contrast could be a sign of New Physics contributions violating the Standard Model universality of leptonic interactions.
The measurement of observables related to the differential decay rate, other than , can shed new light on the observed anomalies, allowing to put complementary constraints on possible New Physics sources [1, 8–13]. However, the only measurement of these observables available to date is a preliminary result for the longitudinal polarization fraction in decays by the Belle experiment [14]which is consistent at 1.4 with the Standard Model prediction [11, 13].
Angular analyses of decays are challenging because finalstate neutrinos can not be reconstructed, implying that the meson rest frame is not precisely determined from the detectable part of the decay. This problem can be mitigated at factories, where the momentum of the meson can be determined from the known centerofmass energy of the collision and the complete reconstruction of the decay of the other meson produced in the interaction. On the contrary, at hadronic colliders the meson momentum is not constrained by the production mechanism since the centerofmass energy of the partonparton collision is unknown.
This article considers the possibility to measure the angular variable distributions of decays by exploiting reconstruction algorithms estimating the meson rest frame only from information related to the detectable finalstate particles, a situation of particular interest for hadron collider experiments like LHCb. The attainable precision on the phase space variables is studied by means of a simulation study set for a forward detector geometry which is detailed in Section 2. It is shown that observables related to the cosine of the polar angle of the meson in the helicity frame, , and the azimuthal angle between the () and () decay planes, , are suitable to be measured in the considered setup. It is shown that and distributions can be extracted using the Plot statistical technique [15] from the template fit selecting decays from background events.
The fully differential decay distribution is reviewed in Section 3 and the observables associated to the aforementioned phase space distributions introduced. These are the longitudinal polarization, the conserving and violating observables related to the angle distributions. The latter are especially interesting being a null test for the Standard Model, since violation in Cabibbofavoured quark transition is strongly suppressed by the CabibboKobayashiMaskawa mechanism.
In Section 4, a method to measure the considered observables while correcting the effect of reconstruction inaccuracies is presented and tested on simulated decays. The decrease in precision due to the use of the reconstruction algorithms is evaluated with respect to ideal measurements in which the phase space distributions are perfectly reconstructed. A discussion on the possible systematic uncertainties associated to the proposed measurements is reported in Section 5. The conclusions of the study are summarized in Section 6.
2. The Decay Reconstruction
2.1. Simulation Configuration
The capability of reconstructing the decay distribution using approximate reconstruction algorithms is studied on simulated semileptonic decays in a detector configuration analogous to the LHCb experiment [16].
Three decay chains are considered: , , and , along with their chargeconjugated decays. The flavour of the meson is determined by the charge of the detectable part of the lepton decay or by that of the pion produced in the decay. The production of mesons from protonproton collisions at a centerofmass energy are simulated using PYTHIA 8.1 [17, 18], their decay to the different final states are simulated by the EVTGEN package [19]. Stable particles are required to be within the nominal LHCb pseudorapidity acceptance , while charged particle momentum cuts and roughly reproducing the LHCb kinematic acceptance (estimated from [16]) have been tried but showed no significant effect on the subsequent studies. A minimum meson flight distance of 3 mm simulates the effect of a displaced vertex trigger requirement. The production and decay vertex positions of the meson have been smeared from their generated values according to Gaussian distributions reproducing the performance of the LHCb VELO detector [20, 21]: for production vertexes the Gaussian widths are 13 and 70 in the transverse and longitudinal directions, respectively, with respect to the beam; for decay vertexes they are 20 and 200 . For decays, a minimum tau lepton flight distance of 1 mm is applied as background rejection cut.
The ROOT package [22] is employed for data handling and graphics.
2.2. Rest Frame Approximate Reconstruction Algorithms
The rest frame reconstruction benefits from the knowledge of the flight direction from its production and decay vertexes, the latter determined by the track combination. Two strategies are considered in this study.
For decays in which a single neutrino is missing, the available information about the decay (the momentum of the detectable part of the decay, the meson flight direction, the and neutrino masses) determines the momentum up to a twofold ambiguity [23]. The two solutions correspond to the forward or backward orientation of the neutrino in the rest frame with respect to the flight direction. If the neutrino is orthogonal to the flight direction a unique, degenerate solution is found. This algorithm will be referred to as “full reconstruction.”
A different momentum approximation can be made assuming that the proper velocity along the beam axis, , of the detectable part of the decay is equal to that of the meson [6]. The magnitude of the momentum in terms of the visible decay system and the angle between flight direction and beam axis is set asThis approach will be referred to as “equal velocity” algorithm and it is applicable also to decays with two or more invisible particles, in which the invariant mass of the unmeasured part of the decay is unknown.
2.3. Resolutions on the Phase Space Variables
The decay is characterized by four degrees of freedom.^{1} Its phase space can be described by the following four kinematic variables: the invariant mass of the system , the cosine of the polar angle of the meson in the helicity frame , the cosine of the polar angle of the lepton in the system helicity frame and the azimuthal angle between the () and () decay planes , see Figure 1. In and helicity frames, the axis is defined by the direction of the and momenta in the rest frame, respectively.
The attainable precision on the four phase space variables is studied computing the resolution defined as the difference between the values measured using the reconstruction algorithms and the true values of the simulated events. Differences of dimensional quantities are divided by the true values.
The rest frame reconstruction for decays is achieved exploiting the full reconstruction algorithm. If a couple of solutions are found, one of the two is selected by random choice, while apparently unphysical configurations, due to experimental uncertainties, in which no momentum solution is available are discarded from the following study, these constituting the 32.7% of the simulated events. Regression techniques based on meson flight direction and magnitude to improve the solution decision [24] have been tried but showed limited improvement. The relative resolution on the momentum magnitude, obtained with the two reconstruction algorithms, is shown in Figure 2. The full reconstruction momentum resolution features a narrow, symmetric distribution peaked at zero, corresponding to events in which the momentum solution corresponding to the true orientation of the neutrino (forward or backward) was chosen, and a broader, asymmetric shape associated to events in which the momentum solution corresponding to the wrong neutrino orientation was assigned. The equal velocity reconstruction presents a more regular but wider distribution. The phase space variables describing the semileptonic decay are computed in the rest frame resulting from the estimated momentum. Their resolutions are reported in Figure 3: the and feature symmetric and unbiased distributions, the distribution is slightly asymmetric but almost unbiased and the relative even if asymmetric peaks at zero. Phase space variable resolutions obtained with the equal velocity algorithm are reported in Figure 4. Their distributions are wider than those resulting from the full reconstruction algorithm, since less information on the decay is employed.
For decays, in which the and vertexes determine the flight direction of the tau lepton, the full reconstruction algorithm is applied sequentially to the tau lepton and meson decays. First, the momentum is estimated from the visible system: if there are two momentum solutions one is chosen randomly. If no solutions are available, the momentum corresponding to the degenerate solution is assigned. Then, the momentum is calculated from the system using the estimated momentum: if there are two momentum solutions one is chosen randomly. If no solutions are available then the other, if any, momentum solution is tried, and the event discarded only if the momentum reconstruction is still impossible. This algorithm tries to retain the maximum information on the decay, however, it rejects 57.7% of the events. The estimated momentum is then used for computing and variables. The relative resolution on the momentum magnitude is shown in Figure 5 along with that obtained using the equal velocity algorithm, the latter being the narrower one. Phase space variables resolutions for full reconstruction algorithm are reported in Figure 6, which are to be compared to those obtained applying the equal velocity algorithm, see Figure 7. Comparing to the muon channel, the distributions are moderately wider, while and resolutions are significantly broader, since they directly depend on the leptonic part of the decay. The distributions are however still unbiased, while the ones are asymmetric and biased, especially for the equal velocity algorithm. Comparing the two algorithms, the distributions are basically equal, while the resolution is better for the full reconstruction one.
For decays, no information on the decay vertex is available and the equal velocity algorithm is applied. The relative resolution on the momentum magnitude is shown in Figure 8 and phase space variables resolutions are reported in Figure 9. The muon momentum is taken as tau lepton momentum for computing and variables. Comparing to the tau lepton hadronic decay channel, the distributions are similar to the more precise resolutions of the full reconstruction algorithm rather than to those obtained with the equal velocity algorithm. Thus, the knowledge of the tau lepton flight direction in the three pion decay mode is not able to add significant information to the decay reconstruction due to the increased ambiguity in the momentum determination.
Summarizing, and resolution distributions have been shown to be symmetric and unbiased for all the decay channels, and the related physical quantities are therefore suitable to be measured even at hadron collider experiments, making use of the presented reconstruction algorithms only. On the contrary, resolution distributions have been found to be biased for lepton decay channels. The measurement of observables depending on would therefore require special care and it is not further considered in this article.
2.4. Extraction of Angular Distributions from the Template Fit Selection
The selection of decays is a challenging task, especially at hadronic colliders. The impossibility of reconstructing all the finalstate particles prevents the direct use of invariant masses as discriminating variables and makes different decays with similar topology but additional unreconstructed particles difficult to distinguish from transitions. In fact, besides discriminating muon from tau lepton decay modes, decays must be separated from decays to , and other higher mass charm meson resonances and decays to double charm resonances in which one has a semileptonic decay. This is usually achieved by means of a template fit to a set of discriminating variables, in which shapes for each decay type are mainly determined from simulation [6, 7].
The extraction of distributions from the fit results can be done straightforwardly by means of the Plot statistical tool [15] only for angular variables independent from the discriminating ones. In this way the distributions are derived using no a priori information about them, but only from the discriminating variables. Distributions which are correlated with the discriminating variables can also be obtained in principle, but since they will depend directly on the construction of the template distributions, their extraction would need a specific statistical treatment and they would be more sensitive to fitrelated systematic uncertainties.
The possibility of deriving angular distributions from a realistic selection is checked by evaluating their correlations, computed as mutual information,^{2} with the set of the three discriminating variables used in [6], in which the detectable part of the leptonic decay, or , is used: the missing mass of the decaythe energy of the system in the rest frame , and , where the rest frame is estimated using the equal velocity algorithm. Correlation plots are presented in Figures 10, 11, and 12 for , and events, respectively. Since the discriminating variables depend on the leptonic part of the decay, correlations for and variables are found to be negligible; for correlations are high for the muon decay mode and small for the tau lepton one, because in the latter case the relationship is blurred by the extra neutrinos coming from the decay.
Detector reconstruction and event selection may introduce additional correlations between discriminating and angular variables, but efficiency corrections are able to subtract these effects. Perevent efficiency corrections are routinely applied in many particle physics analyses, usually obtained from highstatistics simulation samples.
Thanks to their small correlations with the discriminating variables, and distributions can be extracted directly from the template fit using the Plot statistical technique, allowing related observable measurements to be performed on “signalonly” and distributions.
3. The Decay Distribution
Maximum information about the decay is obtained from the fully differential decay distribution [9]in which the dependence on the angular variables , and has been made explicit. The decay is described by twelve angular coefficient functions , dependent on couplings, hadronic form factors and ; NF is a dependent normalization term. The angular coefficients are labelled according to the helicity combinations on which they depend: longitudinal (), transverse (), or mixed ().
The conjugate decay distribution follows from the application of the transformation to (4): the angles are now defined with respect to and antiparticles, and the inversion of the momenta correspond to a transformation and ,Angular terms proportional to and are sensitive to violation, being produced in the interference between amplitudes having different violating weak phases. The associated coefficients, , , and , are practically zero in the Standard Model [9]; therefore a nonzero measurement of these quantities would be a clear sign of beyond the Standard Model physics.
Due to the experimentally available limited statistics, it is useful to integrate the fully differential decay distribution described by (4) to obtain observables retaining specific parts of the decay information. An overview of interesting observables defined for the decay distribution can be found in [9, 12, 25]; the following section will focus on observables constructed from and variables, the most suitable quantities to be measured according to the simulation study presented in Section 2.
3.1. Integrated Distributions and Observables
According to the study detailed in Section 2.3, the best resolution is attained on the polar angle of the meson in the helicity frame, . The singlydifferential distribution over , obtained integrating the complete decay distribution described by (4) over all but the variable, isin which and represent the integrated longitudinal and transverse polarization fractions of the meson, satisfying ; the distribution takes the form of a secondorder polynomial in depending on one single observable ,The longitudinal polarization fraction is sensitive to scalar and tensor New Physics contributions to the quark transition effective Hamiltonian, rather than to vector or axialvector terms [9, 11]. Its ability to constrain New Physics contribution has been recently considered in [13, 26, 27].
Observables derived from dependent decay distributions are especially interesting being clean probes for New Physics violation. Trigonometric functions of the angle can be expressed in terms of the unit vectors orthogonal to the and decay planes in the meson rest frame,asso that observables which are coefficients of or can be extracted as tripleproduct asymmetries. This feature allows violating observables to be extracted by counting rather than by angular fits and will be exploited further on.
The singlydifferential distribution over is obtained by integrating (4)The violating observable is sensitive to vector and axial vector New Physics contributions but not to pseudoscalar ones [9]. It depends linearly on , while for the conjugated decay depends on , changing sign under transformation. The corresponding violating observable can be thus defined asExploiting the odd parity of the term, the observable can be isolated from the distribution described by (10) by defining the tripleproduct asymmetryThe sum of asymmetries measured for the two conjugated decays still represent a violating observable.
Terms proportional to in the full decay distribution are multiplied by and integrate to zero under . The tripleproduct asymmetry defined asis zero even in presence of New Physics, being this angular dependence related to the spin structure of the decay, in which the meson has spin one. The measurement of is therefore a useful crosscheck for the tripleproduct asymmetry measurement, allowing to assess possible biases or contamination from events in which the comes from a spin zero resonance decay, like the , or from a nonresonant system [25].
Observables related to the terms of the decay distribution can be extracted from the dependent angular distribution defined asThe violating observable is sensitive to all vector, axialvector and pseudoscalar couplings [9]. It depends linearly on , while for the conjugated decay depends on , not changing sign under transformation. The corresponding violating observable is thereforeStarting from the distribution reported in (14), a tripleproduct asymmetry equivalent to the observable can be defined asThe difference between asymmetries measured for the two conjugated decays represents a violation observable.
4. Measurement Method for Decay Distribution Observables
The nonnegligible width of the resolution on the angular variables, studied in Section 2.3, must be taken into account when measuring the corresponding observables, which can be biased from their actual value. In Section 4.1, it is shown how the longitudinal polarization can be extracted from maximumlikelihood fits to simulated events by parametrizing the detector response in and as a function of via a polynomial expansion. This way, the nonnegligible experimental resolution effect is subtracted, and the measured values are found compatible with the generated ones. The loss of sensitivity due to the experimental resolution is evaluated. Maximumlikelihood fits have been performed using the ROOFIT package [28].
The same method is then applied for the extraction of and observables, Section 4.2, but found to be successful only for decays, due to the too large uncertainties associated to the angle reconstruction for tau lepton decay modes.
Section 4.3 deals with tripleproduct asymmetries, which can be measured just by counting. In this case, the simulation is used to determine the proportionality factor between the violating observables and the associated reconstructed tripleproduct asymmetry, allowing to correct for the experimental resolution and to quantify the associated loss in precision.
4.1. Longitudinal Polarization
As a first step, a perevent weight is assigned to simulated decays in order to obtain a flat distribution in the generated values, for correcting the distortion due the applied geometry and selection requirements. Different longitudinal polarizations are generated by applying another perevent, polarization dependent, weight such that the generated distribution reproduces (7) for each value. Both weights are normalized in such a way that for each value the mean of the weights is one.
The two perevent weights are multiplied together, assuming the detector efficiency correction is independent of . This assumption has been checked to be valid for the presented simulation study. In a realcase analysis, the generation of events with varying longitudinal polarization should be done before applying the detector reconstruction, so that detector efficiency effects can be taken into account as a function of .
Simulated events are then divided in two samples: a test sample reproducing reconstructed decays with different longitudinal polarizations, and a second used to derive a Legendre polynomial expansion in and . This expansion is used as fit model to model to extract from a maximumlikelihood fit of the test sample. The orthogonality and completeness of Legendre polynomials is exploited to expand the reconstructed decay distribution in and asin which the coefficients are determined asand is the product of the two perevent weights applied. Given the simple dependencies, quadratic in and linear in , only Legendre polynomials up to the second order are sufficient to approximate the decay distribution. The use of a simple parametrization makes the maximumlikelihood fit of the decay distribution fast and robust.
The test samples contain, by choice, ten thousand events per decay mode, while the other samples are five times larger than the test one. This is equivalent to assume that, in a real measurement, the statistics of the simulation sample employed to derive the polynomial expansion is larger enough with respect to the data sample.
The sensitivity to the longitudinal polarization is studied by fitting the test samples using directly the angular distribution equation (7) or the polynomial expansions equation (17) for the three considered decays. The measured polarizations are reported in Table 1. Ideal measurements are simulated by fitting the angular distribution distribution described by (7) to a toy sample generated from the same distribution for varying values, with the same number of events of the test samples. These correspond to measurements made by a detector with perfect resolution, taken as reference to evaluate the decrease in precision due to the reconstruction algorithms employed. Results of these ideal measurements are reported in the last row of Table 1.

Longitudinal polarizations extracted using the true angular distributions are clearly biased towards values for which the distribution is flatter (it is uniform for ). Polynomial expansions allow to correctly measure the generated values within the uncertainties resulting from the the maximumlikelihood fit. The precision for different values with respect to the ideal case decreases by a factor 1.4–1.9 for the muon mode and a factor 1.5–2 for the decay. The precision is therefore similar for muon and tau lepton decay modes, as expected since the variable does not directly depend on the leptonic part of the decay. The exploitation of the tau lepton decay vertex information in the decay reconstruction does not increase the precision on , rather, a larger uncertainty is observed for this mode.
According to this simulation study, the polarization fraction of decays is measurable with the sole use of the employed reconstruction algorithm, with a maximum penalty in sensitivity of a factor 2. This permits an additional search for New Physics in decays complementary to the already measured ratio.
4.2. An Attempt to Directly Measure the and Observables
A simulation study analogous to the one set for the measurement is performed to check the possibility to simultaneously measure the and observables related to the distribution reported in (10). This case is more difficult partly because of the larger resolution on the angle, especially for the tau lepton decay mode, partly because this angular distribution is characterized by fast oscillations ( and terms) more sensitive to reconstruction inaccuracies.
A first study is carried out assuming that the distribution has a simpler form,in which the oscillations are wider. As explained in Section 3, this angular dependence is absent from the actual distribution, so that the angular coefficients , do not correspond to angular observables. They are introduced with the purpose of testing the extraction method already applied to . The fit model is derived from the reconstructed decay distribution by means of a polynomial expansion in , and : Legendre polynomials are used for and , while a Fourier series^{3} up to , terms is employed for the angle,in which the coefficients are determined asand is the product of the two perevent weights applied.
The measured and values using the distribution equation (19) and the polynomial expansions equation (20) are reported in Tables 2 and 3, respectively. Only results in which one of the two observables is zero are shown, since negligible differences in the observables extraction are seen when both and have nonzero values. Ideal measurements are also simulated as done for . Only decays are considered for the tau lepton decay mode.


The polynomial expansions recover the generated values within uncertainties, with a precision on decreased by a factor 2–2.2 for the muon mode and 4.4–4.9 for the tau lepton mode, and a precision on decreased by a factor 1.8–2 for the muon mode and 5.5–6.3 for the tau lepton mode. As a result, the polynomial expansion method proves to be effective but the decrease in precision for the tau lepton decay mode is important to note.
The simulation study is repeated for and using the distribution described by (10) and an analogous polynomial expansion. Unfortunately, the two observables are measurable only for the muon decay mode, the results of which are shown in Tables 4 and 5, with precisions on and observables decreased by a factor 2.9–3.2 and 2.6–2.7, respectively. The measurement is not possible on the tau lepton decay mode because the large uncertainty in the reconstruction completely flattens the angle distribution.


The application of the polynomial expansion method is in principle effective for measuring angle related observables. In practice it is successful only for the decay mode, where and observables can be measured; for tau lepton decay modes the extraction is not possible due to both the larger resolution on the angle and the form of the expected decay distributions. The method has not been attempted for and measurement because its application is complicated by the combined fit to and variables and the need for negativevalued fitting functions (following from the angular distribution described by (14)), which prevent the use of the standard maximumlikelihood fitting technique.
An alternative method for the measurement of violating observables, relying on counting rather than fitting, is explored in the next section.
4.3. TripleProduct Asymmetries
In Section 3.1 it was shown that violating observables related to angle decay distributions can be extracted by defining suitable tripleproduct asymmetries (TPAs). The imperfect reconstruction of the angle leads to an effective dilution of the asymmetries, but this experimental effect can still be subtracted exploiting simulated events, and in a simpler way than for decay angular distribution fits. Moreover, since the angle distribution is unbiased, a measured nonzero value for violating TPAs, even if not corrected for the experimental dilution, would anyway represent an observation of New Physics violation.
The subtraction of reconstruction effects consists in determining the relation between reconstructed TPAs and generated violating observables. The linear function allows to infer from the measured with an uncertainty given by error propagation,in which represent the loss in sensitivity to with respect to the uncertainty on the TPA.
The simulation study is set as follows. Simulated events are weighted to reproduce one of the angle decay distributions at generationlevel, as a function of the violating observables. TPAs are built from the reconstructed value of the angle; for the distribution reported in (14) the dependence is included to take into account uncertainties in the sign determination. Three values for the corresponding conserving quantities have been considered, but it is shown that they have no impact on the TPAs measurement. In fact, and terms still integrate to zero when computing asymmetries using reconstructed angles, since the angle resolution distribution is not biased. The linear relation between reconstructed asymmetries and generated violating observables allows to correct for the dilution effects and to determine the decrease in precision from the inverse of the slope of the straight line.
The study is carried out for , defined in equation (12), from the distribution equation (10), , defined in (16), from the distribution equation (14) and , defined in (13), from the distribution equation (19). The relations for the three decay modes are reported in Figure 13. They are the same for different values within uncertainties. From TPA definitions follow that for perfect reconstruction the factor is for and , one for . The decrease in precision from perfect reconstruction is summarized in Table 6 for the different asymmetries and decays.
