REVIEW 3 major objections 6 minor 18 references
Study of the sensitivity to CKM angle $\gamma$ under simultaneous determination from multiple $B$ meson decay modes
T0 review · 3 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read A shared ratio stabilizes and sharpens CKM-angle fits
desk verdict A legitimate, well-scoped MC validation of the simultaneous gamma-fitting reparameterization; the precision gain is small, the statistical robustness claim is real but still conditional on signal-only toys. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the factorization $z^m_\pm = \xi^m z_\pm$ (equation 6 of the paper), combined with the $\eta(a,b,\kappa)=|a|^2+|b|^2+2\kappa\,\mathrm{Re}(a^*b)$ function used to write each decay mode's rate. The factorization treats every mode as a scaled copy of one reference pair $z_\pm$; each complex $\xi^m$ absorbs the mode-specific amplitude ratio and strong phase, leaving a single $\gamma$ shared by all modes. The parameter count drops because $\xi^m$ is $\gamma$-independent and can be profiled as a nuisance, and because the Cartesian real and imaginary parts of $z_\pm$ and $\xi^m$ are near-Gaussian, the simultaneous fit's error estimates have coverage close to unity.
What would settle it
Run the same comparison with pseudo-experiments that include a realistic efficiency map over the $D\to K_S^0\pi^+\pi^-$ phase space, background contributions, and momentum resolution. If the simultaneous fit's pull widths deviate from one or its outlier rate rises to match the standalone combination, the claimed statistical advantage is not generic.
Extended reading notes
Core claim
The paper's central claim is that the simultaneous parameterization $z^m_\pm = \xi^m z_\pm$ carries all the information on $\gamma$ from multiple $B$ decays in a smaller and better-behaved parameter space. Here $z_\pm = r e^{i\delta} e^{\pm i\gamma}$ describe the reference mode $B^\pm\to DK^\pm$, and $\xi^m$ encodes how each other mode's $D^0$--$\bar D^0$ admixture differs from that reference; by construction $\xi^m$ does not depend on $\gamma$. A fit in these Cartesian variables uses $4+2(N-1)$ parameters instead of $4N$, and the real and imaginary parts are expected to be approximately Gaussian, so uncertainties propagate more honestly. The Monte Carlo study finds that the simultaneous fit is unbiased with correct coverage for every parameter, that adding decay modes decreases the expected uncertainty on $\gamma$ from about $8.73^\circ$ to $7.89^\circ$ at $3\,\text{fb}^{-1}$ for the five modes studied, and that the per-experiment average gain over the standalone combination is about $0.3^\circ$, with underflow/overflow outliers falling from $3.5\%$ to $0.1\%$.
Load-bearing premise
The conclusions depend on the assumption that excluding experimental effects—efficiency variation, backgrounds, momentum resolution, and cross-feed—and using a simplified D-decay amplitude model leaves the relative statistical behavior of the two methods unchanged.
Editorial extensions
If this is right
- The number of independent parameters in a global $\gamma$ fit drops from $4N$ to $4+2(N-1)$, freeing degrees of freedom and simplifying covariance handling.
- For low-yield or low-interference modes such as $B^\pm\to DK^{*\pm}$ and $B^\pm\to D\pi^\pm$, standalone Cartesian fits develop biases up to 40% of the statistical uncertainty; the simultaneous fit removes these biases and restores pull widths near unity.
- Adding each extra decay mode yields a small but consistent reduction in the expected uncertainty on $\gamma$, with the five-mode simultaneous fit giving about $7.9^\circ$ at $3\,\text{fb}^{-1}$ and $3.5^\circ$ at $9\,\text{fb}^{-1}$, compared with $8.5^\circ$ and $3.6^\circ$ for the standalone combination.
- Because $B^\pm\to D\pi^\pm$ is high-statistics but has tiny interference, the paper finds it adds no precision to $\gamma$ but is statistically safe to include, making it usable as an efficiency-calibration channel in real data.
- The same parameterization extends to time-dependent decays such as $B_s^0\to D_s^\mp K^\pm$ through $\lambda_\pm = \xi^m z_\pm e^{\pm i\phi_s}$, so mixing-based and time-integrated $\gamma$ inputs can enter one fit.
Reading between the lines
- Implicit in the paper but not shown: the shared-$z_\pm$ structure would let systematic uncertainties be correlated across decay modes during the fit itself, rather than combined mode-by-mode after each separate extraction.
- A natural extension would be to fit $B^\pm\to D\pi^\pm$ as a pure calibration mode while floating $\xi^{D\pi}$; a fitted value that disagrees with external input would expose an amplitude-model or efficiency-model error that the $\gamma$ fit would otherwise inherit.
- Because $\xi^m$ is independent of $\gamma$, the same parameterization could be recycled in global CKM fits as a Gaussian prior basis, avoiding the need to convert each mode's polar coordinates to Cartesian at combination time.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a strategy for a simultaneous determination of the CKM angle γ from several B meson decay modes, based on a reparameterization of the per-mode D-meson admixture coefficients as z_m± = ξ_m z±, where ξ_m does not depend on γ. This reduces the number of independent parameters from 4N to 4+2(N−1) and is applied to the GLW, ADS, BPGGSZ, and time-dependent methods. Sensitivity studies are performed with 2000 signal-only pseudo-experiments for five B decay modes at luminosities corresponding to 3 fb^-1 and 9 fb^-1, comparing the simultaneous approach with a combination of standalone fits. The paper reports a modest statistical precision gain and, more importantly, improved Gaussian behaviour and error coverage for the simultaneous method.
Significance. If correct, the proposed parameterization is a useful practical tool for future LHCb-style γ combinations: it reduces the fit dimensionality, provides a natural way to include correlations between decay modes, and appears to stabilize fits in low-yield or small-r modes. The algebraic derivation in Sec. 2 is clean, the parameter counting is explicit, and the Monte Carlo study is self-contained in the sense that fitted parameters are validated against known input values. The paper is also honest about its main limitation, namely that the pseudo-experiments are signal-only and use a simplified D-decay amplitude model. However, because the quantitative gain is modest and conditional on idealized assumptions, the principal value of the paper is methodological rather than a definitive demonstration of a large precision improvement.
major comments (3)
- [Sec. 4, combination of standalone measurements] The baseline against which the simultaneous method is compared is not fully defined. After the standalone per-mode fits to z_m±, the text says only that 'the fitted values and full covariance matrix for each pseudo-experiment are used to determine a value for the CP-violating weak phase γ'. It is not specified how the five standalone results are combined: whether by a second-stage χ² fit in (γ, r_m, δ_m), a weighted average of per-mode estimators of γ, or some other procedure. Since the central comparison is 'simultaneous vs combination of standalone measurements', this omission makes Fig. 1 and Table 3 non-reproducible and potentially dependent on the unspecified combination algorithm. Please specify the standalone combination procedure and, if feasible, demonstrate that the qualitative conclusions are insensitive to that choice.
- [Sec. 4, paragraph on Gaussian behaviour] The main claimed advantage is improved statistical behaviour: the text states that for simultaneous fits 'the results exhibit unbiased Gaussian behaviour and the uncertainty estimates provide the appropriate coverage', while the standalone fits do not, but no quantitative pull or coverage information is provided for z±, ξ_m, or γ. The tail analysis in Fig. 2 is informative, but it applies only to a subset of events. Please add a table or figure reporting pull means, pull widths, and outlier fractions for all fitted parameters (and for γ) for both methods, so that the central robustness claim is directly testable.
- [Sec. 4, Table 3 and Fig. 1] The quantitative precision comparison is weakened by the acknowledged unreliability of the standalone uncertainties: the text states that 'the quoted values for the expected precision on γ from the standalone fits actually fluctuate around their unknown true values'. Consequently, the apparent gains in Table 3 (e.g., 8.47° vs 7.89° for Run 1 with all five modes) are hard to interpret as genuine improvements. I recommend presenting the ensemble width (RMS or fitted Gaussian width) of the distribution of fitted γ values as the primary metric for both methods, as is already done in the left panel of Fig. 2, and/or restricting the standalone comparison to modes whose standalone error estimates are reliable.
minor comments (6)
- [Abstract and Conclusions] The claim that the simultaneous analysis results in 'smaller uncertainties and improved statistical behaviour' should be explicitly qualified as applying to signal-only pseudo-experiments, since the study does not include backgrounds, efficiency variations, or systematic effects.
- [Sec. 4, amplitude-model cross-check] The sentence 'This has been cross-checked with a slightly more sophisticated model [17,18] and gives identical results' is not supported by any quantitative comparison. Please report at least one summary statistic (e.g., the resulting uncertainty on γ for both amplitude models) so that the claimed insensitivity can be assessed.
- [Sec. 4, input choices] The choice r_Dπ = 0.005 and δ_Dπ = 300° for B±→Dπ± is said to follow from loose existing constraints, but no reference or sensitivity scan is given. Since the expected contribution of this mode depends directly on these inputs, please provide a justification or a short study of the dependence of the conclusions on them.
- [Tables 1-3 captions] Please state explicitly in the captions that the quoted uncertainties are ensemble standard deviations of fitted values across the pseudo-experiments. For Table 3, clarify whether the numbers are per-experiment uncertainties or widths of the distribution of fitted γ values.
- [Fig. 2 caption] The label 'nToys' in the left panel is not defined in the caption. In addition, the underflow/overflow percentages quoted in the text (0.1% and 3.5%) are not visible in the figure; please indicate them in the caption or in the panel.
- [Eq. (29)] Please verify the normalization convention in Eq. (29): the expression appears to contain both a leading factor 1/2 and a division by 2 in the numerator, which may be a typographical issue.
Circularity Check
No significant circularity: the reparameterization is derived in-paper and the sensitivity study is a self-contained Monte Carlo closure test against generated truth.
full rationale
The derivation chain is self-contained. The simultaneous formalism rests on the algebraic identity z_m± = z_m e^{±iγ} (Eq. 4) and the scaling z_m± = ξ_m z± (Eq. 6), both derived explicitly in the paper rather than imported as unexplained inputs. The parameter-count reduction from 4N to 4+2(N−1) is an exact counting consequence of this reparameterization, not a fitted result. The central sensitivity claims are obtained from 2000 signal-only pseudo-experiments in which γ = 70° and hadronic parameters are fixed input values, and the fitted z± and ξ_m are compared with the known generated values through pulls and coverage checks. Thus the "predictions" are Monte Carlo validation of the estimator, not fitted physical constants presented as predictions. The only self-citation, Ref. [1], identifies the prior methodology, but the relevant equations and the MC study are reproduced and evaluated in this paper, so the citation is not load-bearing. The acknowledged omission of experimental effects (efficiency, backgrounds, resolution) is a stated limitation of the sensitivity study, not a circular step. No quoted reduction of a claimed result to its own input was found.
Assumptions & free parameters
free parameters (2)
- rDπ and δDπ =
0.005 and 300 degrees
- γ test value =
70 degrees
assumptions (4)
- domain assumption The D decay amplitude model of Ref [16] is an adequate description of D0→KS0 π+π− phase space for sensitivity studies.
- domain assumption Experimental effects (efficiency, backgrounds, resolution, cross-feed) can be neglected with small impact.
- domain assumption No direct CP violation in D decays for ADS observables.
- standard math Cauchy-Schwarz inequality for |Xαβ| ≤ 1 and standard algebra.
Cite this review
Pith. "Pith review of Study of the sensitivity to CKM angle $\gamma$ under simultaneous determination from multiple $B$ meson decay modes." pith.science (2026). https://pith.science/paper/AG5PMUSA
@misc{pith2026190900600,
author = {Pith},
title = {Pith review of: Study of the sensitivity to CKM angle $\gamma$ under simultaneous determination from multiple $B$ meson decay modes},
year = {2026},
howpublished = {\url{https://pith.science/paper/AG5PMUSA}},
note = {Machine review of arXiv:1909.00600}
}
abstract
Several methods exist to measure $CP$ violation observables related to the CKM unitarity triangle angle $\gamma$ using $B$ meson decays. These observables are different for every $B$ meson decay considered, although the information they contain on $\gamma$ is encoded in a similar way for all of them. This paper describes a strategy for a simultaneous measurement of $\gamma$ using several $B$ meson decays that takes into account possible correlations between them. Sensitivity studies demonstrate that the simultaneous analysis of several $B$ meson decay modes results in smaller uncertainties and improved statistical behaviour compared to a combination of standalone measurements.
Figures
Reference graph
Works this paper leans on
-
[1]
J. Garra Tic´ o, A strategy for a simultaneous measurement of CP violation parameters related to the CKM angle γ in multiple B meson decay channels, arXiv:1804.05597
-
[2]
Particle Data Group, P. A. Zyla et al., Review of particle physics, to be published in Prog. Theor. Exp. Phys. 2020 (2020) 083C01
work page 2020
-
[3]
Heavy Flavor Averaging Group, Y. Amhis et al., Averages of b-hadron, c-hadron, and τ- lepton properties as of 2018 , arXiv:1909.12524, updated results and plots available at https://hflav.web.cern.ch
arXiv 2018
-
[4]
T. Gershon and V. V. Gligorov, CP violation in the B system, Rept. Prog. Phys. 80 (2017) 046201, arXiv:1607.06746
work page Pith review arXiv 2017
-
[5]
M. Gronau and D. London, How to determine all the angles of the unitarity Triangle from B0 d→DKS and B0 s→Dφ, Phys. Lett. B 253 (1991) 483
work page 1991
-
[6]
M. Gronau and D. Wyler, On determining a weak phase from CP Asymmetries in charged B Decays, Phys. Lett. B 265 (1991) 172
work page 1991
- [7]
- [8]
Show all 18 references
-
[9]
A. Giri, Y. Grossman, A. Soffer, and J. Zupan, Determining γ using B±→ DK ± with multibody D decays, Phys. Rev. D 68 (2003) 054018
2003
-
[10]
A. Bondar. Proceedings of BINP special analysis meeting on dalitz analysis , 24-26 Sep. 2002, unpublished
2002
-
[11]
Fleischer, New strategies to obtain insights into CP violation through Bs→ D± sK ∓, D∗± s K ∓ and Bd→D±π∓, D∗±π∓ decays, Nucl
R. Fleischer, New strategies to obtain insights into CP violation through Bs→ D± sK ∓, D∗± s K ∓ and Bd→D±π∓, D∗±π∓ decays, Nucl. Phys. B 671 (2003) 459
2003
-
[12]
Tanabashi et al
Particle Data Group, M. Tanabashi et al. , Review of particle physics , Phys. Rev. D 98 (2018) 030001
2018
-
[13]
Aaij et al., Measurement of the CKM angle γ from a combination of LHCb results , JHEP 12 (2016) 087
LHCb collaboration, R. Aaij et al., Measurement of the CKM angle γ from a combination of LHCb results , JHEP 12 (2016) 087
2016
-
[14]
Aaij et al., Update of the LHCb combination of the CKM angle γ, LHCb-CONF-2018-002 (2018)
LHCb collaboration, R. Aaij et al., Update of the LHCb combination of the CKM angle γ, LHCb-CONF-2018-002 (2018)
2018
-
[15]
Poluektov et al
Belle Collaboration, A. Poluektov et al. , Measurement of φ3 with dalitz plot analysis of B±→D(∗)K ± decay, Phys. Rev. D 70 (2004) 072003
2004
-
[16]
del Amo Sanchez et al
BABAR Collaboration, P. del Amo Sanchez et al. , Measurement of D0−D 0 mixing pa- rameters using D0→K0 Sπ+π− and D0→K0 SK+K − decays, Phys. Rev. Lett. 105 (2010) 081803
2010
-
[17]
Aaij et al
LHCb collaboration, R. Aaij et al. , Measurement of CP violation and constraints on the CKM angle γ in B±→DK ± with D→K0 Sπ+π− decays, Nucl. Phys. B 888 (2014) 169
2014
-
[18]
Garra Tic´ o,The cfit fitting package , http://www.github.com/cfit
J. Garra Tic´ o,The cfit fitting package , http://www.github.com/cfit. 11
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.