REVIEW 3 major objections 6 minor 1 cited by
Assessing the accuracy of the GENIE event generator with electron scattering data: Reduced cross section and nuclear transparency
T0 review · 3 major / 6 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read This paper shows that electron reduced cross sections can benchmark neutrino generators, and finds GENIE overestimates carbon data below missing momentum 80 MeV/c.
desk verdict A useful and honest GENIE benchmark against semi-exclusive (e,e'p) data with a visible low-pm discrepancy, though the approximate event-to-(pm, εm) mapping and thin quantitative details mean the size of the effect should be treated with caution. 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 reduced cross section $\sigma_{\mathrm{red}}(p_m) = (d\sigma/dp_m)/(K\,\sigma_{\ell N})$, the measured or simulated fivefold cross section divided by phase space and by the elementary electron- or neutrino-nucleon cross section; in the nonrelativistic plane-wave limit it becomes the bound-nucleon momentum distribution, and with final-state interactions it is the distorted spectral function. Because this object is nearly independent of the lepton probe, electron data can be used to test the nuclear parts of a neutrino event generator in isolation. The practical machinery is the mapping of generator events onto $(p_m,\varepsilon_m)$: GENIE's nuclear models only provide one-dimensional $d\sigma/dp_m$ with fixed binding energies, so the paper selects $1\mu1p$ events, defines $\tilde{\varepsilon}_m = \omega - T_p$, and uses the Jacobian $J(\theta_p)$ of Eq. (22) to convert the lepton-proton kinematics into missing momentum and missing energy.
What would settle it
Run the same five GENIE model configurations but reconstruct missing energy with the full recoil term, $\varepsilon_m = \omega - T_p - \varepsilon_B$, instead of $\tilde{\varepsilon}_m = \omega - T_p$, and repeat the comparison; if the low-$p_m$ overestimate disappears or moves to another kinematic region, the central claim fails as stated. A complementary check is to compare generator predictions against shell-tagged $(e,e'p)$ data selected in the $1p_{3/2}$ ($10<\varepsilon_m<25$ MeV) and $1s_{1/2}$ ($30<\varepsilon_m<50$ MeV) missing-energy windows, which GENIE's one-dimensional models cannot currently reproduce.
Extended reading notes
Core claim
The paper's central claim is that the reduced cross section $\sigma_{\mathrm{red}}(p_m)$ of semi-exclusive $(l,l'p)$ scattering is, up to small Coulomb corrections, the same object for electron and neutrino probes: after dividing out phase space and the elementary lepton-nucleon cross section, what remains is a distorted nuclear spectral function determined by the target and by the ejected nucleon's interaction with the residual nucleus. Using this identity, the paper generates $10^8$ charged-current quasielastic neutrino events per GENIE model configuration, keeps events that end as one muon plus one proton after the intranuclear cascade, maps them to $(p_m,\varepsilon_m)$ using the approximate missing-energy definition $\tilde{\varepsilon}_m = \omega - T_p$ and the Jacobian of Eq. (22), and compares the resulting $d\sigma/dp_m$ with four published $^{12}\mathrm{C}(e,e'p)$ datasets. On that comparison the paper reports persistent disagreements: at missing momenta below about 80 MeV/$c$, the GENIE models based on a local Fermi gas momentum distribution (G18_10a and G21_11) and also the spectral-function-based sf1d and effsf models overestimate the data, while at higher $p_m$ they tend to underestimate; the G18_02a Fermi-gas model shows the opposite pattern at low beam energy. A relativistic distorted-wave impulse approximation (RDWIA) treatment of the same kinematics is found to describe the measured reduced cross sections well, and GENIE's nuclear transparency for carbon is found to be roughly consistent with the transparency data, though the paper concludes that both GENIE and RDWIA need additional strength beyond single-nucleon knockout.
Load-bearing premise
The comparison assumes that GENIE events selected as one muon plus one proton and mapped to missing momentum/energy via the approximate formula $\tilde{\varepsilon}_m = \omega - T_p$ and the Jacobian of Eq. (22) reproduce the experimental $(e,e'p)$ reduced cross section; if that kinematic mapping distorts the low-$p_m$ region, the reported overestimate would be an artifact rather than a property of the nuclear models.
Editorial extensions
If this is right
- Neutrino energy reconstruction in experiments that rely on GENIE will inherit a systematic error in the low-missing-momentum region, where the ejected proton carries little recoil momentum.
- The electron-scattering datasets used here provide an out-of-sample test, since they were not used to tune the GENIE configurations; the reported disagreement therefore counts against those nuclear models rather than being a fit artifact.
- RDWIA-style unfactorized calculations, which include the distortion of the outgoing nucleon wave function, reproduce the data better and are a natural starting point for improving event-generator nuclear models.
- Because GENIE's nuclear transparency predictions agree roughly with data, the failure appears concentrated in the ground-state momentum distribution or quasielastic vertex rather than in the whole intranuclear-cascade description.
- The same reduced-cross-section benchmark can be applied to other targets and to other generators, turning precise electron data into a routine validation step for neutrino interaction models.
Reading between the lines
- If the reported low-momentum overestimate is real, its likely origin is that the LFG and effective spectral-function distributions put too much probability at small initial nucleon momenta; testing the generator's pre-FSI momentum distribution against measured spectral functions would separate that cause from distortions introduced by the cascade.
- The paper's comparison discards shell-resolved missing-energy information because GENIE's models are one-dimensional; implementing two-dimensional spectral functions in generators would make the low-$p_m$ comparison much sharper and might resolve whether part of the discrepancy is an artifact of the kinematic mapping.
- A practical extension would be to quantify how much the low-$p_m$ overestimate shifts the reconstructed neutrino energy in a simulated oscillation analysis; that number would translate this physics disagreement directly into a systematic uncertainty.
- Repeating the benchmark on a heavier target such as argon would show whether the failure is specific to carbon or generic to the momentum-distribution parametrizations used in current generators.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes using semi-exclusive electron-carbon scattering data—reduced cross sections and nuclear transparency—as a benchmark for the nuclear models implemented in the GENIE neutrino event generator. The authors simulate charged-current quasielastic neutrino-carbon events with five GENIE model configurations (RFG, LFG, SuSAv2-MEC, spectral function, and effective spectral function), convert the final-state kinematics into a reduced cross section as a function of missing momentum via an approximate missing-energy definition, and compare the results with (e,e'p) data from Tokyo, Saclay, SLAC, and JLab. They report that GENIE models based on the local Fermi gas overestimate the measured reduced cross sections for missing momenta below about 80 MeV/c, while RDWIA calculations agree better. They also compare GENIE nuclear transparency predictions with electron-scattering transparency data and find satisfactory agreement, with some shortcomings attributed to missing multinucleon strength.
Significance. If the comparison is valid, the paper offers a novel, externally grounded test of the nuclear models in a widely used neutrino event generator, using electron-scattering data that were not used to tune those models. The low-missing-momentum region is directly relevant to neutrino energy reconstruction in oscillation experiments, so the claimed discrepancy would be of practical importance. The transparency comparison adds a second, independent validation channel. The paper also demonstrates a general methodology for converting generator events into spectral-function-like observables. However, the central claim depends entirely on the event-to-(pm, epsilon_m) mapping and on matching the missing-energy windows of the experimental data, and the quantitative evidence is reported without a complete goodness-of-fit treatment.
major comments (3)
- [III, Eqs. (17)-(24)] The conversion from GENIE events to the reduced cross section is not sufficiently justified. The paper states in Sec. III that the GENIE nuclear models provide only one-dimensional dsigma/dpm distributions with fixed binding energies, so they cannot produce the shell-resolved d²sigma/dpm/depsilon_m. Nevertheless, Eq. (24) defines dsigma/dpm = Deltaepsilon_m d²sigma/dpm/depsilon_m, which is only correct if d²sigma/dpm/depsilon_m is flat over the window Deltaepsilon_m; for a fixed-binding model the epsilon_m distribution is a narrow peak, not a flat density. Please show explicitly how the event counts are binned and integrated over epsilon_m (for example, by counting events in the full Deltaepsilon_m window and dividing by Delta pm), and state which missing-energy window the electron data actually correspond to. If the data are the shell-separated bins 10-25 MeV and 30-50 MeV mentioned in Sec. III while the GENIE curves are integrated over 80 MeV, the comparison is not apples-to-apples and the reported low-pm overestimate could be an artifact of the window mismatch.
- [V, chi-squared paragraph] The quantitative support for the central claim is incomplete. After Fig. 7 the paper reports chi2/DOF approximately 12.6 for the effsf model with SLAC data and approximately 17.5 for the sf1d model with JLab data, but it never gives the number of degrees of freedom, the definition of chi2, the error treatment (statistical errors of the data, normalization uncertainties, or any uncertainty bands on the GENIE predictions), or the pm range used in the comparison. Without this information the phrase "persistent disagreements" is not quantitatively established; please provide a defined goodness-of-fit statistic with all components and show the sensitivity of the conclusion to the treatment of data normalization uncertainties.
- [III, event selection (1p+Nn vs 1p+0n)] The choice of the 1p+Nn event set as the analogue of the proton-only final state in (e,e'p) needs justification. Figure 2 shows that for the G21_11 and spectral-function models the missing-energy distributions of 1p+Nn and 1p+0n events differ noticeably, and the text notes that inelastic FSI changes proton energies and pushes events to higher missing energy. Since the electron data select a knocked-out proton with no associated neutron emission, including events with additional neutrons in the GENIE comparison may bias sigma_red(pm), and the paper does not quantify this bias. Please show the effect of the two selections on the final reduced cross sections, especially at pm below 80 MeV/c, or justify why the 1p+Nn selection is equivalent to the experimental proton-only selection for this observable.
minor comments (6)
- [Abstract] The sentence "The approach presented in this paper provide a great opportunity" should read "provides a great opportunity".
- [Fig. 7 caption] The caption says "Same as Fig. 3 but for the sf1d (SF) model calculation," but the text and figure content describe the effsf model; the caption should be corrected.
- [Introduction, section numbering] The Introduction states that results are presented in Sec. IV, but the results actually appear in Sec. V; the cross-reference should be corrected.
- [Eq. (22)] "Jacobin" should be "Jacobian," and the sign convention in the second term of the Jacobian expression should be checked for consistency with the definition of pm in Eq. (20).
- [Introduction] "Femilab" should be "Fermilab."
- [Fig. 9 caption] "transparensy" should be "transparency."
Circularity Check
Minor fitted-input overlap for the sf1d comparison; the central LFG-model benchmark is external and non-circular.
-
fitted input called prediction
[Sec. IV (model descriptions) and Sec. V (Fig. 6, sf1d comparison)]
"The mean field piece has been fitted to the (e,e’p) scattering data ... In Fig. 6 we compare the GENIE sf1d model calculations with the data for reduced cross sections of electron scattering on carbon."
The input to the sf1d model, its mean-field spectral function, was previously fitted to (e,e'p) scattering data. Since the paper's Eq. (25) identifies the reduced cross section with the (distorted) spectral function, the Fig. 6 comparison of GENIE sf1d predictions with 12C(e,e'p) reduced-cross-section data is in part a comparison of a fitted quantity with its training data. This is a partial reduction by construction for the sf1d configuration only. The paper's central conclusion about the LFG-based G18_10a and G21_11 models does not rest on this fitted input, so the circularity is minor.
full rationale
The core benchmark is external and non-circular: GENIE v3.04.00 predictions for five model configurations are generated from the generator and compared with published 12C(e,e'p) reduced-cross-section data (Saclay, Tokyo, SLAC, JLab), none of which are fitted in this paper. The headline low-pm overestimate concerns the LFG-based G18_10a and G21_11 models, whose GENIE implementations are not tuned to these data. The RDWIA comparisons cited from Refs. [23-26] are prior parameter-free calculations by the same authors; they were not fitted to the data used here and are externally falsifiable, so invoking them is legitimate independent support rather than circularity. The kinematic mapping in Eqs. (17)-(25), including the approximate missing-energy definition, the Jacobian, and the flatness assumption behind Eq. (24), is an approximation that could bias the low-pm region, but it is a validation/correctness concern rather than a circular definition. The only mild circular element is the sf1d model comparison, since that model's mean-field spectral function was previously fitted to (e,e'p) data; this does not affect the LFG-based central conclusion. Overall circularity is minor, with a score of 2.
Assumptions & free parameters
free parameters (4)
- Fermi momentum and binding energy for 12C in GENIE RFG (G18_02a) =
GENIE defaults, not stated in paper
- Local Fermi gas and RPA parameters (G18_10a, G21_11a) =
Not stated in paper
- Spectral function mean-field normalization (sf1d) =
Fitted to (e,e'p) data in Refs. [46,47]
- Effective spectral function transverse enhancement (effsf) =
Not stated in paper
assumptions (4)
- domain assumption The factorized form d5sigma = K x sigma_lN x sigma_red (Eq. 14) is accurate enough for the kinematics used.
- domain assumption Electron and neutrino reduced cross sections for the same nucleus agree up to Coulomb corrections.
- domain assumption Coulomb corrections can be represented by an effective momentum transfer qeff and are about 10% at low beam energy.
- domain assumption The INTRANUKE hA cascade models proton final-state interactions well enough for reduced cross section and transparency comparisons.
Cite this review
Pith. "Pith review of Assessing the accuracy of the GENIE event generator with electron scattering data: Reduced cross section and nuclear transparency." pith.science (2026). https://pith.science/paper/C5HWEVZG
@misc{pith2026250619355,
author = {Pith},
title = {Pith review of: Assessing the accuracy of the GENIE event generator with electron scattering data: Reduced cross section and nuclear transparency},
year = {2026},
howpublished = {\url{https://pith.science/paper/C5HWEVZG}},
note = {Machine review of arXiv:2506.19355}
}
abstract
The reduced cross sections of the semi-exclusive $(l,l'p)$ lepton scattering process can be identified with distorted nuclear spectral functions. Irrespective of the type of interaction the distorted spectral function is determined mainly by intrinsic properties of the target and the ejected nucleon interaction with residual nucleus. Thus the reduced cross sections of the neutrino and electron scattering on nuclei as functions of nucleon missing momentum and energy are similar up to Coulomb corrections. Here we demonstrate the utility of this approach by benchmarking the GENIE neutrino event generator against a broad set of semi-exclusive electron scattering data for carbon target. We observe persistent disagreements between the generator predictions and data in the range of nucleon missing momenta less than 80 MeV/c. The approach presented in this paper provide a great opportunity to test better the accuracy of nuclear models of quasi-elastic neutrino-nucleus scattering, employed in neutrino event generators.
Figures
Figures from the paper (7 more)
Forward citations
Cited by 1 Pith paper
-
Explicitly on-shell currents in relativistic mean field models
The 'off-shell ambiguity' in relativistic mean-field one-body currents is spurious when currents are defined via free-nucleon operators; the correct mapping (Eq. 26) is parametrization-independent and requires time-li...
Reference graph
Works this paper leans on
-
[1]
M. A. Acero et al. , (NOvA Collaboration), Phys. Rev. Lett. 123, 151803 (2019)
work page 2019
- [2]
- [3]
- [4]
-
[5]
M. Antonello et al. (MicroBooNE, LAr1-ND, ICARUS-W A104 Collaboration), arXiv:1503.01520 [physics.ins-det]
-
[6]
A. Picklesimer, J. W. Van Orden, S. J. Wallace, Phys. Rev. C 32 , 1312 (1985)
work page 1985
-
[7]
J. M. Udias, P. Sarriguren, E. Moya de Guerra, E. Garrido, and J. A. Caballero, Phys. Rev. C 51 , 3246 (1995)
work page 1995
- [8]
Show all 52 references
-
[9]
Gonzalez-Jimenez, M
R. Gonzalez-Jimenez, M. B. Barbaro, J. A. Caballero, T. W . Donnelly, N. Jachowicz, G. D. Megias, K. Niewczas, A. Nikolakopoulos, J. M. Udias, Ph ys. Rev. C 101 , 015503 (2020)
2020
-
[10]
A. V. Butkevich and S. V. Luchuk, Phys. Rev. C 102 024602 (2020)
2020
-
[11]
G. D. Megias, J. E. Amaro, M. B. Barbaro, J. A. Caballero, T. W. Donnelly, Phys. Rev. D 94, 013012 (2016)
2016
-
[12]
G. D. Megias, J. E. Amaro, M. B. Barbaro, J. A. Caballero, T. W. Donnelly, and I. R. Simo, Phys. Rev. D 94 , 093004 (2016)
2016
-
[13]
A. M. Ankowski and A. Friedland, Phys. Rev. D 102 , 053001 (2020)
2020
-
[14]
Papadopoulou et al
A. Papadopoulou et al. (e4v Collaboration) Phys. Rev. D 103 , 113003 (2021)
2021
-
[15]
Khachatryan, A
M. Khachatryan, A. Papapdopoulou, A. Ashkenazi, F. Hau enstein, L. B. Weinstein, O. Hen, E. Piasetzky (CLAS and e4v Collaboration), Nature 599, 565 (2021)
2021
-
[16]
Avarez-Ruso, et al
L. Avarez-Ruso, et al. (GENIE Collaboration), Eur. Phys. J. Spec. Top. 230, 4449 (2021)
2021
-
[17]
Zmuda, K
J. Zmuda, K. Graczyk, C. Juszczak, and J. Sobczyk, Acta. Phys. Pol. B46, 2329 (2015). 24
2015
-
[18]
R. D. Banerjee, A. M. Ankowski, K. M. Graczyk, B. E. Kowal , H. Prasad, and J. T. Sobczyk, Phys. Rev. D 109 , 073004 (2024)
2024
-
[19]
Isaacson, W
J. Isaacson, W. J. Jay, A. Lovato, P. A. N. Machado, and N. Rocco, Phys. Rev. D 107, 033007 (2023)
2023
-
[20]
Abe, Phys
S. Abe, Phys. Rev. D 111 , 033006 (2025)
2025
-
[21]
Dytman, Y
S. Dytman, Y. Hayato, R. Raboanary, J. T. Sobczyk, J. Ten a Vidal and N. Vololoniaina, Phys. Rev. D 104 , 053006 (2021)
2021
-
[22]
Niewczas, and J
K. Niewczas, and J. T. Sobczyk, Phys. Rev. C 100 , 015505 (2019)
2019
-
[23]
A. V. Butkevich and S. A. Kulagin, Phys. Rev. C 76 , 045502 (2007)
2007
-
[24]
A. V. Butkevich, Phys. Rev. C 80 , 014610 (2009)
2009
-
[25]
A. V. Butkevich, Phys. Rev. C 85 , 065501 (2012)
2012
-
[26]
A. V. Butkevich, Phys. Rev. C 109 , 045502 (2024)
2024
-
[27]
Andreopoulos, A
C. Andreopoulos, A. Bell, D. Bhattacharya, F. Cavanna, J. Dobson, S. Dytman, H. Gallagher, P. Guzowski, R. Hatcher, P. Kehayias et al , Nucl. Instrum. Meth. A614, 87 (2010)
2010
-
[28]
The GENIE Event Generator, http://genie-mc.org
-
[29]
J. A. Caballero, T. W. Donnelly, E. Moya de Guerra, and J. M. Udias, Nucl. Phys. A 632 , 323 (1998)
1998
-
[30]
Nikolakopoulos, R
A. Nikolakopoulos, R. Gonzalez-Jimenez, N. Jachowicz , K. Niewczas, F. Sanchez, J. M. Udias, Phys. Rev. C 105 054603 (2022)
2022
-
[31]
Nakamura, and N
K. Nakamura, and N. Izutsu, Nucl. Phys. A 259 , 301 (1976)
1976
-
[32]
J. J. Kelly, Phys. Rev. C 71 , 064610 (2005)
2005
-
[33]
Mougey, M
J. Mougey, M. Bernheim, A. Bussiere, A. Gillebert, P. X. Ho, M. Priou, D. Royer,I. Sick, and G. Wagner, Nucl. Phys. A 262 , 461 (1976)
1976
-
[34]
N. C. R. Makins, R. Ent, M. S. Chapman, J. O. Hansen, K. Lee , Milner R.G. et al. , Phys. Rev. Lett. 72, 1986 (1994)
1994
-
[35]
Dutta et al
D. Dutta et al. , Phys. Rev. C 68 , 064603 (2003)
2003
-
[36]
Andreopoulos, C
C. Andreopoulos, C. Barry, S. Dytman, H. Gallagher, T. G olan, R. Hatcher, G. Perdue, J. Yarba, arXiv:1510.05494
-
[37]
Alvarez-Russo et al
L. Alvarez-Russo et al. , (GENIE Collaboration), Eur. Phys. J. Spec. Top. 230, 4449 (2021)
2021
-
[38]
Tena-Vidal et al
J. Tena-Vidal et al. , (GENIE Collaboration), Phys. Rev. D 106 , 112001 (2022). 25
2022
-
[39]
S. J. Mashnik, A. J. Sierk, K. K. Gudima, and M. I. Baznat, J. Phys. Conf. Ser. 41, 340 (2006)
2006
-
[40]
Bodek and J
A. Bodek and J. L. Ritchie, Phys. Rev. D 23 , 1070 (1981)
1981
-
[41]
Llewellyn Smith, Phys
C. Llewellyn Smith, Phys. Rep. 3, 261 (1972)
1972
-
[42]
Nieves, I
J. Nieves, I. Ruiz Simo, and M. J. Vicente Vacas, Phys. Le tt. B 707 , 72 (2012)
2012
-
[43]
Nieves, J
J. Nieves, J. E. Amaro, and M. Valverde, Phys. Rev. C 70 , 055503 (2004).Phys. Rev. D106, 112001 (2022)
2004
-
[44]
R. Gran, J. Nieves, F. Sanchez, and M. J. Vicente Vacas, P hys. Rev. D88, 113007 (2013)
2013
-
[45]
Dolan, G
S. Dolan, G. D. Megias, and S. Bolognesi, Phys. Rev. C 101 , 033003 (2020)
2020
-
[46]
Benhar, A
O. Benhar, A. Fabrocini, S. Fantoni, Nucl. Phys. A 505 , 267 (1989)
1989
-
[47]
Benhar, A
O. Benhar, A. Fabrocini, S. Fantoni, Nucl. Phys. A 579 , 493 (1994)
1994
-
[48]
Bodek, H
A. Bodek, H. S. Budd, M. E. Christy, Eur. Phys. J. C 71 , 1726 (2011)
2011
-
[49]
Bodek, H
A. Bodek, H. S. Budd, M. E. Christy, B. Coopersmith, Eur. Phys. J. C 74 , 3091 (2014)
2014
-
[50]
McKean, R
J. McKean, R. Gonzalez-Jimenez, M. Kabirnezhad, J. M. U dias, arXiv:2502.10629
-
[51]
Garino, et al
G. Garino, et al. , Phys. Rev. C 45 , 780 (1992)
1992
-
[52]
Rohe et al
D. Rohe et al. , Phys. Rev. C 72 , 054602 (2005)
2005
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.