{"id":"30d7c9db-ab26-4ee3-b68f-58d0d304c5fd","arxiv_id":"1908.09499","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"First measurement of the CKM angle phi3 from B±->D(K_S0 pi+ pi- pi0)K± decays at Belle, giving phi3 = (5.7 +10.2 -8.8 ±3.5 ±5.7) degrees.","lead":"Physicists at the Belle experiment measured the CKM angle phi3 using a new decay mode, B± -> D(K_S0 pi+ pi- pi0) K±, for the first time. The result is imprecise but consistent with the world average, and it adds a new channel for future precision tests of the Standard Model.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The extraction of phi3 relies on externally measured strong-phase inputs whose systematic uncertainties may not transfer cleanly to the Belle environment; a bias in c_i/s_i would propagate directly into the quoted central value.","rationale":"The reader's weakest_assumption correctly identifies the external strong-phase inputs as the most load-bearing step. I agree that the paper internally is a well-executed measurement. However, the lack of any data-driven cross-check of the c_i/s_i inputs in the Belle environment is a specific, testable gap. The paper itself acknowledges the bin-1 resolution issue as significant enough to require a scale factor, and the c_i/s_i uncertainties are quoted as only statistical and systematic in ref. [13] without a full covariance transfer. A conditional acceptance is appropriate: the central claim is plausible, but it carries a hidden assumption about the external inputs that can and should be tested with the same data set. The verdict should not be REJECT because this is a standard practice in the field and the quoted uncertainties are honestly assigned, but the specific check would either confirm or weaken the claim.","tokens_in":23582,"tokens_out":1511,"duration_ms":13643,"concrete_test":"Perform a consistency check where the Belle B±→DK± data are fitted to extract phi3 using the c_i/s_i inputs from ref. [13], and compare the resulting phi3 with the value obtained after varying the c_i/s_i by their full correlated uncertainties (including the covariance matrix, not just individual ±1σ values). More specifically, recompute x± and y± using Eq. (2.5)-(2.6) with c_1/s_1 replaced by values obtained after applying the same Gaussian-smearing scale factors (1.13 and 1.09) to the CLEO-c bin boundaries, and check whether the resulting phi3 shifts by more than the quoted ±5.7° third uncertainty. If the shift is smaller than the uncertainty, the concern is resolved; if larger, the central claim is weakened.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is a first model-independent phi3 measurement with quoted uncertainties (5.7 +10.2/-8.8 ±3.5 ±5.7) degrees. The weakest load-bearing assumption is that the external strong-phase inputs c_i and s_i from CLEO-c (ref [13], Tables 1 and 2) are unbiased for the Belle analysis, and that their quoted uncertainties fully cover experimental differences between the psi(3770) and Upsilon(4S) environments. The analysis relies on these inputs without independent validation: the binning scheme is taken from ref. [13], and the same nine-bin partition is assumed to be self-conjugate and viable for the four-body phase space (Section 3, Table 1). A bias in c_i/s_i would propagate directly through Eqs. (2.5) and (2.6) into x± and y± and hence into phi3. The only bin-specific mitigation is the Gaussian smearing applied to bin 1 (Sections 5.2 and 9), and the paper does not provide a test of whether the c_i/s_i values are consistent with the Belle data themselves. If the strong-phase inputs are off by a correlated amount (e.g., due to the bin-1 migration or the efficiency profile used in ref. [13]), the central value of phi3 could shift by more than the nominal third uncertainty. The paper recognizes this limitation only through the quoted systematic uncertainty, not through a cross-check against the data.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports the first model-independent measurement of the CKM angle φ3 using B±→D(K_S^0 π^+π^-π^0)K± decays, based on the full Belle data set of 772×10^6 BB events. The analysis follows the binned Giri-Grossman-Soffer-Zupan method: the five-dimensional D→K_S^0π^+π^-π^0 phase space is divided into nine exclusive bins, the c_i/s_i strong-phase inputs are taken from an external CLEO-c measurement (ref. [13]), the K_i and \\bar{K}_i fractions are measured from a Belle D*±−tagged sample, and the parameters (x±,y±) are extracted from a simultaneous fit to 36 B±→DK± and B±→Dπ± samples. The quoted results are φ3=(5.7^{+10.2}_{-8.8}±3.5±5.7)° and rB=0.323±0.147±0.023±0.051, with a 95% confidence interval for φ3 of (-29.7,109.5)°, consistent with the current world average. The paper also includes a combination with previous Belle multibody Dalitz analyses, giving φ3=(74^{+13}_{-14})°.","tokens_in":23833,"tokens_out":13833,"duration_ms":129056,"significance":"If the result holds, this paper establishes a new D-decay mode for model-independent φ3 measurements, introducing a nine-bin partition of a four-body phase space and demonstrating that the mode is viable despite the large statistical uncertainty. The analysis is technically careful: it addresses efficiency and migration corrections, uses a B→Dπ control sample for validation, performs a simultaneous fit over all bins and charge modes, lists a comprehensive systematic table, and uses Feldman-Cousins intervals with an explicit check of the Gaussian assumption. The paper also provides concrete reproducibility material, including bin definitions, efficiency tables, migration matrices, and fit projections. The measurement is statistics-limited and the strong-phase input contributes a sizable systematic component, so the result is best viewed as a proof-of-principle for a future Belle II/BESIII program; the authors state this interpretation clearly. I saw no circularity: the external c_i/s_i inputs and the flavour-tag fractions are independent of the target parameter φ3.","major_comments":[{"comment":"As printed, the first term in Eq. (2.5) is K_i and the second is \\bar{K}_i. For B−→DK−, the favoured amplitude produces \\bar{D}^0, so the leading term should be \\bar{K}_i and the suppressed term K_i. If the overline on the first K_i was lost in typesetting, this is a local rendering problem; however, if the formula was actually used with this ordering, the interpretation of x− and y− would be altered and the derivation in Section 2 would be inconsistent. Please confirm that the overlines in Eqs. (2.5) and (2.6) are correct in the manuscript and that the fit uses the correct ordering.","section":"Eq. (2.5)"}],"minor_comments":[{"comment":"The statement that the K_i and \\bar{K}_i values are \"in reasonable agreement\" with ref. [13] is not quantified. Since the authors note a deviation larger than 3σ in bin 9, which contains only 1.2% of the data, please state whether replacing K_9 and \\bar{K}_9 with the values from ref. [13] changes (x±,y±) by a negligible amount or include the numerical effect in the systematics discussion.","section":"Section 7"},{"comment":"The caption says the combination is shown by the \"solid blue curve,\" but three curves are drawn in the figure; please identify the curves unambiguously (e.g., blue for the combination, green for D→K_S^0π^+π^-, brown for D→K_S^0π^+π^-π^0).","section":"Figure 10"},{"comment":"The resolution scale factors 1.13±0.02 and 1.09±0.02 are introduced for the ω-mass smearing in the B and D* samples, respectively; it would be helpful to state explicitly how these factors are varied in the \"m_{πππ0} resolution\" row of Table 9.","section":"Section 5.2"},{"comment":"The footnote about the statistical uncertainties of the constrained s_i values marked with an asterisk is useful; adding a similar note about the treatment of correlated uncertainties for the c_i values would make the propagation of the external inputs easier for the reader to verify.","section":"Section 3 / Table 2"}],"recommendation":"minor_revision","confidential_remarks":"This is a solid Belle measurement with a clear methodology and an honest assessment of the large statistical and strong-phase uncertainties. The main point to check in revision is the apparent typographical issue in Eq. (2.5) concerning the placement of K_i and \\bar{K}_i; this is likely a rendering artifact, but it should be made unambiguous before publication. The skeptic's concern about the transferability of the CLEO-c strong-phase inputs is adequately addressed by the propagated third uncertainty and the internal K_i/\\bar{K}_i comparison; I do not regard it as a blocking issue."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: this is a legitimate first measurement, but the physics payoff is mainly methodological. The result phi3 = 5.7 +10.2/-8.8 +/-3.5 +/-5.7 degrees is consistent with the world average and the 95% interval is roughly (-30, 110) degrees, so it does not sharpen anything on its own. What it does is show that D -> K_S pi+ pi- pi0 is a viable channel for model-independent gamma/phi3 extraction, and it gives Belle II a place to start.\n\nThe paper does the job carefully. The GGSZ formalism is standard, and the authors adapt it to a five-dimensional phase space binned into nine exclusive regions. They use the CLEO-c strong-phase inputs from ref. [13], measure K_i and Kbar_i from their own D* sample, and fit the B+->DK+ and B+->Dpi+ yields in a simultaneous 36-sample fit. The systematic table is thorough: efficiency, migration, resolution, PDF shapes, PID, fit bias. The B->Dpi control channel, the Feldman-Cousins intervals, and the check of the Gaussian likelihood assumption are all appropriate. The combination with the earlier Belle K_S pi pi and K*0 modes moves the Belle average from 78 degrees to 74 degrees, which is reasonable and consistent.\n\nThe soft spot is the one the authors themselves flag as the third uncertainty. The values of c_i and s_i come from CLEO-c, measured at the psi(3770), and the paper assumes they transfer to the Belle environment up to the corrections in sections 5.2 and 9. There is no data-level cross-check that these external inputs are consistent with the Belle D* sample, and the binning is inherited from ref. [13] rather than optimized for this measurement. A correlated bias in the strong phases would not show up in the quoted uncertainty and would shift phi3 directly. I don't consider this a fundamental flaw: it is external-input risk inherent to the model-independent program, and the authors treat it honestly by separating the uncertainty. But a referee should ask whether a closure test using the B->Dpi sample in the (x+, y+, x-, y-) plane, where CP violation is expected to be negligible, was considered. The paper notes the LHCb r_Dpi constraint, so the idea is there.\n\nMinor criticism: no public code or data, but that is normal for Belle analyses. The large statistical uncertainty will limit the paper's near-term impact, but that is the nature of a first measurement.\n\nBottom line: this deserves a serious referee. I would accept it, with a request for a bit more discussion of the strong-phase transfer and the Dpi control sample. For anyone planning Belle II gamma measurements, this is a useful reference.","headline":"A legitimate first measurement that opens a new D four-body channel for gamma/phi3, with uncertainties too large to move the world average but a thorough experimental job that deserves a serious referee.","tokens_in":25349,"tokens_out":3169,"would_cite":true,"duration_ms":34313,"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":"This paper reports the first model-independent measurement of the CKM angle $\\phi_3$ using $B^{\\pm}\\to D(K_S^0\\pi^+\\pi^-\\pi^0)K^{\\pm}$ decays, finding $\\phi_3 = (5.7^{+10.2}_{-8.8}\\pm 3.5 \\pm 5.7)^{\\circ}$, consistent with the world…","keywords":["CKM angle phi3","gamma","B± → D K± decays","model-independent","strong-phase difference","D → KS0 pi+ pi- pi0","CP violation","Belle experiment"],"falsifier":"A measurement of the same mode with a much larger data set and independently determined strong-phase inputs — for example, by re-extracting $(c_i,s_i)$ from an order-of-magnitude larger charm-threshold sample and seeing whether the binned $B^{\\pm}$ yields still converge to the same $\\phi_3$ — would settle the question; a shift larger than the quoted uncertainties would show the external inputs do not transfer.","tokens_in":23396,"feed_emoji":"⚛️","tokens_out":9149,"duration_ms":87222,"temperature":0.7,"pith_summary":"This paper reports the first model-independent measurement of the CKM unitarity-triangle angle $\\phi_3$ using the four-body decay chain $B^{\\pm}\\to D(K_S^0\\pi^+\\pi^-\\pi^0)K^{\\pm}$, where $D$ denotes a $D^0$ or $\\overline{D}^{0}$ meson. The analysis uses the full Belle data set of $772\\times 10^{6}$ $B\\overline{B}$ events and external strong-phase inputs measured at the charm threshold. It obtains $\\phi_3 = (5.7^{+10.2}_{-8.8}\\pm 3.5 \\pm 5.7)^{\\circ}$ and $r_B = 0.323 \\pm 0.147 \\pm 0.023 \\pm 0.051$, with a $95\\%$ confidence interval $(-29.7, 109.5)^{\\circ}$ that is consistent with the world average. A sympathetic reader should take this as evidence that the four-body $D\\to K_S^0\\pi^+\\pi^-\\pi^0$ mode is a viable addition to the suite of $B\\to DK$ decays used to pin down $\\phi_3$, with the dominant uncertainties coming from statistics and from external strong-phase knowledge.","feed_headline":"First phi3 measurement from four-body D decays","feed_subtitle":"Full Belle sample yields phi3 = (5.7 +10.2/-8.8 ±3.5 ±5.7) degrees, consistent with the world average.","key_machinery":"The engine is the binned partial-rate identity built from Eqs. (2.5) and (2.6): for each phase-space bin $i$, the $B^{\\pm}$ partial widths are proportional to $K_i + r_B^2 \\overline{K}_i + 2\\sqrt{K_i\\overline{K}_i}(c_i x_{\\pm} - s_i y_{\\pm})$, where $K_i$ and $\\overline{K}_i$ are flavour-tagged $D$ fractions, $c_i$ and $s_i$ are the external amplitude-weighted strong-phase averages, and $(x_{\\pm},y_{\\pm})$ encode $r_B$, $\\delta_B$, and $\\phi_3$. The analysis determines $K_i$ from $D^{*\\pm}$-tagged events, extracts $(x_{\\pm},y_{\\pm})$ from a simultaneous likelihood fit to $B^{\\pm}\\to D K^{\\pm}$ and the control mode $B^{\\pm}\\to D\\pi^{\\pm}$, corrects for efficiency and bin-migration using a migration matrix, and finally converts $z=(x_+,y_+,x_-,y_-)$ into $(\\phi_3,r_B,\\delta_B)$ with a frequentist confidence-interval construction. The machinery works because the binning keeps the strong-phase variation within each bin small enough that the external $c_i$, $s_i$ inputs capture the interference.","core_discovery":"The central claim is that $\\phi_3$ can be extracted from $B^{\\pm}\\to D(K_S^0\\pi^+\\pi^-\\pi^0)K^{\\pm}$ without resorting to a model of the charm-meson decay amplitude. The paper divides the five-dimensional $D$ phase space into nine exclusive bins chosen around resonances such as $\\omega$, $K^{*}$, and $\\rho$, and uses externally measured bin-averaged strong-phase parameters $c_i$, $s_i$ as input. From the binned $B^{\\pm}\\to D K^{\\pm}$ yields it determines the Cartesian parameters $x_{\\pm} = r_B\\cos(\\delta_B \\pm \\phi_3)$ and $y_{\\pm} = r_B\\sin(\\delta_B \\pm \\phi_3)$, then converts them into $\\phi_3$, $r_B$, and $\\delta_B$ through a frequentist treatment. The quoted result $\\phi_3 = (5.7^{+10.2}_{-8.8}\\pm 3.5 \\pm 5.7)^{\\circ}$ is consistent with the world average within two standard deviations, and the paper notes a local likelihood minimum near $\\phi_3 = 75^{\\circ}$, $\\delta_B = 155^{\\circ}$, reflecting the known twofold ambiguity.","pith_inferences":["If the same binned machinery were run with a finer binning informed by an amplitude model, a shift in $\\phi_3$ beyond the quoted total uncertainty would indicate that strong-phase variation inside the current nine bins is not fully captured by the bin-averaged inputs.","Because this analysis deliberately leaves $B^{\\pm}\\to D\\pi^{\\pm}$ out of the $\\phi_3$ extraction, a natural extension is to check whether its measured $(x_{\\pm},y_{\\pm})$, currently consistent with zero CP violation, would shift the central value when included; that check is not performed here.","A cross-experiment transfer check would be to compare the $c_i,s_i$ inputs obtained from two independent charm-threshold data sets; consistency within uncertainties would validate the approach for other four-body channels, and any discrepancy would set the floor on the external uncertainty."],"forward_implications":["The four-body final state $D\\to K_S^0\\pi^+\\pi^-\\pi^0$ becomes a usable, model-independent channel for $\\phi_3$; its large branching fraction and resonance-rich phase space make it competitive with the three-body $K_S^0\\pi^+\\pi^-$ mode.","Combining this result with Belle's earlier model-independent $D\\to K_S^0\\pi^+\\pi^-$ measurements shifts the combined central value to $\\phi_3 = (74^{+13}_{-14})^{\\circ}$, a modest improvement over $(78^{+14}_{-15})^{\\circ}$ without it.","The quoted $95\\%$ confidence interval $(-29.7, 109.5)^{\\circ}$ contains the world average, so the measurement adds a tree-level $\\phi_3$ constraint that is consistent with the standard model and offers no hint of new physics.","The three quoted uncertainties — statistical, experimental systematic, and external strong-phase — show where future gains lie: larger $B$ samples and better charm-threshold inputs both directly reduce the total error."],"supporting_citations":[{"why":"Supplies the nine-bin strong-phase averages $c_i$ and $s_i$ used as external inputs in Eqs. (2.5) and (2.6).","marker":"[13]"},{"why":"Introduces the model-independent binned method for determining $\\gamma/\\phi_3$ from $B^{\\pm}\\to DK^{\\pm}$ with multibody $D$ decays.","marker":"[8]"},{"why":"Previous Belle measurement of $\\phi_3$ with $D\\to K_S^0\\pi^+\\pi^-$ whose analysis procedure and frequentist treatment this paper extends.","marker":"[11]"},{"why":"Provides the confidence-interval ordering used to convert $(x_{\\pm},y_{\\pm})$ into confidence intervals on $\\phi_3$, $r_B$, and $\\delta_B$.","marker":"[37]"},{"why":"Supplies the world-average values used for consistency comparison and for the confidence-contour cross marks.","marker":"[36]"},{"why":"Prior Belle model-independent analysis of $B^0\\to D K^{*0}$ combined with this result in the multibody $D$ average.","marker":"[38]"}],"fun_headline_variants":["First phi3 from four-body D decays at Belle","Model-independent phi3 measured with KS pi pi pi0","Belle measures phi3 with four-body D K decays","Full Belle sample gives phi3 from 4-body D"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole analysis assumes that the measurements of how the D meson's decay amplitude changes in phase across the nine regions, made with charm-meson data at a different collision energy, are still correct when used to interpret the B-meson data, with any slight mismatch in mass resolution absorbed by one extra smearing adjustment.","fun_headline_variants_meta":{"raw":{"variants":["First phi3 from four-body D decays at Belle","Model-independent phi3 measured with KS pi pi pi0","Belle measures phi3 with four-body D K decays","Full Belle sample gives phi3 from 4-body D"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000821,"raw_usage":{"total_tokens":3676,"prompt_tokens":1111,"completion_tokens":2565,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":727,"completion_tokens_details":{"reasoning_tokens":2499}},"tokens_in":727,"tokens_out":2565,"duration_ms":19905,"temperature":1.0,"reasoning_tokens":2499,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:10:12.037402+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A measurement of the same mode with a much larger data set and independently determined strong-phase inputs — for example, by re-extracting $(c_i,s_i)$ from an order-of-magnitude larger charm-threshold sample and seeing whether the binned $B^{\\pm}$ yields still converge to the same $\\phi_3$ — would settle the question; a shift larger than the quoted uncertainties would show the external inputs do not transfer.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the confidence-interval ordering used to convert $(x_{\\pm},y_{\\pm})$ into confidence intervals on $\\phi_3$, $r_B$, and $\\delta_B$."},{"cited_title":"Negishi et al.,First model-independent Dalitz analysis of B0→DK∗0, D→K 0 Sπ+π− decay, Prog","cited_arxiv_id":null,"evidence_quote":"Prior Belle model-independent analysis of $B^0\\to D K^{*0}$ combined with this result in the multibody $D$ average."}],"review_version":1}