Exploring Neutrino-Nucleus Interactions in the GeV Regime using MINERvA

With the advance of particle accelerator and detector technologies, the neutrino physics landscape is rapidly expanding. As neutrino oscillation experiments enter the intensity and precision frontiers, neutrino-nucleus interaction measurements are providing crucial input. MINERvA is an experiment at Fermilab dedicated to the study of neutrino-nucleus interactions in the regime of incident neutrino energies from one to few GeV. The experiment recorded neutrino and antineutrino scattering data with the NuMI beamline from 2009 to 2019 using the Low-Energy and Medium-Energy beams that peak at 3 GeV and 6 GeV, respectively. This article reviews the broad spectrum of interesting nuclear and particle physics that MINERvA investigations have illuminated. The newfound, detailed knowledge of neutrino interactions with nuclear targets thereby obtained is proving essential to continued progress in the neutrino physics sector.


Introduction
Neutrinos with energies of a few GeV are involved in many different ways among phenomena that present opportunities to probe fundamental aspects of physical reality. Neutrinos produced in accelerators play a central role in precision measurements of the oscillation parameters such as the Dirac CP-violating phase that may be present in the neutrino flavor mixing matrix [1,2,3,4]. Measurement of a non-zero Dirac phase could unlock the mystery of the matter-antimatter asymmetry of the Universe. Neutrino beams serve as potential sources of beyond-Standard-Model (BSM) particles, such as light dark matter and heavy neutral leptons [5,6]. On the other hand, neutrinos could impede the discovery of such new forms of matter by mimicking their BSM signatures. This possibility exists because some neutrino SM processes in detector materials have aspects that are poorly known. Atmospheric neutrinos that are born in cosmic-ray-induced hadronic cascades in the upper atmosphere propagate through the Earth [7], presenting complications as well as opportunities for new physics searches. Atmospheric neutrinos oscillate, and their oscillations undergo highly interesting alterations due to propagation through a matter field. However, these highly penetrating particles also create background to rare-event searches (such as proton decay [8,9]) in deep underground experiments. Understanding how a neutrino interacts with a nucleus is essential for exploiting these opportunities. Given that current GeV-neutrino sources (accelerators or atmospheric) are not monoenergetic, these energy-sensitive interactions are convolved with the neutrino flux, causing major systematic uncertainties in precision measurements.
Neutrino-nucleus (ν-A) interactions arise not only from the primary nucleon-level interaction, but also from the effect that the nuclear environment exerts on the initialstate nucleons and the final-state particles. Since a theory of the complete nuclear response in neutrino-nucleus interactions in the few-GeV regime of incident neutrino energy is yet to be developed [10], comprehensive ν-A measurements are needed to guide and benchmark the development of models. MINERvA (Main INjector Ex-peRiment for ν-A) at Fermilab is a dedicated experiment to illuminate the interplay between hadronic and nuclear degrees of freedom in ν-A interactions and to measure aspects of intranuclear dynamics that are prerequisites for precision neutrino oscillation measurements.
MINERvA received the NuMI (Neutrinos-from-the-Main-Injector) beam at a distance of 1 km from the target of the 120-GeV primary proton beam at Fermilab. In the Low-Energy (LE) beam configurations operated between 2009 and 2012, both the ν µ andν µ fluxes peak at ∼3 GeV, while in the Medium-Energy (ME) configurations used between 2013 and 2019, the fluxes peak at ∼6 GeV. In both beams, there is a high-energy component that extends beyond 50 GeV. The data collected by MIN-ERvA correspond to 4.0 (1.7) and 12.1 (12.4) times 10 20 protons on target (POT) for the LE and ME ν µ (ν µ ) configurations, respectively. In the Sections below, the neutrino interaction physics investigated by MINERvA is reviewed, with the main focus being the techniques developed and measurements reported that are based on the Low-Energy data set.

MINERvA Experiment and Flux Predictions
The MINERvA detector [11] utilizes extruded plastic scintillator as its tracking medium. Most of the active mass is located in its central polystyrene target whose 5.4-t fiducial volume serves as a charged-particle tracker. The upstream section of the detector consists of a series of passive targets (helium, carbon, water, iron, and lead) interleaved with tracking planes. Sampling calorimeters surround both passive and active target regions. Muons produced by ν-A charged-current (CC) interactions in the tracker or in the upstream targets exit the downstream end of the tracker. These muons may then enter and propagate through the magnetized MINOS near detector [12] located 2 m downstream, allowing their trajectories to be momentumanalyzed.
The downstream tracking of muons in the MINOS near detector provides an 8 % muon momentum resolution (at 5 GeV/c). This complements the three-dimensional tracking and energy-loss measurement of final-state particles that is afforded by the MINERvA tracker. In the tracker, the momentum resolution for protons is 2 % at 1 GeV/c with a 450-MeV/c tracking threshold. The hadronic energy response of the detector was calibrated using test-beam measurements [13], with pion calorimetric energy resolution in the range 20-30 %. Furthermore, the tracker is large compared to the interaction length of the neutrons (approximately 10 cm at 20 MeV) produced in ν-A interactions. The 1.5-MeV detection threshold for measuring energy deposit allows interacting neutrons to be registered for a time-of-flight measurement [14]. The neutron timing resolution is 4.5 ns where the hit resolution alone is 3 ns from electronics effects.
The NuMI beam fluxes used by MINERvA are modeled with GEANT4 predictions that are adjusted to match hadron production data [15]. The large fiducial volume of the MINERvA tracker and the intense beam fluxes make it possible to use neutrino scattering on atomic electrons, ( -) ν µ e − → ( -) ν µ e − , to further constrain the flux predictions. In the data from the LE (ME) NuMI beam configurations, 135 (810) neutrino-electron scattering events were identified [16,17]. The constraints thereby provided reduce the ν µ flux normalization in the a priori prediction by 6 % for LE, and 10 % for ME (Fig. 1). Moreover, the uncertainty at the flux peak is reduced from 9 to 6 % for LE, and from 8 to 4 % for ME. In addition, the "low-ν" method has been used to constrain the flux shape [18,19]. The latter method exploits the minimal neutrino energy dependence of the inclusive charged-current cross section at low hadronic recoil energy.

Incoherent Neutrino-Nucleus Interactions
In the few-GeV regime, neutrino-nucleus cross sections are dominated by incoherent processes in which the constituent nucleons can be ejected, possibly accompanied by pions and other mesons. These primary processes are quasielastic (QE), baryon resonance production (RES) including non-resonant background, and deep inelastic scattering (DIS), in ascending order of excitation by incident neutrinos of increasing lab-frame energies. The initial energy and momentum of the struck nucleon and its correlation with other nucleons (including long-range correlations and two-particletwo-hole excitations, or 2p2h [20,21,22,23,24]) contribute to the initial-state conditions. The primary hadronic final-state particles propagate inside the remnant nucleus and may participate in final-state intranuclear interactions (FSI). The latter interactions may further excite the remaining system, causing nucleon emission or even spallation. Because pions can be absorbed or created during FSI, there is no unique experimental signature for a given primary process. As a result, the particle content and energy budget of a ν-A interaction varies with the initial and final states together with the primary reaction.
By restricting the final-state topology, MINERvA can examine exclusive and semiinclusive reactions such as mesonless (i.e., quasielastic-like), pion, and kaon production; the experiment can study inclusive scattering as well. The exclusive-channel studies must take into account pion absorption through FSI, which enables resonance production and DIS reactions to be present in quasielastic-like topologies.
The MINERvA data enable the elucidation of one-particle-one-hole mechanisms that are generally used in ν-A scattering models, and they also allow examination of the significant 2p2h contributions to the quasielastic-like process. Concerning the latter, MINERvA has identified a 2p2h-like enhancement that simultaneously describes both neutrino and antineutrino scattering data at the kinematic region between quasielastic and resonance production. These unmodeled additional event rates are likely an admixture of all three reactions and their precise nature is still under study.
Decoupling the primary process and medium effects is challenging, especially in a wide-band neutrino beam where the neutrino energy (E ν ) is unknown. In first order, reaction E ν -dependence comes from the primary interactions and the final-state momenta and angles depend strongly on E ν . The nuclear response, on the other hand, affects these elementary distributions as a perturbation; it depends on E ν through the medium coupling to the primary initial-and final-state hadrons. With a neutrino beam where the neutrino direction is precisely known, the kinematics projected onto the transverse plane to the neutrino will have less dependence on E ν .
In a charged-current measurement, certain final-state correlations of the lepton and the hadronic system, such as the transverse kinematic imbalance (TKI) [25,26], avoid or cancel the primary-level dependence on, for example, E ν and axial form factors, and are directly sensitive to the nuclear response with minimal dependence on the neutrino energy. Because exclusivity can be achieved in the transverse plane to the neutrino direction, TKI can probe the hidden dynamics inside the target nucleus. MINERvA has systematically explored the potential of TKI to identify the medium properties and interaction dynamics in exclusive processes.
A review of MINERvA's measurements of incoherent interactions is presented in this section. These measurements encompass both exclusive and elastic reactions to inclusive and inelastic processes.

First Measurements of Quasielastic and Quasielastic-Like Cross Sections
The charged-current quasielastic interaction (CCQE) is an important channel in neutrino oscillation experiments. This interaction gives rise to the majority of events in T2K [1] and to a sizable fraction of the events recorded by NOvA [2]. In quasielastic scattering, the neutrino energy can in principle be inferred from the outgoing lepton kinematics. However, this relatively simple interaction incurs significant effects from the nuclear environment.
Candidate CCQE interactions of neutrinos and antineutrinos on nucleons, ν µ n → µ − p andν µ p → µ + n, respectively, were extracted from the Low-Energy data in the first two cross-section measurements reported by MINERvA [27,28]. After subtracting the non-QE background processes, namely resonance production and DIS, the fluxintegrated differential cross sections in four-momentum transfer squared, Q 2 , were compared to model predictions by GENIE [29] and NuWro [30] (Fig. 2). In the neutrino generator predictions, the shape of the Q 2 distribution is parameterized by the axial vector mass, M A , plus model representations of the nuclear state (relativistic     Fermi gas, or RFG [31], and Spectral Function, or SF [32]). Both the neutrino and antineutrino data sets were found to favor a parametric enhancement in the magnetic form factor (Transverse Enhancement Model, or TEM [33]), in addition to the atthat-time-standard choice of RFG and the world-average M A value of 0.99 GeV/c 2 . Becasue TEM is extracted from a fit to electron scattering data to describe the contributions from two-nucleon knock-out processes, these results suggest possible contributions from 2p2h in neutrino and antineutrino scattering. This interpretation is further supported by the observed pattern of energy deposits near the interaction vertices (vertex energy) in both measurements.
MINERvA also measured the ν e CCQE cross section in the Low-Energy neutrino data set [34]. In the Standard Model, lepton couplings are universal. However, because the final-state lepton mass is different in ν e versus ν µ CC events, the nuclei respond to slightly different phase spaces-a difference that cannot be ignored in oscillation experiments [35]. As the CC-induced electrons are reconstructed in the tracker which is charge-insensitive, a flux-averaged cross section including limited contributions from ν e is obtained. By comparing to the ν µ counterpart [28], it has been shown that the two CCQE cross sections are consistent (Fig. 3). Moreover, GENIE describes the ν e -to-ν µ cross-section ratio within the experimental error which is of order 10-20%.
In MINERvA quasielastic measurements, the inelastic background arising from primary pion production that is followed by absorption during FSI must be subtracted. In subsequent MINERvA studies of the quasielastic dynamics, the inelastic pionless events are instead considered as part of the "QE-like" signal definition in order to mitigate the model uncertainties for pion absorption. In MINERvA's first QE-like cross-section measurement [36], the Q 2 distribution is extracted. In addition to the change of signal definition, the Q 2 calculation is completely hadron-based, using only the momentum of the leading final-state proton above the tracking threshold. This was the first MINERvA measurement to use reconstructed proton tracks. The distribution is compared to model predictions in order to test different hypotheses concerning initial-state correlations and the data are found to disfavor the modeled correlations (Fig. 4). On the other hand, proton kinematics such as the proton-based Q 2 are sensitive to FSI whose strength depends on the size of the target nucleus due to the intranuclear energy loss. This motivated further quasielastic-like cross-section measurements by MINERvA using the detector's upstream target planes (carbon, iron, and lead) together with the active tracker (CH) [42]. The extracted cross-section ratios of iron and lead to CH are shown in Fig. 5. The data are found to be reasonably well described by the A-dependent FSI prescriptions used by the GENIE and NuWro event generators.

Initial-State Correlations and 2p2h-Like Enhancement
The nuclear response are often described theoretically using the energy transfer, q 0 (also called ω or ν in the literature), and the three-momentum transfer, q 3 or | q|, from the lepton to the target nucleus, whereby q 2 3 − q 2 0 = Q 2 . However, the estimation of q 0 experimentally introduces uncertainties due to missing energy from nucleon unbinding and neutrons in the final state. To avoid systematic errors that may enter in this way,  MINERvA introduced a new observable called "available energy", E avail , as a proxy for q 0 in the "low-recoil" (meaning q 3 < 0.8 GeV) analyses using the Low-Energy ν µ andν µ CC inclusive samples [43,44]. Available energy is defined as the sum of the energies from final-state particles that can be calorimetrically measured in the scintillator and excludes removal energy and neutron energy. An estimate of q 0 is still used with muon kinematics to calculate E ν and then q 3 , but its model uncertainties are subdominant for q 3 . The measured distributions for reconstructed E avail of ν µ scattering [43] are shown in Fig. 6 (I) for two regions of the reconstructed q 3 . They are compared to the default simulation based on GENIE with modifications to pion production [45,46]. The simulation is further improved by taking into account initial-state correlations: the collective long-range medium effect in quasielastic events calculated with a Random Phase Approximation (RPA) approach [37], and the Valencia QE-like 2p2h model [39,40,23,41]. Inclusion of RPA brings the event rates at low E avail into better agreement with the data, though beyond-Fermi-gas models may have a similar effect [47,48]. The Valencia 2p2h model helps reduce the model deficit at the dip region between quasielastic and resonance production. The best data-model agreement is achieved by separately scaling up the 2p2h event rates in regions of q 0 and q 3 : across all q 0 -q 3 regions, an enhancement by 50 % is required to enable the model to match the data, and in the dip region an enhancement of a factor of 2 is needed. Since this ad hoc enhancement is based on the Valencia 2p2h model, it was first interpreted as a correction to the modeled 2p2h mechanism. The GENIE model that has evolved to this end is denoted as MnvGENIE in subsequent MINERvA publications.
The correspondingν µ measurement [44] is shown in Fig. 6 (II), where the aforementioned default simulations and MnvGENIE are compared in (II.a) and (II.b), respectively. Importantly, the 2p2h enhancement in MnvGENIE is only tuned to the ν µ measurement. Consequently, its predictive power forν µ , coming from the same underlying q 0 -q 3 -dependence as for ν µ , is intriguing. It suggests that the enhancement is at the level of structure functions rather than of FSI or Fermi motion. However, given its empirical nature, this enhancement should be understood as 2p2h-like; the analysis does not currently rule out a quasielastic or resonant enhancement in addition to 2p2h.
The inclusion of RPA and 2p2h in the predictions has made a significant improvement in the description of the MINERvA data. The evolution of modeling improve-  ment from GENIE to MnvGENIE is illustrated by a reanalysis of the vertex energy in a ν µ CCQE-like sample [49]. As is shown in Fig. 7, the distributions of the reconstructed vertex energies in events with and without proton tracks are described by GENIE after the modifications with RPA, Valencia 2p2h, and the 2p2h-like enhancement are taken into account.

State-of-the-Art Quasielastic-Like Measurements
In the few-GeV range of the neutrino energy, the quasielastic cross section is nearly constant as a function of E ν , while the phase space for resonance production and DIS is opening up. In addition to extracting the CCQE-like cross sections in Q 2 and E ν which can be calculated with the quasielastic hypothesis (that is invalid for the non-QE components), MINERvA measured the CC muon transverse (p T ) and longitudinal momenta (p || ) using the Low-Energyν µ [50], ν µ [49], and Medium-Energy ν µ [51] data sets. These particular muon momentum projections respectively approximate the true Q 2 and E ν . As shown in Fig. 8 for the Low-Energy ν µ results, the different interaction contributions in the MnvGENIE predictions are relatively stable across the p || bins. This is in contrast to the inclusive measurement discussed in Sect. 3.6 below where the p || -(and therefore, E ν -) dependence of the DIS processes is evident. In both the QE-like and the inclusive measurements, while MnvGENIE describes the data in most bins, there is an overall model deficit at large p || .
By measuring the transverse kinematic imbalance which cancel out the primary interaction kinematics, the cross section dependence on the incoming neutrino energy is lessened and also the initial-and final-state effects can be directly probed [26]. Using the ν µ CCQE-like events in the tracker with the Low-Energy data [52], the direction and magnitude of the transverse momentum imbalance (δ p T ) between the muon and the leading proton, δα T and δp T , respectively, are calculated. The angle δα T has the most sensitivity to FSI and to the unaccounted-for momentum carried by missing particles such as absorbed pions or the correlated nucleon of the proton from 2p2h. Figure 9 shows that, within the uncertainties, MnvGENIE describes the data. Here, FSI are classified into three categories: 1. A flat distribution for events which do not experience FSI (both the "no-FSI" and "p-FSI non-interacting" categories in Fig. 9) reflecting the isotropy of the Fermi motion; 2. The deceleration region (δα T → 180 • ) for energy-dissipating processes-decelerating FSI, pion absorption, and 2p2h; 3. The acceleration region (δα T → 0 • ) for accelerating FSI if such a mechanism exists. Interestingly, GENIE did predict FSI acceleration for protons, roughly half of them singularly occupying the acceleration region, while the other half falls into the deceleration region due to the transverse projection. By assuming a carbon-11 target remnant, the transverse momentum imbalance magnitude δp T is promoted to the three-dimensional momentum imbalance, p n , following from an additional constraint by energy conservation [53]. For the FSI-noninteracting events, p n can be interpreted as the momentum of the struck neutron in the CCQE initial state. The location of the Fermi-motion peak in data is well captured by the MnvGENIE prediction. However, its accelerating FSI component causes the predicted peak shape to deviate from data. This component was identified as the elastic component of the GENIE v2.8 hA FSI model and has been removed in later versions of GENIE [54,55]. The measured cross section is further compared to NuWro predictions with alternative nuclear states: local Fermi gas (LFG) and Spectral Function. The latter model better describes the Fermi motion peak, but neither model provides enough strength in the transition region between the quasielastic peak and the non-QE tail [52].
An enhanced sensitivity to another initial-state condition, the binding energy, is achieved by further projecting δ p T onto the lepton scattering plane, which defines the δp Ty variable [55]. Various generator implementations of the interaction energy on carbon are compared (Fig. 10) and the data favor approximate corrections to GENIE [56].  Fig. 10. Mean δpTy calculated from the νµ CCQE-like differential cross section in δpTy within -0.20 to 0.20 GeV/c. The p-value is the probability, assuming normal distribution, that the observed result would have been produced by the test model. Uopt and V eff are the optical and Coulomb potentials experienced by the CCQE proton and muon, respectively, as proposed in Ref. [56]. NEUT [57] and GiBUU [58] predictions are also compared. Figure  from Ref. [55].

Charged-Current Pion Production
Neutrino-induced pion production is an important channel in neutrino oscillation experiments because it accounts for a large part of the signal in NOvA and DUNE [3] far detector event samples, and is both a low-statistics signal and background in the far detector samples of T2K. The process can proceed through baryon resonance production or through non-resonant interaction. It can also occur through a CC coherent interaction with a nucleus-see discussions in Sect. 4. Because of the additional final-state hadrons, incoherent pion production may be more affected by the nuclear environment than the quasielastic process, and the difficulty with reconstructing pions means that the process is harder to measure. One would like to use the lessons learned about nuclear effects in quasielastic-like scattering to (overt) pion production, but the effects may not be the same. The nuclear responses to the pion production fall into three categories according to their relation with those of the quasielastic-like processes: 1. Nuclear responses that are generic but might be quantitatively different in pion production are Fermi motion, binding energy, initial-state correlations, and Pauli blocking [59,60]; 2. Pion absorption and charge exchange are nuclear responses that migrate primary pion production channels among each other and into quasielastic-like topology [61,62]; 3. Some nuclear responses are specific to pion production such as ∆-resonance inmedium modifications [63,64,65].
Moreover, in the energy region of MINERvA, the contribution of higher resonances could be important along with the nuclear medium effects in the pion production processes. MINERvA has measured the following processes on its plastic scintillator tracker with the Low-Energy beam: ν µ CC π + production (with limited contributions from π − ) [66,67], ν µ CC single π 0 production [68,69], -ν µ CC single π 0 production [70,67], and -ν µ CC single π − production [71].
The nuclear effects described by RPA and 2p2h are crucial to the description of quasielastic-like processes in MINERvA (Sect. 3.2). Since the RPA effect creates a suppression of the QE-like cross section at low Q 2 , measuring Q 2 in pion production could provide relevant information. However, as shown by the results in Fig. 11, suppression at low Q 2 ranges from being nonexistent (ν µ π − [71]), to mild or insignificant (ν µ π + [66,67] andν µ π 0 [70,67]), and to fairly pronounced (ν µ π 0 [68]). A combined fit of the various underlying mechanisms using a subset of these measurements can  be found in Ref. [72]. On the other hand, while 2p2h in pion production has not been incorporated into predictions, current models describe the MINERvA data with sufficient strength without it. The effect of initial-state correlations on pion production is not well understood.
Because of the granularity of the tracker and the event statistics obtained for ν µ charged-current proton-π 0 final states with the Low-Energy beam, MINERvA is able to investigate the transverse kinematic imbalances that may arise in pion production [69]. In this way, MINERvA can probe the initial state and FSI in parallel to the CCQE-like measurement discussed in Sect. 3.3 (Fig. 12). The TKI between the muon and the p-π 0 hadronic system is calculated by reconstructing the π 0 momentum and combining it with the proton momentum [73]. As with ν µ CCQE, this channel also has an initial-state neutron and as expected, the Fermi motion peaks in Fig. 12 (b) from both channels are consistent. The consistency in the p n tail size and in the trend of δα T in Fig. 12 is purely coincidental. As a specific example, if nature had less pion absorption, the QE-like p n tail and the δα T deceleration region would fall but the overall π 0 production would increase. Currently, for the same Fermi motion peak, generator model predictions do not describe both channels simultaneously [69], clearly illustrating the challenges inherent to consistent modeling of few-GeV ν-A interactions.

Kaon Production
Charged and neutral kaons are produced in both charged-current and neutral-current interactions at the beam energies of the MINERvA exposures. Kaon production is important to measure because the neutral-current (NC) channels yielding K + mesons constitute a background to proton decay searches [8,9]. This background arises from interactions of atmospheric neutrinos that are incident on proton decay detectors. The production mechanism is associated production with another strange meson or with a hyperon. Below the threshold of associated particle production, kaon production would take place through ∆S = 1 currents (S is the strangeness quantum number). At the energies of present interest, single kaon production may be important [74]. A finalstate K + with energy less then ∼600 MeV can stop inside the MINERvA tracker and then decay at rest. Once a K + decay chain is identified through the decay-associated time delay, the event topology, and the energy loss, the K + initial momentum can be calculated from track range. For neutral-current interactions, the final-state lepton does not identify the vertex location, consequently the 100-MeV tracking threshold sets a lower limit on the measured kaon kinetic energy. For CC interactions on the other hand, the final-state muon determines the primary vertex as the K + starting point, hence the kaon kinetic energy threshold can be lower since a kaon track needs not be identified. Figure 13 shows the measured charged-current [75] and neutralcurrent [76] cross sections. The analysis showed that the GENIE prediction was in agreement with the neutral-current measurement, while the NEUT prediction, which was used by Super-K [8], predicted a lower cross section by about 20 %.

Inelastic Reactions
In quasielastic and resonance production, the dynamic degrees of freedom are baryons and mesons. The invariant mass of the hadronic system, W , is of the order of 1-2 GeV/c 2 . In the high-W high-Q 2 region where DIS dominates, the target nucleon breaks up in the reaction and the QCD dynamics in the nuclear environment can be studied. This kinematic region is only accessible in inclusive measurements where any hadronic final state is allowed.
MINERvA measures CC inclusive cross sections using the muon kinematics. These measurements require one momentum-analyzed muon per event; no requirement is placed on the final-state hadronic system [77]. Figure 14 shows the p || -evolution of the inclusive p T spectra. In contrast to the CCQE-like measurement of Fig. 8 of Sect. 3.3, there are significant contributions from (GENIE) DIS in all p || bins and its relative contribution increases with p || . As in the CCQE-like case, there is an overall deficit of MnvGENIE at large p || .
Charged-current inclusive cross sections have also been measured in the passive target region as a function of the reconstructed Bjorken-x. These measurements are carried out for iron and lead targets; they are based on both muon kinematics and calorimetric hadronic energy [78]. The ratios of ν µ -Fe and ν µ -Pb cross section to the cross section in the tracker (CH) are shown in Fig. 15. The default GENIE simulations show a strong deficit in the elastic region (Bjorken-x ∼ 1) which might be due to the unmodeled 2p2h contributions.
By restricting the sample to the DIS region, W > 2 GeV/c 2 and Q 2 > 1 GeV 2 , these Bjorken-x-dependent cross section ratios can be compared to additional DIS models [79] (Fig. 16). As is shown in Fig. 16 (b) for Pb, the data at Bjorken-x < 0.1 suggests possible shadowing effects beyond those predicted. Because these models are only tuned to charged-lepton scattering data, which are by definition insensitive to the axial-vector current, these data also reflect neutrino-specific effects in the DIS region.

Coherent Interactions
Coherent pion production is a relatively rare process that is worthy of measurement because the neutral-current channel contributes a small but poorly constrained background to electron neutrino appearance measurements. Charged-current coherent pion production, ( -) ν µ + A → µ ∓ + π ± + A, can be considered as the inverse process of the pion decay permitted in a nuclear environment, in analogy to pair production in QED. The momentum transfer from the pion to the nucleus-denoted by t ≡ (P − P ) 2 where P and P are the final and initial 4-momenta of the nucleus-is small and leaves the nucleus intact, and the process has a characteristic exponential fall-off of the cross section with increasing |t|. By reconstructing t, MINERvA measured this process with both ν µ andν µ beams [82]. This measurement has provided a critical validation of the model implementation in generators (Fig. 17). Further model comparisons in a reanalysis [83] show preference for the Berger-Sehgal model [84] over the Rein-Sehgal model [85,86]   nuclear model has been used to study coherent pion production (see, for example, Refs. [87,88]).
The t-measurement was augmented with K + tagging (Sect. 3.5) in order to search for the analogous kaon production, ν µ + A → µ − + K + + A. MINERvA found 6 candidates with 1.77 predicted background events. This observation comprises evidence at 3.0 σ that CC coherent K + production does indeed occur [89].
A process that is the analog of neutrino NC coherent π 0 production on nuclei, is neutrino diffractive π 0 production on hydrogen: ( -) ν + H → ( -) ν + π 0 + H. This reaction relies solely on vacuum-quantum-number (Pomeron) exchange, hence diffractive [90]. In the past, diffractive π 0 production on hydrogen was not included in neutrino event generators. It has been identified by MINERvA as the cause of a data excess observed in neutral-current events containing electromagnetic showers in the tracker [91].

Conclusions and Outlook
With the NuMI Low-Energy data, MINERvA has investigated a wide variety of neutrino interactions in the GeV region of incident neutrino and antineutrino energies, including elastic scattering on electrons as well as neutrino-nucleus incoherent and coherent scattering processes. The reported measurements include utilization of neutrino scattering observations to constrain the flux, precision measurements of model parameters and model validation, and discoveries of novel processes. The results show neutrino-nucleus interactions to involve complex phenomena which challenge many of the current theoretical descriptions. Incremental improvements to the models have been essential to progress with MINERvA data and interpretation. The experiment's Medium-Energy data, analyses of which are currently underway, will enable comparisons of interaction channels on a range of nuclei at new levels of statistical precision, and expansion to kinematic phase space that has heretofore not been accessible. These data, consolidated by a preservation campaign [92,93], will not only facilitate future precision measurements of neutrino oscillations, but also further extend the knowledge of electroweak phenomena that lie at the intersection of nuclear and particle physics.