Recognition: no theorem link
Measurement of branching fractions of D^+_sto K⁰_SK⁰_S π^+π⁰ and D^+_sto K⁰_S K^+π⁰π⁰
Pith reviewed 2026-05-12 04:16 UTC · model grok-4.3
The pith
D_s^+ decays observed with branching fractions 4.08 and 3.32 x 10^{-3}
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The hadronic decays D_s^+ to K^0_S K^0_S pi^+ pi^0 and D_s^+ to K^0_S K^+ pi^0 pi^0 are observed, with branching fractions determined to be (4.08 plus or minus 0.46 statistical plus or minus 0.45 systematic) times 10^{-3} and (3.32 plus or minus 0.64 statistical plus or minus 0.31 systematic) times 10^{-3} respectively.
What carries the argument
Reconstruction of final-state particles including K0S from pi^+ pi^- pairs and pi0 from photon pairs, followed by invariant-mass fits to extract signal yields and efficiency corrections from simulation.
If this is right
- These modes now contribute known fractions to the total decay width of the D_s^+ meson.
- The measured rates can be added to global averages of D_s^+ branching fractions used by other experiments.
- The results serve as benchmarks for Monte Carlo generators modeling multi-body charm decays.
- The same reconstruction methods can be applied to search for analogous decays in other charmed mesons.
Where Pith is reading between the lines
- The ratio of the two branching fractions could test expectations from isospin symmetry or Cabibbo-favored versus suppressed amplitudes.
- Combining these measurements with the D_s^+ lifetime would tighten the uncertainty on the total decay width.
- Dalitz-plot or amplitude analyses of the same final states could reveal resonant substructure not addressed here.
Load-bearing premise
Monte Carlo simulations accurately reproduce detector efficiencies, background shapes, and signal extraction without significant bias from fit models or data selection criteria.
What would settle it
Finding no significant excess above background in the D_s^+ candidate mass distribution for either reconstructed final state would invalidate the claimed observations and branching fractions.
Figures
read the original abstract
By analyzing $e^+e^-$ collision data corresponding to an integrated luminosity of 7.33~fb$^{-1}$ collected with the BESIII detector at center-of-mass energies ranging from 4.128 to 4.226~GeV, we report the observations of the hadronic decays $D^+_s\to K^0_SK^0_S\pi^+\pi^0$ and $D^+_s\to K^0_S K^+\pi^0\pi^0$. Their decay branching fractions are determined to be ${\mathcal B}(D^+_s\to K^0_SK^0_S \pi^+\pi^0)=(4.08\pm0.46_{\rm stat}\pm0.45_{\rm syst})\times 10^{-3}$ and ${\mathcal B}(D^+_s\to K^0_S K^+\pi^0\pi^0)=(3.32\pm0.64_{\rm stat}\pm0.31_{\rm syst})\times 10^{-3}$, where the first uncertainties are statistical and the second are systematic.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports the first observations of the decays D_s^+ → K_S^0 K_S^0 π^+ π^0 and D_s^+ → K_S^0 K^+ π^0 π^0 in 7.33 fb^{-1} of e^+e^- collision data collected by BESIII at √s = 4.128–4.226 GeV. Branching fractions are extracted via standard single- or double-tag techniques with MC-derived efficiencies, yielding B(D_s^+ → K_S^0 K_S^0 π^+ π^0) = (4.08 ± 0.46_stat ± 0.45_syst) × 10^{-3} and B(D_s^+ → K_S^0 K^+ π^0 π^0) = (3.32 ± 0.64_stat ± 0.31_syst) × 10^{-3}.
Significance. These absolute branching-fraction measurements add new experimental input on multi-body hadronic D_s decays, which can constrain models of charm decay dynamics and improve predictions for related processes. The direct use of collision data with quoted statistical and systematic uncertainties, without reliance on parameter-free derivations or self-referential fits, is a methodological strength.
major comments (1)
- The central results depend on Monte Carlo modeling of detector efficiencies and background shapes for these four-body final states. A dedicated section or appendix should quantify the agreement between data and MC in control samples (e.g., sideband or tag-side distributions) to demonstrate that any residual mismatch does not bias the extracted signal yields at the level of the quoted uncertainties.
minor comments (2)
- The abstract states the observations but does not quote the statistical significance or raw event yields; adding these numbers would help readers assess the strength of the claims immediately.
- Clarify in the text whether the analysis employs a single-tag or double-tag method and how the normalization mode is chosen, as this choice directly affects the absolute branching-fraction scale.
Simulated Author's Rebuttal
We thank the referee for the positive evaluation of our manuscript and for the constructive comment on the validation of Monte Carlo modeling. We address the point below and have revised the paper accordingly.
read point-by-point responses
-
Referee: The central results depend on Monte Carlo modeling of detector efficiencies and background shapes for these four-body final states. A dedicated section or appendix should quantify the agreement between data and MC in control samples (e.g., sideband or tag-side distributions) to demonstrate that any residual mismatch does not bias the extracted signal yields at the level of the quoted uncertainties.
Authors: We agree that explicit validation of the Monte Carlo description is important for four-body final states. In the revised manuscript we have added a new Appendix A that presents direct data-MC comparisons for the relevant kinematic variables (invariant masses of K_S^0 pairs, pion momenta, and missing-mass distributions) in both the signal region and the sideband regions, as well as for the single-tag D_s^- candidates. The comparisons demonstrate agreement within the available statistics; any small residual differences are already incorporated into the systematic uncertainties quoted in the paper. This addition confirms that modeling uncertainties do not bias the extracted yields beyond the reported errors. revision: yes
Circularity Check
No significant circularity
full rationale
This is a direct experimental measurement of two hadronic branching fractions from e+e- collision data at BESIII. The reported values are obtained by counting signal yields in data, dividing by reconstruction efficiencies derived from Monte Carlo, and normalizing to the known D_s^+ production cross section and integrated luminosity. No derivation chain reduces a claimed prediction to a fitted parameter by construction, no self-citation supplies a load-bearing uniqueness theorem, and no ansatz is smuggled in. The analysis follows standard single- or double-tag procedures with quoted statistical and systematic uncertainties; the central results remain independent of the paper's own inputs once external benchmarks (luminosity, cross sections, detector response) are accepted.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
P. U. E. Onyisi et al. (CLEO Collaboration), Phys. Rev. D 88, 032009 (2013)
work page 2013
- [2]
-
[3]
Y. K. Hsiao and C. Q. Geng, Phys. Lett. B 770, 348 (2017)
work page 2017
-
[4]
Z. Rui, Y. Li and H. n. Li, JHEP 05, 082 (2021)
work page 2021
-
[5]
Y. K. Hsiao, Phys. Lett. B 845, 138158 (2023)
work page 2023
-
[6]
M. Ablikim et al . (BESIII Collaboration), Phys. Rev. D 104, 032011 (2021)
work page 2021
-
[7]
M. Ablikim et al . (BESIII Collaboration), Phys. Rev. D 103, 092006 (2021)
work page 2021
- [8]
-
[9]
M. Ablikim et al. (BESIII Collaboration), Nucl. Instrum. Meth. A 614, 345 (2010)
work page 2010
-
[10]
C. H. Yu et al. , Proceedings of IPAC2016, Busan, Korea, 2016
work page 2016
- [11]
-
[12]
M. Ablikim et al. (BESIII Collaboration), Chin. Phys. C 39, 093001 (2015) ; Chin. Phys. C 46, 113002 (2022) . These articles described the integrated luminosity measurement for data taken at √s = 4.189, 4.199, 4.209, 4.219, and 4.226 GeV. The integrated luminosity values for the other data samples have been obtained by a similar procedure
work page 2015
-
[13]
S. Agostinelli et al. (GEANT4 Collaboration), Nucl. Instrum. Meth. A 506, 250 (2003)
work page 2003
- [14]
-
[15]
D. J. Lange, Nucl. Instrum. Meth. A 462, 152 (2001) ; R. G. Ping, Chin. Phys. C 32, 599 (2008)
work page 2001
-
[16]
J. C. Chen, G. S. Huang, X. R. Qi, D. H. Zhang and Y. S. Zhu, Phys. Rev. D 62, 034003 (2000) ; R. L. Yang, R. G. Ping and H. Chen, Chin. Phys. Lett. 31, 061301 (2014)
work page 2000
-
[17]
E. Barberio, B. van Eijk and Z. Was, Comput. Phys. Commun. 66, 115 (1991)
work page 1991
-
[18]
R. M. Baltrusaitis et al. (MARKIII collaboration), Phys. Rev. Lett. 56, 2140 (1986). J. Adler et al. (MARKIII collaboration), Phys. Rev. Lett. 60, 89 (1988)
work page 1986
-
[19]
B. C. Ke, J. Koponen, H. B. Li and Y. Zheng, Ann. Rev. Nucl. Part. Sci. 73, 285 (2023)
work page 2023
-
[20]
M. Ablikim et al. (BESIII Collaboration), Phys. Rev. D 108, 112001 (2023)
work page 2023
-
[21]
M. Ablikim et al . (BESIII Collaboration), Phys. Rev. D 83, 112005 (2011)
work page 2011
-
[23]
M. Ablikim et al. (BESIII Collaboration), Phys. Rev. D 104, 032001 (2021)
work page 2021
-
[24]
M. Ablikim et al. (BESIII Collaboration), Phys. Rev. D 92, 112008 (2015)
work page 2015
-
[25]
M. Ablikim et al. (BESIII Collaboration), Phys. Rev. D 108, 092003 (2023)
work page 2023
-
[26]
M. Ablikim et al. (BESIII Collaboration), Phys. Rev. D 109, 032011 (2024)
work page 2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.