{"id":"9d7233af-af98-4394-8af2-2052ab5aafa7","arxiv_id":"2511.12007","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Applying the Hofstadter-Herman plot to 1994 SLAC data suggests the 4% normalization applied to the 1.6 GeV spectrometer subset moves the data away from modern global fits.","lead":"This paper re-plots 1994 SLAC proton-scattering data with a 1960 visualization method to compare two versions of the data: one with the original 4% normalization correction and one without. The plots suggest the uncorrected data agree better with modern JLab global fits, so the 4% correction may be wrong.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 4%-normalization conclusion rests on an unquantified, possibly circular comparison to the JLab global fit; whether unnormalized SLAC data are preferred is not yet demonstrated.","rationale":"The reader's weakest-assumption analysis matches my own: the physical conclusion is benchmark-dependent and visually asserted. I add one concrete observation from the manuscript itself: Tables I and II show the 1.6 and 8 GeV spectrometers disagree by ~5% at identical kinematics (Q2=1.75 GeV), and the 4% factor was the original cross-calibration. Thus 'uncorrected data' are not a neutral alternative; preferring them declares the 8 GeV normalization wrong, which requires more than a visual curve comparison. The paper's methodological contribution—reproducing the SLAC plot and demonstrating that HH bands expose low-epsilon normalization effects—remains valuable and independently checkable. The correct disposition is unchanged: conditional on adding a quantitative chi-square/pull comparison and demonstrating that the chosen global fit is independent (or re-fitting with SLAC excluded), the physical claim would be supported. I therefore leave the reader's CONDITIONAL verdict in place.","tokens_in":8501,"tokens_out":8301,"duration_ms":78282,"concrete_test":"Take the JLab global fit from ref. [13] App. A (including its parameter covariance). Fix the form factors to this fit and predict dσ/dΩ at each SLAC kinematics. Compute χ² for the full dataset with the 1.6 GeV data normalized (×0.96) and unnormalized (×1.00), using the published total errors and accounting for normalization correlations. Also re-fit with all 1994 SLAC data excluded from the global fit inputs (checking ref. [13]'s data list) and repeat. If Δχ² between normalized and unnormalized is <4, or disappears when SLAC is excluded from the fit, the conclusion that the 4% correction should be removed is not supported; if Δχ² clearly favors unnormalized and survives SLAC exclusion, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's physical conclusion (Section IV, Figs. 5–6) is that the 1994 SLAC 1.6 GeV data should not have been normalized by 4%. This rests wholly on the assertion that the unnormalized data 'align more closely' with the JLab Global Fit of ref. [13, App. A]. That comparison is load-bearing but not demonstrative: (1) no chi-square, pull, or other quantitative agreement statistic is computed; in (G_E^2, G_M^2) space the allowed bands are long, steep strips (slope ~ -tau/epsilon), so a visual impression of closeness can be dominated by the projection rather than by actual statistical preference; (2) the fit curve is plotted without its uncertainty, so it cannot be determined whether both the normalized and unnormalized cases are equally compatible once fit errors are included; (3) independence is not established—if the global fit includes the normalized Andivahis cross sections, the comparison is at least partly circular and cannot adjudicate the normalization; (4) the 4% factor was not an arbitrary correction: at overlapping kinematics (e.g., Q2=1.75 GeV in Tables I and II), the 1.6 GeV cross section is ~5% higher than the 8 GeV value, and the original analysis applied the normalization to reconcile the two spectrometers. Preferring the unnormalized data therefore means implicitly declaring the 8 GeV absolute normalization wrong, which is not justified anywhere except by the same unquantified global-fit comparison. The visualization method itself is unaffected, but the specific physical claim is unsupported without a quantitative, independent benchmark.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper revisits the classic 1994 SLAC elastic electron-proton scattering data of Andivahis et al. and compares the conventional Rosenbluth separation visualization with the Hofstadter-Herman (HH) method, which plots allowed bands in (G_E^2, G_M^2) space. The authors reproduce the original SLAC reduced cross-section plot, generate HH plots with and without the 4% normalization correction applied to the 1.6 GeV spectrometer data, and visually compare them with the Jefferson Lab Global Fit of Ref. [13]. From this visual comparison they conclude that the unnormalized data align more closely with the global fit, suggesting that the original 4% normalization shifted the measurements away from the region favored by modern extractions. They recommend the HH visualization as a routine diagnostic cross-check at future facilities such as the EIC.","tokens_in":8913,"tokens_out":5104,"duration_ms":46816,"significance":"If the central physical claim were quantitatively established, the paper would be a useful cautionary note for global form-factor analyses and for normalization handling at future facilities. The paper's pedagogical value is real: it demonstrates a compact visualization that can make systematic tensions easier to perceive, and it includes tabulated data and reproduces the original SLAC plot. However, the conclusion about the 4% normalization is currently supported only by an unquantified visual impression, and the external benchmark may be correlated with the very data in question. The visualization method itself is an algebraic rearrangement of the same cross-section formula and adds no new statistical information; its value is heuristic. The significance is therefore conditional on a quantitative, non-circular demonstration, which the present manuscript does not provide.","major_comments":[{"comment":"The central conclusion that the unnormalized SLAC data 'more accurately reflect the underlying form factor behavior' rests entirely on visual comparison of the HH bands to the JLab Global Fit curves. No chi-square, pull, or other quantitative agreement measure is computed. The global fit curves are drawn without uncertainty bands, so it is unknown whether both the normalized and unnormalized datasets are statistically compatible with the fit once fit errors are included. Moreover, the bands in (G_E^2, G_M^2) space are long, steep strips with slope approximately -tau/epsilon (see Eq. 9); visual proximity can be dominated by this projection rather than by genuine statistical preference. The paper should compute a quantitative compatibility metric, e.g., chi-square or pull with full uncertainties, for both datasets relative to the global fit.","section":"Section IV, Figs. 5-6"},{"comment":"The comparison to the JLab Global Fit is used as an external benchmark, but the paper never establishes that this fit is independent of the SLAC Andivahis dataset. If the global fit in Ref. [13] includes the 4%-normalized SLAC cross sections, then the comparison is at least partly circular: the fit has already absorbed the normalization, and using it to adjudicate the normalization is not valid. If the fit does not include the SLAC data, that should be stated explicitly. This point is load-bearing because the entire physical conclusion depends on the benchmark being independent; without this assurance, the conclusion collapses.","section":"Section IV, Figs. 5-6 and Ref. [13]"},{"comment":"The original 4% normalization was not an arbitrary correction: at overlapping kinematics, e.g., Q^2 = 1.75 GeV^2, the 1.6 GeV spectrometer cross section (1.514 x 10^-1 in Table II) is about 5% higher than the 8 GeV value (1.440 x 10^-1 in Table I). The original analysis applied the normalization to reconcile the two spectrometers. Preferring the unnormalized data therefore implicitly declares the 8 GeV absolute normalization to be wrong. The paper provides no independent evidence for this, for example from a relative normalization check, a radiative-correction check, or a comparison to a truly independent dataset. This alternative interpretation must be addressed before the claim is accepted.","section":"Section III, Tables I and II, and Section IV"}],"minor_comments":[{"comment":"The manuscript would benefit from a data/code availability statement. Since the analysis is purely computational and the authors mention using Python with NumPy and Matplotlib, making the code available would greatly increase reproducibility and would facilitate the quantitative checks recommended above.","section":"General"},{"comment":"The notation in Eq. (2) is confusing: the left-hand side is written as sigma_R/G_D^2(Q^2) while the right-hand side includes the factor (1+Q^2/0.71)^4. The numerical constants (5.18, 0.71) and their units should be defined explicitly. Also, the phrase 'with numerical factor 0.71 from fitting to existing 1994 data' is imprecise; the dipole mass parameter is 0.71 GeV^2.","section":"Eq. (2)"},{"comment":"The table headings and beam-energy columns are confusing. Table I is labeled '8 GeV spectrometer' but contains rows with E = 9.800 GeV, and Table II is labeled '1.6 GeV spectrometer' but contains rows with E up to 5.507 GeV. Please clarify the spectrometer naming convention (e.g., maximum scattered momentum versus incident beam energy) or fix the table entries.","section":"Tables I and II"},{"comment":"The legends in Figs. 5 and 6 are redundant and confusing: 'Normalized Q2' appears in the legend box as well as in the axis label, and 'JLab Global Fit 1' and 'JLab Global Fit 2' are listed twice. Please clean up the legend formatting.","section":"Figs. 5-6"},{"comment":"The claim that the HH method 'reveals previously obscured regions of form factor parameter space' is overstrong. Since Eq. (8) is an algebraic rearrangement of the same cross-section formula, the method cannot reveal information not already present in the data; it can only present it in a different geometric form. A wording such as 'makes more apparent' would be more accurate.","section":"Abstract and Section IV"},{"comment":"There are several typographical issues, e.g., 'elasticepscattering' in the introduction and 'T A' in the author affiliation line. A careful proofreading pass is needed.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is a short, clearly written diagnostic note, but the central physical claim is not currently supported. The lack of any quantitative comparison and the unresolved question of the global fit's independence from the SLAC data are substantial. I would encourage the editor to require the authors to address these points before further consideration. The visual method has merit as a pedagogical tool, but the paper should be framed accordingly unless the quantitative analysis is added."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe short version: the paper shows that the Hofstadter-Herman band plot, applied to the old SLAC data, makes a normalization tension visible that Rosenbluth plots hide. That is a genuinely useful demonstration. But the physical takeaway — that the 4% normalization applied to the 1.6 GeV spectrometer data was wrong — rests on a purely visual comparison to a global fit whose relation to the SLAC data is never established. That part is not demonstrated.\n\nWhat is new: taking the 1960 band-plot method into Sachs space and applying it to Andivahis et al. works. At Q2 = 1.75 GeV, the Rosenbluth plot looks like a clean line, while the band plot shows the 1.6 GeV and 8 GeV data intersecting in a way that shifts visibly under the 4% correction. The paper also reproduces the original SLAC plot with high fidelity, which is a good sanity check. As a diagnostic cross-check for future experiments, the tool is worth having.\n\nThe soft spots are in the interpretation. First, the comparison to the JLab Global Fit is purely visual — no chi-square, no pull, no uncertainty band on the fit curve. In (GE^2, GM^2) space the allowed bands are long, steep strips, so “closer” can be dominated by projection rather than statistical preference. Second, and more important, the 4% factor was not arbitrary. At overlapping kinematics (e.g., Q2 = 1.75 GeV), the 1.6 GeV cross section is about 5% higher than the 8 GeV value. The original normalization was how the authors reconciled the two spectrometers. Preferring the unnormalized data therefore means declaring the 8 GeV absolute normalization wrong by roughly that amount. That is a stronger claim than the paper makes, and the paper offers no independent justification for it. Third, we never learn whether the global fit includes the normalized SLAC cross sections. If it does, the comparison is at least partly circular.\n\nSo: the tool is good and should be seen. The specific 4% conclusion needs either a quantitative treatment — chi-square with fit uncertainties, or a fit that explicitly excludes the SLAC points — or a softer framing: “this discrepancy merits checking,” not “the uncorrected data better reflect the underlying physics.”\n\nWho gets value from this: anyone analyzing cross-section data for form factors, and people thinking about systematic checks for EIC-era experiments. It is a short, readable methods paper with an honest core. I would send it to peer review, but the authors should expect a request to either quantify the comparison or reframe the conclusion.\n\nRecommendation: send to peer review, with the expectation of revisions focusing on the strength of the physical claim.","headline":"Useful diagnostic demo, but the 4% normalization claim needs quantitative support.","tokens_in":9393,"tokens_out":2566,"would_cite":true,"duration_ms":26104,"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":"Replotting 1994 proton-scattering data in the Hofstadter-Herman plane suggests the original 4% normalization correction moved measurements away from modern global-fit values, so the uncorrected points may be closer to the true form factors.","keywords":["proton electromagnetic form factors","Rosenbluth separation","Hofstadter-Herman visualization","elastic electron-proton scattering","normalization systematics","form-factor extraction","global fit comparison","1994 electron-scattering data"],"falsifier":"A decisive test would be a two-version refit of the world elastic-scattering data, once with the 1.6 GeV spectrometer points scaled by 4% and once unscaled, using a global fit constructed without those points; if the scaled version yields an equal or better chi-square per point and the unscaled band intersections fall within the global-fit uncertainty contours, the paper's suggestion is falsified. Even before refitting, plotting the global fit's uncertainty band inside the Hofstadter-Herman panels would show whether the normalized data's intersections actually lie outside it at any Q^2.","tokens_in":8436,"feed_emoji":"⚛️","tokens_out":7430,"duration_ms":62506,"temperature":0.7,"pith_summary":"The paper revisits the 1994 Rosenbluth-separation dataset for elastic electron-proton scattering and re-expresses each measurement as a band in (G_E^2, G_M^2) space, the Hofstadter-Herman representation. It claims that when multiple measurements at the same Q^2 are overlapped, the band intersections show directly where the consistent form-factor pair lies, and that this geometry exposes a systematic problem invisible in standard Rosenbluth plots: the 4% normalization applied in the original analysis to the 1.6 GeV spectrometer subset pulls the data away from the region favored by a modern global fit. On this visual evidence, the paper concludes that the uncorrected measurements may reflect the underlying form-factor behavior better than the published corrected set. The practical point is diagnostic: a cheap geometrical cross-check can flag normalization and calibration shifts before they are absorbed into global fits, which matters for forthcoming precision measurements.","feed_headline":"A 4% correction to 1994 proton data may point the wrong way","feed_subtitle":"Hofstadter-Herman bands show corrected data missing the global-fit region that line fits hid.","key_machinery":"The Hofstadter-Herman visualization in the Sachs basis. Starting from the reduced cross section, the paper rewrites the relation as G_E^2 = sigma_R(1+tau) - (tau/epsilon) G_M^2. For a fixed Q^2 and a single measurement, scanning G_M^2 traces a straight band in (G_E^2, G_M^2) whose width is set by the cross-section uncertainty; measurements taken at different beam energies and angles give bands of different slopes, and their overlap region identifies the consistent form-factor pair. This replaces the conventional two-step Rosenbluth extraction—fit a line to reduced cross sections, then read slope and intercept—with a direct geometric intersection, so systematic shifts such as a 4% normalizati","core_discovery":"The central claim is a diagnostic demonstration rather than a new measurement: in the 1994 elastic-scattering dataset, the 4% normalization correction applied to the 1.6 GeV spectrometer data is not supported by a modern global fit. In the conventional reduced-cross-section presentation, normalized and unnormalized data are nearly indistinguishable. In the Hofstadter-Herman plane, however, the two treatments produce visibly different band intersections, and the global-fit curve passes through the unnormalized intersections at several Q^2 values while missing the normalized ones. The paper takes this as evidence that the original normalization step shifted the data away from, rather than towa","pith_inferences":["The visual comparison could be converted into a quantitative diagnostic: compute the distance between the normalized and unnormalized band intersections and the global-fit curve, with the fit's own uncertainty included, turning 'closer' into a number.","The same band-intersection logic is transferable to any two-parameter extraction in which a measured observable depends linearly on two unknowns, such as separating longitudinal and transverse structure functions or neutron form factors from quasielastic data.","A stronger test of the paper's suggestion would refit the world elastic-scattering data with the 1.6 GeV subset scaled and unscaled, using a fit that is constructed without those points, and compare goodness of fit—this would separate the diagnostic claim from the choice of benchmark.","If the conclusion is later confirmed quantitatively, historical normalization procedures from other 1960s-1990s electron-scattering measurements may deserve similar reexamination, which could slightly shift averaged form-factor and proton-radius values."],"forward_implications":["If the claim is right, the 4% normalization applied in the original analysis should be revisited, since it may have pushed the 1994 data away from modern global-fit values.","The Hofstadter-Herman plot is a low-cost per-Q^2 diagnostic that can flag normalization and calibration shifts before they propagate into global fits.","Showing data as overlap bands lets a reviewer see at a glance whether a fit curve is consistent with the raw measurements at each Q^2, complementing chi-square summaries.","Future electron-ion collider elastic-scattering measurements can use this geometric check to decide, at the data-analysis stage, whether a normalization adjustment is actually justified.","If the uncorrected points are indeed closer to the true form factors, global fits that include the published corrected values carry a small systematic bias that a re-fit with either treatment could quantify."],"fun_headline_variants":["1994 proton data correction may be off target","New plot method questions 4% correction to old proton data","Hofstadter-Herman plot flags wrong-way correction in 1994 data","Visualization exposes questionable normalization in classic proton data","Did a 1994 proton data fix point the wrong way?"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The argument rests on treating the modern global fit as a trustworthy benchmark for G_E^2 and G_M^2 at these Q^2 values; the fit's uncertainty is not shown and it is not established whether the fit is independent of the (normalized) 1994 data, and the conclusion that the unnormalized data are more accurate is based purely on visual band-curve comparisons without a quantitative agreement measure.","fun_headline_variants_meta":{"raw":{"variants":["1994 proton data correction may be off target","New plot method questions 4% correction to old proton data","Hofstadter-Herman plot flags wrong-way correction in 1994 data","Visualization exposes questionable normalization in classic proton data","Did a 1994 proton data fix point the wrong way?"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000581,"raw_usage":{"total_tokens":2564,"prompt_tokens":725,"completion_tokens":1839,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":469,"completion_tokens_details":{"reasoning_tokens":1765}},"tokens_in":469,"tokens_out":1839,"duration_ms":10933,"temperature":1.0,"reasoning_tokens":1765,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T22:06:01.779216+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test would be a two-version refit of the world elastic-scattering data, once with the 1.6 GeV spectrometer points scaled by 4% and once unscaled, using a global fit constructed without those points; if the scaled version yields an equal or better chi-square per point and the unscaled band intersections fall within the global-fit uncertainty contours, the paper's suggestion is falsified. Even before refitting, plotting the global fit's uncertainty band inside the Hofstadter-Herman panels would show whether the normalized data's intersections actually lie outside it at any Q^2.","supporting_citations":[],"review_version":1}