{"id":"2cbb3e16-27ab-4ba2-807d-6554bbf01197","arxiv_id":"2411.15423","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A new event-based likelihood analysis of gamma p to pi+ pi- p CLAS data gives cross sections, spin-density matrix elements, and N rho(770) branching fractions for many N* and Delta* resonances.","lead":"This paper presents an event-based likelihood analysis of CLAS data on photon-proton collisions producing two charged pions and a proton. It extracts new branching fractions for nucleon and Delta resonances decaying into a rho meson and a nucleon, with direct implications for hadron spectroscopy.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The N rho branching fractions in Table IV depend on the BnGa assumption that no coherent non-resonant N pi pi amplitude is present; the paper shows no test of this, and the SAPHIR comparison indicates the cross sections shift when such a background is included.","rationale":"The reader's weakest_assumption identifies the BnGa isobar-model assumption that non-resonant N pi pi production is absent, which is exactly the point on which the central claim depends. I agree with that assessment and with the CONDITIONAL verdict. The paper has real strengths: it is based on a large CLAS data set, checks four event topologies for consistency, reports pull distributions, and gives systematic uncertainties from a spread of fits. These support the qualitative findings of strong rho0 and Delta++ production. However, the quantitative branching fractions are only as secure as the model assumption that all unpredicted strength can be absorbed into isobar and exchange amplitudes. The SAPHIR comparison in Section IV.C is direct evidence that adding a phase-space background changes the N rho cross section considerably, so the omission of such a background in BnGa is not innocuous. The lack of a sensitivity test, plus the reliance on unpublished formulas in Ref. [6], means the paper should not be accepted as the last word on the branching fractions. Since the reader already reached CONDITIONAL, no verdict change is needed; the concrete test above would either lift or confirm the remaining caveat.","tokens_in":26208,"tokens_out":4408,"duration_ms":42426,"concrete_test":"Augment the BnGa fit used for the event-based sample with a coherent non-resonant N pi pi amplitude, e.g., a smooth energy-dependent background in the lowest partial waves with floating magnitude and phase, and refit exactly the same data and weighting scheme. Compare every entry in Table IV (total and spin-separated branching fractions) with the published values. If any entry shifts by more than its quoted uncertainty, the 'uniquely determined' claim is falsified and the branching fractions should be presented as model-dependent; if no entry shifts, the assumption is validated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claims, the N rho branching fractions in Table IV, are extracted from a BnGa coupled-channel fit in which the reaction is, by construction, decomposed into quasi-two-body isobar channels plus t-channel exchange. Section IV.C states this explicitly: 'Non-resonant production of N pi pi is not needed.' The abstract and Summary go further, calling the branching ratios 'uniquely determined by the new CLAS data.' For that uniqueness claim to hold, the fit must be able to distinguish a coherent non-resonant N pi pi amplitude from modified resonance couplings, and the paper provides no such test. The SAPHIR comparison in Section IV.C shows why this is load-bearing: SAPHIR fitted the same reaction with a large phase-space (non-resonant) contribution and obtained a markedly smaller N rho excitation function, while BnGa obtains a larger one by omitting that background. The difference is attributed to model choices, not to the data. Since amplitudes are complex and unitarized, a background term can interfere with and partially substitute for resonance contributions, shifting the extracted branching fractions. The condition is also unverifiable as published: the explicit BnGa formulas are deferred to Ref. [6] (in preparation), so the model is not fully specified. These facts do not show the result is wrong, but they show the central claim rests on an untested model assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a high-statistics measurement of the reaction γp → π+π−p using CLAS g11a data, with about 400 million events processed through four event topologies and a subsample of nearly 1.8 million events used in an event-based likelihood fit. The analysis reports total and differential cross sections for the isobars pρ0(770), Δ(1232)++π−, and Δ(1232)0π+, spin-density matrix elements, and, from the Bonn–Gatchina coupled-channel analysis, branching ratios of N* and Δ* resonances into Nρ(770) (Table IV). The authors state that this is the first extraction of Nρ branching ratios from an event-based likelihood fit to γp → π+π−p and that the new CLAS data uniquely determine these branching fractions.","tokens_in":26520,"tokens_out":4300,"duration_ms":38029,"significance":"If the central claims hold, the paper provides a comprehensive resonance-region dataset for γp → π+π−p and the first event-based, multi-topology likelihood extraction of Nρ branching fractions, which would be a valuable step for baryon spectroscopy. The internal consistency checks are genuine strengths: the four event topologies agree after acceptance correction (Fig. 7), the kinematic-fit pulls are Gaussian with unit width (Fig. 4), the three-pion background is simulated and quantified (Section II.D), and systematic uncertainties are estimated from the spread of PWA solutions rather than from a single fit (Table II, Section IV.C). The cross sections and spin-density matrix elements will be useful references for future coupled-channel analyses even if the model-dependent branching ratios shift.","major_comments":[{"comment":"The central quantitative claim, that the Nρ branching ratios in Table IV are 'uniquely determined by the new CLAS data' (Section VI), rests on the untested BnGa assumption that the reaction contains no coherent non-resonant Nππ amplitude. Section IV.C states 'Non-resonant production of Nππ is not needed,' but no test is shown in which such an amplitude is added and its effect on Table IV is quantified. The SAPHIR comparison in Fig. 12 makes this load-bearing: SAPHIR fitted the same reaction with a large phase-space contribution and obtained a markedly smaller Nρ excitation function, and the difference is attributed to model choices, not to the data. Because amplitudes are complex and can interfere, the omitted term could be partially absorbed by modified resonance couplings, shifting the extracted branching fractions. Please add a concrete test (e.g., fits with a model non-resonant Nππ amplitude or an equivalent phase-space/Deck term) and report the resulting variations in Table IV, or restrict the claims accordingly.","section":"Section IV.C / Table IV"},{"comment":"The model that produces Table IV is not fully specified in this manuscript. Section III.B refers to the explicit D-matrix formulae as 'given elsewhere [6],' and Ref. [6] is in preparation, while the polarization data of Ref. [5] are also in preparation. A reader cannot reproduce the branching fractions or assess their model dependence. Please either include the relevant amplitude definitions and data set descriptions in an appendix or supplementary material, or state clearly which parts of the result are contingent on the unpublished analysis.","section":"Section III.B / Refs. [5,6]"},{"comment":"The acceptance correction for the event-based sample is computed with the GSIM detector simulation and the JM05 event generator, and Section III.C asserts that 'the details of the reaction model are irrelevant to the likelihood fit.' This is not demonstrated: if JM05 misrepresents the kinematic distributions in the 1.6–2.4 GeV region, the reconstructed-versus-generated weights in Eq. (5) are biased. The agreement among the four topologies (Fig. 7) is a necessary check but not sufficient, since all topologies share the same generator. Please quantify the acceptance-model dependence, for example by reweighting Monte Carlo events with the final BnGa solution and comparing acceptances to those obtained with JM05, and add the resulting uncertainty to Table II.","section":"Sections III.C and II.C"}],"minor_comments":[{"comment":"The abstract says branching ratios are obtained 'from an event based likelihood fit,' but Table IV is a product of the full BnGa coupled-channel fit that also includes the full data set and other channels; please clarify the specific role of the likelihood sample.","section":"Abstract / Section VI"},{"comment":"The meaning of σi as a differential cross section 'calculated for the reconstructed data events' normalized by the Monte Carlo sum needs a precise definition (binning, phase-space density, normalization) to be reproducible.","section":"Eq. (5)"},{"comment":"The caption ends with the stray string 'begindocument/before'; this appears to be a leftover from the manuscript preparation and should be removed.","section":"Figure 11 caption"},{"comment":"The name 'Kroll–Rudermann' is misspelled; the standard form is Kroll–Ruderman.","section":"Sections II.C and IV.C"},{"comment":"The relationship between the first line (this work, with uncertainties) and the second line (RPP ranges) should be stated explicitly in the caption; several rows have no RPP entry, which is not explained.","section":"Table IV"},{"comment":"The text alternates between 'data are' and 'data is' (e.g., 'This data is called the full data set'); please harmonize.","section":"Section III.C"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely of interest to the journal's audience, and the experimental material appears solid. The main issue is the gap between the promotional language ('uniquely determined') and the untested isobar-model assumption underlying Table IV. A revision that adds a model-dependence test, specifies the model sufficiently, and softens the uniqueness claim would be appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"One thing to know: this is a solid, careful paper whose cross sections and spin-density matrix elements are probably right, but the headline branching ratios in Table IV are BnGa model outputs, not unique data extractions, and the Summary oversells them.\n\nWhat's genuinely new: the first event-based likelihood fit to γp→π+π−p using four event topologies, with a ~1.8M-event subsample and 400M events as distributions. That yields new differential cross sections and SDMEs for pρ0, Δ++π−, Δ0π+, plus Nρ branching fractions for a long list of N* and Δ* resonances. The analysis is workmanlike: pulls are Gaussian, the four topologies agree after acceptance correction, background contributions are quantified (three-pion contamination 1–6%, CL<0.1 cut correction factor 1/0.72), and systematic uncertainties come from an honest spread of PWA solutions with varied weights and added resonances. The qualitative findings—ρ0 dominates at high W, Δ++ at low W through the Kroll–Ruderman mechanism, Δ0 suppressed—are credible and consistent with the earlier literature.\n\nSoft spots, in proportion. The main concern is the model dependence of the Nρ branching ratios. The fit assumes no coherent non-resonant Nππ amplitude; Section IV.C says so explicitly. The SAPHIR comparison in that same section shows that including a phase-space non-resonant term substantially lowers the Nρ excitation function. Since resonance couplings and a background amplitude can interfere and trade strength, the branching fractions are conditional on that assumption, and the paper runs no test. The Summary phrase \"uniquely determined by the new CLAS data\" is not supported; at best they are determined within the BnGa framework. Second, the BnGa formulas are deferred to Ref. [6] (in preparation), so the model is not checkable from this paper alone. Third, the acceptance correction uses the JM05 generator with GSIM; the authors argue the generator details don't matter, but they don't demonstrate that for the specific observables. These are standard PWA caveats; they don't make the results wrong, but they mean the branching ratios should not be quoted as data-driven values.\n\nWho benefits: hadron spectroscopists and anyone comparing quark-model or lattice decay predictions. The cross sections and SDMEs are useful measurements regardless of the PWA. The branching ratios need to be consumed with the model caveat in mind.\n\nRecommendation: yes, send to peer review. A serious referee should push for tempering the uniqueness language and either releasing the data/code or at least making the BnGa inputs and a non-resonant amplitude test available. The paper is worth publishing after that.","headline":"Solid CLAS cross sections and spin-density matrix elements; the Nρ branching fractions are BnGa-model outputs that the Summary overstates as uniquely determined by the data.","tokens_in":27080,"tokens_out":2739,"would_cite":false,"duration_ms":26333,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["13.60.Le","14.20.Gk","25.20.Lj"],"model":"deepseek-v4-flash","headline":"Event-based fit fixes N rho(770) branching fractions of baryon resonances.","keywords":["gamma p to pi+ pi- p","nucleon resonances","Delta resonances","rho(770) photoproduction","event-based likelihood fit","BnGa coupled-channel analysis","spin-density matrix elements","CLAS"],"falsifier":"A fit to the same CLAS data with an explicit non-resonant $N\\pi\\pi$ three-body amplitude added, or an independent partial-wave analysis of high-statistics data that does not impose the isobar decomposition, would settle the point: if the goodness of fit improves substantially and the extracted $N\\rho$ branching fractions in Table IV change by more than their quoted uncertainties, the central claim fails. A direct measurement of the $\\gamma p \\to \\pi^+\\pi^-p$ spin-density matrix elements over a wider angular range could also expose missing amplitudes.","tokens_in":25977,"feed_emoji":"⚛️","tokens_out":5072,"duration_ms":43050,"temperature":0.7,"pith_summary":"The paper reports a new measurement of the reaction $\\gamma p \\to \\pi^+\\pi^-p$ with the CLAS detector at Jefferson Lab, covering final-state invariant masses from 1.6 to 2.4 GeV. For the first time, this reaction is fitted event-by-event in five-dimensional phase space inside the Bonn-Gatchina coupled-channel analysis, rather than through binned histograms. The analysis produces total and differential cross sections and spin-density matrix elements for the isobars $p\\rho^0(770)$, $\\Delta(1232)^{++}\\pi^-$, and $\\Delta(1232)^{0}\\pi^+$. It also yields branching fractions for decays of most known $N^*$ and $\\Delta^*$ resonances into $N\\rho(770)$, values that the authors argue are uniquely determined by the new data. The result matters because these decay rates test the quark structure of baryon excitations and have previously been poorly known or contradictory.","feed_headline":"Event-based fit fixes N rho(770) branching fractions","feed_subtitle":"New CLAS data on gamma p -> pi+ pi- p determine N* and Delta* decay modes for the first time.","key_machinery":"The carrying tool is the Bonn-Gatchina coupled-channel partial-wave amplitude, fitted to the data with an event-based maximum-likelihood term. The amplitude uses a D-matrix based on dispersion relations with a one-step subtraction, K-matrix poles for resonances, and Blatt-Weisskopf form factors, and it describes the reaction as a sum of quasi-two-body isobar channels - $p\\rho^0(770)$, $\\Delta(1232)^{++}\\pi^-$, $\\Delta(1232)^{0}\\pi^+$, and smaller contributions - with interference between all amplitudes. Non-resonant $N\\pi\\pi$ production is deliberately not included. The event-based likelihood exploits all correlations in the five-dimensional phase space, whereas the full 400-million-event data set enters only as binned mass and angular distributions; the two are fitted jointly to the BnGa database of pion- and photo-induced reactions.","core_discovery":"The central claim is that an event-based likelihood fit of the CLAS g11a data on $\\gamma p \\to \\pi^+\\pi^-p$, combined with the BnGa coupled-channel amplitude, determines the $N\\rho(770)$ branching fractions of most known nucleon and $\\Delta$ resonances. The fit uses nearly two million events in four detection topologies treated together, so the five-dimensional correlations of the three-body final state are preserved. The resulting Table IV is the first set of $N\\rho$ branching ratios obtained from an event-based likelihood for this reaction. The authors find in particular that the cross section for $\\gamma p \\to N\\rho^0(770)$ is largely diffractive above $E_\\gamma \\sim 1.4$ GeV, while the strong $\\Delta(1232)^{++}\\pi^-$ production is driven by the Kroll-Ruderman mechanism, and that resonance contributions show up as deviations in the spin-density matrix elements and in backward-angle intensity.","pith_inferences":["If the isobar assumption is relaxed by adding an explicit non-resonant $N\\pi\\pi$ three-body amplitude, the $N\\rho$ branching fractions would likely shift because Table IV is obtained by integration over the $\\rho$ line shape; a direct three-body term could absorb some of the intensity currently assigned to $N\\rho$.","The method suggests a template for event-based partial-wave extraction in other reactions with large event samples, for example double-pion electroproduction, where the same BnGa machinery is applied to virtual photons.","The near-constant $N\\rho$ cross section above $E_\\gamma \\sim 1.3$ GeV, interpreted as diffractive, implies that the resonance contributions to $N\\rho$ are best isolated at backward angles or through the spin-density matrix elements; a dedicated backward-angle measurement with higher statistics could test the $N\\rho$ branching fractions without relying on the forward-dominated fit.","A testable extension is to compare the $N\\rho$ branching fractions in Table IV with values from an independent analysis of $\\pi^- p \\to \\pi^+\\pi^-n$ data in the same mass region, since the coupled-channel fit ties the photoproduction amplitudes to the pion-induced sector."],"forward_implications":["If correct, Table IV provides the first set of $N\\rho(770)$ branching fractions for resonances such as $N(1520)3/2^-$, $N(1675)5/2^-$, $N(1720)3/2^+$, $\\Delta(1620)1/2^-$, and $\\Delta(1920)3/2^+$, replacing or sharpening the wide ranges quoted in the Review of Particle Physics.","The measured total cross section and isobar excitation functions become a benchmark for future models of two-pion photoproduction in the resonance region.","The spin-density matrix elements for $\\rho^0(770)$ and $\\Delta(1232)^{++}$ decays provide new constraints on the exchange mechanisms (Pomeron exchange, pion exchange, Kroll-Ruderman) and on the size of resonance contributions.","The joint treatment of four detection topologies demonstrates that acceptance biases can be controlled at the level needed for event-based extraction, opening the same treatment for other three-body photoproduction channels.","The full binned data set and the event-based sample are not fully consistent at low photon energies, and the authors fold the spread of different fits into the systematic uncertainty rather than into the quoted central values."],"supporting_citations":[{"why":"Supplies the earlier CLAS analysis of $\\gamma p \\to \\pi^+\\pi^-p$, whose full data set and published cross sections this paper extends with the event-based fit.","marker":"[4]"},{"why":"New double-polarization data on $\\gamma p \\to \\pi^+\\pi^-p$ included in the BnGa database; the event-based sample is checked for consistency against them.","marker":"[5]"},{"why":"GSIM detector simulation used to compute the acceptance correction for the Monte Carlo events.","marker":"[70]"},{"why":"JM05 event generator used to simulate $\\gamma p \\to \\pi^+\\pi^-p$ events for acceptance and background studies.","marker":"[71]"},{"why":"JM05 reaction model from which the event generator is derived.","marker":"[72]"},{"why":"Defines how branching ratios of overlapping resonances are computed, the method behind Table IV.","marker":"[84]"},{"why":"Review of Particle Physics values used for comparison of the resulting branching ratios.","marker":"[85]"}],"fun_headline_variants":["First event-based N rho branching fractions from CLAS","Event-based fit extracts N* and Delta rho decay modes","Two-pion photoproduction yields N rho branching ratios","CLAS data determine N rho(770) branching fractions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The extraction rests on the isobar-model assumption that $\\gamma p \\to \\pi^+\\pi^-p$ is fully described by quasi-two-body channels ($p\\rho^0(770)$, $\\Delta(1232)\\pi$, plus a few smaller isobars) with no direct non-resonant $N\\pi\\pi$ production; if a sizable three-body amplitude exists, the $N\\rho(770)$ branching fractions and cross sections would shift.","fun_headline_variants_meta":{"raw":{"variants":["First event-based N rho branching fractions from CLAS","Event-based fit extracts N* and Delta rho decay modes","Two-pion photoproduction yields N rho branching ratios","CLAS data determine N rho(770) branching fractions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000315,"raw_usage":{"total_tokens":1717,"prompt_tokens":811,"completion_tokens":906,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":427,"completion_tokens_details":{"reasoning_tokens":845}},"tokens_in":427,"tokens_out":906,"duration_ms":8931,"temperature":1.0,"reasoning_tokens":845,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:19:00.211864+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A fit to the same CLAS data with an explicit non-resonant $N\\pi\\pi$ three-body amplitude added, or an independent partial-wave analysis of high-statistics data that does not impose the isobar decomposition, would settle the point: if the goodness of fit improves substantially and the extracted $N\\rho$ branching fractions in Table IV change by more than their quoted uncertainties, the central claim fails. A direct measurement of the $\\gamma p \\to \\pi^+\\pi^-p$ spin-density matrix elements over a wider angular range could also expose missing amplitudes.","supporting_citations":[{"cited_title":"CLAS GEANT Simulation,","cited_arxiv_id":null,"evidence_quote":"GSIM detector simulation used to compute the acceptance correction for the Monte Carlo events."},{"cited_title":"Cross Section of the reaction γp → pπ+π− from the G11A run period","cited_arxiv_id":null,"evidence_quote":"JM05 event generator used to simulate $\\gamma p \\to \\pi^+\\pi^-p$ events for acceptance and background studies."},{"cited_title":"A Phenomenological description of π−∆++ photo- and electroproduction in nucleon reso- nance region,","cited_arxiv_id":null,"evidence_quote":"JM05 reaction model from which the event generator is derived."},{"cited_title":"Photon dissociation model for vector meson photoproduction,","cited_arxiv_id":null,"evidence_quote":"Defines how branching ratios of overlapping resonances are computed, the method behind Table IV."},{"cited_title":"Elastic and proton dissociative ρ0 photoproduction at HERA,","cited_arxiv_id":null,"evidence_quote":"Review of Particle Physics values used for comparison of the resulting branching ratios."}],"review_version":1}