Improved measurement of the branching fractions of the inclusive decays $D^+ \to K_S^0X $ and $D^0 \to K_S^0X $

TL;DR

This paper reports improved measurements of the inclusive decay branching fractions of D+ and D0 mesons to K_S^0X using BESIII data, achieving higher precision than previous results.

Contribution

The study provides more precise measurements of D+ and D0 inclusive decay branching fractions to K_S^0X, refining previous world average values.

Findings

Branching fraction for D+ → K_S^0X is (32.78 ± 0.13 ± 0.27)%

Branching fraction for D0 → K_S^0X is (20.54 ± 0.12 ± 0.18)%

Results are consistent with previous measurements but with improved accuracy.

Abstract

By analyzing 2.93 fb$^{-1}$ of $e^+e^-$ collision data taken at the center-of-mass energy of 3.773 GeV with the BESIII detector, the branching fractions of the inclusive decays $D^+\to K^0_S X$ and $D^0\to K^0_S X$ are measured to be $(32.78\pm 0.13\pm 0.27)\%$ and $(20.54\pm 0.12\pm 0.18)\%$, respectively, where the first uncertainties are statistical and the second are systematic. These results are consistent with the world averages of previous measurements, but with improved precision.

Tables

Table 1. Table 1: Measurements of S measured ± subscript superscript 𝑆 plus-or-minus measured S^{\pm}_{\rm measured} , M measured ± subscript superscript 𝑀 plus-or-minus measured M^{\pm}_{\rm measured} , f CP + subscript 𝑓 limit-from CP f_{\rm CP+} and the correction factor for D 0 → K S 0 ​ X → superscript 𝐷 0 subscript superscript 𝐾 0 𝑆 𝑋 D^{0}\to K^{0}_{S}X due to quantum correlation effect.

CP tag mode	${\overline{D}}^0\to K^{+}{K}^{-}$	${\overline{D}}^0\to K_S^0{\pi }^0$
$S_{\mathrm{measured}}^{\pm }$	$57779\pm 287$	$70512\pm 311$
$M_{\mathrm{measured}}^{\pm }$	$4760\pm 81$	$4068\pm 78$
$f_{\mathrm{CP}+}$	$0.413\pm 0.010$
Correction factor	$1.0204\pm 0.0024$

Table 2. Table 2: Numbers of single-tag D − ​ ( D ¯ 0 ) superscript 𝐷 superscript ¯ 𝐷 0 D^{-}\,(\bar{D}^{0}) mesons ( N tag subscript 𝑁 tag N_{\rm tag} ), fitted double-tag yields in the K S 0 subscript superscript 𝐾 0 𝑆 K^{0}_{S} signal and sideband regions ( N K S 0 ​ sig subscript 𝑁 subscript superscript 𝐾 0 𝑆 sig N_{K^{0}_{S}\rm sig} and N K S 0 ​ sid subscript 𝑁 subscript superscript 𝐾 0 𝑆 sid N_{K^{0}_{S}\rm sid} ), net numbers of signal events ( N DT net superscript subscript 𝑁 DT net N_{\rm DT}^{\rm net} ), efficiencies of finding D + ( 0 ) → K S 0 ​ X → superscript 𝐷 0 superscript subscript 𝐾 𝑆 0 𝑋 D^{+(0)}\to K_{S}^{0}X in the presence of the single-tag D − ​ ( D ¯ 0 ) superscript 𝐷 superscript ¯ 𝐷 0 D^{-}\,(\bar{D}^{0}) meson ( ϵ sig subscript italic-ϵ sig \epsilon_{\rm sig} ), and measured branching fractions of D + ( 0 ) → K S 0 ​ X → superscript 𝐷 0 superscript subscript 𝐾 𝑆 0 𝑋 D^{+(0)}\to K_{S}^{0}X ( ℬ sig subscript ℬ sig \mathcal{B}_{\rm sig} ). The efficiency values include the branching fraction of K S 0 → π + ​ π − → subscript superscript 𝐾 0 𝑆 superscript 𝜋 superscript 𝜋 K^{0}_{S}\to\pi^{+}\pi^{-} as well as the correction factors described in the text. The uncertainties are statistical only.

Signal mode	$N_{\mathrm{tag}}$	$N_{K_S^0\mathrm{sig}}$	$N_{K_S^0\mathrm{sb}}$	$N_{\mathrm{DT}}^{\mathrm{net}}$	${ϵ}_{\mathrm{sig}}$ (%)	${ℬ}_{\mathrm{sig}}$ (%)
$D^{+}\to K_S^{0}X$	$798935\pm 1011$	$102186\pm 365$	$13981\pm 140$	$95195\pm 372$	$36.35\pm 0.08$	$32.78\pm 0.13$
$D^0\to K_S^{0}X$	$529227\pm 761$	$42166\pm 215$	$8688\pm 97$	$37822\pm 220$	$34.81\pm 0.10$	$20.54\pm 0.12$

Table 3. Table 3: Relative systematic uncertainties (in %) in the branching fraction measurements. A “…” denotes that there is no systematic uncertainty.

Source	$D^{+}\to K_S^{0}X$	$D^0\to K_S^{0}X$
$N_{\mathrm{tag}}$	0.50	0.50
$K_S^0$ reconstruction	0.53	0.53
$K_S^0$ multiplicity	0.20	0.28
$K_S^0$ sideband	0.15	0.10
Quantum-correlation effect	…	0.24
MC statistics	0.25	0.28
Quoted branching fraction	0.07	0.07
Scale factor $f_{\mathrm{co}}$	0.02	0.05
Total	0.82	0.87

Table 4. Table 4: Comparisons of the branching fractions of D + ( 0 ) → K S 0 ​ X → superscript 𝐷 0 superscript subscript 𝐾 𝑆 0 𝑋 D^{+(0)}\to K_{S}^{0}X measured in this study and Mark-III, BES, and PDG. For the values from BES, Mark-III, and PDG, we take half of the branching fractions of D → K 0 / K ¯ 0 ​ X → 𝐷 superscript 𝐾 0 superscript ¯ 𝐾 0 𝑋 D\to K^{0}/\bar{K}^{0}X . ℬ exclusive sum subscript superscript ℬ sum exclusive {\mathcal{B}}^{\rm sum}_{\rm exclusive} denotes the total branching fractions summing over all known exclusive D + ( 0 ) superscript 𝐷 0 D^{+(0)} decays involving a K S 0 subscript superscript 𝐾 0 𝑆 K^{0}_{S} .

Decay mode	Mark-III (%) Mark3-KSX	BES (%) BESII-KSX	PDG (%) pdg2020	This study (%)	${ℬ}_{\mathrm{exclusive}}^{\mathrm{sum}}$ (%)
$D^{+}\to K_S^{0}X$	$30.60\pm 3.25\pm 2.15$	$30.25\pm 2.75\pm 1.65$	$30.5\pm 2.5$	$32.78\pm 0.13\pm 0.27$	$31.68\pm 0.32$
$D^0\to K_S^{0}X$	$22.75\pm 2.50\pm 1.60$	$23.80\pm 2.40\pm 1.50$	$23.5\pm 2.0$	$20.54\pm 0.12\pm 0.18$	$18.16\pm 0.72$

Table 5. Table 5: Initial and final states contributing to the inclusive decay D + → K S 0 ​ X → superscript 𝐷 superscript subscript 𝐾 𝑆 0 𝑋 D^{+}\to K_{S}^{0}X and the corresponding branching fractions.

Initial state	$ℬ$ (%)	Final state	$ℬ$ (%)
$K_S^0{\pi }^{+}$	$1.49\pm 0.03$ ref18	$K_S^0{\pi }^{+}$	$1.49\pm 0.03$
${\overline{K}}^0{e}^{+}{\nu }_e$	$8.73\pm 0.10$ pdg2020	$K_S^0{e}^{+}{\nu }_e$	$4.37\pm 0.05$
${\overline{K}}^{\ast 0}{e}^{+}{\nu }_e$	$5.40\pm 0.10$ pdg2020	$K_S^0{\pi }^0{e}^{+}{\nu }_e$	$0.89\pm 0.02$
${\overline{K}}^0{\mu }^{+}{\nu }_{\mu }$	$8.76\pm 0.19$ pdg2020	$K_S^0{\mu }^{+}{\nu }_{\mu }$	$4.380\pm 0.095$
${\overline{K}}^{\ast 0}{\mu }^{+}{\nu }_{\mu }$	$5.27\pm 0.15$ pdg2020	$K_S^0{\pi }^0{\mu }^{+}{\nu }_{\mu }$	$0.88\pm 0.02$
$K_S^0{K}_S^0{K}^{+}$	$0.254\pm 0.013$ pdg2020	$K_S^0{K}_S^0{K}^{+}$	$0.254\pm 0.013$
$K_S^0{K}_L^0{K}^{+}$	$0.508\pm 0.026$ pdg2020	$K_S^0{K}_L^0{K}^{+}$	$0.508\pm 0.026$
$K_S^0{K}_S^0{\pi }^{+}$	$0.270\pm 0.013$ pdg2020	$K_S^0{K}_S^0{\pi }^{+}$	$0.270\pm 0.013$
${\overline{K}}^0{\pi }^{+}{\pi }^{+}{\pi }^{-}$	$6.20\pm 0.18$ pdg2020	$K_S^0{\pi }^{+}{\pi }^{+}{\pi }^{-}$	$3.10\pm 0.09$
${\overline{K}}^0{K}^{+}{K}^{-}{\pi }^{+}$	$0.048\pm 0.010$ pdg2020	$K_S^0{K}^{+}{K}^{-}{\pi }^{+}$	$0.024\pm 0.005$
$K_S^0{\pi }^{+}{\pi }^0$	$7.36\pm 0.21$ pdg2020	$K_S^0{\pi }^{+}{\pi }^0$	$7.36\pm 0.21$
${\overline{K}}^{0}K_{}^{\ast }{}_{}{}^{+}$	$1.738\pm 0.182$ bes3-kskpi0	$K_S^{0}K_{}^{\ast }{}_{}{}^{+}$	$0.869\pm 0.091$
${\overline{K}}_1^0\left(1270\right)e^{+}{\nu }_e$	$0.231\pm 0.030$ bes3-Dp-K1ev	$K_S^0{\pi }^0{\pi }^0{e}^{+}{\nu }_e$	$0.0055\pm 0.0007$
${\overline{K}}_1^0\left(1270\right)e^{+}{\nu }_e$	$0.231\pm 0.030$ bes3-Dp-K1ev	$K_S^0{\pi }^{+}{\pi }^{-}{e}^{+}{\nu }_e$	$0.039\pm 0.005$
${\overline{K}}_1^0\left(1270\right){\mu }^{+}{\nu }_{\mu }$	$0.231\pm 0.030$ bes3-Dp-K1ev	$K_S^0{\pi }^{+}{\pi }^{-}{\mu }^{+}{\nu }_{\mu }$	$0.039\pm 0.005$
${\overline{K}}_1^0\left(1270\right){\mu }^{+}{\nu }_{\mu }$	$0.231\pm 0.030$ bes3-Dp-K1ev	$K_S^0{\pi }^0{\pi }^0{\mu }^{+}{\nu }_{\mu }$	$0.0055\pm 0.0007$
$K_S^0{\pi }^{+}\eta$	$1.309\pm 0.048$ bes3-etaX	$K_S^0{\pi }^{+}\eta$	$1.309\pm 0.048$
$K_S^0{K}^{+}\eta$	$0.0185\pm 0.0050$ bes3-etaX	$K_S^0{K}^{+}\eta$	$0.0185\pm 0.0050$
$K_S^0{\pi }^{+}{\pi }^0\eta$	$0.12\pm 0.03$ bes3-etaX	$K_S^0{\pi }^{+}{\pi }^0\eta$	$0.12\pm 0.03$
$K_S^0{K}^{+}{\pi }^{+}{\pi }^{-}$	$0.189\pm 0.013$ bes3-D-KKpipi	$K_S^0{K}^{+}{\pi }^{+}{\pi }^{-}$	$0.189\pm 0.013$
$K_S^0{K}^{-}{\pi }^{+}{\pi }^{+}$	$0.227\pm 0.013$ bes3-D-KKpipi	$K_S^0{K}^{-}{\pi }^{+}{\pi }^{+}$	$0.227\pm 0.013$
$K_S^0{K}^{+}{\pi }^0{\pi }^0$	$0.058\pm 0.013$ bes3-D-KKpipi	$K_S^0{K}^{+}{\pi }^0{\pi }^0$	$0.058\pm 0.013$
$K_S^0{K}_S^0{\pi }^{+}{\pi }^0$	$0.134\pm 0.021$ bes3-D-KKpipi	$K_S^0{K}_S^0{\pi }^{+}{\pi }^0$	$0.134\pm 0.021$
$K_S^0{K}_L^0{\pi }^{+}{\pi }^0$	$0.268\pm 0.042$ bes3-D-KKpipi	$K_S^0{K}_L^0{\pi }^{+}{\pi }^0$	$0.268\pm 0.042$
${\overline{K}}^0{K}^{+}$	$0.604\pm 0.024$ k0k	$K_S^0{K}^{+}$	$0.302\pm 0.012$
$K_S^0{K}^{+}{\pi }^0$	$0.507\pm 0.003$ k0k	$K_S^0{K}^{+}{\pi }^0$	$0.507\pm 0.003$
$\varphi {\pi }^{+}$	$0.570\pm 0.015$ phipi	$K_S^0{K}_L^0{\pi }^{+}$	$0.193\pm 0.006$
${\overline{K}}^0{\pi }^{+}{\eta }^{\prime }$	$0.380\pm 0.042$ k0pieta	$K_S^0{\pi }^{+}{\eta }^{\prime }$	$0.190\pm 0.021$
$K_S^0\omega {\pi }^{+}$	$0.707\pm 0.050$ ksompi	$K_S^0\omega {\pi }^{+}$	$0.707\pm 0.050$
$K_S^0{\pi }^{+}{\pi }^0{\pi }^0$	$2.2516\pm 0.1068$ kspipi0pi0	$K_S^0{\pi }^{+}{\pi }^0{\pi }^0$	$2.2516\pm 0.1068$
$K_S^0{\pi }^{+}{\pi }^{+}{\pi }^{-}{\pi }^0$	$0.4001\pm 0.0550$ kspipi0pi0	$K_S^0{\pi }^{+}{\pi }^{+}{\pi }^{-}{\pi }^0$	$0.4001\pm 0.0550$
$K_S^0{\pi }^{+}{\pi }^0{\pi }^0{\pi }^0$	$0.126\pm 0.050$ kspipi0pi0	$K_S^0{\pi }^{+}{\pi }^0{\pi }^0{\pi }^0$	$0.126\pm 0.050$
$\varphi {\pi }^{+}{\pi }^0$	$0.579\pm 0.055$ phipipi0	$K_S^0{K}_L^0{\pi }^{+}{\pi }^0$	$0.196\pm 0.019$
Sum	…	…	$31.68\pm 0.32$

Table 6. Table 6: Initial and final states contributing to the inclusive decay D 0 → K S 0 ​ X → superscript 𝐷 0 superscript subscript 𝐾 𝑆 0 𝑋 D^{0}\to K_{S}^{0}X and the corresponding branching fractions.

Initial state	$ℬ$ (%)	Final state	$ℬ$ (%)
$K_S^0\omega$	$1.11\pm 0.12$ pdg2020	$K_S^0\omega$	$1.11\pm 0.12$
$K_S^0{\eta }^{\prime }{\pi }^0$	$0.252\pm 0.027$ pdg2020	$K_S^0{\eta }^{\prime }{\pi }^0$	$0.252\pm 0.027$
$K^{-}{a}_0^{+}$	$0.059\pm 0.018$ pdg2020	$K_S^0{K}^{+}{K}^{-}$	$0.003\pm 0.001$
$K_S^0\varphi$	$0.444\pm 0.015$ pdg2020	$K_S^0\varphi$	$0.444\pm 0.015$
$\varphi K_L^0$	$0.4440\pm 0.0304$ pdg2020	$K_S^0{K}_L^0{K}_L^0$	$0.151\pm 0.010$
${\overline{K}}^{\ast 0}{\pi }^0$	$1.82\pm 0.07$ pdg2020	$K_S^0{\pi }^0{\pi }^0$	$0.303\pm 0.012$
$\varphi \gamma$	$0.00282\pm 0.00019$ pdg2020	$K_S^0{K}_L^0\gamma$	$0.0010\pm 0.0001$
$K_S^0{\pi }^{+}{\pi }^{-}$	$2.80\pm 0.18$ pdg2020	$K_S^0{\pi }^{+}{\pi }^{-}$	$2.80\pm 0.18$
$K_S^0{K}^{+}{\pi }^{-}$	$0.217\pm 0.034$ pdg2020	$K_S^0{K}^{+}{\pi }^{-}$	$0.217\pm 0.034$
$K_S^0{K}^{-}{\pi }^{+}$	$0.33\pm 0.05$ pdg2020	$K_S^0{K}^{-}{\pi }^{+}$	$0.33\pm 0.05$
$K_S^0{K}_S^0{K}^{-}{\pi }^{+}$	$0.059\pm 0.013$ pdg2020	$K_S^0{K}_S^0{K}^{-}{\pi }^{+}$	$0.059\pm 0.013$
$K_S^0{K}_L^0{K}^{-}{\pi }^{+}$	$0.1180\pm 0.0013$ pdg2020	$K_S^0{K}_L^0{K}^{-}{\pi }^{+}$	$0.1180\pm 0.0013$
$K_S^0{K}_S^0{K}^{+}{\pi }^{-}$	$0.059\pm 0.013$ pdg2020	$K_S^0{K}_S^0{K}^{+}{\pi }^{-}$	$0.059\pm 0.013$
$K_S^0{K}_L^0{K}^{+}{\pi }^{-}$	$0.1180\pm 0.0013$ pdg2020	$K_S^0{K}_L^0{K}^{+}{\pi }^{-}$	$0.1180\pm 0.0013$
$K_S^0{\rho }^0{\pi }^{+}{\pi }^{-}$	$0.11\pm 0.07$ pdg2020	$K_S^0{\pi }^{+}{\pi }^{-}{\pi }^{+}{\pi }^{-}$	$0.11\pm 0.07$
$K^{\ast -}{\pi }^{+}{\pi }^{+}{\pi }^{-}$	$0.10\pm 0.07$ pdg2020	$K_S^0{\pi }^{+}{\pi }^{-}{\pi }^{+}{\pi }^{-}$	$0.033\pm 0.023$
$K_S^0{\pi }^{+}{\pi }^{-}{\pi }^0$	$4.09\pm 0.60$ pdg2020	$K_S^0{\pi }^{+}{\pi }^{-}{\pi }^0$	$4.09\pm 0.60$
$\varphi {\pi }^{+}{\pi }^{-}$	$0.02\pm 0.01$ pdg2020	$K_S^0{K}_L^0{\pi }^{+}{\pi }^{-}$	$0.0068\pm 0.0034$
${\overline{K}}^{\ast 0}\gamma$	$0.042\pm 0.007$ pdg2020	$K_S^0{\pi }^0\gamma$	$0.007\pm 0.001$
$K_S^0{\pi }^0\eta$	$1.006\pm 0.046$ bes3-etaX	$K_S^0{\pi }^0\eta$	$1.006\pm 0.046$
$K_S^0{K}_S^0\eta$	$0.0130\pm 0.0062$ bes3-etaX	$K_S^0{K}_S^0\eta$	$0.0130\pm 0.0062$
$K_S^0{K}_L^0\eta$	$0.0270\pm 0.0124$ bes3-etaX	$K_S^0{K}_L^0\eta$	$0.0270\pm 0.0124$
$K_S^0{\pi }^{+}{\pi }^{-}\eta$	$0.280\pm 0.021$ bes3-etaX	$K_S^0{\pi }^{+}{\pi }^{-}\eta$	$0.280\pm 0.021$
$K_S^0{\pi }^0{\pi }^0\eta$	$0.176\pm 0.026$ bes3-etaX	$K_S^0{\pi }^0{\pi }^0\eta$	$0.176\pm 0.026$
$K_S^0{K}^{-}{\pi }^{+}{\pi }^0$	$0.132\pm 0.016$ bes3-D-KKpipi	$K_S^0{K}^{-}{\pi }^{+}{\pi }^0$	$0.132\pm 0.016$
$K_S^0{K}_S^0{\pi }^{+}{\pi }^{-}$	$0.053\pm 0.009$ bes3-D-KKpipi	$K_S^0{K}_S^0{\pi }^{+}{\pi }^{-}$	$0.053\pm 0.009$
$K_S^0{K}_L^0{\pi }^{+}{\pi }^{-}$	$0.106\pm 0.018$ bes3-D-KKpipi	$K_S^0{K}_L^0{\pi }^{+}{\pi }^{-}$	$0.106\pm 0.018$
$K_S^0{K}^{-}{\pi }^{+}{\pi }^0$	$0.132\pm 0.016$ bes3-D-KKpipi	$K_S^0{K}^{-}{\pi }^{+}{\pi }^0$	$0.132\pm 0.016$
$K_S^0{K}^{+}{\pi }^{-}{\pi }^0$	$0.065\pm 0.008$ bes3-D-KKpipi	$K_S^0{K}^{+}{\pi }^{-}{\pi }^0$	$0.065\pm 0.008$
$\varphi {\pi }^0$	$0.117\pm 0.004$ phipi	$K_S^0{K}_L^0{\pi }^0$	$0.040\pm 0.001$
$\varphi \eta$	$0.018\pm 0.005$ phipi	$K_S^0{K}_L^0\eta$	$0.006\pm 0.001$
$K_S^0{\pi }^0\omega$	$0.848\pm 0.055$ ksompi	$K_S^0{\pi }^0\omega$	$0.848\pm 0.055$
$K_S^0{\pi }^0{\pi }^0{\pi }^0$	$0.591\pm 0.042$ kspipi0pi0	$K_S^0{\pi }^0{\pi }^0{\pi }^0$	$0.591\pm 0.042$
$K_S^0{\pi }^{+}{\pi }^{-}{\pi }^0{\pi }^0$	$0.327\pm 0.045$ kspipi0pi0	$K_S^0{\pi }^{+}{\pi }^{-}{\pi }^0{\pi }^0$	$0.327\pm 0.045$
$K_S^0{\eta }^{\prime }$	$0.949\pm 0.016$ ref31	$K_S^0{\eta }^{\prime }$	$0.949\pm 0.016$
$K_S^0{\pi }^0$	$1.240\pm 0.022$ ref31	$K_S^0{\pi }^0$	$1.240\pm 0.022$
$K_{}^{\ast }{}_{}{}^{-}{e}^{+}{\nu }_e$	$2.03\pm 0.66$ k*enu	$K_S^0{\pi }^{-}{e}^{+}{\nu }_e$	$0.68\pm 0.22$
${\overline{K}}^0{\pi }^{-}{e}^{+}{\nu }_e$	$0.158\pm 0.034$ k*enu	$K_S^0{\pi }^{-}{e}^{+}{\nu }_e$	$0.079\pm 0.017$
$K_{}^{\ast }{}_{}{}^{-}{\mu }^{+}{\nu }_{\mu }$	$1.91\pm 0.62$ k*enu	$K_S^0{\pi }^{-}{\mu }^{+}{\nu }_{\mu }$	$0.63\pm 0.20$
$K_1^{}{}_{}{}^{-}\left(1270\right)e^{+}{\nu }_e$	$0.109\pm 0.018$ k1enu	$K_S^0{\pi }^0{\pi }^{-}{e}^{+}{\nu }_e$	$0.026\pm 0.004$
$K_1^{}{}_{}{}^{-}\left(1270\right){\mu }^{+}{\nu }_{\mu }$	$0.109\pm 0.018$ k1enu	$K_S^0{\pi }^0{\pi }^{-}{\mu }^{+}{\nu }_{\mu }$	$0.026\pm 0.004$
$K_S^0{K}_S^0$	$0.0170\pm 0.0015$ ksks	$K_S^0{K}_S^0$	$0.0170\pm 0.0015$
$K_S^0{K}_S^0{K}_S^0$	$0.072\pm 0.007$ ksks	$K_S^0{K}_S^0{K}_S^0$	$0.072\pm 0.007$
$K_S^0{K}_L^0{K}_L^0$	$0.216\pm 0.021$ ksks	$K_S^0{K}_L^0{K}_L^0$	$0.216\pm 0.021$
$K_S^0{K}_S^0{K}_L^0$	$0.216\pm 0.021$ ksks	$K_S^0{K}_S^0{K}_L^0$	$0.216\pm 0.021$
Sum	…	…	$18.16\pm 0.72$

Equations

$$N_{\mathrm{tag}}=2N_{D\bar{D}}{\mathcal{B}}_{\mathrm{tag}}{ϵ}_{\mathrm{tag}},$$

$$N_{\mathrm{DT}}=2N_{D\bar{D}}{\mathcal{B}}_{\mathrm{tag}}{\mathcal{B}}_{\mathrm{sig}}{ϵ}_{\mathrm{DT}},$$

$${\mathcal{B}}_{\mathrm{sig}}=\frac{N_{\mathrm{DT}}}{N_{\mathrm{tag}}{ϵ}_{\mathrm{sig}}},$$

$${\mathcal{B}}_{\mathrm{sig}}^{\mathrm{cor}}=\frac{1}{1-C_f(2f_{\mathrm{CP}+}-1)}{\mathcal{B}}_{\mathrm{sig}},$$

$$C_f=\frac{2r\cos \delta }{1+r^2}.$$

$$f_{\mathrm{CP}\pm }=\frac{N^{\pm }}{N^{+}+N^{-}},$$

$$N^{\pm }=\frac{M_{\mathrm{measured}}^{\pm }}{S^{\pm }},$$

$$S_{\pm }=\frac{S_{\mathrm{measured}}^{\pm }}{1-{\eta }_{\pm }{y}_D},$$

$$N_{\mathrm{DT}}^{\mathrm{net}}=N_{K_S^0\mathrm{sig}}-f_{\mathrm{co}}\cdot N_{K_S^0\mathrm{sb}}.$$

Full text

I Introduction

Weak decays of charmed mesons offer an ideal laboratory for improving our understanding of both strong and weak interactions. The $D^{+}$ and $D^{0}$ mesons decay dominantly into final states involving $K^{-}$ and $\bar{K}^{0}$ via Cabibbo-favored processes. The branching fractions of the inclusive decays $D^{+(0)}\to K^{0}/\bar{K}^{0}X$, where $X$ denotes any possible particle combinations, were first measured by Mark-III Mark3-KSX in 1991 and later by BES BESII-KSX in 2006. Assuming $\mathcal{B}(K^{0}/\bar{K}^{0}\to K^{0}_{S})=0.5$ and considering the world average branching fractions of $D^{+(0)}\to K^{0}/\bar{K}^{0}X$, the branching fractions of $D^{+(0)}\to K^{0}_{S}X$ are expected to be $(30.5\pm 2.5)\%$ and $(23.5\pm 2.0)\%$, respectively pdg2020 . In recent years, the precision of the branching fractions of known decay modes has been significantly improved, and many new channels have been observed. Following this progress, the sum of the branching fractions of the exclusive $D^{+}$ and $D^{0}$ decays involving $K^{0}_{S}$ are found to be $(31.68\pm 0.32)$% and $(18.16\pm 0.72)$%, as summarized in Appendix A. In contrast, however, the branching fractions of the inclusive decays $D^{+}\to K^{0}_{S}X$ and $D^{0}\to K^{0}_{S}X$ suffer from large uncertainties due to the limited samples sizes used to perform the measurements. Improved determinations of these inclusive decay branching fractions will help to quantify and guide the search for the missing decay modes involving a $K^{0}_{S}$ meson.

In this paper, we report the improved measurements of the branching fractions of $D^{+}\to K^{0}_{S}X$ and $D^{0}\to K^{0}_{S}X$, obtained from the analysis of $2.93\ \text{fb}^{-1}$ of $e^{+}e^{-}$ collision data lum_bes3 taken at the center-of-mass energy $\sqrt{s}$ of 3.773 GeV with the BESIII detector. The charge-conjugated processes are implied throughout the discussion.

II BESIII detector and Monte Carlo simulation

The BESIII detector BESIII records symmetric $e^{+}e^{-}$ collisions provided by the BEPCII storage ring Yu:IPAC2016-TUYA01 , which operates in the center-of-mass energy range from 2.0 to 4.95 GeV, with a peak luminosity of $1\times 10^{33}$ cm*-2s-1* achieved at $\sqrt{s}=3.773$ GeV. BESIII has collected large data samples in this energy region Ablikim:2019hff . The cylindrical core of the BESIII detector covers 93% of the full solid angle and consists of a helium-based multilayer drift chamber (MDC), a plastic scintillator time-of-flight system (TOF), and a CsI(Tl) electromagnetic calorimeter (EMC), which are all enclosed in a superconducting solenoidal magnet providing a 1.0 T magnetic field. detvis The solenoid is supported by an octagonal flux-return yoke with resistive plate counter muon identification modules interleaved with steel. The charged-particle momentum resolution at $1~{}{\rm GeV}/c$ is $0.5\%$, and the specific energy loss (${\rm d}E/{\rm d}x$) resolution is $6\%$ for electrons from Bhabha scattering. The EMC measures photon energies with a resolution of $2.5\%$ ($5\%$) at $1$ GeV in the barrel (end-cap) region. The time resolution in the TOF barrel region is 68 ps, while that in the end-cap region is 110 ps.

Simulated samples produced with the geant4-based geant4 Monte Carlo (MC) package, which includes the geometric description of the BESIII detector and the detector response, are used to determine the detection efficiency and to estimate the backgrounds. The simulation includes the beam-energy spread and initial-state radiation in the $e^{+}e^{-}$ annihilations modeled with the generator kkmc kkmc . The inclusive MC samples consist of the production of $D\bar{D}$ pairs with consideration of quantum coherence for all neutral $D$ modes, the non-$D\bar{D}$ decays of the $\psi(3770)$, the initial-state radiation production of the $J/\psi$ and $\psi(3686)$ states, and the continuum processes. The known decay modes are modeled with evtgen evtgen using the branching fractions taken from Particle Data Group (PDG) pdg2020 , and the remaining unknown decays of the charmonium states are modeled by lundcharm lundcharm . Final-state radiation from charged final-state particles is incorporated using the photos package photos .

The three-body decays $D^{+}\to K_{S}^{0}\pi^{+}\pi^{0}$ and $D^{+}\to K_{S}^{0}\pi^{+}\eta$ are modeled with a modified data-driven generator BODY3 evtgen , which was developed to simulate different intermediate states in data for a given three-body final state. The Dalitz plot from data, corrected for backgrounds and efficiencies, is taken as input for the BODY3 generator. The decay $D^{+}\to K_{S}^{0}\pi^{+}\pi^{+}\pi^{-}$ is modeled with the partial-wave analysis model from Ref. bes3-kspipipi . Other decays reported in recent BESIII measurements bes3-D-KKpipi ; ksompi ; lanxing are modeled by constructing a cocktail of possible intermediate resonances.

III Analysis Method

At $\sqrt{s}=3.773$ GeV, $D^{0}\bar{D}^{0}$ or $D^{+}D^{-}$ pairs are produced without accompanying hadron(s), thereby offering a clean environment in which to investigate hadronic $D$ decays with the double-tag (DT) method Li:2021iwf . To measure the branching fractions of the hadronic decays of $D^{+}$($D^{0}$) mesons, as a first step the single-tag $D^{-}$($\bar{D}^{0}$) mesons are selected by using the hadronic decay modes $D^{-}\to K^{+}\pi^{-}\pi^{-}$ and $\bar{D}^{0}\to K^{+}\pi^{-}$. The number of the expected single-tag $D^{-}$ ($\bar{D}^{0}$) mesons is given by

[TABLE]

where $N_{{D\bar{D}}}$ is the total number of $D^{+}D^{-}$ ($D^{0}\bar{D}^{0}$) pairs, ${\mathcal{B}}_{\rm tag}$ is the branching fraction of $D^{-}$ ($\bar{D}^{0}$) decay into the tag mode, and $\epsilon_{\rm tag}$ is the efficiency of reconstructing the single-tag $D^{-}$ ($\bar{D}^{0}$) mesons.

The double-tag sample is comprised of events in which the signal decay $D^{+(0)}\to K_{S}^{0}X$ can be selected in the presence of the single-tag $D^{-}$ ($\bar{D}^{0}$) mesons. The number of the expected double-tag events is given by

[TABLE]

where ${\mathcal{B}}_{{\rm sig}}$ is the branching fraction of $D^{+(0)}\to K_{S}^{0}X$, and $\epsilon_{{\rm DT}}$ is the efficiency of detecting the double-tag events. The branching fraction of the signal decay is then

[TABLE]

where $\epsilon_{{\rm sig}}=\frac{\epsilon_{{\rm DT}}}{\epsilon_{\rm tag}}$ is the effective efficiency of reconstructing $D^{+(0)}\to K_{S}^{0}X$.

Quantum correlations between the $D^{0}\bar{D^{0}}$ pair modify the double-tag rates. The measured branching fraction of $D^{0}\to K^{0}_{S}X$ tagged by $\bar{D}^{0}\to K^{+}\pi^{-}$ must be corrected to

[TABLE]

where $f_{CP+}$ is the fractional $CP$-even content of the decay $D^{0}\to K^{0}_{S}X$, and $C_{f}$ is the strong-phase factor

[TABLE]

Here $r=0.0587\pm 0.0002$ A1 is the ratio between doubly-Cabibbo-suppressed and Cabibbo-favored amplitudes for $D^{0}/\bar{D}^{0}\to K^{\pm}\pi^{\mp}$, and $\delta=(191.7_{-3.8}^{+3.6})^{\circ}$ A1 is the strong-phase difference between two amplitudes refcp3 ; refcp4 . These inputs yield $C_{f}=(-11.5\pm 0.2)\%$ for $D^{0}/\bar{D}^{0}\to K^{\pm}\pi^{\mp}$.

The $CP$-even fraction $f_{CP+}$ of $D^{0}\to K^{0}_{S}X$ decays is not known, but can be measured from the $D^{0}\bar{D^{0}}$ data set. Following Refs. refcp1 ; refcp2 , the $CP$-even fraction is calculated by

[TABLE]

[TABLE]

and

[TABLE]

where $N^{\pm}$ is the ratio between double-tag and single-tag yields with $CP$-odd(even) tags, $M_{\rm measured}^{\pm}$ is the number of double-tag candidates for a signal channel versus CP-odd(even) tags, $S^{\pm}_{\rm measured}$ is the number of single-tag candidates for $CP$-even(odd) decay modes, $S^{\pm}$ is the number of the corrected single-tag candidates for $CP$-even(odd) decay modes. Finally, $\eta_{\pm}=\pm 1$ for $CP$-even(odd) mode and $y_{D}=(6.47\pm 0.24)\times 10^{-3}$ is the mixing parameter taken from the latest average from HFLAV A1 .

IV Yields of single-tag $D^{-}\,(\bar{D}^{0})$ mesons

Charged kaons and charged pions are selected with the same selection criteria as those used in our previous works bes3-pimuv ; bes3-etaetapi ; bes3-etaX . All charged tracks, apart from those from $K^{0}_{S}$ decays, are required to have a polar angle $\theta$ with respect to the MDC symmetry axis within the MDC acceptance $|\rm{cos\theta}|<0.93$, a distance of closest approach with respect to the interaction point along the beam direction $|V_{z}|<$ 10 cm, and a distance of closest approach with respect to the interaction point in the plane transverse to the beam direction $|V_{xy}|<$ 1 cm. Particle identification (PID) for charged kaons and pions is performed by exploiting combined d$E/$d$x$ and TOF information. The confidence levels for pion and kaon hypotheses ($CL_{\pi}$ and $CL_{K}$) are calculated. Charged tracks satisfying $CL_{K}>CL_{\pi}$ and $CL_{\pi}>CL_{K}$ are assigned as charged kaons and pions, respectively.

When selecting $\bar{D}^{0}\to K^{+}\pi^{-}$ candidates, contamination from cosmic rays and Bhabha events is suppressed with the following requirements deltakpi . First, the two charged tracks must have a TOF time difference less than 5 ns and they must be inconsistent with being a muon-antimuon or an electron–positron pair. Second, in the event there must be at least one EMC shower with an energy greater than 50 MeV or at least one additional charged track detected in the MDC.

The tagged $D^{-}\,(\bar{D}^{0})$ mesons are identified using the energy difference $\Delta E\equiv E_{\bar{D}}-E_{\rm beam}$ and the beam-constrained mass $M_{\rm BC}\equiv\sqrt{E^{2}_{\rm beam}-|\vec{p}_{\bar{D}}|^{2}}$. Here, $E_{\rm beam}$ is the beam energy, and $\vec{p}_{\bar{D}}$ and $E_{\bar{D}}$ are the momentum and energy of the $D^{-}\,(\bar{D}^{0})$ candidate in the rest frame of $e^{+}e^{-}$ system, respectively. For each tag mode, if there are multiple candidates in an event, only the one with the smallest $|\Delta E|$ is kept. The tagged $\bar{D}$ candidates are required to be within $\Delta E\in(-25,25)$ MeV.

To determine the yields of single-tag $D^{-}\,(\bar{D}^{0})$ mesons for individual tag modes, binned maximum-likelihood fits are performed on the corresponding $M_{\rm BC}$ distributions of the accepted single-tag candidates bes3-pimuv ; bes3-etaetapi ; bes3-etaX . In the fits, the $D^{-}\,(\bar{D}^{0})$ signal is modeled by an MC-simulated shape convolved with a double-Gaussian function describing the resolution difference between data and MC simulation. The combinatorial background is described by an ARGUS function ARGUS . The best fit results for the $M_{\rm BC}$ distributions are shown in Fig. 1. The yields of the single-tag $D^{-}$ and $\bar{D}^{0}$ mesons are $798935\pm 1011$ and $529227\pm 761$ respectively, where the uncertainties are statistical only. The efficiencies of reconstructing the single-tag $D^{-}\,(\bar{D}^{0})$ mesons are estimated to be $(51.90\pm 0.08)\%$ and $(65.60\pm 0.09)\%$ respectively, by analyzing the inclusive MC sample.

V Yields of double-tag $D^{+(0)}$ mesons

The $K^{0}_{S}$ mesons are reconstructed in the $K^{0}_{S}\to\pi^{+}\pi^{-}$ decay mode. The charged pions making up the $K^{0}_{S}$ candidates are required to satisfy $|V_{z}|<$ 20 cm and $|\!\cos\theta|<$ 0.93. To improve the resolution of the $K^{0}_{S}$ invariant mass, a primary vertex and a $K^{0}_{S}$ decay vertex fit are performed on the $\pi^{+}\pi^{-}$ pairs. To suppress the combinatorial background, we require $L/\sigma_{L}>$ 2, where $\sigma_{L}$ is the uncertainty of the flight distance $L$. The invariant mass of the chosen charged pions, $M_{\pi^{+}\pi^{-}}$, is required to be within $|M_{\pi^{+}\pi^{-}}-M_{K_{S}^{0}}|<$ 12 MeV/$c^{2}$, which corresponds to about $\pm 3\sigma$ around the known $K^{0}_{S}$ mass ($M_{K_{S}^{0}}$) pdg2020 . To estimate the combinatorial $\pi^{+}\pi^{-}$ background, the $K^{0}_{S}$ sideband region is defined by $20<|M_{\pi^{+}\pi^{-}}-M_{K_{S}^{0}}|<$ 44 MeV$/c^{2}$, which corresponds to about $\pm 5\sigma$ to $\pm 11\sigma$ away from the known $K^{0}_{S}$ mass. The definitions of the $K^{0}_{S}$ signal and sideband regions are indicated in Fig 2. The background is mainly from the $D\bar{D}$ and $q\bar{q}$ ($q=u,d,s$) continuum process. Especially, events with wrong tag and one real $K^{0}_{S}$ meson will form a peak in the $M_{\pi^{+}\pi^{-}}$ distribution but do not form a peak in the $M_{\rm BC}$ distribution. There are some discrepancy in resolution between data and MC simulation and this effect will be taken into account in the systematic uncertainty of $K_{S}^{0}$ reconstruction.

The distributions of the recoil mass, momentum, and $\cos\theta$ of the $K^{0}_{S}$ candidates from the selected $D^{+}\to K_{S}^{0}X$ and $D^{0}\to K_{S}^{0}X$ events are shown in Fig. 3. The experimental and the MC distributions are in good agreement.

The decay yields of the signal are determined from the fits to the $M_{\rm BC}$ distributions of the accepted double-tag events in data. However, combinatorial background events with $M_{\pi^{+}\pi^{-}}$ in the $K^{0}_{S}$ signal region may form peaking background in the $M_{\rm BC}$ distribution. So, the net signal yields are determined by

[TABLE]

where $N_{K^{0}_{S}\rm sig}$ and $N_{K^{0}_{S}\rm sb}$ are the fitted double-tag yields in the $K^{0}_{S}$ signal and sideband regions, respectively; the sideband scaling factor $f_{\rm co}$ is taken as 0.5 under the assumption that the background contribution in the $M_{\pi^{+}\pi^{-}}$ spectrum is flat. Figure 4 shows the $M_{\rm BC}$ distributions of the accepted double-tag events in data with $M_{\pi^{+}\pi^{-}}$ in the $K^{0}_{S}$ signal and sideband regions. We perform similar fits to these spectra as for the fits to the single-tag candidates. In these fits, the parameters of the smeared double-Gaussian resolution functions are fixed at the values obtained from the fits to the single-tag $D^{-}\,(\bar{D}^{0})$ candidates. The fit results are also shown in Fig. 4. From these fits, we obtain the yields of $N_{K^{0}_{S}\rm sig}$ and $N_{K^{0}_{S}\rm sb}$, as summarized in Table 2.

The efficiencies of reconstructing the inclusive decays $D^{+(0)}\to K_{S}^{0}X$ are determined by analyzing the inclusive MC sample ($\epsilon_{\rm sig}$). Inserting the numbers of $N_{\rm DT}^{\rm net}$, $N_{\rm tag}$, and $\epsilon_{\rm sig}$ into Eq. (3), we obtain the branching fractions of $D^{+}\to K_{S}^{0}X$ and $D^{0}\to K_{S}^{0}X$.

In order to account for the effect of quantum correlations on the branching fraction of $D^{0}\to K_{S}^{0}X$, we apply the correction given in Eq. 4. This procedure requires knowledge of the $CP$-even fraction of these decays, which is measured in the data by using the $CP$-even tag $\bar{D}^{0}\to K^{+}K^{-}$ and the CP-odd tag $\bar{D}^{0}\to K^{0}_{S}\pi^{0}$. The yields of these tags and the resulting value for $f_{CP+}$ are given in Table 1. The correction factor is determined to be $1.0204\pm 0.0024$.

The net efficiencies, including all factors discussed above, and the final values for the branching fractions are summarized in Table 2.

VI Systematic uncertainties

With the double-tag method, the branching-fraction measurements have a greatly reduced sensitivity to systematic bias arising from the selection criteria on the tag side.

The systematic uncertainties in the yields of the single-tag $D^{-}$ or $\bar{D}^{0}$ mesons have been estimated in Refs. bes3-pimuv to be 0.50%, by examining the relative change between data and MC simulation yields when varying the fit range, signal shape, and end point of the ARGUS function.

The data/MC difference of the $K^{0}_{S}$ reconstruction is determined to be $(1.01\pm 0.53)\%$ by using control samples of $J/\psi\to K^{0}_{S}K^{\pm}\pi^{\mp}$ and $\phi K^{0}_{S}K^{\pm}\pi^{\mp}$ decays kso_reconstruction . After correcting the MC efficiencies to match those in data, we assign 0.53% as the systematic uncertainty per $K_{S}^{0}$ decay.

Two methods are used to estimate the systematic uncertainty associated with the $K_{S}^{0}$ multiplicity. Firstly, we vary the components of the decays containing one, two, and three $K^{0}_{S}$ mesons by $\pm 1\sigma$ according to their known branching fractions. The changes in the detection efficiencies for $D^{+}\to K^{0}_{S}X$ and $D^{0}\to K^{0}_{S}X$, 0.05% and 0.07% respectively, are assigned as the uncertainties. Secondly, we measure the detection efficiencies for all known exclusive decays and determine the change of the re-weighted detection efficiency by varying each of decays within $\pm 1\sigma$ according to the known branching fractions. By adding in quadrature all the changes in the detection efficiencies, we obtain 0.20% for $D^{+}\to K^{0}_{S}X$ and 0.28% for $D^{0}\to K^{0}_{S}X$, which are assigned as the uncertainties. We assign the larger effect coming from the two methods as the corresponding systematic uncertainty for each signal decay.

To estimate the systematic uncertainty due to the choice of the $K^{0}_{S}$ sideband region, we examine the effects on the branching fractions by shifting the $K^{0}_{S}$ sideband region by $\pm 2$ or $\pm 4$ MeV/$c^{2}$. If the difference of the re-measured branching fraction is greater than the statistical uncertainty after considering the signal correlations, it is assigned as the systematic uncertainty; otherwise, it is neglected. Following this procedure, we assign 0.15% for $D^{+}\to K_{S}^{0}X$ and 0.10% for $D^{0}\to K_{S}^{0}X$ as the corresponding systematic uncertainties.

The uncertainty on the correction factor for quantum correlations is propagated to the branching fraction of $D^{0}\to K_{S}^{0}X$ decays. The uncertainties due to the limited signal double-tag MC samples for $D^{+}\to K_{S}^{0}X$ and $D^{0}\to K_{S}^{0}X$ are estimated to be 0.25% and 0.28%, respectively. The uncertainty of the branching fraction of $K^{0}_{S}\to\pi^{+}\pi^{-}$ is 0.07%, taken from the PDG pdg2020 .

The uncertainty arising from the scale factor $f_{\rm co}$ used in Eq. (6) is assessed with the alternative scale factors derived from the inclusive MC sample, which are 0.549 and 0.534 for $D^{+}\to K_{S}^{0}X$ and $D^{0}\to K_{S}^{0}X$, respectively. The changes of the re-measured branching fractions, 0.02% and 0.05%, are assigned as the corresponding systematic uncertainties.

The systematic uncertainties are summarized in Table 3. Assuming that all these components are independent, they are summed in quadrature to give the totals 0.82% and 0.87% for the branching fractions of $D^{+}\to K_{S}^{0}X$ and $D^{0}\to K_{S}^{0}X$, respectively.

VII Summary

By analyzing 2.93 fb*-1* of $e^{+}e^{-}$ collision data taken at the center-of-mass energy of 3.773 GeV with the BESIII detector, the branching fractions of $D^{+}\to K^{0}_{S}X$ and $D^{0}\to K^{0}_{S}X$ are measured to be $(32.78\pm 0.13\pm 0.27)\%$ and $(20.54\pm 0.12\pm 0.18)\%$, respectively, where the first uncertainties are statistical and the second are systematic. These results are consistent with previous measurements by Mark-III Mark3-KSX and BES BESII-KSX within uncertainties, as summarized in Table 4. Compared with the average values from the PDG, the precision of the branching fractions of $D^{+}\to K^{0}_{S}X$ and $D^{0}\to K^{0}_{S}X$ is improved by factors of 9.0 and 8.1, respectively. Summing over the branching fractions of the known $D^{+(0)}$ decay modes containing $K^{0}_{S}$, we obtain ${\mathcal{B}}^{\rm sum}_{\rm exclusive}(D^{+}\to K^{0}_{S}X)=(31.68\pm 0.32)\%$ and ${\mathcal{B}}^{\rm sum}_{\rm exclusive}(D^{0}\to K^{0}_{S}X)=(18.16\pm 0.72)\%$. The differences between the inclusive and exclusive decay branching fractions are $(1.10\pm 0.41)\%$ and $(2.38\pm 0.75)\%$ for $D^{+}$ and $D^{0}$ decays, respectively. These results indicate that there may be some missing decay modes involving $K_{S}^{0}$ for both $D^{+}$ and $D^{0}$ yet to be observed.

VIII Acknowledgement

The BESIII collaboration thanks the staff of BEPCII and the IHEP computing center for their strong support. This work is supported in part by National Key R&D Program of China under Contracts Nos. 2020YFA0406400, 2020YFA0406300; National Natural Science Foundation of China (NSFC) under Contracts Nos. 11635010, 11735014, 11835012, 11935015, 11935016, 11935018, 11961141012, 12022510, 12025502, 12035009, 12035013, 12192260, 12192261, 12192262, 12192263, 12192264, 12192265; the Chinese Academy of Sciences (CAS) Large-Scale Scientific Facility Program; Joint Large-Scale Scientific Facility Funds of the NSFC and CAS under Contract No. U1832207, U1932102; the CAS Center for Excellence in Particle Physics (CCEPP); 100 Talents Program of CAS; The Institute of Nuclear and Particle Physics (INPAC) and Shanghai Key Laboratory for Particle Physics and Cosmology; ERC under Contract No. 758462; European Union’s Horizon 2020 research and innovation programme under Marie Sklodowska-Curie grant agreement under Contract No. 894790; German Research Foundation DFG under Contracts Nos. 443159800, Collaborative Research Center CRC 1044, GRK 2149; Istituto Nazionale di Fisica Nucleare, Italy; Ministry of Development of Turkey under Contract No. DPT2006K-120470; National Science and Technology fund; National Science Research and Innovation Fund (NSRF) via the Program Management Unit for Human Resources & Institutional Development, Research and Innovation under Contract No. B16F640076; STFC (United Kingdom); Suranaree University of Technology (SUT), Thailand Science Research and Innovation (TSRI), and National Science Research and Innovation Fund (NSRF) under Contract No. 160355; The Royal Society, UK under Contracts Nos. DH140054, DH160214; The Swedish Research Council; U. S. Department of Energy under Contract No. DE-FG02-05ER41374.

Appendix A Branching fractions of the known exclusive $D^{+(0)}$ decays involving $K_{S}^{0}$

Tables 5 and 6 summarize the branching fractions of the known exclusive $D^{+}$ and $D^{0}$ decays involving $K_{S}^{0}$, respectively.