Neutrino interaction measurements with the MicroBooNE and ArgoNeuT liquid argon time projection chambers

Precise modeling of neutrino interactions on argon is crucial for the success of future experiments such as the Deep Underground Neutrino Experiment (DUNE) and the Short-Baseline Neutrino (SBN) program, which will use liquid argon time projection chamber (LArTPC) technology. Argon is a large nucleus, and nuclear effects—both on the initial and final-state particles in the interaction—are expected to be large in neutrino–argon interactions. Therefore, measurements of neutrino scattering cross sections on argon will be of particular importance to future DUNE and SBN oscillation measurements. This article presents a review of neutrino–argon interaction measurements from the MicroBooNE and ArgoNeuT collaborations, using two LArTPC detectors that have collected data in the NuMI and Booster Neutrino Beams at Fermilab. Measurements are presented of charged-current muon neutrino scattering in the inclusive channel, the ‘0π\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi $$\end{document}’ channel (in which no pions but some number of protons may be produced), and single pion production (including production of both charged and neutral pions). Measurements of electron neutrino scattering are presented in the form of νe+ν¯e\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu _e+\bar{\nu }_e$$\end{document} inclusive scattering cross sections.


Introduction
One of the largest challenges faced by accelerator-based neutrino oscillation experiments is understanding the interactions of neutrinos with the target nucleus constituting the detectors. Neutrino-nucleus interactions are inherently challenging to model, and experiments must rely on data to tune models for use in oscillation analyses [1,2]. Even for a simple interaction, such as quasielastic scattering off a single nucleon in a nucleus, models must correctly consider the initial state motion and binding energy of that nucleon, as well as its interactions with the residual nucleus as it makes its way from an interaction point out of the nucleus (these are known as final-state interactions, or FSIs). Recent electron scattering measurements have shown the importance of correlations between nucleons inside a nucleus, further complicating the picture. For these reasons, many models have historically relied heavily on data [2]. Ab initio calculations that consider all nuclear effects are in development (e.g. [3]), which will reduce reliance on data, but nevertheless it is important that these calculations are compared to data. a e-mail: kirsty.duffy@physics.ox.ac.uk (corresponding author) A significant amount of neutrino interaction measurements have been made with carbon targets over the past decades, as well as some measurements on other similar sized nuclei; generators and models have been tuned to this data to sufficient precision for current oscillation measurements. However, for the future Deep Underground Neutrino Experiment (DUNE) [4] and the Short Baseline Neutrino (SBN) program [5] it is critical that our understanding of neutrino-nucleus interactions is extended to argon targets, since both will use liquid argon technology as the primary neutrino detection method. As argon is a significantly larger nucleus than carbon, nuclear effects generally have a larger impact on the kinematics of initial and final state particles so even where models can successfully describe data on carbon, it is important to verify that they scale to argon correctly. In addition to this, future oscillation measurements will require further reduced neutrino interaction uncertainties, requiring more precise measurements on argon than exist for carbon.
This article presents a review of neutrino interaction measurements by two experiments-ArgoNeuT and MicroBooNE-which together form the majority of all neutrino-argon interaction measurements globally. The energy ranges of these measurements are of particular relevance to DUNE and SBN: the ArgoNeuT measurements sample the primary DUNE oscillation Due to the need for higher precision, modeling of neutrino-nucleus interactions has been developed considerably over the past decade. Therefore the models used for developing analyses, as well as the models compared to, have some differences throughout this article. Both the experiments discussed used the GENIE monte carlo (MC) generator [6,7] as their primary simulation. The ArgoNeuT experiment used GENIE v2. 6 and GENIE v2.8 for most analyses. The majority of MicroBooNE analyses use GENIE v2.12 as the primary simulation. The largest difference between these versions is the addition of meson-exchange current (MEC) interactions (also known as 2p2h or multinucleon interactions) in v2.12. Data are additionally compared to GENIE v3.0 [8] predictions, which include a number of theoretical updates, as well as NuWro [9], GiBUU [10], and NEUT [1,11]. Care must be taken when interpreting comparisons of data to predictions, as the newer data are generally compared to newer simulations.
In addition to being the first measurements of various interaction types with an argon target, the use of liquid argon TPC technology means these measurements take advantage of the full 3D imaging, precise tracking, calorimetric energy reconstruction, particle identification, and low energy thresholds available. This has led to a number of measurements of exclusive topologies with detailed characterization of the full final state, leading to more stringent tests of models.

LArTPCs in action: ArgoNeuT and MicroBooNE
A comprehensive introduction to liquid argon time projection chamber (LArTPC) detector technology may be found in [12,13]. Several features of this technology make these detectors a versatile tool for the study of neutrino interactions; LArTPCs are effectively fullyactive tracking calorimeters. They bestow tridimensional imaging, tracking, calorimetry, and-in many cases-a self-triggering mechanism provided by scintillation light in liquid argon. Leveraging both the highresolution tracking and calorimetry, excellent particle identification capabilities have been achieved in LArT-PCs over a large kinematic phase space.
This article presents a review of neutrino interaction measurements in two LArTPC detectors: ArgoNeuT [14] and MicroBooNE [15]. ArgoNeuT was the first LArTPC detector in a low-energy neutrino beam: it was designed as a test experiment, but obtained physics results. ArgoNeuT's physics run lasted from September 2009 through February 2010, and featured a ∼ 0.25 ton active volume LArTPC, with dimensions of 47.5 cm (drift) × 40.0 cm (height) × 90.0 cm (length) exposed to Fermilab's Neutrinos from the Main Injector (NuMI) beamline [16]. The TPC had three anode planes, all with 4 mm wire pitch, with two planes instrumented for electronic readout. ArgoNeuT did not have a light collection system, and instead relied upon timing signals from the accelerator complex to trigger on NuMI beam spills. The detector was located underground in the NuMI tunnel at Fermilab, immediately upstream of the MINOS Near Detector (MINOS-ND) [17] which served as a spectrometer for muons exiting ArgoNeuT. The magnetized MINOS-ND detector allowed ArgoNeuT to make measurements separating muon neutrino and antineutrino interactions, and the underground location means that the cosmic background was negligible. Sitting on-axis with respect to the NuMI beamline, a data sample of 1.335 × 10 20 protons on target (POT) was collected in the "low-energy" NuMI configuration, composed of 8.5 × 10 18 POT in neutrino mode ( E ν = 4.3 GeV) and 1.25 × 10 20 POT in antineutrino mode ( E ν = 9.6 GeV, Eν = 3.6 GeV).
MicroBooNE [15] is the first detector of the SBN program to be fully operational, starting its first neutrino data-taking run in October 2015. The SBN program is comprised of three surface LArTPCs positioned on-axis in Fermilab's Booster Neutrino Beam (BNB) [18]: SBND (110 m from the decay pipe), MicroBooNE (470 m), and ICARUS (600 m) [5]. MicroBooNE is an 85 ton active volume LArTPC with TPC dimensions of 2.6 m (drift) × 2.3 m (height) × 10.4 m (length). The TPC has three anode planes, with 3 mm wire pitch, all instrumented for electronic readout (a total of 8256 wires in the detector readout). MicroBooNEs light collection system consists of 32 photomultiplier tubes (PMTs), each of 200 mm diameter, which register scintillation light produced in neutrino interactions to provide the start time for electron drift (allowing three-dimensional track reconstruction). Optical signals are also used to discriminate between TPC activity caused by neutrino interactions and cosmic rays. Because MicroBooNE sits on the surface, cosmic ray backgrounds are significant; the impact of these events is estimated using data collected when the beam is not running. In December 2017, a cosmic ray tagger [19] was installed, although the measurements presented in this article use only data collected prior to that date. The on-axis BNB neutrino flux has an average energy of ∼ 800 MeV, and is used for the measurements of muon neutrino interactions presented in this article. Micro-BooNE also detects neutrinos produced by Fermilab's NuMI beam [16]. The MicroBooNE detector is situated close to 8 • off axis with respect to the NuMI beamline, leading to a neutrino flux in the "medium energy" neutrino-mode configuration with a mean energy of 905 MeV and a composition of roughly 96% ν μ +ν μ and 4% ν e +ν e [20]. This enables high-statistics measurements of electron neutrino interactions. To date, Micro-BooNE has collected an exposure of 1.52 × 10 21 POT in the BNB flux and 2.33 × 10 21 POT in NuMI, although only subsamples of these data sets have currently been analysed (exposures for individual measurements are given in the following sections).

Measurements of ν μ and ν μ charged-current inclusive scattering
A charged-current (CC) inclusive event selection can be relatively simple, if it requires only the outgoing muon (μ + in the case ofν μ interactions) to be reconstructed. These event selections allow for overall tests of neutrino interaction models, as well as forming possible foundations for studies of more complex event topologies (for example, the MicroBooNE ν μ CC inclusive selection in Sect. 3.2 is used as a preselection for the ν μ CC0πNp measurement in Sect. 4.2.1). In terms of interaction modes, the inclusive signal definition includes many different types of neutrino interaction, with the exact mixture depending on the composition and energy of the neutrino beam flux being studied. The on-axis BNB neutrino flux at MicroBooNE is relatively low-energy, with a mean neutrino energy of 0.8 GeV. Therefore, the dominant interaction processes in MicroBooNE's CC inclusive measurement will be quasielastic (QE) and 2p2h scattering. ArgoNeuT, on the other hand, uses the on-axis NuMI low-energy beam, with a mean neutrino energy of 4.3 GeV (ν μ in the neutrino-mode beam), 9.6 GeV (ν μ in the antineutrino-mode beam), or 3.6 GeV (ν μ in the antineutrino-mode beam). Due to this higher neutrino energy, ArgoNeuT expects to see significant contributions from neutrino-induced resonance production (RES) and deep inelastic scattering (DIS) in the CC inclusive selection, in addition to quasielastic and 2p2h scattering.

ArgoNeuT ν μ andν μ CC inclusive cross-section measurement
References [21,22] report the first measurements of the inclusive muon neutrino and antineutrino charged current cross sections on argon using the ArgoNeuT detector in the NuMI beamline.
The initial muon neutrino measurement was made using a data sample of 8.5 × 10 18 POT, collected during a very brief ∼ 2 week run of ArgoNeuT when the NuMI beam was in the low-energy neutrino-mode configuration. Event candidates are required to have a candidate neutrino interaction vertex inside of a fiducial volume (3 cm < x vtx < 44 cm, −16 cm < y vtx < 16 cm, 6 cm < z vtx < 86 cm), with one reconstructed track originating from the vertex matched to a negatively-charged track starting within the first 20 cm of the MINOS-ND. Figure 1 shows the flux-integrated differential CC inclusive cross section measured in this analysis as a function of reconstructed muon momentum and scattering angle, compared against a prediction made using GENIE v2.6. In general, good agreement with the model was found, except potentially in the lowest bins of both muon momentum and scattering angle. The small statistical size of the data sample makes it difficult to draw further conclusions.
A subsequent measurement of both muon neutrino and antineutrino CC inclusive cross sections was made using a larger 1.25 × 10 20 POT data sample recorded by ArgoNeuT over several months when the NuMI beam was in the low-energy antineutrino-mode configuration. The sizeable neutrino component of the antineutrino-mode beam, combined with the capability of the MINOS-ND to perform sign selection of muons emanating from ArgoNeuT, enabled this comprehensive suite of measurements. The same fiducial volume definition and muon matching requirements were utilized in this analysis, with the addition that both negatively and positively charged MINOS-ND tracks were allowed, to differentiate between neutrino and antineutrino interactions respectively. Figure 2 shows the measured flux-integrated differential CC inclusive cross section for neutrinos (Left) and antineutrinos (Right) as a function of reconstructed muon momentum (Bottom) and scattering angle (Top), compared against predictions made using GENIE v2.8.0 and NUWRO version 11m (with a spectral function and "TE" meson exchange currents). As in the case of the neutrinomode analysis described above, though this time with larger statistics, generally good agreement is observed between the measurement and the models across the momentum and scattering angle distributions.

MicroBooNE ν μ CC inclusive cross-section measurement
Reference [23] reports the first measurement of the double-differential (in muon momentum and scattering angle) muon neutrino charged current inclusive cross section on argon, using the MicroBooNE detector. This measurement uses an exposure of 1.6 × 10 20 POT from the on-axis BNB, and selects 27,200 candidate events with an estimated ν μ CC interaction purity of 50.4%; this corresponds to the largest sample of neutrino interactions on argon published to date. The main background to this measurement (29.1% of all selected events) comes from cosmic rays, and is estimated using data.
The initial stage of the selection consists of cuts to remove backgrounds from cosmic rays using the timing of the interaction, correlations between optical and TPC information, and event topology. Muons are then identified using calorimetric information in the form of a truncated mean of the deposited charge per unit length (dQ/dx) and track length, to reject neutralcurrent interactions. The muon momentum is reconstructed using multiple Coulomb scattering [24], which allows a measurement of muon momentum for muons either contained in, or exiting the detector. This means that this measurement is able to probe the outgoing muon phase space with full acceptance in both angle and momentum.
The double-differential ν μ CC cross section is reported as a function of reconstructed kinematic variables, and obtained using a "forward folding" approach [23,25]. 1 Figure 3 shows the measurement compared to a number of neutrino interaction generator predictions. The   [22] largest disagreements between data and the predictions are in the high-momentum bins in the two most forward-going muon angular bins (0.86 ≤ cos θ μ ≤ 94 and 0.94 ≤ cos θ μ ≤ 1). In this region the estimated backgrounds are low, and the uncertainty on the measurement is small. χ 2 values are calculated between data and the predictions, taking into account the full covariance matrix with off-diagonal elements; the lowest χ 2 is found for the GENIE v3 prediction, χ 2 = 103.9/42 bins. GENIE v3 predicts a lower cross section in the forward direction than GENIE v2.12.10 (the model used for background simulation and efficiency estimation in this analysis); this is driven predominantly by adopting the local Fermi gas nuclear initial state model, and (to a lesser extent) the inclusion of the lishing matrices that can be used to apply detector resolution effects to predictions. random phase approximation (RPA) correction. These effects have the largest impact for low neutrino energies and heavy nuclear targets, so it is expected that these data may be more sensitive to these effects than measurements using lighter target nuclei. Further investigation of these effects will be pursued in two ways: through measurements of exclusive final states, aiming to detangle the impacts of different interaction processes and nuclear effects; and through improved measurement of the inclusive cross section with reduced systematic uncertainties (in particular, the uncertainties on detector modeling for all MicroBooNE analyses presented in this article are large, and have been the focus of significant development efforts by the collaboration [26][27][28][29]).

Discussion
While the MicroBooNE and ArgoNeuT measurements of the inclusive CC ν μ cross section on argon are performed at different neutrino energies-leading to a different underlying composition of the interaction channels-some similarities in the results arise. Overall, both experiments find good agreement with the generators. Interestingly, some tension is present in the forward-angle region of the MicroBooNE measurement (Fig. 3). However, while indications of a similar disagreement are seen in the earlier ArgoNeuT measurement ( Fig. 1, left side), the same tendency is not seen in the subsequent measurement (Fig. 2). It should also be noted that the measurements are compared to different generators (in particular, the GENIE versions used in the ArgoNeuT measurement do not include 2p2h models), so it is difficult to compare data-simulation agreement across these data sets.

Focus on protons: measurements of ν μ charged-current 0π scattering
Charged-current neutrino interactions that do not produce pions in the final state (commonly called CC0π interactions) are of particular interest to many neutrino oscillation experiments [4,5,30,31] because of their relatively simple topology and because they are the dominant interaction type at < 1 GeV neutrino energies. The CC0π interaction topology is expected to be composed largely of CCQE and CC2p2h interactions, as well as CCRES interactions that produce a pion which is subsequently reabsorbed in the nucleus through finalstate interactions. Both ArgoNeut and MicroBooNE have published measurements of CC0π cross section, with particular focus on measuring the protons produced in such interactions, making use of the LArTPC technology's low thresholds for reconstructing and identifying protons.

ν μ andν μ CC0π cross section
Results from the analysis of CC0π muon neutrino and antineutrino events in argon collected by the ArgoNeuT experiment on the NuMI beam at Fermilab are presented in Ref. [32]. Measurements of ν μ andν μ CC0π cross sections are reported for different proton multiplicities, using an exposure of 1.25 × 10 20 POT in the low-energy antineutrino-mode NuMI beam configuration. An analysis method based on the reconstruction Two-dimensional views of one of the ArgoNeuT "hammer events", with a forward going muon and a backto-back proton pair. Figure from [33] of exclusive topologies, fully exploiting the LArTPC detector's capabilities, is used to analyze the events and study nuclear effects in neutrino interactions on argon nuclei. CC0π events are characterized by the presence at the vertex of a leading μ − or μ + track, which may be accompanied by one or more highly ionizing tracks (although a proton candidate is not required), and the absence of pion tracks. With a proton reconstruction threshold in ArgoNeuT of 21 MeV kinetic energy, multiple protons accompanying the leading muon and the presence of vertex activity are clearly visible and measured.

Two-proton CC0π interactions
Reference [33] reports a measurement of chargedcurrent interactions in the ArgoNeuT detector which produced no pions and exactly two protons in the final state, using the full ArgoNeuT data set (8.5 × 10 18 POT in the low-energy neutrino-mode NuMI beam configuration, and 1.25 × 10 20 POT in the low-energy antineutrino-mode configuration). The kinematics of both protons in these events were studied, showing interesting characteristics. A fraction of the sample detected with ArgoNeuT is found with the two protons in a strictly back-to-back high-momenta configuration in the final state ("hammer events", see Fig. 4). Another equivalent fraction is found to be compatible with a reconstructed back-to-back configuration of a np pair in the initial state inside the nucleus. The detection of these two subsamples suggests that mechanisms directly involving nucleon-nucleon short range correlated pairs in the nucleus are active and can be efficiently explored in neutrino-argon interactions with LArTPC technology. However, the event statistics from ArgoNeuT are very limited and cannot provide definitive conclusions.

MicroBooNE measurements of CC0π scattering
MicroBooNE has reported two measurements of CC0π interactions in the BNB neutrino flux: a measurement of interactions that produce one muon, one or more protons, and no pions in the final state [25], and a more exclusive measurement of CCQE-like scattering [34]. Both measurements focus on similar topologies (including muons and protons, and excluding pions), but have different signal definitions and selections, leading to different expected proportions of CCQE, CC2p2h and CCRES interactions in the selected events, and complementary constraints on the interaction models. Both measurements are made in detector-measurable quantities such as individual particles' momenta or directions in order to minimise model dependence. Both measurements have a proton detection threshold of 300 MeV/c (47 MeV kinetic energy), which is driven by a combination of reconstruction algorithm failures at low momenta, and particle identification (PID) failures for short tracks. In both measurements, the PID uses exclusively the (vertical) collection plane, limiting its efficiency at higher proton angles; a limitation that will be eliminated in future measurements by using all three wire planes for PID.

Charged-current pionless interactions with protons
Reference [25] reports a measurement of charged current scattering with protons and no pions in the final state (referred to as "CC0πNp". The signal definition includes all events that produce protons within the detector's acceptance and don't produce pions. Though primarily composed of quasi-elastic interactions, this signal definition is sensitive to the presence of 2p2h interactions and resonant events with absorbed pions, as well as any other potential interactions that produce protons and not pions. This measurement used a data sample corresponding to 1.6 × 10 20 POT. The event selection begins with an inclusive muon neutrino charged-current selection from the analysis described in Sect. 3, which removes most backgrounds and identifies a muon candidate. All non-muon clusters are then required to be proton-like (using a particle identification score PID prot defined by comparing the reconstructed dE/dx to the expected energy deposition of a stopping proton [25]). Secondary proton candidates with fewer than five hits on the collection plane are assumed to be protons, as this particle identification method becomes inaccurate with fewer than five hits.
Although the signal definition allows for any number of protons, the measurement focuses on the most energetic (leading) proton. The measurement follows the same forward-folding strategy as [23], whereby predictions must be smeared to compare to the data, and measurements are made of five kinematic variables: the momentum and scattering angle of both the muon and leading proton, as well as the opening angle between those two particles.  Figure 5 shows the measurement of the muon direction compared to a number of generator predictions. Echoing the observations of the CC inclusive measurement in Sect. 3.2, the data in this more exclusive CC0π signal definition also show a suppression at the most forward muon directions. Figure 6 shows comparisons of the same generators across the other four variables. With the exception of the muon momentum distribution, the other variables show reasonable agreement with all models, although tensions can also be seen at the lowest proton momentum.

Charged-current quasielastic-like scattering
Reference [34] presents a more exclusive measurement of CCQE-like charged-current muon neutrino scattering on argon with a muon and a single proton in the final state (referred to as "CC1p0π"). The CCQE interaction process can be reasonably well-approximated as a two-body interaction, and these interactions are the focus of many accelerator-based oscillation studies. For this reason, precise measurements of CCQE processes would be valuable inputs to oscillation measurements. However, FSIs make it impossible to measure CCQE processes directly: the proton produced in a CCQE interaction could be reabsorbed in the nucleus, or additional particles (nucleons or pions) could be produced as a result of proton reinteractions. Other processes can also imitate the CCQE signal of one muon and one proton due to FSIs, for example if an additional nucleon or pion is absorbed in the nucleus before exiting. For this reason, Ref. [34] defines 'CCQE-like' events as the signal for this measurement, in which there is one outgoing lepton and exactly one outgoing proton with momentum greater than the reconstruction threshold of 300 MeV/c (kinetic energy of 47 MeV). Additional cuts are placed on the signal phase space to further enhance the contribution of CCQE interactions (|Δφ μp − 180 • | < 35 • , and |p μ T + p p T | < 350 MeV/c, where p T is the reconstructed momentum transverse to the beam) and to keep only regions of phase space where the selection efficiency is high and well-understood (0.3 < p p < 1.0 GeV/C, cos θ p > 0. 15 This measurement was performed using the on-axis BNB neutrino flux, with an exposure of 4.59 × 10 19 POT. CCQE-like candidate events are selected by identifying pairs of tracks, of which one is classified as a muon candidate and the other as a proton candidate using the same PID prot as in Sect. 4.2.1. Further restrictions are placed on the reconstructed kinematics of the muon and proton, to match the signal phase space. Figure 7 shows the flux-integrated single-differential CCQE-like cross section as a function of the cosine of the measured muon scattering angle. In general the measurement finds good agreement between the data and predictions from most of the generators, except in the most forward-going bin (0.8 < cos θ μ < 0.95); the overprediction at forward-going muon angles shown in Sects. 3.2 and 4.2.1 is also seen in this more exclusive measurement. The fact that the overprediction grows with increasing exclusivity is indicative that this effect is due to nuclear effects impacting quasielasticlike interactions, either in true quasielastic interactions, or in other interactions that mimic the signature. Figure 8 shows the flux-integrated single-differential cross sections as a function of muon momentum, proton momentum, and cosine of the proton scattering angle. Results are shown both for the full phase-space of cos θ μ and for cos θ μ < 0.8. The data-prediction agreement is improved across all variables when excluding the most forward-going muon bin.

Discussion
LArTPC technology allows the reconstruction of neutrino interactions producing protons in the final state with never-before-accessible low energy thresholds, in high statistics datasets. This is a powerful line of research when testing nuclear models in neutrino generators, as the multiplicity and kinematics of protons in the final state can be impacted by several nuclear mechanisms, such as short range nuclear correlations and pion nuclear absorption, which have large theoretical uncertainties. The ArgoNeuT measurements already highlight how proton multiplicity is a thorny territory for neutrino generator predictions.
The MicroBooNE measurements show tension in the lepton forward angle, in which the models overpredict the measured data. A similar tension was also observed in the CC inclusive measurement presented in Sect. 3.2, and this is seen to be more substantial in the more exclusive measurements (with the largest tension in the most exclusive CCQE-like measurement). The high cos θ μ phase space corresponds to low momentumtransfer events, so this may indicate a mismodeling of signal or background events in that phase space (the more exclusive measurements rely more on simulation  [25] to subtract backgrounds). These results justify further investigation into the behaviour of cross sections at low four-momentum transfer, and future measurements may learn more about this region by reducing the proton detection threshold.
ArgoNeuT's two-proton channel measurement is particularly interesting for the implications of exploring short-range correlations in neutrino scattering, and this will be extensively studied in a future dedicated search with MicroBooNE. Further investigation of CC0π topologies will also include analysis of transverse imbalance variables and other derived quantities that are more sensitive to nuclear models.

ν μ charged pion production
Charged-current pion production is an important process in few-GeV neutrino-nucleus interactions, where the underlying interactions mainly proceed through the excitation of nucleon resonances, deep inelastic scattering or coherent scattering of the neutrino on the entire nucleus. A precise measurement of the observed charged pion rates and understanding of the underlying physics processes is crucial for future neutrino oscillation measurements, particularly at DUNE [4] where single pion production is expected to be the dominant interaction process.

ArgoNeuT measurement of charged-current single pion production
Reference [35] reports the first CC1π ± differential cross section measurement on argon using the ArgoNeuT detector. The signal is defined to be a charged-current ν μ orν μ interaction in the detector, with only one charged pion above 100 MeV/c momentum exiting the target nucleus. The measurements are based on data taken with the NuMI beam line operating in the lowenergy antineutrino-mode configuration, corresponding to 1.25 × 10 20 POT. The event selection requires that there is at least one track in ArgoNeuT that is matched to a track in the MINOS detector. The matched track is identified as the muon track. A boosted decision tree (BDT) is trained using the calorimetric and topological information in the TPC to select events with a single charged pion in the final state. The selection efficiency is approximately 30% for both the ν μ andν μ CC1π ± samples. The purity is 36% (56%) for the selected ν μ (ν μ ) CC1π ± sample. Figure 9 shows the ArgoNeuT ν μ andν μ CC1π ± differential cross sections compared to four neutrino generators: GENIE, NuWro, GiBUU and NEUT. The GiBUU predictions are in a reasonable agreement with the measured cross sections. The NuWro, NEUT and GENIE predictions are higher than the measured cross section reported in the paper.

ArgoNeuT measurement of coherent charged pion production
Neutrinos can produce single pion final states by coherently scattering from the entire nucleus. Both neutral current and charged current coherent-scattering processes are possible. In these interactions, the squared four-momentum transfer to the target nucleus, |t|, is small so the nucleus remains unchanged. Reference [36] reports the CC coherent pion production from muon neutrinos and antineutrinos on argon using the ArgoNeuT detector: (1) where the low |t| condition entails that the pions and muons are forward going with respect to the incoming neutrino direction. The measurements are based on data taken with the NuMI beam line operating in the low-energy antineutrino-mode configuration, corresponding to 1.25 × 10 20 POT. The event selection requires exactly two tracks consistent with being minimum-ionizing particles (MIPs), with one track matched to a track in the MINOS detector as the muon candidate. The lack of any particles other than the muon and the pion emerging from the vertex is further reinforced by rejecting events with high charge ArgoNeuT νμ andνμ CC1π ± differential cross sections compared to GENIE, NuWro, GiBUU and NEUT. Thick error bars refers to statistical errors while thin error bars refers to statistical and systematic errors summed together. Figure from [35] deposition around the vertex. A BDT is trained using the calorimetric and kinematic information to select events with low |t|. The selection efficiency is approximately 20% for both the ν μ andν μ coherent charged pion samples. The purity is close to 100%. A comparison between the ArgoNeuT's CC coherent pion cross section measurements, existing data, and the predictions from GENIE and NuWro is shown in Fig. 10. The antineutrino measurement agrees well with the Rein-Seghal model [37] while the measured neutrino cross section deviates by ∼ 1.2σ.

ν μ neutral pion production
Neutral pions produced in neutrino interactions generally decay to two photons inside the detector. Photons, like electrons, produce electromagnetic showers in LArTPC detectors; while photon-induced showers can be distinguished from electron-induced showers using characteristics such as distance between the neutrino interaction vertex and the shower start and the energy deposited per unit length, dE /dx (as discussed in Sect. 6), photons from π 0 decay may still form important backgrounds to ν e appearance oscillation searches due to reconstruction failures. Therefore, it is important that these interactions be understood. Because of the well-known mass of the π 0 , neutral pions can also form a "standard candle" for tests of the shower energy resolution in LArTPC reconstruction.

Fig. 10
ArgoNeuT's CC coherent pion cross section measurements compared to model predictions and measurements from the SKAT [38], CHARM [39] and MIN-ERvA [40] experiments. These measurements are scaled to Argon assuming the A 1/3 dependence from the Rein-Seghal model. Figure from

ArgoNeuT measurement of neutral-current neutral pion production
Reference [41] reports ArgoNeuT's measurement of the cross-section for neutral current (NC) π 0 production in ν μ -argon interactions, using an exposure of 1.25 × 10 20 POT in the low-energy antineutrino-mode NuMI beam Fig. 11 The NC π 0 cross-section measured in ArgoNeuT compared with predictions from GENIE (solid) and NuWro (dashed), Left: for neutrino (red) and antineutrinos (blue) separately, and Right: for a combined νμ +νμ measurement. Figures from [41] configuration. The event selection requires the presence of two reconstructed electromagnetic shower candidates consistent with a π 0 → γγ decay. Consistency with the γ hypothesis is obtained by requiring that at least 75% of the first 4 cm of each shower candidate presents a dE/dx value greater than 3.5 MeV/cm. The selection further requires the absence of additional reconstructed tracks to reduce the contribution of charged-current interactions. Candidate interactions are required to originate within a fiducial volume that is the same as that defined for the CC inclusive analyses described in Sect. 3.1. Given the very small size of the ArgoNeuT detector almost none of the candidate showers are fully contained; thus, the analysis develops a series of model-driven corrections leveraging the angular and deposited energy information reconstructed in the contained part of the shower to estimate the true energy of each shower. The analysis is further complicated by the impossibility to distinguish between neutrino-and antineutrino-induced NC π 0 production-both sizeable contributions in the antineutrino-mode NuMI beam. As such, the final result-shown in Fig. 11-reports both a combined ν μ +ν μ cross-section (right) as well as separate neutrino and antineutrino cross-sections that assumed a fractional composition in the data that matched that of the GENIE (v2.8.0) model.

MicroBooNE measurement of charged-current neutral pion production
While neutral-current π 0 production forms the majority of backgrounds to ν μ →ν e oscillation searches, the spectrum and rate of π 0 production can be sculpted by FSI in the nucleus, so a good understanding of the impact of this is needed to estimate this background. To first order, assuming isospin symmetry, the interaction of a neutral pion with the residual nucleus is the same for both the neutral-current interaction mode an the charged-current interaction mode, so charged-current interactions that produce a π 0 are a useful way of con-straining the impact of FSI on this important background. Reference [42] reports the first measurement of ν μ CC single π 0 production on argon, using the Micro-BooNE detector and a BNB exposure of 1.62 × 10 20 POT. The measurement is performed using an inclusive final-state topology of at least one photon from a π 0 decay exiting the nucleus, in addition to a muon, and any number of other particles. Events are selected based on the presence of reconstructed showers, which are assumed to be produced by photons, and a reconstructed muon track. Tracklike topologies are rejected, to reduce cosmic ray contamination, although this leads to a reduced showerreconstruction efficiency for photons with kinetic energies below 50 MeV (which can appear more track-like). This has an impact on the selection efficiency, because the energy distribution for sub-leading photons produced by the decays of π 0 s at BNB neutrino energies peaks below 50 MeV. When requiring two showers to be reconstructed, the selection has an estimated purity of 64% and efficiency of 6%; when requiring at least one reconstructed shower, the selection has an estimated purity of 56% and efficiency of 16%. The result, calculated using the single-photon selection and shown in Fig. 12, is a flux-integrated total cross section for semiinclusive CC single π 0 production on argon. A slight deficit (1.2σ) is seen relative to the GENIE prediction for argon. The measurement is also compared to measurements of exclusive CC single π 0 production in the MiniBooNE [43] and Argonne National Laboratory (ANL) bubble chamber [44] experiments; the reasonable agreement of all three measurements with model predictions indicates that the nuclear scalings implemented in these models are applicable, within uncertainties, for neutrino-argon scattering. Total flux integrated νμ CC single π 0 production cross section measured in the ANL bubble chamber [44], MiniBooNE [43], and MicroBooNE [42]. These are compared to the GENIE v2.12.00 prediction with the Rein-Sehgal model (RS) [45], the RS model with FSI removed (No FSI), an alternative GENIE prediction using the Berger-Sehgal model (BS) [46], and the NuWro v18.02.01 default prediction. Figure from [42] producing pions in the final state in terms of signal definition. ArgoNeuT explored the CC1π ± channel, CC coherent charged-pion production, and the NC production of π 0 , while MicroBooNE has focused on the production of π 0 in CC interactions. The ArgoNeuT result for the CC1π ± channel is particularly interesting as it presents tension with almost all generators in the low momentum and low angle region for both ν μ andν μ . When exploring coherent pion production, ArgoNeuT measurements for the antineutrino data are in good agreement with the Rein-Seghal model for neutrino excitation of baryon resonances but present a slight deviation in the neutrino data. Both coherent and resonant CC1π ± channels will be explored further in MicroBooNE with dedicated high-statistics searches. Hints of tension can also be found in the MicroBooNE's measurement of the total flux integrated ν μ CC single π 0 production cross section, where the generators' prediction tends to be higher than the measured value (although within estimated uncertainties).

ν e andν e scattering
The measurement of electron neutrinos appearing in a muon-neutrino beam is essential to measure the neutrino mass ordering [47] and to search for CP violation in the neutrino sector [4,48] in long-baseline experiments. In short-baseline experiments, the analogous signature has the potential to reveal the existence of sterile neutrinos [5].
Generally, accelerator neutrino oscillation experiments use ν μ interactions to constrain the ν e flux and cross section models [18,49], mitigating the lack of pre-cise electron-neutrino cross section measurements. Any uncertainty on the ν e /ν μ cross section ratio weakens the reliability of such a procedure. In principle under the assumption of lepton universality, the only difference between ν e and ν μ interactions is in the final state lepton mass. However, that different mass also leads to differences in the cross section through radiative corrections, as well as through changes to form factors. The impact of radiative corrections is expected to be approximately 10% [50]. Additionally, nuclear effects in general depend on the energy and momentum transferred to the nucleus, which will be different due to the lepton mass, and therefore uncertainties in those nuclear effects can increase the uncertainty on the ν e /ν μ cross section ratio. Recent theoretical calculations predict large differences between the ν e and ν μ cross sections, particularly for forward-going leptons in the sub-GeV range [51]. Independent direct measurements of electron-neutrino cross sections on the target medium of choice for future neutrino oscillation experiments are therefore of paramount importance to further inform our understanding of different flavor neutrino interactions and perform oscillation searches to the desired precision.

ArgoNeuT measurement of inclusive CC ν e and ν e cross section
Reference [52] reports measurement of the combined ν e +ν e flux-averaged charged-current inclusive cross section performed using the on-axis flux from the NuMI beam in the low-energy antineutrino-mode configuration, corresponding to 1.25 × 10 20 POT. In the small ArgoNeuT detector, showers are rarely completely contained; this prohibits the use of total charge for complete shower characterization and energy reconstruction. Instead, charge ratios characterizing the shape and evolution of each candidate electron shower are used to select electrons. A BDT is constructed using the charge ratios as well as the shower angle and vertex dE /dx. The event selection is expected to have a 10.5% efficiency with a 78.9% purity based on studies of simulation. A flux-averaged ν e + ν e total cross section and a differential cross section as a function of lepton angle are extracted using 13 ν e and ν e events identified with fully-automated selection and reconstruction. A total cross section of (1.04 ± 0.38 (stat.) +0. 15 −0.23 (syst.)) × 10 −36 cm 2 on argon ( E νe = 10.5 GeV and Eν e = 4.3 GeV) was extracted, consistent with the GENIE v2.12.10 expectation. Figure 13 shows the ArgoNeuT ν e + ν e CC differential cross section for electron/positron angle with respect to the neutrino beam compared to the GENIE prediction. These are the first measurements of electron neutrino scattering cross sections on argon, which provide a good validation of the shower reconstruction technique and the modeling of electron neutrino cross sections using the same target nucleus and over the same energy range that will be used by DUNE. Fig. 13 The ArgoNeuT νe + νe CC differential cross section for electron/positron angle with respect to the neutrino beam compared to the GENIE prediction. Figure from [52]

MicroBooNE measurement of inclusive CC ν e andν e cross section
Reference [20] reports measurement of the combined ν e +ν e flux-averaged charged-current inclusive cross section performed using the off-axis flux from the NuMI beam in the medium energy neutrino mode, as measured in the MicroBooNE detector. The integrated neutrino flux above 250 MeV has an average of 0.905 GeV. The measurement selects 214 candidate ν e +ν e interactions with an estimated exposure of 2.4 × 10 20 protons on target. Given the estimated purity of 38.6%, this implies the observation of 80 ν e +ν e interactions in argon.
The analysis follows a staged approach for background removal. The initial stages of the selection focus on reducing cosmic ray contamination, by a combined strategy that leverages the timing of the interaction, correlations between optical and TPC information, and quality requirements on the reconstructed interaction vertex. The identification of a reconstructed electromagnetic shower abates backgrounds from ν μ CC interactions. The characteristics of the leading shower such as dE/dx, Molière angle, and distance from the reconstructed neutrino interaction vertex are leveraged to select ν e CC inclusive interactions. Indeed, the analysis includes a demonstration of a fully-automated dE/dxbased particle discrimination technique for electronand photon-induced showers in a LArTPC neutrino detector, as well as a study of electron and photon separation. The example dE/dx distribution for electromagnetic showers reported by the analysis can be found in Fig. 14.
While the selection accomplishes an overall decrease in the cosmic ray contamination of a factor of 10 5 compared to the initial sample, the main background for this first ν e analysis in MicroBooNE remains cosmic ray contamination (42% of all selected events). Significantly higher purity is expected in underground detectors, as well as with next-generation reconstruction algorithms where cosmic rays are rejected via improved light and charge matching, topology reconstruction and CRT cosmic ray tagging. The leading uncertainties for the measurements are the detector modeling (23%), the modeling of hadron production for the flux prediction (22%) and statistical uncertainty (22%); all are expected to reduce with future iterations.

Discussion
MicroBooNE and ArgoNeuT both measured the inclusive CC ν e +ν e cross section on argon analyzing neutrinos from the NuMI beam. Despite using the same beam, there are significant differences in the measured neutrino flux. The average flux energy is much higher for the ArgoNeuT measurement given its position onaxis to the NuMI beam. As a consequence, the dominant interaction channels predicted in for the inclusive channel are DIS and RES for ArgoNeuT, compared to predominantly CCQE for MicroBooNE. The neutrino to antineutrino flux ratio is also different given the configurations of the NuMI beam for the data employed: antineutrino mode for ArgoNeuT and neutrino mode for MicroBooNE. Both experiments find similar agreement to the theory predictions for the flux-averaged cross section: the agreement between the GENIE v2.12.2 used in MicroBooNE and the GENIE v2.12.10c used in ArgoNeuT is equivalent. Following the first ArgoNeuT measurement of differential electron neutrino cross section in lepton angle, a more detailed exploration of inclusive and exclusive electron neutrino interaction channels is expected from MicroBooNE's higher statistics data sample.

New analysis techniques
Future LArTPC measurements will continue to develop analysis techniques that take advantage of the detectors' low thresholds and high-resolution reconstruction. Future measurements from MicroBooNE will benefit from improved cross-section, detector, and cosmic ray modeling, and associated reductions in systematic uncertainties (particularly those associated with the detector model). This will allow more precise measurements of neutrino interactions, including detailed measurements of hadronic final states and variables not   [58] previously measured on argon such as those associated with Transverse Kinematic Imbalance [53][54][55][56].
Extending the LArTPC's range of physics sensitivity down to the sub-MeV level is crucial for future studies of supernova, solar, and beam neutrino interactions. Furthermore, including charge deposits from deexcitation photons in GeV accelerator neutrino events in LArTPCs is expected to improve neutrino energy resolution by about 25% [57]. ArgoNeuT has developed new techniques for low-energy LArTPC reconstruction for this purpose [58]. The total number, reconstructed energies, and positions of isolated low-energy depositions in the TPC volume were compared to those from simulations using the FLUKA Monte Carlo generator [59][60][61]. These measured features are consistent with energy depositions from photons produced by de-excitation of the neutrino's target nucleus and by inelastic scattering of primary neutrons produced by neutrino-argon interactions (see Fig. 15). Similar studies of low-energy reconstruction are underway in MicroBooNE.
Studies of low-energy new physics scenarios, such as millicharged particles [62], light mediators [63], and inelastic scatterings with small splittings [64], could also invaluably profit from such low-energy reconstruction. A search for millicharged particles, a simple extension of the standard model, has been performed with the ArgoNeuT detector [65].

New detectors
As well as new analysis techniques, new detectors are under construction that will add to the landscape of neutrino-argon cross section measurements. In the next few years, SBND will collect data in the same beamline as MicroBooNE, but will collect an order of magnitude more interactions. Expected statistics at SBND are 1.5 million CC-inclusive events, as well as 12,000 electron neutrino interactions [5]. This huge increase in statis-tics will allow multi-dimensional cross section measurements of a range of exclusive final states, including rare final states such as those including kaons and hyperons. Additionally, the ICARUS detector will make measurements of two GeV-scale neutrino beams over the next few years (the BNB on-axis and the NuMI beam off-axis), which will also include approximately 10,000 electron neutrino interactions in the NuMI dataset [5]. Further in the future, the DUNE near detector complex will include a liquid argon detector and a gaseousargon detector in a wide-band beam peaked at 4 GeV, expected to collect 30 million and 1.6 million ν μ CCinclusive events respectively per year in neutrino-mode beam running [66]. These data sets will be crucial to reach the precise understanding of neutrino-argon interactions needed for DUNE's oscillation program.

Conclusions
This article has summarized existing measurements from the ArgoNeuT and MicroBooNE collaborations, demonstrating the potential for LArTPCs to provide new insight into neutrino-nucleus interactions. The capabilities of this detector technology to study finalstate particle properties spanning a wide kinematic phase space with exquisite resolution opens up new avenues of inquiry, while the use of a nuclear target like argon requires improvements in modeling and analysis techniques to make precise measurements. Future highstatistics data samples collected by the SBN program, and the DUNE Near Detector, will be able to expand on the measurement techniques described in this article and further exploit the LArTPC technology to make exciting measurements of neutrino-nucleus interactions with improved precision.