{"id":"a59bc524-825d-4dcc-b375-ab2adea6b56d","arxiv_id":"1909.00600","paper_version":8,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A simultaneous fit of five B meson decay modes using the z±, ξm reparameterization yields robust Gaussian behavior and a small improvement in the statistical uncertainty on the CKM angle gamma compared with combining standalone fits.","lead":"This paper tests a way of measuring the CKM angle gamma by fitting several B meson decay modes at once, instead of combining separate measurements. The simultaneous method gives slightly smaller statistical errors and much more stable fits for modes that are individually hard to fit.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed advantage of simultaneous fit rests on simplified signal-only MC; realistic cross-check needed.","rationale":"The reader's weakest assumption and my concern overlap: the MC is signal-only and simplified. I am not objecting to the formalism itself; Sec. 2's derivation that z± = xi^m z± (Eq. 6) is internally plausible and the parameter-count reduction is real. The paper is transparent in Sec. 4 about omitted experimental effects. But because the paper's headline is about statistical behaviour in realistic LHCb-like conditions, the claim is stronger than the evidence. The cross-check to a more sophisticated amplitude model is mentioned but not quantified, and no code is provided for reproduction. Therefore the appropriate verdict is conditional: accept the methodology but require a realistic cross-check before relying on the quantitative precision and robustness claims. This is consistent with the reader's ACCEPT, but I would add a condition.","tokens_in":779,"tokens_out":1387,"duration_ms":129454,"concrete_test":"Run a new ensemble of 2000 pseudo-experiments at 3 fb^-1 including a binned Dalitz-plot efficiency map, a sideband-defined background component, and momentum resolution smearing, with the same five modes and world-average parameters. If simultaneous pulls remain near unity and the standalone-vs-simultaneous uncertainty difference is consistent with the current result, the claim is supported; if pull widths or bias change materially, the conclusion needs qualification.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that simultaneous analysis improves statistical behavior and reduces gamma uncertainty (Abstract, Sec. 4, Fig. 1, Table 3). The only evidence is signal-only pseudo-experiments with 2000 toys, using a simplified D->Ks pi pi amplitude model [16]. The authors explicitly state that experimental effects are not considered and expect small impact. This limitation is acknowledged, not an internal inconsistency. However, the main claimed advantage is statistical robustness near the Gaussian regime: standalone fits for B->Dpi and B->DK* are biased due to small r or low yields, while the simultaneous xi parameterization resolves this. Whether that behavior survives realistic efficiency, backgrounds, and resolution is untested. Also, the asserted cross-check with Refs. [17,18] giving identical results is not shown. Thus the central statistical advantage is conditional on unmodeled experimental conditions.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":9990,"tokens_out":12228,"duration_ms":109557,"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":[{"comment":"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.","section":"Sec. 4, combination of standalone measurements"},{"comment":"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.","section":"Sec. 4, paragraph on Gaussian behaviour"},{"comment":"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.","section":"Sec. 4, Table 3 and Fig. 1"}],"minor_comments":[{"comment":"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.","section":"Abstract and Conclusions"},{"comment":"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.","section":"Sec. 4, amplitude-model cross-check"},{"comment":"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.","section":"Sec. 4, input choices"},{"comment":"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.","section":"Tables 1-3 captions"},{"comment":"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.","section":"Fig. 2 caption"},{"comment":"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.","section":"Eq. (29)"}],"recommendation":"major_revision","confidential_remarks":"This is a methodological paper that builds directly on the authors' own preprint [1], and the incremental contribution here is the sensitivity study. The editor may wish to ensure that the relationship to [1] is clear to readers and that the standalone-combination procedure (Sec. 4) is fully documented before publication. No concerns about data integrity are raised."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The useful thing here is not the formalism—the xi reparameterization is from the previous paper by the same group—but the first multi-mode Monte Carlo test of it. The paper shows, with 2000 signal-only pseudo-experiments across five B decay modes, that the simultaneous fit gives pull widths near unity and lower outlier rates than the standalone fits, while the gain in sigma(gamma) is modest, about 0.3 degrees on average. That is a genuine and honest result, and the authors are appropriately restrained in their conclusions.\n\nThe algebra in Sec. 2 is clean, and the setup is described in enough detail to reproduce the study in principle. The authors also state their limitations plainly: no efficiency variation, no backgrounds, no resolution, no systematic correlations. The claim that the inclusion of those effects will have a small impact is an expectation, not a demonstrated result, and the asserted cross-check against the more sophisticated model in Refs. [17,18] is not shown. That is the softest spot. I would like to see either the cross-check plot or a clear statement of what was compared.\n\nThe comparison to the standalone baseline is also a little uneven, though not in a way that undermines the main point. The standalone fits for B±→Dπ± and B±→DK*± are known to be unreliable at these sample sizes, so the improved pull behavior of the simultaneous fit is partly a rescue of a poor baseline. The paper acknowledges this. Still, the baseline is the standard practice in the field, so the comparison is relevant.\n\nOther soft spots are minor. The paper does not ship the code, which would help, but the method is described well enough that an interested group could implement it. The citation pattern is fine; Ref. [1] is the source of the formalism and is legitimate, and the other references are standard.\n\nWho is this for? Anyone planning a future LHCb gamma combination, or working on statistical methods for Dalitz-based gamma measurements. It is not a breakthrough, but it is a solid, useful sensitivity study that does what it claims within stated limits. I would send it to a serious referee. The main thing I would ask the referee to push on is the unshown cross-check and the quantitative impact of the omitted experimental effects.","headline":"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.","tokens_in":10491,"tokens_out":1341,"would_cite":true,"duration_ms":15750,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A shared ratio stabilizes and sharpens CKM-angle fits","keywords":["CKM angle gamma","CP violation","B meson decays","simultaneous fit","reparameterization","Dalitz analysis","statistical sensitivity","pseudo-experiments"],"falsifier":"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.","tokens_in":9688,"feed_emoji":"⚛️","tokens_out":8791,"duration_ms":83462,"temperature":0.7,"pith_summary":"This paper proposes and tests a way to measure the CKM angle $\\gamma$ (a CP-violating phase of the quark-mixing matrix) from several $B$ meson decays at once, instead of fitting each decay mode separately and then combining the results. The key move is to write each mode's pair of complex admixture coefficients as $z^m_\\pm = \\xi^m z_\\pm$, so that all modes share the same $\\gamma$-dependent pair $z_\\pm$ and differ only through $\\gamma$-independent nuisance parameters $\\xi^m$. This cuts the number of free parameters from $4N$ to $4+2(N-1)$ and lets the fit share information across modes. In signal-only pseudo-experiments the simultaneous fit gives Gaussian pulls with widths near one and far fewer outliers, while the standalone combination shows biases as large as 40% of the statistical uncertainty in the low-rate modes. If the method holds in data, the payoff is a modest precision gain and, more importantly, reliable uncertainty estimates for future $\\gamma$ combinations.","feed_headline":"A shared ratio stabilizes and sharpens CKM-angle fits","feed_subtitle":"Sharing one amplitude ratio across five B-meson decay modes cuts fit parameters and cures the largest error-estimation failures.","key_machinery":"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.","core_discovery":"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\\%$.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Introduces the simultaneous measurement strategy that this paper develops and tests.","marker":"[1]"},{"why":"Supplies the simplified D-decay amplitude model used to generate signal-only pseudo-experiments.","marker":"[16]"},{"why":"Provides the more sophisticated D-decay model used to cross-check the simplified model and the experimental context for $D\\to K_S^0\\pi^+\\pi^-$.","marker":"[17]"},{"why":"Define the GLW observables for decays to CP eigenstates, which the simultaneous formalism recovers as ratios of $\\eta$ functions.","marker":"[5,6]"},{"why":"Define the ADS observables for favoured and suppressed final states, written in the paper as ratios of $\\eta$ functions.","marker":"[7,8]"},{"why":"Define the direct BPGGSZ approach in which $z_\\pm$ and $\\xi^m$ are fitted from multibody self-conjugate $D$ decays.","marker":"[9,10,15]"},{"why":"Defines the time-dependent $\\lambda_\\pm$ parameterization that the simultaneous fit extends to $B_s^0\\to D_s^\\mp K^\\pm$.","marker":"[11]"},{"why":"Provide the world-average hadronic parameters $r_m$ and $\\delta_m$ used to generate the pseudo-experiments.","marker":"[3,12]"}],"fun_headline_variants":["Shared ratio across B modes cuts CKM gamma uncertainty","Simultaneous B-decay fits sharpen CKM angle measurement","Joint analysis of five B decays improves CKM gamma errors","One shared ratio, five B modes, tighter CKM constraints","Multi-mode B fit reduces CKM gamma uncertainty"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Shared ratio across B modes cuts CKM gamma uncertainty","Simultaneous B-decay fits sharpen CKM angle measurement","Joint analysis of five B decays improves CKM gamma errors","One shared ratio, five B modes, tighter CKM constraints","Multi-mode B fit reduces CKM gamma uncertainty"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00068,"raw_usage":{"total_tokens":3081,"prompt_tokens":926,"completion_tokens":2155,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":542,"completion_tokens_details":{"reasoning_tokens":2075}},"tokens_in":542,"tokens_out":2155,"duration_ms":18015,"temperature":1.0,"reasoning_tokens":2075,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:41:39.998685+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"A strategy for a simultaneous measurement of $CP$ violation parameters related to the $C\\!K\\!M$ angle $\\gamma$ in multiple $B$ meson decay channels","cited_arxiv_id":"1804.05597","evidence_quote":"Introduces the simultaneous measurement strategy that this paper develops and tests."},{"cited_title":"del Amo Sanchez et al","cited_arxiv_id":null,"evidence_quote":"Supplies the simplified D-decay amplitude model used to generate signal-only pseudo-experiments."},{"cited_title":"Aaij et al","cited_arxiv_id":null,"evidence_quote":"Provides the more sophisticated D-decay model used to cross-check the simplified model and the experimental context for $D\\to K_S^0\\pi^+\\pi^-$."},{"cited_title":"Fleischer, New strategies to obtain insights into CP violation through Bs→ D± sK ∓, D∗± s K ∓ and Bd→D±π∓, D∗±π∓ decays, Nucl","cited_arxiv_id":null,"evidence_quote":"Defines the time-dependent $\\lambda_\\pm$ parameterization that the simultaneous fit extends to $B_s^0\\to D_s^\\mp K^\\pm$."}],"review_version":1}