Feasibility study of measuring $b\to s\gamma$ photon polarisation in $D^0\rightarrow K_1(1270)^- e^+\nu_e$ at STCF

We report a sensitive study of measuring $b\to s\gamma$ photon polarisation in $D^{0}\to K_1(1270)^-e^+\nu_e$ with an integrated luminosity of $\mathcal L$ = 1 ab$^{-1}$ at a center-of-mass energy of 3.773 GeV at future Super Tau Charm Facility. More than 61,000 signals of $D^{0}\to K_1(1270)^-e^+\nu_e$ are expected. Based on a fast simulation software package, the statistical sensitivity for the ratio of up-down asymmetry is estimated to be $1.5\times 10^{-2}$ by performing a two-dimensional angular analysis in $D^{0}\to K_1(1270)^-e^+\nu_e$. Combining with measurements of up-down asymmetry in $B\to K_1\gamma$, the photon polarisation in $b\to s\gamma$ can be determined model-independently.


Introduction
The new physics (NP) and related phenomena beyond the Standard Model (SM) could be explored by indirect searches in b → sγ processes. The photon emitted from the electroweak penguin loop in b → sγ transitions is predominantly polarised in SM. New sources of chirality breaking can modify the b → sγ transition strongly as suggested in several theories beyond the SM [1][2][3]. A representative example is the leftright symmetric model (LRSM) [4,5], in which the photon can acquire a significant right-handed component. An observation of right-handed photon helicity would be a clear indication for NP [4].
The effective Hamiltonian of b → sγ is where C 7L and C 7R are the Wilson coefficients for left-and right-handed photons, respectively. In the SM, the chiral structure of W ± couplings to quarks leads to a dominant polarisation photon and a suppressed right-handed configuration, a e-mail: E-mail: yulanf@whu.edu.cn b e-mail: E-mail: sunl@whu.edu.cn (corresponding author) c e-mail: E-mail: zxrong@ustc.edu.cn (corresponding author) and the photon from radiativeB (B) decays is predominantly left-(right-) handed, i.e., |C 7L | 2 |C 7R | 2 (|C 7L | 2 |C 7R | 2 ). Various methods have been proposed to determine the photon polarisation of the b → sγ. The first method [6] suggests that the CP asymmetries which depend on the photon helicity could be measured in time dependent asymmetry in the charged and neutral B(t) →X CP s/d γ decays. The second method [7] to determine the photon polarisation is based on b → s l + l − transition where the dilepton pair originates from a virtual photon. The third method in Λ b → Λ γ also could be used to measure the photon polarisation directly. The forward-backward asymmetry defined in [8,9] is proportional to the photon polarisation.
Measuring the photon polarisation in radiative B decays into K resonance states, K res (→ Kππ), is proposed in [10,11]. The photon polarisation parameter λ γ could be described by an up-down asymmetry (A UD ) of the photon momentum relative to the Kππ decay plane in K res rest frame. The photon polarisation in B → K res γ is given in terms of Wilson coefficients [11]: with λ γ −1 for b → sγ and λ γ +1 forb →sγ. The integrated up-down asymmetry which is proportional to photon polarisation parameter λ γ for the radiative process proceeding through a single resonance K res is defined [10,11] where θ K is defined as the relative angle between the normal direction n of the K res decay plane and the opposite flight di-arXiv:2107.06118v3 [hep-ex] 3 Nov 2021 rection of the photon in the K res rest frame, and J denotes the K res → Kππ decay amplitude [10]. In charm sector, radiative D 0 − decays into CP eigenstate are expected to determine the photon polarization by means of the charm meson's finite width difference [12]. Recently, LHCb collaboration reported the direct observation of the photon polarisation with a significance of 5.2σ in B + → K + π − π + γ decay [13]. In the Kππ mass interval [1.1,1.3] GeV/c 2 which is dominated by K 1 (1270), A UD is extracted to be (6.9 ± 1.7) ×10 −2 . However, the currently limited knowledge of the structure of the Kππ mass spectrum, which includes interfering kaon resonances, prevents the translation of a measured asymmetry into an actual value for λ γ .
To solve this dilemma, three methods are proposed [14][15][16]. In Ref. [14], by using B → J/ψK 1 → J/ψKππ channel, the hadronic information of Kππ can be determined. Along the lines of the method known to B → K 1 (→ Kππ)γ decays, the extraction of photon polarization in D ( s) → K 1 (→ Kππ)γ decays is introduced in the Kππ system [15]. A novel method is proposed in Ref. [16] to determine the photon helicity in b → sγ by combining the B → K 1 γ and semileptonic decay D → K 1 l + ν l (l = µ + , e + ) model-independently. A ratio of up-down asymmetries is introduced in [16]. Ref. [16] introduces two angles θ K and θ l in D 0 → K − 1 e + ν e shown in Figure 1 and the ratio of up-down asymmetry A UD is defined as [16] Here the definition of the normal direction of K 1 decay Fig. 1 The kinematics for D 0 → K − 1 (K − π + π − )e + ν e . The relative angle between the normal direction of K − 1 decay plane and the opposite of the D 0 flight direction in the K − 1 rest frame is denoted as θ K , where the normal direction of K − 1 decay plane is defined as p π,slow × p π,fast in which p π,slow and p π,fast corresponding to the momenta of the lower and higher momentum pions, respectively. The θ l is introduced as the relative angle between the flight direction of e + in the e + ν e rest frame and e + ν e in the D 0 rest frame. [16] plane is the same as in B → K res γ in LHCb [13].
Then photon helicity parameter of b → sγ could be extracted by [16] So the photon polarisation of b → sγ could be determined model-independently by the combination of A UD in D 0 → K − 1 (→ K − π + π − )e + ν e and A UD in B + → K + 1 (→ K + π − π + )γ. Experimentally, the semileptonic decay of D 0 → K 1 (1270) − e + ν e has been observed for the first time with a statistical significance greater than 10σ by using 2.93 fb −1 of e + e − collision data at √ s = 3.773 GeV by BESIII [17]. About 109 signals are observed, and the measured branching fraction is where the first and second uncertainties are the statistical and systematic uncertainties, respectively, and the third uncertainty is the external uncertainty from the assumed branching fractions (BFs) of K 1 subdecays. Still, the statistics of the current BESIII data set are insufficient to measure the ratio of up-down asymmetry in D 0 → K 1 (1270) − e + ν e . A much larger data sample with similarly low background level is urgently needed for performing the angular analysis in D 0 → K 1 (1270) − e + ν e , which calls for the construction of a next generation e + e − collider operating at the τ-charm energy region with much higher luminosity.
The Super Tau Charm Facility (STCF) is a scientific project proposed in China for high energy physics frontier [18]. The STCF plans to produce charmed hadron pairs near the charm threshold which allow for exclusive reconstruction of their decay products with well-determined kinematics. Such samples at the threshold allow for a double-tag technique [19] to be employed where the full events can be reconstructed and provide a unique environment to measure A UD in D 0 → K 1 (1270) − e + ν e with very low background level.
In this work, we present a feasibility study of a ratio of up-down asymmetry in D 0 → K 1 (1270) − e + ν e at STCF. Throughout this paper, charged conjugated modes are always implied. This paper is organised as follows: In Sect. 2, detector concept for STCF is introduced as well as the Monte Carlo (MC) samples used in this feasibility study. In Sect. 3, the event selection and analysis method are described. The optimisation of detector response is elaborated in Sect. 4 and the results are presented in Sect. 5. Finally, we conclude in Sect. 6.

Detector and MC simulation
The proposed STCF is a symmetric electron-positron beam collider designed to provide e + e − interactions at a centerof-mass (c.m.) energy √ s from 2.0 to 7.0 GeV. The peaking luminosity is expected to be 0.5 × 10 35 cm −2 s −1 at √ s = 4.0 GeV, and the integrated luminosity per year is 1 ab −1 . Such an environment will be an important low-background playground to test the SM and probe possible new physics beyond the SM. The STCF detector is a general purpose detector designed for e + e − collider which includes a tracking system composed of the inner and outer trackers, a particle identification (PID) system with 3σ charged K/π separation up to 2 GeV/c, and an electromagnetic calorimeter (EMC) with an excellent energy resolution and a 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, the expected detection efficiency and resolution can be found in [18,20,21].
Currently, the STCF detector and the corresponding offline software system are under active development. A reliable fast simulation tool for STCF has been developed [21], which takes the most common event generators as input to perform a fast and realistic simulation. The simulation includes resolution and efficiency responses for tracking of final state particles, PID system and kinematic fit related variables. Besides, the fast simulation also provide some functions for adjusting performance of each sub-system which can be used to optimise the detector design according to physical requirement.
This study uses MC simulated samples corresponding to 1 ab −1 of integrated luminosity at √ s = 3.773 GeV. The simulation includes the beam-energy spread and initial-state radiation (ISR) in the e + e − annihilations modeled with the generator KKMC [22]. The inclusive MC samples consist of the production of the DD pairs, the non-DD decays of the ψ(3770), the ISR production of the J/ψ and ψ(3686) states, and the continuum process incorporated in KKMC [22]. The known decay modes are modeled with EVTGEN [23] using BFs taken from the Particle Data Group [24], and the remaining unknown decays from the charmonium states with LUNDCHARM [25]. Final-state radiation (FSR) from charged final-state particles is incorporated with the PHOTOS package [26].
Included in the inclusive DD MC sample, the D 0 → K 1 (1270) − e + ν e decay is generated with the ISGW2 model [27] with BF comes from Ref. [24] and K 1 (1270) − meson is allowed to decay into all intermediate processes that result in a K − π + π − final state. The resonance shape of the K 1 (1270) − meson is parameterised by a relativistic Breit-Wigner function. The mass and width of K 1 (1270) − meson are fixed at the known values as shown in Table 1, and the BFs of K 1 (1270) subdecays measured by Belle [28] are input to generate the signal MC events, since they give better consistency [17] between data and MC simulation than those reported in [24]. Table 1 Mass, width [24] and ratios of subdecays of K 1 (1270) − (Fit2) [28] used in this analysis.

Event section and anasysis
The feasibility study employs the e + e − → ψ(3770) → D 0D0 decay chain. TheD 0 mesons are reconstructed by three channels with low background level,D 0 → K + π − , K + π − π 0 and K + π − π + π − . These inclusively selected events are referred to as single-tag (ST)D 0 mesons. In the presence of the ST D 0 mesons, candidates for D 0 → K 1 (1270) − e + ν e are selected to form double-tag (DT) events. Each charged track is required to satisfy the vertex requirement and detector acceptance in fast simulation. The combined confidence levels under the positron, pion and kaon hypotheses (CL e , CL π and CL K , respectively) are calculated. Kaon (pion) candidates are required to satisfy CL K > CL π (CL π > CL K ). Positron candidates are required to satisfy CL e / (CL e + CL K + CL π ) > 0.8. To reduce the background from hadrons and muons, the positron candidate is further required to have a deposit energy in the EMC greater than 0.8 times its momentum in the MDC. The π 0 meson is reconstructed via π 0 → γγ decay. The γγ combination with an invariant mass in the range (0.115, 0.150) GeV/c 2 are regarded as a π 0 candidates, and a kinematic fit by constraining the γγ invariant mass to the π 0 nominal mass [24] is performed to improve the mass resolution. The STD 0 mesons are identified by the energy difference ∆ E ≡ ED0 − E beam and the beam-constrained mass M BC ≡ E 2 beam − | pD0| 2 , where E beam is the beam energy, and ED0 and pD0 are the total energy and momentum of the STD 0 in the e + e − rest frame. If there are multiple combinations in an event, the combination with the smallest ∆ E is chosen for each tag mode. The combinatorial backgrounds in the M BC distributions are suppressed by requiring ∆ E within (-29, 27), (-69, 38) and (-31, 28) MeV forD 0 → K + π − , K + π − π 0 and K + π − π + π − , respectively, which correspond to about 3.5σ away from the fitted peak. Particles recoiling against the STD 0 mesons candidates are used to reconstruct candidates for D 0 → K 1 (1270) − e + ν e decay, where the K 1 (1270) − meson is reconstructed using its dominant decay K 1 (1270) − → K − π + π − . It is required that there are only four good unused charged tracks available for this selection. The charge of the lepton candidate is required to be the same as that of the charged kaon of the tag side. The other three charged tracks are identified as a kaon and two pions, based on the same PID criteria used in the ST selection. The kaon candidate must have charge opposite to that of the positron.
The main peaking background comes from misidentifying a pion to a positron, and additional criteria as in [17] are used to improve the π/e separation.
Information concerning the undetectable neutrino inferred by the kinematic quantity M 2 miss ≡ E 2 miss − | p miss | 2 , where E miss and p miss are the missing energy and momentum of the signal candidate, respectively, calculated by E miss ≡ E beam − ∑ j E j and p miss ≡ − pD0 − ∑ j p j in the e + e − center-of-mass frame. The index j sums over the K − , π + , π − and e + of the signal candidate, and E j and p j are the energy and momentum of the j-th particle, respectively. To partially recover the energy lost to the FSR and bremsstrahlung, the fourmomenta of photon(s) within 5 • of the initial positron direction are added to the positron four-momentum. Figure 2 shows the distribution of M K − π + π − vs. M 2 miss of the accepted D 0 → K − π + π − e + ν e candidate events in the MC sample after combining all tag modes. A clear signal, which concentrates around the K 1 (1270) − nominal mass in the M K − π + π − distribution and around zero in the M 2 miss distribution, can be seen. The selection efficiency of signal candidates with the ST modesD 0 → K + π − , K + π − π 0 and K + π − π − π + are 12.11%, 6.93% and 6.25%, respectively. In order to determine the angular distributions of cos θ K and cos θ l , a two-dimensional (2-D) fit to M 2 miss and M K − π + π − is performed to extract the signal yield in each angle bin. The 2-D fit projections to the M 2 miss and M K − π + π − distributions are shown in Figure 3. In the fit, the 2-D signal shape is described by the MC-simulated shape extracted from the signal MC events while the 2-D background shape is modeled by those derived from the inclusive MC sample. The smooth 2-D probability density functions of signal and background are modeled by using RooNDKeysPdf [30,31]. The reconstructed efficiencies of signal candidates in each cos θ K and cos θ l interval are shown in Figure 4. The signal reconstruction efficiency shows a clear trend of increasing monotonically with cos θ l , which is due to strong correlation between cos θ l and electron momentum, and D 0 candidates with lower momentum electrons are less likely to satisfy electron tracking and PID requirements.
The signal yields in each cos θ K and cos θ l interval corrected by the signal reconstruction efficiency are fitted with a polynomial function [16] f (cos θ K , cos θ l ;A UD , d + , d − ) = where the d ± are the angular coefficients, defined as: Fig. 4 The signal reconstruction efficiencies (in percentage) in bins of cos θ K and cos θ l . In each bin j, the signal reconstruction efficiency is where the B i ST denotes the known BF of each tag mode i, the ε i j DT represents the DT efficiency in each bin j for tag mode i.
The coefficients c ± and c 0 correspond to the nonperturbative amplitudes for D decays into K 1 with transverse and longitudinal polarisations, respectively. The ratio of up-down asymmetry A UD can be extracted directly. Besides, the fraction of longitudinal polarisation |c 0 | 2 |c 0 | 2 +|c + | 2 +|c − | 2 can be derived from the fitted d ± values. Form factor calculations based on different approaches such as covariant light-front quark model(LFQM) and light-cone QCD sum rules(LCSR) obtain the different results significantly. [32,33].
The possibility of some events migrating from a bin to its neighbor caused by the detection resolution is considered by calculating the full width at half maximum(FWHM) of cosθ l and cosθ K , respectively. The value of FWHM is 0.115 and 0.05, which indicates that the bin migration effects can be ignored, due to the larger bin width of 0.5.
A 2-D χ 2 fit to the cos θ l and cos θ K distributions allows to extract A UD , and the fit projections are shown in Figure 5. The statistical sensitivity of A UD based on 1 ab −1 MC sample is thus determined to be in the order of 1.8×10 −2 . As a cross-check, A UD is also determined with a counting method according to Eq. 4 and the corresponding result is compat- ible with the angular fit method. However, the angular fit method yields a more precise result on A UD and is taken as the nominal result.

Optimization of detector response
The main loss of the signal efficiency comes from the effects of charged tracking selection, neutral selection and identification of electron at low momentum. These effects correspond to the sub-detectors of the tracking system, the EMC and the PID system. By studying the DT efficiencies or signalto-background ratios for this process with variation of the sub-detector's responses, the requirement of detector design can be optimised accordingly. With the help of fast simulation software package, three kinds of detector responses are studied as introduced below: a.Tracking efficiency The tracking efficiency in fast simulation is characterised by two dimensions: transverse momentum P T and polar angle cosθ , which are correlated with the level of track bending and 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 technique in the tracking system design at STCF, or with advanced track finding algorithm, the efficiency is expected to be improved for lowmomentum tracks.
Benefiting from the flexible approach to change the response of charged track, the detection efficiency is scaled with a factor from 1.1 to 1.5 in the fast simulation. The figure-of-merit, defined by DT efficiency for characterising the performance of tracking efficiency is shown in Figure 6(a). From Figure 6(a), it is found that the DT efficiency can be significantly improved with the given scale factors. The resolution of momentum and position can also be optimised in fast simulation with proper functions. Insignificant improvement is found among optimisation of absolute σ xy from 30 µm to 150 µm and absolute σ z from 500 µm to 2500 µm yet, where σ xy and σ z are the resolution of the tracking system in the xy plane and z direction, respectively. This can be understood since the main source that affects the momentum resolution comes from electromagnetic multiple scattering on the material in the detector, instead of the position resolution. Therefore, material with low atomic number Z is required in the tracking system.
b.Detection efficiency for photon In this analysis, π 0 s are selected as part of the tag modeD 0 → K + π − π 0 and the π 0 selection also helps to suppress the main background of D 0 → K − π + π − π + π 0 in signal side. The figure-of-merit, defined by S √ S+B , to characterise the effect of photon detection efficiency on signal significance, in which S denotes the expected signal yield of D 0 → K 1 (1270) − (→ K − π + π − )e + ν e versus the scale factor of photon detection efficiency scanned from 1.1 to 1.5 is shown in Figure 6(b). c.π/e identification Misidentification from a pion to electron in the momentum smaller than 0.6 GeV/c forms the main peaking background D 0 → K − π + π − π + and D 0 → K − π + π − π + π 0 . As the fast simulation provides the function for optimising the π/e identification which allows to vary the misidentification rate for π/e, the π/e misidentification rate at 0.2 GeV/c is scanned from 5.7% to 0.64%, shown in Figure 6(c). And S √ S+B defined before is used to characterise the effect of misidentification of π + to e + on the signal significance.
In summary, three sets of optimization factors for different sub-detector responses are calculated separately: compared with our fast simulation with default settings, the DT efficiency is improved by ∼27% if the reconstructed efficiency for charged track is scaled by the factor of 1.1, and the value of S √ S+B is improved by 4% if the photon detection efficiency is scaled by a factor of 1.1, or 7% if the π/e misidentification rate is lowered by half to 3.2%, as reasonable assumptions in real case scenarios. With the above three factors applied altogether, the DT efficiency is improved by a factor of 33%. The corresponding angular 2-D χ 2 fit based on updated efficiency-corrected signal yields in different angular bins is performed. From the fit, the statistical uncertainty of the ratio of up-down asymmetry is extracted to be 1.5×10 −2 , that is, improved by 17% compared with the no optimisation scenario.

Statistical analysis
With the above selection criteria and optimisation procedure, the 2-D simultaneous fit to M 2 miss vs. M Kππ in the different interval of cos θ K vs. cos θ l is performed, the semileptonic decay signal yields produced are used for fitting the angular distribution. Therefore, the sensitivity of the ratio of up-down asymmetry in D 0 → K 1 (1270) − e + ν e with an integrated luminosity of 1 ab −1 is extracted as 1.5×10 −2 .
Besides, the selection efficiency for this process at √ s = 3.773 GeV where the cross section for e + e − → D 0D0 = 3.6 nb [34] is studied by a large MC sample, with a negligible error.
Eq. (2) indicates that the Wilson coefficients can be constrained by measuring the uncertainty of the photon polarisation parameter λ γ . Combining the uncertainty of A UD measurement [13] with the uncertainty of A UD measurement in this analysis, the sensitivity of λ γ can be determined using Eq. (5). Thus, the Wilson coefficients can be translated by the sensitivity of A UD . Figure 7 depicts the dependency of Wilson coefficients on ratio of A UD , using the A UD measured in the Kππ mass range of (1.1,1.3) GeV/c 2 [13] as the input, shown in the blue solid line. Considering the uncertainty of A UD , the corresponding constraints shown in the green parts. The photon polarisation parameter λ γ is predicted to be λ γ -1 forb →sγ in SM, which translated to A UD (9.2 ± 2.3)×10 −2 shown in the red and black solid line in Figure 7. For the systematic uncertainty on A UD , possible sources include the electron tracking and PID efficiencies as functions of electron momentum which cannot cancel out in the cos θ l distribution due to strong correlation between cos θ l and electron momentum as mentioned before. With the current binning scheme as shown in Fig. 4, the possibility of some events migrating from an angular bin to its neighbor because of the detector resolution effects on cos θ K and cos θ l is expected to be small and the related systematic uncertainty should be manageable. Moreover, as in the BESIII analysis [17], the signal and background shape modeling would affect the signal yields considerably in different angular bins, due to imprecise knowledge on the K 1 (1270) line shape, and background events such as D 0 → K − π + π − π + π 0 .
Our simulation does not include non-K 1 (1270) − sources of K − π + π − in the D 0 → K − π + π − e + ν e decay, which are estimated to be at least one order of magnitude lower than our signal decay of K 1 (1270) − [17]. We expect the systematic effect of the non-K 1 (1270) − sources on A UD to be small, although detailed studies on the K 1 (1400) − contribution are needed when more data become available.

Summary and prospect
In this work, the statistical sensitivity of a ratio of up-down asymmetry D 0 → K 1 (1270) − e + ν e with an integrated luminosity of L = 1 ab −1 at √ s = 3.773 GeV and the optimised efficiency with the fast simulation, is determined to be 1.5×10 −2 by performing an angular analysis. The hadronic effects in K 1 → Kππ can be quantified by A UD , therefore, combined with the measured up-down asymmetry A UD in B + → K + 1 (→ K + π − π + )γ [13], the photon polarisation in b → sγ can be measured to probe the new physics.