Loading [MathJax]/jax/element/mml/optable/BasicLatin.js

Feasibility study of CP violation in τKSπντ decays at the Super Tau Charm Facility

Cited by

1. Achasov, M., Ai, X.C., An, L.P. et al. STCF conceptual design report (Volume 1): Physics & detector[J]. Frontiers of Physics, 2024, 19(1): 14701. doi: 10.1007/s11467-023-1333-z
2. Chen, F.-Z., Li, X.-Q. CP asymmetries in τ→KSπντ decays[J]. International Journal of Modern Physics A, 2024. doi: 10.1142/S0217751X24420107
3. Pich, A.. Tau physics opportunities at the super tau-charm facility[J]. International Journal of Modern Physics A, 2024. doi: 10.1142/S0217751X24420028
4. Li, H., Luo, T., Shi, X. et al. Feasibility study of Ds+→τ+ντ decay and test of lepton flavor universality with leptonic Ds+ decays at STCF[J]. European Physical Journal C, 2022, 82(4): 310. doi: 10.1140/epjc/s10052-022-10265-5
5. Li, Z.-Y., Dai, A.-X., Xie, J.-J. Electromagnetic Form Factors of Λ Hyperon in the Vector Meson Dominance Model and a Possible Explanation of the Near-Threshold Enhancement of the Reaction[J]. Chinese Physics Letters, 2022, 39(1): 011201. doi: 10.1088/0256-307X/39/1/011201
6. Chen, F.-Z., Li, X.-Q., Peng, S.-C. et al. CP asymmetry in the angular distributions of τ → KS πντ decays. Part II. General effective field theory analysis[J]. Journal of High Energy Physics, 2022, 2022(1): 108. doi: 10.1007/JHEP01(2022)108
7. Zhou, X.. Experimental program at super tau-charm facility[J]. Proceedings of Science, 2021.

Figures(7) / Tables(3)

Get Citation
Haoyu Sang, Xiaodong Shi, Xiaorong Zhou, Xianwei Kang and Jianbei Liu. Feasibility study of CP violation in τKSπντ decays at the Super Tau Charm Facility[J]. Chinese Physics C. doi: 10.1088/1674-1137/abeb07
Haoyu Sang, Xiaodong Shi, Xiaorong Zhou, Xianwei Kang and Jianbei Liu. Feasibility study of CP violation in τKSπντ decays at the Super Tau Charm Facility[J]. Chinese Physics C.  doi: 10.1088/1674-1137/abeb07 shu
Milestone
Received: 2020-12-15
Article Metric

Article Views(1288)
PDF Downloads(27)
Cited by(8)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Feasibility study of CP violation in τKSπντ decays at the Super Tau Charm Facility

    Corresponding author: Xianwei Kang, kxw198710@126.com
  • 1. State Key Laboratory of Particle Detection and Electronics, Hefei 230026, China
  • 2. University of Science and Technology of China, Hefei 230026, China
  • 3. Key Laboratory of Beam Technology of Ministry of Education, College of Nuclear Science and Technology, Beijing Normal University, Beijing 100875, China
  • 4. Beijing Radiation Center, Beijing 100875, China

Abstract: We report a feasibility study for violation in τKSπντ decays at a Super Tau Charm Facility (STCF). With an expected luminosity of 1 ab1 collected by STCF per year at a center-of-mass energy of 4.26 GeV, the statistical sensitivity for CP violation is determined to be of order 9.7×104 by measuring the decay-rate difference between τ+KSπ+ˉντ and τKSπντ. The analysis is performed using a reliable fast simulation software package, which can describe the detector responses properly and vary the responses flexibly for further optimization. Moreover, the energy-dependent efficiencies for reconstructing τKSπντ are presented. The expected CP sensitivity is proportional to 1/L in the energy region from 4.0 to 5.0 GeV. The sensitivity of CP violation is of order 3.1×104 with 10 ab1 integrated luminosity, which is equivalent to ten years' data taking in this energy region at STCF.

    HTML

    I.   INTRODUCTION
    • The Cabibbo-Kobayashi-Maskawa (CKM) matrix can accommodate CP violation (CPV) in a natural way with a complex phase in the quark mixing sector [1, 2]. Experimentally, CPV has been observed in the meson sector, first in the K meson [3-5], subsequently in the B meson [6-11] and most recently in charm meson decays [12]. All results to date are consistent with the predictions of the CKM mechanism in the Standard Model (SM).

      However, the origin of CPV has remained an unsolved problem, as it is not clear whether the CKM mechanism is the unique source. The matter-antimatter asymmetry of the universe [13, 14] indicates that there must be non-SM CPV sources, and such additional sources are indeed incorporated in many extensions of the SM [15, 16]. Currently, CPV has not been observed in the lepton sector, and any significant observation would be a clear indication of physics beyond the SM. In an intriguing scenario [17], baryogenesis is argued to be mainly driven by leptogenesis, and then CPV is required in leptodynamics. Thus, exploration of CPV in the lepton sector provides a different and complementary landscape, at least. In the lepton sector, τ decay is a good place to seek for CPV either within or beyond the SM [18-22], since it has abundant hadronic decay channels with sizable branching ratios. Furthermore, half of the matrix element is purely electroweak and the hadrons involved are generated by quark-antiquark pairs, such that they are factorized, avoiding the complications of pure hadronic decay.

      In this paper, we will discuss the specific channel τKSπντ and concentrate on CPV and related New Physics (NP) searches. This channel is also of great importance for the measurements of form factors and extraction of the |Vus| element in the CKM matrix. It also has rich intermediate resonances which provide an important playground for non-perturbative QCD studies. Experimentally, a rich source of τ leptons comes from electron-positron colliders, and the physics programs related to τ decays are usually performed at a B-factory or Z-factory. There are a lot of prospects for the potential of τ decays at future such factories [23]. However, τ decays are rarely investigated in the tau-charm energy region, due to the low statistics available and the difficulty in reconstruction. A feasibility study of τ decays in the tau-charm region will be of great importance to complete the global picture of whole τ physics. The proposed STCF [24] in China is designed to have a peak luminosity of >0.5×1035 cm2s1 at s=4 GeV. It is a symmetric electron-positron beam collider designed to provide e+e interactions at c.m. energy s from 2.0 to 7.0 GeV. The τ leptons are selected via the process e+eτ+τ, where the cross-section of e+eτ+τ peaks at 3.5 nb, around s=4.26 GeV [25]. The number of τ pairs produced at STCF is estimated to be 3.5×109, with an expected integrated luminosity of 1 ab1 per year.

      This paper is organized as follows. In Sec. II, we elaborate the physical significance of investigating the channel τKSπντ. In Sec. III, the detector concept for STCF is introduced, as well as the Monte Carlo (MC) samples used for this study. The event selection for the signal process and detector optimization are presented in Sec. IV and Sec. V respectively. Section VI discusses the results and Sec. VII gives our conclusion.

    II.   PHYSICS PROGRAM OF τKSπντ
    • Within the SM, there is no direct CPV in hadronic τ decays at tree level in weak interactions. However, the well-measured CPV in KLπl±ν (l=e,μ) produces a difference in Γ(τ+KLπ+ˉντ) vs. Γ(τKLπντ) due to the K0ˉK0 oscillation. The same asymmetry also appears in Γ(τ+KSπ+ˉντ) vs. Γ(τKSπντ), and is calculated to be [26, 27]:

      ACP(τKSπντ)=Γ(τ+KSπ+ˉντ)Γ(τKSπντ)Γ(τ+KSπ+ˉντ)+Γ(τKSπντ)=(0.33±0.01)%.

      (1)

      Experimentally, the BaBar experiment has found evidence for CPV in τ decays [28]:

      ACP(τKSπντ[0π0])=(0.36±0.23±0.11)%,

      (2)

      and there is a 2.8 standard deviation difference from the theoretical prediction in Eq. (1). It should be noted that the K0ˉK0 oscillation leads to the same value of CPV for inclusion of any π0s in τ decays from the SM prediction.

      Motivated by the above CPV disagreement, various NP scenarios have been proposed, e.g. introducing a non-standard tensor interaction [29, 30]. In Ref. [30], a very large angular-weighted CPV is predicted additionally, with such a tensor interaction. However, those conclusions are questioned in Ref. [31], where the authors point out the tensor contributions used in Refs. [29, 30] are overestimated due to the lack of the well-known Watson theorem of final-state interactions, resulting in very large values of the imaginary parts of the Wilson coefficients to explain the CP anomaly; those large values of Wilson coefficients contradict the bound from the neutron electric dipole moment (EDM) and DˉD mixing (unless NP occurs below the electroweak breaking scale). The standpoint of Ref. [31] is adopted later in Refs. [32, 33], and especially, Ref. [33] finds that once the combined constraints from the branching ratio and decay spectrum are imposed, it is hard to explain the experimental values of CPV. But again, Ref. [31] is questioned by the authors of Ref. [34] – they point out that the tensor operator used in Ref. [31] corresponds only to a specific way of imposing gauge invariance, and instead, introducing a new type of tensor operator could account for the CP anomaly while evading bounds from the neutron EDM and keeping the extraction of |Vus| from exclusive τ decay unaffected. In brief, BaBar's CPV observables lead to a series of theoretical discussions, and those theoretical discussions are in a state of confusion. Therefore, a new measurement on the experimental side is essential to approach any conclusive statement on the NP signal. We stress the importance of precision here.

      Moreover, as discussed in Ref. [19], probing ACP(τ[Kπ]ντ), ACP(τ[K2π]ντ) and ACP(τ[K3π]ντ) separately can help to establish the existence of new dynamics, since the global CPV as expressed in Eq. (2) is often much reduced, and one needs to understand the basis of the observed data on τKSπντ[0π0] vs. τ[Kπ]ντ, τ[K2π]ντ and τ[K3π]ντ. In the SM, the same CPV holds for decays including any number of π0 mesons, due to the sole source from K0ˉK0 oscillation. However, the situation will change in the presence of NP. For example, the two-Higgs doublet model has been used to show its influence on the different observable sets for the channels τKππ0ντ and τKπντ [35, 36]. Certainly, a more complicated impact of resonances appears in the multibody channels, and more data is needed in this case. We also note that in four- and five-body decays, one can also access CPV by T-odd observables [37-39].

      Apart from the above BaBar measurement, the CLEO [40] and Belle [41] experiments have also focused on the CPV that could arise from a charged scalar boson exchange [42]. This type of CPV can be detected by measuring the τ± decay angular distributions. The CPV is found to be compatible with zero with a precision of O(103) [41]. However, current experimental sensitivity cannot make a conclusion on the CPV from τ decay, due to the large uncertainty. Therefore, a higher-precision result is required to hunt for signals of NP. At Belle II, CPV is studied in the angular distribution of τKSπντ decays, and the precision is expected to be 70 times improved with a 50 ab1 data sample [43].

      In order to capture a faint NP signal, or instead, to constrain a NP model, accurate knowledge of the form factor is also an essential input [30, 32]. The process τKsπντ decay in the SM is described by vector and scalar form factors. The vector form factor receives mainly the contribution of K(892), while K(1410) is needed in the higher-energy region tail. For the scalar form factor, there is no clear dominance of a single resonance. Currently the most sophisticated theoretical description of the form factors is obtained by the model-independent dispersive representation, imposing constraints from chiral symmetry and their asymptotic QCD behavior [44-46]. However, other descriptions are available [36, 47]. In Ref. [36], the form factors are calculated up to one-loop level using the chiral Lagrangian. The resulting shape differs from the dispersive one in Ref. [44]. Such a difference may be attributed to the perturbative vs. non-perturbative treatment. Reference [47] aims to achieve a satisfactory description of experimental data in a phenomenological way, and then a superposition of Breit-Wigner functions is used. As already commented in Ref. [31], the resulting phase does not vanish at threshold and also violates the Watson theorem before the inelasticity sets in. In short, the shapes of the P-wave form factor agree, while the S-wave shapes can differ dramatically in different model calculations. However, from a pragmatic point of view, all these variants of form factor could describe the data for the total mass spectrum almost equally well. A partial-wave analysis will certainly help to pin down this issue, which again needs more high-quality data. This point has been noted by the Belle collaboration [47]: in future, more data on the invariant mass spectrum of Ksπ combined with angular analysis will elucidate the nature of the scalar form factor and check various theoretical approaches. At STCF, high statistics of τKSπντ can put strong constraints on the vector and scalar form factors. Meanwhile, the precision in the determination of the form factor parameters can be improved.

      On the other hand, |Vus| can be extracted from the measurement of the τKSπντ decay. Indeed, the decay rate of τKSπντ can be expressed as [48]:

      Γ=G2FNC2K,τSEW(|Vus|fK0π+(0))2IK(1+δEM+δSU(2))2,

      (3)

      where N is a normalization coefficient (N=m3τ/(48π3)), GF is the Fermi constant, CK,τ is the Clebsch-Gordan coefficient, IK is the phase space integrals, SEW is the electroweak short-distance, δEM is the electromagnetic long-distance, δSU(2) is the isospin breaking corrections, and fK0π+(0) is the form factor at zero momentum transfer. According to the above equation, the measurement of |Vus| requires: i) accurate measurement of Γ; ii) precise calculation of IK; iii) good knowledge of the radiative corrections, which are SEW, δEM and δSU(2); and iv) a determination of the value of fK0π+(0). In general, with higher precision for τKSπντ branching fraction measurements, the error on |Vus| will be reduced [49].

      In short, the physical importance of measuring τKSπντ is discussed from three aspects: i) probing the CPV of this process and testing several NP scenarios; ii) the resonance parameters (mass and width) can also be determined, which is itself an important topic of the hadron spectrum, assisting us to understand more of the strong impact of intermediate resonances and their interference; we also notice that a large strong phase will be rendered in the resonance region, which together with the weak phase is a requisite for generating large CPV; and iii) the |Vus| element can also be determined, which is worth doing, at least as another cross check for the extraction from exclusive decay modes. For a review see Ref. [21]. As has been stressed, to make a conclusive remark the precision in the measurement is an essential ingredient. We indeed need to enter the era of high precision. A detailed analysis of points ii) and iii) will be presented in a future study. In this paper, we will investigate the reconstruction strategy for τ decays at STCF, and will focus on the sensitivity estimation of CPV by measuring the decay rate difference between τ+KSπ+ˉντ and τKSπντ. Moreover, to make the best use of the large statistics, it is necessary to have a compatibly sophisticated detector at STCF to precisely detect and measure particles, where the process τKSπντ can serve as a benchmark physics process to offer constraints to the detector design.

    III.   DETECTOR AND MC SIMULATION
    • The STCF detector is a general purpose detector designed for an e+e collider. It includes a tracking system composed of inner and outer trackers, a particle identification (PID) system with 3σ charged K/π separation up to 2 GeV/c, an electromagnetic calorimeter (EMC) with excellent energy resolution and good time resolution, a super-conducting solenoid, and a muon detector (MUD) that provides good charged π/μ separation. The detailed conceptual design for each sub-detector can be found in Ref. [50]. Currently, the STCF detector and the corresponding offline software system are in the research and development (R&D) phase, and it is necessary to have a reliable simulation tool which can access the physics reaches. A fast simulation tool for STCF has been developed [50], which takes the most common event generator as input to perform a fast and realistic simulation. The simulation includes resolution and efficiency responses for the tracking of final state particles, the PID system and kinematic fit related variables. The fast simulation also provides a flexible interface for adjusting the performance of each sub-system, which can be used to optimize the detector design according to physical requirements.

      The process τKsπντ is studied for 1 ab1 integrated luminosity at s=4.26 GeV, where the cross-section of e+eτ+τ is the largest. Unless specified, the charge conjugate decays are always implied throughout the analysis. The MC events for e+el+l(l=e,μ) and e+eγγ are generated with BABAYAGA [51, 52], and hadronic production processes (e+eqˉq)(q=u,d,s,c) are generated with LUNDARLW [53]. The e+eτ+τ process is generated with KKMC [54], which implements TAUOLA to describe τ decays inclusively. All the MC samples described above are generated according to the expected amount of 1 ab1 integrated luminosity and no CPV is assigned. Passage of the particles through the detector is simulated by the fast simulation software [50].

      The signal process τKSπντ is generated with vector and scalar configurations of KSπ, and the parameterized spectrum is described by:

      dΓds1s(1sm2τ)2(1+2sm2τ)P(s)×{P2(s)|FV|2+3(m2KSm2π)2|FS|24s(1+2sm2τ)},

      (4)

      where s is the squared invariant mass of KSπ, and mτ,mKS and mπ are the masses of τ,KS and charged π from PDG [55]. P(s) is the momentum of KS in the (KSπ) c.m. frame, given by:

      P(s)=(s(mKS+mπ)2)(s(mKSmπ)2)2s.

      (5)

      FS and FV are the scalar and vector form factors to parameterize the amplitudes of K0(800), K(892) and K(1410):

      FS=aK0(800)BWK0(800),

      (6)

      FV=BWK(892)+aK(1410)BWK(1410)1+aK(1410),

      (7)

      where BW denotes the Breit-Wigner function, and aK0(800) and aK(1410) are complex coefficients for the fractions of the K0(800) and K(1410) resonances as presented in Ref. [47].

    IV.   EVENT SELECTION AND ANALYSIS
    • The signal events are selected with τ+ decay to leptons, τ+l+νlˉντ(l=e,μ), denoted as the tag side, and τKSπντ with KSπ+π, denoted as the signal side. At least four charged tracks are required in one event after passing fast simulation, with only the requirement of acceptance applied, i.e. |cosθ|<0.93, where θ is the polar angle with respect to the beam direction. The KS candidates are selected from pairs of oppositely charged tracks, which satisfy a vertex-constrained fit to a common point. The two charged tracks with minimum χ2 of vertex fit are assumed to be pions produced from KS. The KS is required to have an invariant mass in the range 0.485<Mπ+π<0.512 GeV/c2. Furthermore, considering the finite decay length of KS [55], the significance of the flight length of KS candidates is required to be larger than 2 and the flight length of KS should be larger than 0.5 cm, as shown in Fig. 1. Background processes with the same final states as the signal process but which contain no KS, such as τπ+ππντ, can be significantly suppressed by the above selections.

      Figure 1.  (color online) Flight length of KS candidates for signal and background, where the solid line in red is the signal process, the square-hatched area in blue is hadronic final states e+eqˉq (hadrons), and the diagonally hatched area in yellow is other τ decay channels in e+eτ+τ, excluding the signal process (ditau), respectively. MC samples are normalized to 1 ab1 integrated luminosity and all other selection criteria have been applied.

      After KS selection, a stricter vertex requirement is performed on other charged tracks than the pions in KS decay, and the number of these tracks satisfying the vertex requirement should equal 2. The lepton candidate on the tag side and the bachelor pion on the signal side will be further selected. For electron candidates, E/p is required to be larger than 0.8, where E is the deposited energy in the EMC and p is the momentum of the charged track. In fast simulation, there are two efficiency curves provided for muon identification. One is from the responses of the reference detector [56], e.g. a mixing requirement with P(μ)>0.001, P(μ)>P(K) and P(μ)>P(e) in the PID system, where P(X) is the probability of a candidate identified as particle X, 0.1<E<0.3 GeV in the EMC, and a layer requirement in the MUD. The efficiency curve is represented by the black dots in Fig. 2. The second efficiency curve is the response from preliminary R&D results for the MUD detector, where three kinds of π/μ separation levels are provided, 3%, 1.7% and 1%, according to different requirements on π/μ responses. The efficiency curve of muon identification of the STCF MUD is also shown in Fig. 2 by the red dots, for the 3% π/μ mis-identification requirement. After the selection of lepton candidates, the remaining track is identified as a pion if it satisfies P(π)>P(K) and P(π)>P(p). Finally, a signal event is required to have one lepton, one pion and one KS.

      Figure 2.  (color online) Detection efficiency of muon and the mis-identification rate of π/μ for two efficiency curves provided by fast simulation.

      To suppress the background from τKSπ(1π0)ντ processes, π0 is reconstructed by selecting two good photons, and the event is rejected if the invariant mass of any two photons is located in the region between 0.12 and 0.15 GeV/c2. The good photon candidates are selected with the efficiency sampled from fast simulation.

      After the above selection, the distribution of KSπ mass spectrum from generic MC is shown in Fig. 3(a). There are huge backgrounds from other τ decay channels in e+eτ+τ, excluding the signal process and the hadronic final states containing multiple π final states. The selection efficiency of the signal and the suppression rate of backgrounds with the above selection criteria are listed in Table 1.

      Figure 3.  (color online) Invariant mass spectra of KSπ with combined e-tag and μ-tag. (a) Distribution without detector optimization or likelihood requirement yL. (b) Distribution with likelihood requirement but no detector optimization. (c) Distribution with likelihood requirement and detector optimization. MC samples are normalized to 1 ab1 integrated luminosity.

      Processεraw(%)εyL(%)εyL&optimization(%)
      e+eτ+τ, τ+l+νlˉντ,τKSπντ22.80±0.0222.22±0.0232.82±0.02
      e+eτ+τ(excluding signal process)0.1820.04030.0422
      hadronic final states0.3040.02110.0179

      Table 1.  Selection efficiencies for signal and background processes, with neither detector optimization nor likelihood requirement εraw, with no detector optimization but applying likelihood requirement εyL, and with both detector optimization and likelihood requirement εyL&optimization. The errors are statistical only. For background processes, the statistical uncertainties are too small and can be neglected.

      To further suppress the background events, a likelihood ratio yL is used, defined by

      yL(x)=LS(x)LS(x)+LB(x),

      (8)

      where LS and LB are the likelihood functions for signal and background events, respectively. x is a set of variables used for likelihood. Each likelihood function LS/B is the product of the probability density function (PDF) of the input variables, defined by

      LS/B=nvark=1PS/B,k(xk),

      (9)

      where PS/B,k is the signal/background PDF of the k-th input variable xk. For this analysis, the set of variables x includes the number of neutral clusters, momenta of decay products of KS, decay length of KS, mass spectrum of KS, χ2 of KS from secondary vertex fit, E/p ratio of the electron, momentum of μ, cosine of the polar angle of μ, and momenta of π and KSπ. The likelihood ratio yL(x) from these variables between signal and backgrounds is shown in Fig. 4, as well as the figure-of-merit. The classifier requirement on yL(x) is determined, by optimizing the figure-of-merit, to be 0.9979. The mass spectrum of KSπ after the above selections is shown in Fig. 3(b), and it is obvious that after likelihood requirements, the background level has been significantly suppressed. The selection efficiencies for signal and background processes with likelihood requirement are also shown in Table 1.

      Figure 4.  (color online) (a) Training and test sample distributions. The x-axis denotes yL and the y-axis, which is set at log scale, denotes the entries. The signal is the process τKSπντ. The backgrounds are e+eτ+τ with τ decays excluding the signal processes and the hadronic final states. (b) The figure-of-merit for the likelihood selection. MC samples are normalized to 1 ab1 integrated luminosity.

      The background events surviving the above selection criteria come from hadronic final states or e+eτ+τ decay products, where e+eτ+τ is mainly from the final processes with one additional π0, or multiple pion final states. At higher c.m. energies, the thrust T can be used to separate the hadronic final states from the signal process, defined as

      Tmax

      (10)

      where \overrightarrow{p}_{i} is the 3-momenta of the final state particles, \overrightarrow{n} is a 3-vector with unit norm. Finding \overrightarrow{n} in any direction of space makes the thrust T maximum. The typical distributions of thrust T at c.m. energy \sqrt{s} = 4.26 and 7.0 GeV are shown in Fig. 5 for both signal and hadronic processes. It is found that the thrust T can distinguish clearly between signal and hadronic background at high c.m. energies such as \sqrt{s} = 7.0 GeV. At lower c.m. energies, it is hard to distinguish the signal from background due to insufficient boost. However, with a good description of the hadronic final states from LUNDARLW, these backgrounds can be described well with MC simulation.

      Figure 5.  (color online) Distribution of thrust T (log scale in y-axis) at c.m.energy (a) \sqrt{s}=4.26 GeV and (b) \sqrt{s}=7.0 GeV. The solid red line shows the signal process and the blue shaded area shows hadronic processes. MC samples are normalized to 1 ab^{-1} integrated luminosity.

    V.   OPTIMIZATION OF DETECTOR RESPONSE
    • After the above selection criteria, the signal process is selected with an efficiency of 22.22%, where the main loss of efficiency comes from the effects of track selection, reconstruction of K_{S} , and particle identification of leptons. These effects correspond to the sub-detectors of the inner tracker, the PID system, and the MUD. By studying the signal-to-background ratios in the selection of the signal process with varying sub-detector responses, the requirements of detector design can be optimized accordingly. With the interfaces provided by fast simulation software, four kinds of detector response can be studied.

      a.Tracking efficiency The tracking efficiency in fast simulation is characterized by two dimensions: transverse momentum P_{T} and polar angle \cos\theta , which are correlated with the level of track bending and the hit positions of tracks in the tracker system. For low-momentum tracks ( P_{T}<0.2 GeV/c), it is difficult to reconstruct efficiently due to stronger electromagnetic multiple scattering, electric field leakage, energy loss, etc. However, with different techniques of inner tracker design at STCF, or with advanced track-finding algorithms, the reconstruction efficiency of low-momentum tracks can be improved.

      In this analysis, the efficiency is scaled in the fast simulation, with a ratio from 1 to 1.4 for low-momentum tracks, which indicates increasing the tracking efficiency by a factor of 10% to 40%. For high-momentum tracks or tracks with efficiency approximately 100%, the tracking efficiency remains the same. The figure-of-merit for each scale factor of low-momentum tracking efficiency is shown in Fig. 6, defined by \dfrac{S}{\sqrt{S+B}} , where S denotes the signal yields of \tau^{-}\to K_{S}\pi^{-}\nu_{\tau} , and B denotes the backgrounds, both normalized to 1 ab ^{-1} integrated luminosity. From Fig. 6, it is found that the efficiency can be significantly improved with an optimization factor of 1.1 to 1.2.

      Figure 6.  The figure-of-merit for optimization of the tracking efficiency for low-momentum tracks.

      In fast simulation, the default tracking efficiency is sampled from the performance of the helium-gas-based cylindrical main drift chamber [50]. There are several alternative options for the tracker technique, especially for the inner tracker, such as a silicon pixel detector [57] or micro-pattern gaseous detector [58]. From the current R&D activities, the tracking efficiency for low-momentum tracks is expected to fulfill the request from our analysis of these techniques at STCF.

      b. Momentum/spatial resolution The momentum and spatial resolution can also be optimized in fast simulation with proper scaling factors. A better resolution will improve the significance of signal with the improvement of resolution in the K_{S}\pi^{-} mass spectrum. The effects of momentum and position resolutions on K_{S} reconstruction efficiency are also studied, and it is found that the spatial resolution is more sensitive than the momentum resolution for K_{S} reconstruction. From the reconstruction method for a charged track, these two resolutions are related to the same sources, which are the position resolution of a single wire in the xy plane and z direction, and multiple scattering. Micro-pattern gaseous detectors can achieve excellent position accuracy [58]. Besides, a MC study on the silicon pixel detector shows that the momentum resolution and tracking efficiency can be significantly improved due to its good spatial resolution [59].

      c. \pi/\mu separation As discussed before, there are two efficiency curves provided for \mu identification. The default muon identification in the fast simulation has a large \pi/\mu mis-identification rate, as shown by the black points in Fig. 2. As the fast simulation provides a function for optimizing the \pi/\mu mis-identification, we can vary the mis-identification rate for \pi/\mu based on the provided efficiency curves. The signal-to-background ratio can be significantly improved with a better \pi/\mu mid-identification rate. The \pi/\mu separation power can also be significantly improved with a hybrid design for MUD [50], and the efficiency curves are shown by the red points in Fig. 2 for the requirement of \pi/\mu mis-identification to be 3%. Three mis-identification rates are provided for the MUD, with different efficiencies for \mu . Table 2 summarizes the selection efficiency for these three efficiency responses. The detector efficiency for the signal process is highest for the \pi/\mu mis-identification at 3%, because a moderate selection on the performance of the MUD is applied with this requirement. Moreover, at STCF, the Ring Imaging CHerenkov detector (RICH) can be used to identify muons from pions at low momentum.

      FR(%)Eff. \mu(%)Eff. bkg(%)Final eff.(%)
      376.80\pm0.020.38532.82\pm0.02
      1.766.79\pm0.020.37531.40\pm0.02
      158.63\pm0.030.37030.20\pm0.02

      Table 2.  Efficiency of selecting \mu and final selection efficiency using a hybrid design for MUD at STCF with different \pi/\mu fake rate (FR) requirements.

      d. Position/energy resolution for photons. In this analysis, processes including \pi^{0} are suppressed by vetoing events with a reconstructed \pi^{0} . The position and energy resolutions for neutral tracks are optimized with the fast simulation. It is found that the current resolutions provided in the fast simulation, 6 mm for position resolution and 2.5% for energy resolution of a 1 GeV photon, can satisfy the physical requirements for this analysis.

      From the above study, a set of optimization factors for sub-detector responses for this analysis is provided: an improvement of tracking efficiency by a factor of 10\% \sim 20\% for low-momentum tracks; a good momentum/spatial resolution for charged tracks; a good muon identification with the mis-identification rate of \pi/\mu to be less than 3%; and a good position/energy resolution for photons. With these optimization factors applied, the selection efficiency for signal process is improved from 22.22% to 32.82%. The invariant mass of K_{S}\pi^{-} after optimization is shown in Fig. 3(c). The background level is significantly suppressed, where the dominant \tau pair decay background events are \tau^{-}\to\pi^{+}\pi^{-}\pi^{-}\nu_{\tau} and \tau^{-}\to K_{S}\pi^{-}\pi^{0}\nu_{\tau} . The optimized selection efficiency and the background levels are summarized in Table 1. The detailed selection efficiencies for each step are listed in Table 3, where the improvement of the selection efficiency compared to the MC without optimization is given.

      Selection\varepsilon (%)\varepsilon_{\rm{rela}} (%)\Delta\varepsilon_{\rm{rela}} (%)
      N_{\rm{charged\_track}}69.75\pm0.0269.757.52
      select K_{S}54.34\pm0.0377.9013.10
      lepton ID37.89\pm0.0387.3119.13
      \pi ID36.64\pm0.0396.702.66
      veto \pi^{0}36.64\pm0.0399.99
      K_{S} flight length33.45\pm0.0291.380.27
      likelihood method32.82\pm0.0298.100.52
      Total32.82\pm0.02- -

      Table 3.  Cut flow for each selection criterion for the signal process, where \varepsilon is the overall efficiency, \varepsilon_{\rm{rela}} is the relative efficiency for each criterion, and \Delta\varepsilon_{\rm{rela}} is the relative improvement of optimized detector response compared to the original.

    VI.   STATISTICAL ANALYSIS
    • With the above selection criteria and optimization procedure, the numbers of signal events from \tau^{-} and \tau^{+} decays with generic MC normalized to 1 ab ^{-1} integrated luminosity are obtained by fitting the K_{S}\pi invariant mass with the RooFit tool, where the signal can be parameterized by the function shown in Eq. (4) and the background is described with simulated background. The efficiency-corrected numbers for \tau^{-}\to K_{S}\pi^{-}\nu_{\tau} and \tau^{+}\to K_{S}\pi^{+}\bar{\nu}_{\tau} are 3681017\pm5034 and 3681127\pm5091 , respectively. This shows a good consistency with the input values. The statistical sensitivity of CPV with decay rate can be calculated, using Eq. (1), to be 9.7\times10^{-4} . Since the cross-section of e^{+}e^{-}\to\tau^{+}\tau^{-} is around 3.5 nb in the energy region from \sqrt{s} = 4.0 to 5.0 GeV, we studied the selection efficiency for the signal process in this energy region, as shown in Fig. 7. The efficiencies in this energy region vary from 32.82% to 32.03%, with a relative difference of less than 3%. Therefore, the statistics of the signal process increase linearly with more data collected. As the event selection is not background-free, the sensitivity of CPV is proportional to 1/\sqrt{L} , as shown in Fig. 7(b), where different sizes of MC samples, from 0.1 ab ^{-1} to 1.0 ab ^{-1} are applied for the study. Finally, we can conclude that, with 10 ab ^{-1} integrated luminosity collected at STCF from \sqrt{s} = 4.0 to 5.0 GeV, the sensitivity of CPV for the process \tau^{-}\to K_{S}\pi^{-}\nu_{\tau} is 3.1\times10^{-4} .

      Figure 7.  (a) Selection efficiencies for the signal process at different c.m. energies from \sqrt{s}=4.0 GeV to 5.0 GeV. (b) CPV sensitivity of \tau^{-}\to K_{S}\pi^{-}\nu_{\tau} for varying integrated luminosity L, following 1/\sqrt{L} extrapolation.

      Though this analysis is performed purely with MC simulation, we need to discuss possible systematic uncertainties which may bias the decay-rate asymmetry measurement between \tau^{-}\to K_{S}\pi^{-}\nu_{\tau} and \tau^{+}\to K_{S}\pi^{+}\bar{\nu}_{\tau} . The sources of uncertainty are estimated as follows. i) The detection asymmetry for charged particles, which can be studied from a control sample of \tau^{-}\to \pi^{+}\pi^{-}\pi^{-}\nu_{\tau} . The difference from the detector can be corrected by comparing the asymmetry between MC simulation and data, and the remaining uncertainty is related to the statistical uncertainty of the control sample. Since the branching fraction of the control sample is over one order of magnitude larger than the signal process, the uncertainty will be significantly smaller than the statistical uncertainty of \tau^{-}\to K_{S}\pi^{-}\nu_{\tau} . ii) Uncertainty related to the selection criteria, which can be studied by varying the selection criteria. By performing the Barlow test [60], any bias from selection will be studied and corrected until the systematic accuracy matches the statistical precision. iii) Uncertainty from the MC generator. In this analysis, no CPV is modeled and the output asymmetry in MC is about 3\times10^{-5} , which can be neglected. iv) Uncertainty from background contamination. In this analysis, there is no decay-rate asymmetry observed in the background, and the uncertainty can be tested by applying different criteria to the background between data and MC simulation. v) The different nuclear-interaction cross sections of the K^{0} and \bar{K}^{0} meson. In Ref. [28], a correction to the asymmetry accounting from K^{0} and \bar{K}^{0} mesons with material in the detector was applied by taking the K^{\pm} nucleon cross-sections as an analogy under the assumption of isospin invariance [61]. We take an uncertainty of 0.01% for this uncertainty, similar to that from Ref. [28], ignoring the material difference at STCF. As a consequence of the above discussion, the systematic uncertainty for a decay-rate asymmetry study of \tau^{-}\to K_{S}\pi^{-}\nu_{\tau} will be at the same level as its statistical uncertainty.

      As introduced before, the sensitivity of CPV in the processes of \tau^{-}\to K_{S}\pi^{-}\nu_{\tau} at the Belle II experiment is expected to be 1/\sqrt{70} of the current precision of the Belle experiment with a 50 ab ^{-1} data sample collected [43], which is better than the sensitivity obtained at STCF with one year of data taking, but comparable with that at STCF with ten years of data taking. The study of \tau^{-}\to K_{S}\pi^{-}\nu_{\tau} at a B-factory benefits from a larger luminosity, and a good separation of the two hemispheres from \tau^{+} and \tau^{-} decay. However, at STCF, the cross section of e^{+}e^{-}\to\tau^{+}\tau^{-} is the largest and the selection efficiency for the signal process is better with high detection/PID efficiency for low-momentum tracks. When the c.m. energy goes up, e.g. at \sqrt{s} = 7 GeV, the two hemispheres from \tau^{+} and \tau^{-} decays can also be separated well, as shown in Fig. 5. Therefore, it is of great potential to study the CPV and relevant physics program of the process \tau^{-}\to K_{S}\pi^{-}\nu_{\tau} at STCF.

    VII.   SUMMARY AND PROSPECTS
    • In this paper, the sensitivity of the decay-rate asymmetry in \tau^{-}\to K_{S}\pi^{-}\nu_{\tau} decays is studied at \sqrt{s} = 4.26 GeV with 1 ab ^{-1} inclusive MC. Benefitting from the largest cross-section of 3.5 nb for e^{+}e^{-}\to\tau^{+}\tau^{-} at \sqrt{s} = 4.26 GeV, and the optimized efficiency with the fast simulation, the statistical sensitivity of CPV is 9.7\times10^{-4} , which is 2.3 times better than that of BaBar [28]. With 10 ab ^{-1} luminosity collected at STCF, the sensitivity is expected to be at a level of 3.1\times10^{-4} , which is comparable to the uncertainty of theoretical prediction in the SM. A new round of measurements of A_{CP} defined by the decay rate difference is important and desirable, c.f. Eq. (2), and will certainly shed light on the existence of the NP signal. Several theoretical models discussed in Sec. II will be reexamined. CPV in the exchange of charged Higgs bosons has also been explored in experiments, and is defined by the difference between angular observables, as a helpful complement. With this precision of 10^{-4} , re-measurement of this channel will be pursued. Meanwhile, the other prospects mentioned (form factors, |V_{us}| extraction etc.) could be examined.

    ACKNOWLEDGEMENTS
    • The authors thank the USTC Supercomputing Center and the Hefei Comprehensive National Science Center for their strong support.

Reference (61)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return