REVIEW 3 major objections 4 minor 66 references
First Measurement of the Charged Current $\overline{\nu}_{\mu}$ Double Differential Cross Section on a Water Target without Pions in the final state
T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read The T2K Collaboration reports the first measurement of the charged-current antineutrino-muon double differential cross section on water with no pions in the final state, with an integrated cross section of (1.11 ± 0.18) × 10^-38 cm² per…
desk verdict A solid first measurement of antineutrino CC0π on water; the unfolding's generator dependence is real but disclosed, and the 'model independent' phrase in Sec. V.D should be toned down before publication. 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 smearing (migration) matrix $S_{ij}$ combined with two independent sets of scale parameters, $c_j$ for water-target signal and $d_j$ for non-water signal, in the likelihood of Eq. (11). The water-out sample constrains $d_j$, so the fitted $c_j$ directly yield the true CC0$\pi$ rate on water in each $p$–$\cos\theta$ bin; a penalty term regularizes neighbouring momentum bins and an L-curve chooses its strength.
What would settle it
Repeat the unfolding with a different event generator in place of the nominal one and see if the 19 unfolded cross-section bins move by more than the quoted uncertainties; if they do, the generator-shape assumption behind the measurement is wrong. A dedicated high-statistics water-in/water-out cycling run would also directly test the linear-subtraction premise by checking whether the fitted non-water scale parameters stay at unity.
Extended reading notes
Core claim
The central claim is that the flux-averaged $\overline{\nu}_\mu$ charged-current no-pion (CC0$\pi$) cross section on oxygen in water, measured in 19 true bins of muon momentum and $\cos\theta$, is consistent at the 1$\sigma$ level with the NEUT, GENIE, and NuWro model predictions over most of the phase space, with three bins at momenta around 0.67–1.0 GeV/c sitting about 2$\sigma$ below NEUT. The integrated cross section over the measured phase space is $\sigma=(1.11\pm0.18)\times10^{-38}\,\mathrm{cm^2}$ per water molecule in the regularized fit. The result is unfolded from detector smearing with a binned-likelihood fit in which water and non-water target event rates are scaled by independent parameters, so the extracted water signal is separated from interactions in surrounding material.
Load-bearing premise
The unfolding assumes the simulated event generator correctly describes the shape of true events inside each momentum-angle bin, and that subtracting the water-out sample linearly removes interactions in non-water material.
Editorial extensions
If this is right
- Antineutrino oscillation analyses at T2K can now use a data-driven CC0$\pi$ cross section on water, reducing the model uncertainty they inherit from neutrino-nucleus scattering.
- The 19-bin double-differential result and its covariance give generator tunes a concrete target: three bins near 0.67–1.0 GeV/c sit about 2$\sigma$ below NEUT, defining where models need adjustment.
- Regularized and unregularized results agree, so the measurement is stable against the choice of unfolding regularization.
- The flux-averaged integrated cross section of $(1.11 \pm 0.18)\times10^{-38}$ cm$^2$ per water molecule provides a normalization benchmark for future water-Cherenkov detectors.
Reading between the lines
- If the three low bins reflect a genuine oxygen nuclear effect rather than an unfolding artifact, retuning generators to reproduce them would change reconstructed antineutrino energies and could shift CP-violation sensitivity at long-baseline experiments.
- A natural extension is to apply the same water-in/water-out subtraction to CC1$\pi$ and neutral-current samples, giving a fuller map of nuclear effects on oxygen with the same detector systematics.
- With more exposure, the flux-averaged result could be re-binned in reconstructed neutrino energy, turning the double-differential cross section into a flux-unfolded measurement that model comparisons would read more directly.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents the first measurement of the charged-current antineutrino-muon cross section on a water target with no pions in the final state (CC-0π), using the T2K off-axis near detector ND280. The analysis selects CC-0π events in the P0D water target, using both water-in and water-out data samples to separate water from non-water contributions. A binned maximum-likelihood fit with penalty terms for flux, detector, and background/FSI systematics is used to unfold the observed reconstructed muon momentum and angle distributions into true p-cosθ bins. The smearing matrix and signal templates are taken from the NEUT Monte Carlo generator. The paper reports an unfolded flux-averaged double-differential cross section in 19 bins, an integrated cross section of σ = (1.11 ± 0.18) × 10^-38 cm² per water molecule in the restricted phase space, and comparisons with NEUT, GENIE, and NuWro predictions. Both regularized and unregularized results are provided, together with a data release containing central values and covariance matrices.
Significance. If the result holds, this is the first differential antineutrino cross-section measurement on a water target and provides a new datum for nuclear-model validation in the energy region relevant to T2K and future long-baseline experiments. The paper has several concrete strengths: the water-in/water-out subtraction is handled by a simultaneous fit with independent scale parameters; the analysis provides both regularized and unregularized unfolded results with full covariance matrices in a public data release; and the fit is validated with Monte Carlo closure tests, including mass-scale checks and fits to pseudo-experiments. The main caveat is that the unfolding is anchored to the NEUT migration matrix and signal shape, so the numerical cross-section values inherit a generator dependence that is discussed but not fully quantified. This is a standard limitation for this class of measurements and does not by itself undermine the central result, but it should be stated more carefully and, if possible, quantified.
major comments (3)
- [Sec. V.D, after Eq. (14)] The sentence immediately after Eq. (14) states that the unregularized results are "fully correct and model independent," but the fitted signal in Eq. (11) is constructed from the NEUT-based smearing matrix S_ij and the NEUT signal templates N_sig,water,MC_j. The unfolded bin contents therefore inherit NEUT's migration and efficiency corrections even in the unregularized case. I recommend either rephrasing this claim to something like "independent of the regularization choice" or adding a quantitative cross-check, for example by re-unfolding with GENIE or NuWro smearing matrices, so that the generator dependence of the central values is estimated rather than asserted away.
- [Sec. VI.B, Table VI and Figs. 9-10] The chi-squared values for NEUT are inconsistent between Table VI and the figure captions. Table VI quotes 29.2 (regularized) and 33.1 (unregularized) for the comparison of data to NEUT, while the captions of Fig. 9 and Fig. 10 quote 22.22 and 25.12 for what is described as the same quantity defined by Eq. (13). Because the model-comparison discussion draws conclusions from these numbers, the discrepancy must be traced and corrected in one place.
- [Sec. V.B and Sec. VI.A] The fit is performed in all 28 true bins, but only 19 bins are used in the final cross-section result. The nine excluded bins still have fitted scale factors c_j that enter the smeared predictions for neighboring reconstructed bins, so a mismodeled signal in those bins could in principle distort the extracted c_j for the reported bins. I ask for an explicit closure test in which the signal is deliberately mis-modeled in the excluded bins while the data are fitted, to show that the 19 reported values are unbiased, or for a short discussion of why the excluded bins' contributions to the used bins are negligible.
minor comments (4)
- [Eq. (8), Sec. V.A] The notation f_i^n is described as "the fraction of antineutrinos in flux energy bin n for reconstructed bin i," but a flux parameter should be an energy-bin scale factor common to all reconstructed bins; please clarify the intended index structure and how the 11 flux parameters enter Eq. (11).
- [Title and Abstract] The title and abstract use "νμ" rather than "ν̄μ" in the phrase "charged-current νμ double differential cross section," although the paper measures antineutrinos; please make the overbar explicit for consistency with the introduction and the body of the paper.
- [Sec. VI.B, Eq. (13)] The quoted chi-squared values are computed with a post-fit covariance matrix and therefore are not a standard goodness-of-fit statistic; the current wording implicitly treats them as such. This should be stated explicitly, both in the text and in the Table VI caption.
- [Sec. VII] The statement that NuWro has the lowest chi-squared should be accompanied by a caveat that the differences between models are not assigned a statistical significance, given the non-standard chi-squared definition and the use of the post-fit covariance.
Circularity Check
No circularity: the cross section is extracted from data by fitting free per-bin scale factors; NEUT provides only the unfolding response, which is disclosed model dependence rather than a circular input.
full rationale
The paper is an experimental measurement, not a derivation of a predicted quantity from its inputs. The extracted double-differential cross section is obtained by fitting per-bin scale factors c_j and d_j in Eq. (11) to water-in and water-out data with a binned likelihood (Eqs. 6-12), then converting the fitted c_j N_sig,water,MC_j terms to a flux-averaged cross section. The NEUT generator provides the smearing matrix S_ij and the signal/background template shapes, which is a disclosed response correction and the usual model dependence of an unfolding analysis, not a circular reduction: the 19 fitted coefficients are free parameters determined by the data, and the paper also reports both unregularized and regularized results and compares them to three independent generators (NEUT, GENIE, NuWro). References to previous T2K analyses ([32], [41], [47], [48]) supply technical methods and prior parameter values, but they do not define the measured cross section or forbid alternatives; no uniqueness theorem or ansatz is imported from the authors' prior work. The statement in Sec. V.D that unregularized results are 'fully correct and model independent' overstates the residual NEUT dependence of the migration matrix, but that is a correctness and model-dependence caveat, not a circular step under the definitions of this review. No self-definitional, fitted-input-called-prediction, self-citation-load-bearing, uniqueness-imported, ansatz-smuggled, or renaming pattern is present.
Assumptions & free parameters
free parameters (6)
- Signal scale factors on water, c_j =
Post-fit values near 1.0, with three bins (6, 7, 11) roughly 0.6-0.8
- Signal scale factors on non-water, d_j =
Post-fit values centered near 0.9
- Flux parameters, f_i (11 energy bins) =
Constrained by prior covariance V_flux_cov; typical uncertainties around 10% near the peak
- Detector response parameters, r_i (76 parameters) =
Nominally 1.0 with priors from detector uncertainties
- Background and FSI model parameters, a (15 parameters) =
Priors from [41], [47], [48], including axial mass, nonresonant background fraction, DIS normalization, and pion FSI
- Regularization strength, p_reg =
1 (chosen from L-curve largest curvature)
assumptions (6)
- domain assumption NEUT 5.3.3 correctly models the true p-cosθ distribution of CC0π events within each kinematic bin.
- domain assumption The predicted antineutrino flux from FLUKA2011, GEANT3, and GCALOR, reweighted with NA61/SHINE thin-target data, is accurate to the quoted uncertainties.
- domain assumption Particle identification using TPC dE/dx and P0D track energy loss correctly separates muons from pions and protons.
- domain assumption Water-in versus water-out subtraction isolates the water target contribution by linear scaling of non-water target material.
- standard math The binned likelihood, chi-square, and L-curve regularization are standard statistical tools that need no new derivation.
- domain assumption Excluding 9 of 28 bins with low efficiency or low statistics does not bias the reported restricted-phase-space cross section.
Cite this review
Pith. "Pith review of First Measurement of the Charged Current $\overline{\nu}_{\mu}$ Double Differential Cross Section on a Water Target without Pions in the final state." pith.science (2026). https://pith.science/paper/TAF4WSZZ
@misc{pith2026190810249,
author = {Pith},
title = {Pith review of: First Measurement of the Charged Current $\overline\nu_\mu$ Double Differential Cross Section on a Water Target without Pions in the final state},
year = {2026},
howpublished = {\url{https://pith.science/paper/TAF4WSZZ}},
note = {Machine review of arXiv:1908.10249}
}
abstract
This paper reports the first differential measurement of the charged-current $\overline{\nu}_{\mu}$ interaction cross section on water with no pions in the final state. The unfolded flux-averaged measurement using the T2K off-axis near detector is given in double differential bins of $\mu^+$ momentum and angle. The integrated cross section in a restricted phase space is $\sigma=\left(1.11\pm0.18\right)\times10^{-38}$ cm$^{2}$ per water molecule. Comparisons with several nuclear models are also presented.
Figures
Figures from the paper (6 more)
Reference graph
Works this paper leans on
-
[28]
R. Gran, J. Nieves, F. Sanchez, and M. V. Vacas, Phys. Rev. D88, 113007 (2011)
work page 2011
-
[1]
Primary ¯νµ and background νµ, νe, and ¯νe beam production in the graphite target and propagation through the following horns and decay volume. The hadronic rates from the beam target were gen- erated by FLUKA2011 which was tuned to the NA61/SHINE measurements and the GEANT3 simulation software predicted the flux and energy spectrum for the different neutri...
-
[2]
The antineutrino and neutrino interactions in the ND280 detector, where the NEUT [44] MC gen- erator (version 5.3.3) is used to calculate the in- teraction cross sections and the final state particle kinematics
-
[3]
The detector response, which used the GEANT4
-
[4]
These first four requirements define the CC-Inc event selection
There is a muon track candidate that is the highest momentum positively charged track, the highest momentum track in the event, and has a TPC dE/dx track measurement consistent with a muon energy loss. These first four requirements define the CC-Inc event selection
-
[5]
and HyperKamiokande [6]. The main difficulty in the description of (anti)neutrino-nucleus interactions derives from the intrinsic nature of the nucleus, where nucleons are bound together and nuclear effects must be taken into account. Many models are currently avail- able, describing different pieces of this complex sce- nario such as relativistic Fermi gas [...
work page 1900
-
[6]
Overall ND280 data quality flags are good such that the detector was operational and stable dur- ing taking data
-
[7]
There is a reconstructed track in the P∅D match- ing a track in the TPC with the start of the track reconstructed in the fiducial volume of the P ∅D water target
Show all 66 references
-
[8]
There is at least one track reconstructed in TPC1
-
[9]
A. M. Ankowski, O. Benhar, and N. Farina, Phys. Rev. D82, 013002 (2010), arXiv:1001.0481 [nucl-th]
2010 arXiv
-
[10]
This cut removes charged current events with a π0
There are no reconstructed P ∅D showers in the event. This cut removes charged current events with a π0
-
[11]
smearing ma- trix
Remaining events are then separated into 3 cate- gories based on the number of µ-like P∅D tracks in the event. (a) Events with only a muon track candidate de- fine the CC-0π selection. (b) Events with a muon track candidate and one µ-like track define the CC-1π selection. (c) Al...
2000
-
[12]
This improves the stability of the fit results
The number of events in each 2-D bin should have reasonable statistics,∼100 events. This improves the stability of the fit results
-
[13]
The selection efficiency should be relatively high to minimize model dependence of the efficiency correction, and event populations should not dif- fer very much between adjacent bins which also improves the stability of the fit results
-
[14]
We expect these choices should reduce regularization complications, which are discussed in later sections, or possibly even the need for regularization
The bin sizes should be fine enough so local detec- tor resolution effects are well represented and the detector resolutions do not change too much from bin to bin, however not too fine such that there are too few events in the bin. We expect these choices should reduce regulariz...
2010
-
[15]
Abe et al
K. Abe et al. (T2K Collaboration), Phys. Rev. Lett. 121, 171802 (2018)
2018
-
[16]
Adamson et al
P. Adamson et al. (NOvA Collaboration), Phys. Rev. Lett. 118, 151802 (2017)
2017
-
[17]
Abe et al
K. Abe et al. (T2K Collaboration), Nucl. Instrum. Meth. A659, 106 (2011), arXiv:1106.1238 [physics.ins- det]
2011
-
[18]
Adamson et al
P. Adamson et al. (NOvA Collaboration), Phys. Rev. Lett. 116, 151806 (2016)
2016
-
[19]
Acciarri et al
R. Acciarri et al. (DUNE Collaboration), (2015), arXiv:1512.06148 [physics.ins-det]
2015
-
[20]
Abe et al
K. Abe et al. (Hyper-Kamiokande Working Group) (2014) arXiv:1412.4673 [physics.ins-det]
2014 arXiv
-
[21]
R. A. Smith and E. J. Moniz, Nucl. Phys. B43, 605 (1972), B101, 547(E) (1975)
1972
-
[22]
Benhar, A
O. Benhar, A. Fabrocini, S. Fantoni, and I. Sick, Nucl. Phys. A579, 193 (1994)
1994
-
[23]
Gallmeister, U
K. Gallmeister, U. Mosel, and J. Weil, Phys. Rev. C94, 035502 (2016), arXiv:1605.09391 [nucl-th]
2016 arXiv
-
[24]
S. K. Singh and E. Oset, Nucl. Phys. A542, 587 (1992)
1992
-
[25]
A. Gil, J. Nieves, and E. Oset, Nucl. Phys. A627, 543 (1997), arXiv:nucl-th/9711009 [nucl-th]
1997 arXiv
-
[26]
Nieves, J
J. Nieves, J. E. Amaro, and M. Valverde, Phys. Rev. C70, 055503 (2004), [Erratum: Phys. Rev.C72,019902(2005)], arXiv:nucl-th/0408005 [nucl- th]
2004 arXiv
-
[27]
Nieves, M
J. Nieves, M. Valverde, and M. J. Vicente Vacas, Phys. Rev. C73, 025504 (2006), arXiv:hep-ph/0511204 [hep- ph]
2006 arXiv
-
[29]
Martini, M
M. Martini, M. Ericson, G. Chanfray, and J. Marteau, Phys. Rev. C80, 065501 (2009), arXiv:0910.2622 [nucl- th]
2009 arXiv
-
[30]
Delorme and M
J. Delorme and M. Ericson, Phys. Lett. 156B, 263 (1985)
1985
-
[31]
Marteau, J
J. Marteau, J. Delorme, and M. Ericson, Neutrino factories based on muon storage rings. Proceedings, ICFA/ECFA Workshop, nuFACT’99, Lyon, France, July 5-9, 1999 , Nucl. Instrum. Meth. A451, 76 (2000)
2000
-
[32]
Martini, M
M. Martini, M. Ericson, G. Chanfray, and J. Marteau, Phys. Rev. C81, 045502 (2010), arXiv:1002.4538 [hep- ph]
2010 arXiv
-
[33]
Nieves, I
J. Nieves, I. Ruiz Simo, and M. J. Vicente Vacas, Phys. Lett. B707, 72 (2012), arXiv:1106.5374 [hep-ph]
2012 arXiv
-
[34]
Martini, M
M. Martini, M. Ericson, and G. Chanfray, Phys. Rev. C84, 055502 (2011), arXiv:1110.0221 [nucl-th]
2011 arXiv
-
[35]
Nieves, I
J. Nieves, I. R. Simo, and M. V. Vacas, Phys. Lett. B 707, 72 (2012)
2012
-
[36]
O. Buss, T. Gaitanos, K. Gallmeister, H. van Hees, M. Kaskulov, O. Lalakulich, A. B. Larionov, T. Leit- ner, J. Weil, and U. Mosel, Phys. Rept. 512, 1 (2012), arXiv:1106.1344 [hep-ph]
2012 arXiv
-
[37]
Zeitnitz and T
C. Zeitnitz and T. A. Gabriel, Proc. of International Conference on Calorimetry in High Energy Physics (1993). 19
1993
-
[38]
G. D. Megias, J. E. Amaro, M. B. Barbaro, J. A. Ca- ballero, and T. W. Donnelly, Phys. Rev. D94, 013012 (2016)
2016
-
[39]
Abe et al
K. Abe et al. (Hyper-Kamiokande Proto- Collaboration), PTEP 2015, 053C02 (2015), arXiv:1502.05199 [hep-ex]
2015 arXiv
-
[40]
Abe et al
K. Abe et al. (T2K Collaboration), Phys. Rev. D96, 052001 (2017)
2017
-
[41]
Abe et al
K. Abe et al. (T2K Collaboration), Phys. Rev. D93, 112012 (2016), arXiv:1602.03652 [hep-ex]
2016 arXiv
-
[42]
Abe et al
K. Abe et al. (T2K Collaboration), Phys. Rev. D97, 012001 (2018)
2018
-
[43]
Abe et al
K. Abe et al. (T2K Collaboration), Phys. Rev. D90, 052010 (2014)
2014
-
[44]
Abe et al
K. Abe et al. (T2K Collaboration), Phys. Rev. D93, 072002 (2016), arXiv:1509.06940 [hep-ex]
2016 arXiv
-
[45]
simulation package to transport the final state particles through the ND280 detector complex. IV. EVENT AND KINEMA TIC SELECTION The event selection for antineutrino interactions is opti- mized to identify the observable charged current events with no charged or neutral pions i...
- [46]
-
[47]
Abe et al
K. Abe et al. (T2K Collaboration), Phys. Rev. D98, 032003 (2018)
2018
-
[48]
Abe et al
K. Abe et al. (T2K Collaboration), Phys. Rev. D87, 012001 (2013)
2013
-
[49]
Bohlen, F
T. Bohlen, F. Cerutti, M. Chin, A. Fasso, A. Ferrari, P. Ortega, A. Mairani, P. Sala, G. Smirnov, and V. Vla- choudis, Nuclear Data Sheets 120, 211 (2014), fLUKA transport code
2014
-
[50]
Ferrari, P
A. Ferrari, P. R. Sala, A. Fasso, and J. Ranft, CERN- 2005-010, SLAC-R-773, INFN-TC-05-11 FLUKA trans- port code
2005
-
[51]
R. Brun, F. Carminati, and S. Giani, CERN-W5013 (1994)
1994
-
[52]
Abgrall et al
N. Abgrall et al. (NA61/SHINE Collaboration), Phys. Rev. C 84, 034604 (2011)
2011
-
[53]
Abgrall et al
N. Abgrall et al. (NA61/SHINE Collaboration), Phys. Rev. C 85, 035210 (2012)
2012
-
[54]
Abgrall et al
N. Abgrall et al. (NA61/SHINE Collaboration), Eur. Phys. J. C 76, 84 (2016)
2016
-
[55]
Abe et al
K. Abe et al. (T2K Collaboration), Phys. Rev. D96, 092006 (2017), arXiv:1707.01048 [hep-ex]
2017 arXiv
-
[56]
Assylbekov et al
S. Assylbekov et al. (T2K ND280 P0D), Nucl. Instrum. Meth. A686, 48 (2012), arXiv:1111.5030 [physics.ins- det]
2012 arXiv
-
[57]
Abgrall et al
N. Abgrall et al. (T2K ND280 TPC), Nucl. Instrum. Meth. A637, 25 (2011), arXiv:1012.0865 [physics.ins- det]
2011 arXiv
-
[58]
Hayato, Acta Phys
Y. Hayato, Acta Phys. Pol. B40, 2477 (2009)
2009
-
[59]
Agostinelli et al
S. Agostinelli et al. (GEANT4), Nucl. Instrum. Meth. A506, 250 (2003)
2003
-
[60]
Baker and R
S. Baker and R. D. Cousins, Nucl. Instrum. Meth. A221, 437 (1984)
1984
-
[61]
Abe et al
K. Abe et al. (T2K Collaboration), Phys. Rev. D96, 011102 (2017), arXiv:1704.06409 [hep-ex]
2017 arXiv
-
[62]
Abe et al
K. Abe et al. (T2K Collaboration), Phys. Rev. D91, 112002 (2015), arXiv:1503.07452 [hep-ex]
2015 arXiv
-
[63]
Calvetti, S
D. Calvetti, S. Morigi, L. Reichel, and F. Sgallari, Jour- nal of Computational and Applied Mathematics 123, 423 (2000)
2000
-
[64]
P. C. Hansen, SIAM Rev. 34, 561580 (1992)
1992
- [65]
-
[66]
http://t2k-experiment.org/results/nd280data-anti- numu-cc-zero-pi-xs-on-water-2019,
“http://t2k-experiment.org/results/nd280data-anti- numu-cc-zero-pi-xs-on-water-2019,”
2019
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.