Measurement of the Absolute Branching Fraction of $D_s^+ \to \tau^+ \nu_{\tau}$ via $\tau^+ \to e^+ \nu_e \bar{\nu}_{\tau}$

Using a dataset of 6.32 fb$^{-1}$ of $e^+ e^-$ annihilation data collected with the BESIII detector at center-of-mass energies between 4178 and 4226 MeV, we have measured the absolute branching fraction of the leptonic decay $D_s^+ \to \tau^+ \nu_{\tau}$ via $\tau^+ \to e^+ \nu_e \bar{\nu}_{\tau}$, and find $\mathcal{B}_{D_s^+ \to \tau^+ \nu_{\tau}}=(5.27\pm0.10\pm0.12)\times10^{-2}$, where the first uncertainty is statistical and the second is systematic. The precision is improved by a factor of 2 compared to the previous best measurement. Combining with $f_{D_s^+}$ from Lattice quantum chromodynamics calculations or the $|V_{cs}|$ from the CKMfitter group, we extract $|V_{cs}|=0.978\pm0.009\pm0.012$ and $f_{D_s^+}= (251.1\pm2.4\pm3.0)$ MeV, respectively. Combining our result with the world averages of $\mathcal{B}_{D_s^+ \to \tau^+ \nu_{\tau}}$ and $\mathcal{B}_{D_s^+ \to \mu^+ \nu_{\mu}}$, we obtain the ratio of the branching fractions $\mathcal{B}_{D_s^+ \to \tau^+ \nu_{\tau}} / \mathcal{B}_{D_s^+ \to \mu^+ \nu_{\mu}} = 9.72\pm0.37$, which is consistent with the standard model prediction of lepton flavor universality.

Using a dataset of 6.32 fb −1 of e + e − annihilation data collected with the BESIII detector at centerof-mass energies between 4178 and 4226 MeV, we have measured the absolute branching fraction of the leptonic decay D + s → τ + ντ via τ + → e + νeντ , and find B D + s →τ + ντ = (5.27 ± 0.10 ± 0.12) × 10 −2 , where the first uncertainty is statistical and the second is systematic. The precision is improved by a factor of 2 compared to the previous best measurement. Combining with f D + s from Lattice quantum chromodynamics calculations or the |Vcs| from the CKMfitter group, we extract |Vcs| = 0.978 ± 0.009 ± 0.012 and f D + s = (251.1 ± 2.4 ± 3.0) MeV, respectively. Combining our result with the world averages of B D + s →τ + ντ and B D + s →µ + νµ , we obtain the ratio of the branching fractions B D + s →τ + ντ /B D + s →µ + νµ = 9.72 ± 0.37, which is consistent with the standard model prediction of lepton flavor universality. Leptonic decays of charged pseudoscalar mesons can provide accurate determinations of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements and a clean setting for tests of the lepton flavor universality (LFU). In the standard model (SM), the partial decay width of D + s → ℓ + ν ℓ (ℓ = e, µ, τ ) to the lowest order is given by [1] where G F is the Fermi coupling constant, f D + s is the decay constant parameterizing strong-interaction effects between the two initial-state quarks, V cs is the c → s CKM matrix element, and µ ℓ is the ratio of the ℓ + lepton mass to the D + s meson mass, m D + s .
From Eq. (1), it is clear that the ratio of the decay widths Γ D + s →τ + ντ /Γ D + s →µ + νµ only depends on µ τ and µ µ , and is equal to 9.75 ± 0.01 [2]. This is consistent with the experimental measurements [2] given the current relatively large uncertainties. In recent years, some hints of LFU violation in semileptonic B decays have been reported by different experiments [3][4][5][6][7]. It is argued that the violation may occur in c → s transitions due to the interference between different amplitudes [8] or the interactions with scalar operators [9]. Therefore, the improved precision on the measurements of D + s → ℓ + ν ℓ will be of great interest in testing LFU.
Given f D + s from Lattice quantum chromodynamics (LQCD) calculations [10], a precision measurement of the branching fraction (BF) of D + s → ℓ + ν ℓ (B D + s →ℓ + ν ℓ ) allows an accurate determination of |V cs | and a stringent test of the unitarity of the CKM matrix. Conversely, given |V cs | from the SM global fit [2], the B D + s →ℓ + ν ℓ value allows one to extract f D + s and verify the theoretical determinations of the decay constant.
Simulated data samples produced with a geant4based [26] Monte Carlo (MC) package, which incorporates the geometric description of the BESIII detector and the detector response, are used to determine detection efficiencies and to estimate backgrounds. The MC simulation includes the beam-energy spread and initialstate radiation (ISR). Moreover, a generic MC sample corresponding to 40 times the integrated luminosity of data is produced, including the production of opencharm final states, ISR production of vector charmonium states, and continuum processes. Details can be found in Ref. [18].
The absolute B D + s →τ + ντ is given by where i indicates the data samples at the six energy points, N i ST (N DT ) and ǫ i ST (ǫ DT ) are the ST (DT) yields and efficiencies, and B τ + →e + νeντ is the world average BF of τ + → e + ν eντ [2]. The ST yields N i ST are obtained for each energy point, while the DT yield N DT is obtained for the combined data sample because of limited statistics.
The ST D − s candidates are reconstructed from the above eleven hadronic decay modes with almost the same selection criteria as those in Ref. [29], the differences are described here. For D − s → K 0 S K − , the K 0 S candidate is not required to have a decay length larger than twice the vertex resolution. For D − s → π + π − π − , the K 0 S candidates are retained in the invariant mass of either π + π − combination. The specific ionization energy loss (dE/dx) in the main drift chamber (MDC) and the time-of-flight (TOF) information are combined and used for particle identification (PID) by forming confidence levels for the charged pion and kaon hypotheses (CL π , CL K ). Kaons are identified by requiring CL K > CL π , while pions except for those from K 0 S decays are demanded to satisfy CL π > CL K . To suppress background from Bhabha events, the relative-probability sum, , must be smaller than 2.0 (0.9) for n > 1 (n = 1), where n is the number of charged tracks for the ST mode.
The recoil mass against the tag In the fits, the D − s signal is modeled by the MC-simulated shape convolved with a Gaussian resolution function. The combinatorial background is described by a first to third order Chebychev function. For D − s → K 0 S K − , the peaking background from D − → K 0 S π − is incorporated with the MC-simulated shape and its yield is left free. The fit results for the data sample at E cm = 4178 MeV in the range of M ST ∈ [1.89, 2.04] GeV/c 2 are shown in Fig. 1 as an example. For each tag mode, the ST yield is obtained within a tag-dependent M ST window corresponding to approximately ±3σ of the D s mass, where σ is the fitted resolution of M ST . Following the same procedure, the ST efficiencies are estimated using the generic MC sample. Table I lists the sum (R) of the ratios of the ST yield over the ST efficiency from the six data samples.
In the presence of the ST D − s candidate, we select the signal candidate for D + s → τ + ν τ with τ + → e + ν eντ . Given the very small variation of detection efficiencies and backgrounds at different energy points, the six data samples are combined for further analysis. The DT candidates are required to have exactly one charged track in addition to the daughters of the tag side and that track must have a charge opposite that of the tag side decay. The track is also required to be associated with , and N Xe + νe DT ), the numbers of DT events after removing backgrounds (NDT), the DT efficiencies (ǫDT), and the obtained B D + s →τ + ντ . Both ǫ i ST in R and ǫDT include the relevant BFs for K 0 S , π 0 , η, η ′ , and ρ 0 decays. The N is fixed in the fit and its uncertainty will be treated later as a systematic uncertainty. For B D + s →τ + ντ , the first, second, and third uncertainties are statistical, tag-mode dependent systematic and tag-mode independent systematic, respectively. The statistical uncertainty from the ST yield in R is included in the second uncertainty of B D + s →τ + ντ . For the other numbers, the uncertainties are statistical only. Mode The pairs of arrows denote the signal regions.
a good electromagnetic calorimeter (EMC) shower, as described in Ref. [29]. To identify the positron, the combined dE/dx, TOF, and EMC information is used to determine a CL e . We assign the track as a positron if it satisfies CL e > 0.1% and CL e /(CL e + CL π + CL K ) > 0.8. The candidate track is further required to have a mo-mentum in the MDC greater than 0.2 GeV/c and a ratio of the energy deposited in the EMC to the momentum greater than 0.8.
To effectively discriminate signal from background, we adopt the variable E tot extra following Refs. [11,12]. It is the total energy of the good EMC showers [31], excluding those associated with the ST selection and those within 5 • of the initial direction of the positron. The latter eliminates energy associated with final state radiation (FSR) from the positron track. We do not exclude soft photons or π 0 s originating directly from the D * s . The E tot extra distributions of the DT candidates for various tag modes in the combined data sample are shown in Fig. 2. Entries in the lowest bin include those events without any extra good EMC showers. We determine the background yield by a fit to the region E tot extra > 0.6 GeV, where the signal is negligible, and extrapolate the backgrounds into the signal region using MC-derived shapes. The signal yield is then determined in the region E tot extra < 0.4 GeV by statistically subtracting the expected backgrounds from the DT events seen in data. This procedure is insensitive to the signal shape, except for the inefficiency introduced by the definition of the signal region.
The backgrounds in the E tot extra distributions can be divided into three categories. The first one is the non-D − s background with an incorrectly reconstructed ST D − s . The second is the D + s → K 0 L e + ν e background, which survives when the K 0 L passes through the detector without decaying or significantly interacting. The third is the D + s → Xe + ν e background, which is dominated by the six semileptonic decays D + s → ηe + ν e , η ′ e + ν e , φe + ν e , f 0 (980)e + ν e , K * (892) 0 e + ν e , and K 0 S e + ν e . The latter two GeV/c 2 ). For the tag modes with neutral daughters, the resolution difference between data and MC simulation (called data-MC difference) has been corrected. The shape of the D + s → K 0 L e + ν e background is modeled by the MC-derived shape corrected by a two dimensional (polar angle and momentum) data-MC difference for the K 0 L detector response. These correction factors are obtained by using a control sample of D 0 → K 0 L π + π − from 2.93 fb −1 of data collected at E cm = 3.773 GeV [32]. The background yield is calculated with B D + s →K 0 e + νe = (3.25 ± 0.38 ± 0.16) × 10 −3 , quoted from our previous work [33]. The peaking background from D − → K 0 S π − is present only for the D − s → K 0 S K − tag mode and its yield is estimated from the M ST fit. The yield of the D + s → Xe + ν e background is left free, with the shape extracted from the MC simulation with the individual BFs for the six background channels fixed as B D + s →ηe + νe = (2.32 ± 0.08)% [34][35][36], B D + s →η ′ e + νe = (0.80 ± 0.07)% [34][35][36], B D + s →φe + νe = (2.37 ± 0.11)% [36][37][38], B D + s →f0(980)e + νe, f0(980)→ππ = (0.30 ± 0.05)% [39,40], B D + s →K * (892) 0 e + νe = (0.21 ± 0.03)% [33,36], and B D + s →K 0 e + νe = (0.34 ± 0.04)% [33,36]. Moreover, the MC-based shapes have been further weighted by the individual ST yields at various energy points. Table I lists the obtained DT yields and DT efficiencies, where the latter are evaluated from the generic MC sample. In the BF measurement, most uncertainties related to the ST selection are cancelled. The remaining systematic uncertainties are divided into two cases. The first case is from tag-mode dependent systematic uncertainties.
Systematic uncertainties in the ST yield are examined by changing the fit range, the signal and background shapes, and the bin size and the background fluctuation of the fitted ST yield. The alternative fit range is chosen as [1.895, 2.035] GeV/c 2 , corresponding to approximately 1σ of M ST reduced from both sides of the nominal range. The nominal signal shapes obtained from the generic MC sample are replaced with those from the signal MC sample. The background shape is changed to a different order of the Chebychev function. The bin size is doubled or halved. For each variation, the efficiency-corrected yields are found to be consistent. The differences in the ratio of the ST yield over the ST efficiency for a given ST mode for all variations, and the background fluctuation of the fitted ST yield, are weighted by the ST yields in various data samples, and added in quadrature. The resulting overall systematic uncertainty of the ST yield is 0.61%.
Tag bias related to the ST selection arises from different event environments (e.g., charged and neutral multiplicities). The ratios of the ST efficiencies from the generic MC sample and the signal MC sample for various tag modes are examined. The difference between the two ST efficiencies is weighted by the ST yields in various data samples. The systematic uncertainty due to tag bias is assigned to be 0.26%.
For the non-D − s background, we replace the distributions from the M ST sidebands with those for the background events in the M ST signal region from the generic MC sample. The change of the measured BF, 0.07%, is assigned as the systematic uncertainty. In the nominal fit, the MC-based correction signal events in the M ST sideband regions are normalized using the ratio of signal yields in the M ST distributions of data and the generic MC sample. Alternative fits are performed 10 4 times with the ratio varied according to a random Gaussian sampling based on its statistical error. The distribution of the relative difference on the DT yield is fitted by a Gaussian function, and the width of 0.07% is taken as the systematic uncertainty. Here, and below where this method is also used, no significant biases in the means are observed. The systematic uncertainty from the D − → K 0 S π − background for the D − s → K 0 S K − tag mode is assigned to be 0.05% with a similar technique. The uncertainty from the limited MC statistics is 0.31%. The quadratic sum of all the tag-mode dependent systematic uncertainties is 0.74%.
The second case is from tag-mode independent systematic uncertainties. The systematic uncertainty associated with the E tot extra < 0.4 GeV signal region is estimated using the DT events of D + s → π + (π 0 )η as the control sample. MC studies have verified that the efficiency for the requirement of E tot extra < 0.4 GeV in the control sample is consistent with that in our signal decay. The systematic uncertainty is assigned to be 1.0% considering the data-MC difference of the efficiencies.
The systematic uncertainty related to extra charged tracks is studied with the DT events of D + s → π + φ(→ K + K − ) and D + s → K +K * (892) 0 (→ K − π + ). The data-MC difference of the efficiencies, 1.4%, is conservatively taken as the corresponding systematic uncertainty.
The systematic uncertainties associated with the e + tracking [40] and PID efficiencies are studied with a control sample of radiative Bhabha events. The two dimensional (momentum and the polar angle) efficiencies in the control sample are re-weighted to match those in our signal decay. The data-MC difference on e + tracking (PID) efficiency, 0.3% (0.3%), is assigned as the relevant systematic uncertainty.
The systematic uncertainty from the shape of D + s → K 0 L e + ν e background is estimated by replacing the corrected shape of E tot extra with the uncorrected one. The change of the measured BF is negligible. We also examine the size of the D + s → K 0 L e + ν e background by sampling the quoted BF [33] 10 4 times based on a Gaussian distribution given by its uncertainty. The width of the relative difference of the fitted DT yield, 1.2%, is assigned as the systematic uncertainty.
The uncertainty coming from the fixed D + s → Xe + ν e background shape is examined by varying the proportion of each of the six main background modes via sampling their BFs [33][34][35][36][37][38][39][40] 10 4 times, based on Gaussian distributions given by their uncertainties. The width of the relative difference of the DT efficiency, weighted by the ST yields in various data samples, 0.1%, is taken as the systematic uncertainty.
The uncertainty due to FSR effect is checked with radiative Bhabha events. The data-MC difference of the efficiencies for the requirement of including FSR photons, 0.5%, is taken as the systematic uncertainty. The uncertainty from the quoted BF for τ + → e + ν eντ is 0.2% [2]. The quadratic sum of all tag-mode independent systematic uncertainties is 2.2%.
Based on our result of f D + s |V cs | and the measured f D + |V cd | in Ref. [41], along with |V cd /V cs | = 0.23259 ± 0.00049 from the SM global fit [2], it yields f D + s /f D + = 1.244 ± 0.017 ± 0.021, which is consistent with the LQCD calculation [10] within 2.4σ. Alternatively, taking f D + s /f D + = 1.1783 ± 0.0016 calculated by LQCD [10] as input, we extract |V cd /V cs | 2 = 0.049 ± 0.001 ± 0.002, which agrees with the value obtained from the CKM fitter within 2.3σ.
In summary, with data samples corresponding to an integrated luminosity of 6.32 fb −1 collected at center-ofmass energies between 4178 and 4226 MeV, we present a precise measurement of the absolute BF for D + s → τ + ν τ . The precision is improved by a factor of two compared to the previous best measurement [18]. Taking inputs from the SM global CKM fit and LQCD separately, the decay constant f D + s and the magnitude of the CKM matrix element |V cs | are also extracted individually; all are the most precise results to date. Combining our result with the world average B D + s →τ + ντ and B D + s →µ + νµ , we test the LFU in τ -µ flavors, and no LFU violation is found with the current precision.