Rejecting cosmic background for exclusive charged current quasi elastic neutrino interaction studies with Liquid Argon TPCs; a case study with the MicroBooNE detector

Cosmic ray (CR) interactions can be a challenging source of background for neutrino oscillation and cross-section measurements in surface detectors. We present methods for CR rejection in measurements of charged-current quasielastic-like (CCQE-like) neutrino interactions, with a muon and a proton in the final state, measured using liquid argon time projection chambers (LArTPCs). Using a sample of cosmic data collected with the MicroBooNE detector, mixed with simulated neutrino scattering events, a set of event selection criteria is developed that produces an event sample with minimal contribution from CR background. Depending on the selection criteria used a purity between 50 and 80% can be achieved with a signal selection efficiency between 50 and 25%, with higher purity coming at the expense of lower efficiency. While using a specific dataset and selection criteria values optimized for the MicroBooNE detector, the concepts presented here are generic and can be adapted for various studies of exclusive νμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu _{\mu }$$\end{document} CCQE interactions in LArTPCs.


Introduction
Liquid argon time projection chambers (LArTPCs) [1] serve as both the target and the detection medium for several operating and planned neutrino oscillation experiments [2,3]. The main advantage of using liquid argon TPCs in neutrino experiments is their fine-grained spatial resolution and precise charge measurements that allow for a low-threshold, threedimensional reconstruction of charged particles and photons produced in neutrino interactions.
The final-state detection capabilities of LArTPCs make them excellent detectors for next generation neutrino oscillation experiments, such as searches for sterile neutrinos in short-baseline oscillations [2], and determination of the neutrino mass ordering and CP violating phase of the neua e-mail: hen@mit.edu trino mixing matrix in long-baseline oscillation measurement analyses [3]. Their superb final-state reconstruction capabilities also make them ideal for studies of neutrinonucleus interaction cross-sections that serve as an important input to many oscillation measurements and searches for new physics. In addition, liquid argon TPCs have been proposed as possible detectors for proton decay, neutrinos from corecollapse supernovae, solar neutrinos, the diffuse supernova neutrino flux [3], and dark matter [4].
The Short Baseline Neutrino (SBN) program at Fermilab is composed of a three-detector complex, all using LArTPCs [2]. MicroBooNE [5] is the first detector commissioned and has been taking data since October 2015. Its primary goal is to address observations of an excess of low-energy electronneutrino like neutrino interactions [6] that could give hints of physics beyond the Standard Model, such as sterile neutrino interactions. Achieving the goals of SBN and the future DUNE experiment [3] requires a precise understanding of neutrino-argon interaction cross sections, as well as detailed knowledge of nuclear final-state interactions (FSI) and other nuclear effects. The excellent particle reconstruction capabilities of LArTPCs allow them to detect neutrino interactions using exclusive final states that allow for the required understanding of the neutrino-nucleus interaction and impact of nuclear effects.
An example of such a process is charged-current quasielastic (CCQE) scattering [7] off a bound neutron in a nucleus. These are events in which the neutrino interacts and knocks out from the nucleus a single intact nucleon. This is one of the simplest lepton-nucleus interactions in the energy regime relevant for most neutrino oscillation experiments. Measurements of such interactions also allow benchmarking theoretical models of the neutrino interaction cross-section, needed for the precision extraction of oscillation parameters. Therefore, CCQE interactions are an appealing neutrinointeraction channel for use in high-precision oscillation measurements [8,9]. The three wire planes, shown on the right side of the cylindrical cryostat in the right-hand image, contain vertical wires and wires ± 60 • to the vertical. The wire planes are the anode A typical working definition of CCQE interactions in LArTPCs requires the reconstruction of a vertex with exactly one muon, one proton, and no additional particles. Such a classification can include contributions from non-CCQE interactions that lead to production of additional particles that are below the detection threshold or are absent from the final state due to nuclear effects, e.g., pion absorption; we refer to these as "CCQE-like" events. A significant challenge in studying such interactions stems from cosmic rays whose interactions in a LArTPC can mimic CCQE-like interaction topologies. This is especially significant for detectors located on the Earth's surface, where such cosmic ray (CR) interactions are vastly more abundant than neutrino interactions.
There are various ways by which cosmogenic background particles can mimic CCQE-like interactions, including: a single CR traversing the detector and being misidentified as two trajectories with a common vertex due to local non-active regions in the detector; cosmogenic neutron interactions in the detector that result in the production of multiple charged particles in the final state; muon decays that result in the production of Michel electrons; and CR scattering off an argon nucleus misidentified as two charged particles with a common vertex. Cosmic rays can also produce electromagnetic activity that can mimic electron-neutrino interactions and thus impact ν e measurements.
This article presents a collection of cosmogenic background rejection methods for studies of muon-proton pairs from CCQE-like neutrino interactions within a LArTPC. We use a mixed sample of cosmic data collected with the MicroBooNE detector and simulated neutrino interactions generated using the GENIE Monte Carlo (MC) event generator [10]. While focusing on muon-proton final states and using CR data collected by the MicroBooNE detector, the methods presented here can be adapted for other stud-ies using LArTPC detectors and different exclusive neutrino interactions, such as charged pion production, e.g. ν μ n → μ − Δ + → μ − nπ + , where only the muon and the pion emerge from the nucleus.
This article is organized as follows. In Sect. 2, we briefly present the MicroBooNE detector [5], describe its trigger and the effect on the accepted event rates, and present the event sample used in this study. In Sect. 3, we present a collection of event selection criteria, denoted as "cuts", for CCQE-like studies that are based on the detector response, and examine their impact on CR data and the simulated neutrino signal. In Sect. 4, we present complementary kinematic cuts based on known properties of CCQE interactions. Finally, in Sect. 5, we present the sequential impact of each cut on the CR background rejection and neutrino-interaction selection signal efficiency and purity.
We note that this work focuses on CR background rejection using real data. The impact of the event selection criteria on the simulated neutrino signal is shown as a reference.

Experimental setup
The MicroBooNE detector is a 170 tonne LArTPC with an active mass of 85 tonnes. It is located 463 m downstream from the production target of the Booster Neutrino Beam (BNB) at FNAL. The BNB energy spectrum ranges up to 2 GeV and peaks around 0.7 GeV. See [5] for details. The detector consists of a rectangular cuboid TPC with dimensions of 256 cm (width) × 232 cm (height) × 1036 cm (length), displayed in Fig. 1. The TPC includes three wire readout planes with coordinate system in whichẑ is along the beam direction, the drift directionx is horizontal, and theŷ direction is vertical. The origin (x = 0, y = 0, z = 0) is chosen at the anode plane, the detector vertical center, and the upstream side of the detector, respectively. The stainless steel cathode operates at −70 kV.
Neutrino interactions occur in the LArTPC when a neutrino from the beam interacts with an argon nucleus. The secondary charged particles produced in the interactions travel through the liquid argon, lose energy through ionization processes, and leave a trail of ionization electrons. These ionization electrons drift to the wire planes in an electric field of 273 V/cm. It takes about 2.3 ms for an ionization electron to travel the maximal drift distance from cathode to anode of 256 cm. Consequently, the triggered time window of the TPC is opened for slightly more than 4.8 ms (about twice the maximal drift time). This large trigger time window provides information on tracks exceeding beyond one full drift window, thus allowing proper reconstruction of cosmic tracks outside the beam time window. The wire planes are biased in such a way that drift electrons merely induce a pulse in the first two (U and V induction) planes and are collected in the plane farthest from the cathode, the collection plane. These signals are used to create three distinct two-dimensional views in terms of wire position and time [11].
Thirty-two photomultipliers tubes (PMTs) are placed outside the TPC volume facing the anode plane to collect the scintillation light from the de-excitation of argon dimers produced along with the ionization electrons. A TPB plate in front of each PMT absorbs the UV argon scintillation light and re-emits in a wavelength range to which the PMTs are sensitive. The PMT signals in coincidence with the BNB beam gate (or an external, off-beam window) is used as a trigger condition and to define the reference time for the event. An "event" in MicroBooNE consists of a continuous readout of the TPC and the associated PMT light during a 4.8 ms drift window.
During data-taking, a software trigger rejects events with light levels below that expected for a neutrino interaction. To enrich the recorded event sample with neutrino induced events in the beam window, we implement an additional trigger in the data-acquisition software during event building which scans the PMT waveforms within the beam and external windows. If the sum of all PMT waveforms exceeds a pre-set threshold equivalent to 6.5 photo-electrons (PEs), in any 100 ns time window within the 1.6 µs beam gate, the event is saved. To open this time window and sum the PMTs signals, a 0.5 PE signal is required to be recorded in at least one individual PMT.
At nominal beam intensity and without a PMT trigger, one in roughly 500 BNB beam spills is expected to contain a neutrino interaction. The PMT trigger provides a richer sample in which approximately 1 in 10 spills will contain a neutrino interaction. Even with this PMT trigger, further suppression of the CR background by one to two orders of magnitude is required for isolating neutrino interactions.
The cosmogenic muon rate in MicroBooNE is estimated to be 5.5 kHz [12], which corresponds to about 20 muons per TPC drift time window of 4.8 ms and constitutes the main source of background to neutrino interaction events.
To study the cosmogenic background, additional drift windows are recorded during periods between neutrino beam spills, or when the Fermilab accelerator is off. This acquisition window, being out of time with the neutrino beam, is referred to as the external trigger and is configured to be superseded by beam triggers in the event of overlapping acquisition windows.

The event sample used for this study
For CR background studies we use externally triggered data, during which the beam is off, as mentioned above. Each real CR data event is mixed with a simulated neutrino interaction event produced by GENIE [10]. The overlay is performed such that each simulation event is combined with a single external triggered event. Propagation of final state particles following the simulated neutrino interactions is simulated using GEANT4 [13]. Signals along the wires from ionization electrons are simulated using the LArSoft simulation [14].
The sample contains about 10 6 cosmic data events and the same number of simulated events, which is equivalent to that expected for about 10 21 protons on target (POT) of data. In MicroBooNE triggered data, fewer than 10% of the triggered events contain a neutrino interaction; the true mixture of CRneutrino data should be one simulated signal for about every ten beam triggered events. The external triggered events do not have the PMT trigger applied, and thus represent an unbiased sample of cosmic data.
We define the CCQE-like signal as events with a muon and a proton in the final state, both fully contained within the TPC, that originate from charged-current ν μ -Ar scattering. We look for events with exactly one proton with momentum greater than 200 MeV/c and any number of protons below 200 MeV/c; We denote these events as CC1 p0π . Our signal definition allows any number of neutrons at any momenta, any number of charged pions with momentum lower than 70 MeV/c (corresponding to a conservative requirement of about 3 cm track length in liquid argon, well above the reconstructability limit of the detector), and no neutral pions, electrons or photons (at any momentum). The minimal proton momentum requirement is due to its stopping range in LAr, which is about 6 mm, corresponding to two 3 mm wire pitches of the TPC.

Track reconstruction and cosmic background rejection prior to vertex building
Track reconstruction consists of three main stages: hit reconstruction, candidate CR track reconstruction (PandoraCosmic), and candidate neutrino-induced track reconstruction (PandoraNu) [15]. The first stage includes reconstructing individual hits on the TPC wires. In the second stage, the PandoraCosmic algorithm attempts to combine hits to construct as many candidate CR tracks as possible, identified by their geometric information, e.g. downward tracks. The hits associated with identified CR tracks are removed from further analysis. In the third step, the PandoraNu algorithm is used for track reconstruction from the remaining set of hits [15]. The impact of the PandoraCosmic removal pass was studied by comparing two simulated samples of cosmic events, with and without the application of the PandoraCosmic removal pass. We found it rejects about 20% of track pairs, induced by CRs, that might imitate neutrino interaction events.

Event selection based on detector observables
The LArTPC response to CCQE-like neutrino interaction is different from its response to pairs of reconstructed tracks induced from CRs, as the former produces a muon and a proton with typical momenta of a few hundred MeV/c primarily traveling horizontally, while the latter consists primarily of vertically traveling high-momentum muons. Consequently, we can reject CR backgrounds using information related to the energy of the tracks produced in CCQE-like neutrino interactions (i.e., CC1 p0π events), the distance between the tracks and the scintillation light produced, the length of the muon and proton tracks, and the location of vertices within a fiducial volume. This section describes these cuts and quantifies their impact on the cosmic data rejection and simulated neutrino signal loss.
The set of cuts presented in this section rely on low-level detection quantities, which we find are well-modeled by the MicroBooNE detector simulation [11]. As a result, these methods represent a robust prescription for isolating μp pairs primarily originating from neutrino-argon interactions.

Close track identification and reconstruction efficiencies
We consider pairs of tracks that are fully contained in the fiducial volume of MicroBooNE, defined by 3 < x < 253 cm, −110 < y < 110 cm, 5 < z < 1031 cm.
(1) To determine the desired proximity between the start points of the two tracks, we study the distribution of the reconstructed three-dimensional distance between the start points of the μ and p tracks in all simulated CC1 p0π events. The resulting distribution is shown in Fig. 2. The resolution for the start point of a track is on the order of 1 cm. We conservatively include vertices in which the muon and proton candidate track start points are closer than 11 cm, to minimize possible detector model dependence.
Out of the CC1 p0π simulated events contained in the active volume of the detector, about 50% have a muon track reconstructed and fully contained in the detector fiducial volume. Requiring the same also for the proton track, reduces this fraction to about 25%. An additional requirement of no other track reconstructed within 11 cm of the vertex, reduces it further to 20%. We denote their selection as the "preselection" stage. However, this preselection alone does not remove any "broken" CR trajectories, where a single particle from a CR is reconstructed as two tracks, typically characterized by a separation of the start points less than 11 cm.
The reconstruction efficiency for cosmic-induced trackpairs, as well as the efficiency in reconstructing an artificial pair of tracks at close proximity, were studied using a set of simulated cosmic events generated using the CORSIKA generator [16]. On average, for each event, there are 0.32 ± 0.02 pairs of contained PandoraNu CR tracks reconstructed with a common start point (closer than 11 cm). About 75% of these pairs originate from mis-reconstructed broken tracks, 5% from misconstruction of the trajectories of two distinct particles with intersecting trajectories, and 20% from misinterpreted muon-electron pairs due to Michel decays of CR muons. These "fake" vertices are reduced by the cuts discussed below.

Energy deposition profile of the proton and muon candidates
Particle identification (PID) can be used to discriminate broken muon tracks from real muon-proton pairs. PID is based on energy deposition. Particle identification in LArTPCs is usually done using calorimetry as measured by the energy deposition profile d E/dx along the track, as a function of its residual range (the distance of each trajectory point from the end of the track). PID is implemented here by comparing the measured d E/dx profile to the Bethe-Bloch expectation for a proton (or a muon) in liquid argon using a χ 2 test, where we normalize χ 2 to the number of hits in the track. We assume two tracks reconstructed at close proximity are a muon-proton pair and label the proton candidate as the tracks with the smaller χ 2 p (i.e. χ 2 compared to the proton expectation). Figure 3 shows the distributions of the χ 2 p of the proton candidate track ((χ 2 p ) p ) vs. χ 2 p of the muon candidate track ((χ 2 p ) μ ) for the mixed cosmic data and simulated signal sample.
As can be seen, muon-proton pairs coming from neutrino interactions populate a region where χ 2 p is very low (< 30), whereas CR pairs have a χ 2 p around 200. Most of the CR events are rejected due to the higher energy loss for protons then muons. Some are also rejected due to the (a) (b) (c) (d) Fig. 3 The distribution of χ 2 p of the proton vs. the muon candidate for two identified tracks in a close-proximity pair in MicroBooNE, for a cosmic data and b CC1 p0π simulated signal. The proton candidate is labeled as the track with the smaller χ 2 p . Also shown are the onedimensional projections of the distributions. A cut of 80 < (χ 2 p ) μ and (χ 2 p ) p < 30 to suppress CR background contribution is depicted in the figure. c The impact of a cut on (χ 2 p ) μ on the number of simulated signal CC1 p0π events lost and the number of CR background events rejected. d The impact of a cut on (χ 2 p ) p on the number of simulated CC1 p0π events lost and CR background events rejected electromagnetic debris, such as delta rays and shower fragments, which produce very little ionization in liquid argon. In addition, there is a population of events for which the calorimetric reconstruction has failed for at least one of the particles. We remove these events with a quality cut requiring (χ 2 p ) p > 0. Requiring 80 < (χ 2 p ) μ and (χ 2 p ) p < 30 suppresses the CR background by about a factor of 20, at the cost of losing about 34% of the CC1 p0π events.

Optical filtering
The data from the 32 PMTs spread in the y-z plane behind the TPC wire planes is also useful for CR rejection, as CRs typically produce less scintillation light, with a larger spread along the tracks, as compared to the more localized light produced in the vicinity of the neutrino interactions vertex. MicroBooNE event reconstruction produces a "flash", if: (1) at least three separate PMT pulses exceed the single PE (SPE) level within 30 ns, and (2) the sum of the three PMT pulses above SPE exceed five PEs. The flash includes the sum of all the PMT pulses integrated during an eight microseconds time window from the time of the first coincidence. Each event can include several such flashes.
We identify the flash associated with the vertex by comparing the expected PE yield in each PMT to the observed one. We use two selection criteria based on this flash: (1) the two-dimensional distance of the reconstructed vertex from the center of the associated flash in the y-z plane, d YZ , and (2) the total number of PEs recorded in this flash, N PE . Figure 4 shows the distribution of d YZ vs N P E for the mixed cosmic data and simulated CC1 p0π signal sample. The events shown have an identified vertex and pass the energy deposition profile cut discussed above. Also shown is a cut on d YZ < 200 cm and N P E > 150. The results of these cuts is a rejection of the cosmic contribution by about a factor of 2, with a signal loss of about 15%. We have adopted

Muon and proton track lengths
The reconstructed track length is used to further distinguish neutrino interactions from cosmic background. Pandora reconstructs segmented lengths of particle tracks, i.e., the accumulated length of multiple straight track segments, which do not necessarily sit on a straight line due to multiple coulomb scattering and other interactions. Figure 5 shows the distributions of the segmented length of the muon candidate, l μ , vs. the proton candidate, l p , after application of the cuts on the d E/dx profile and optical filtering described above. We expect that the muon tracks will be longer than the proton track with no cuts applied; this is indeed predicted to be the case in about 78% of the simulated CC1 p0π events when we identify the proton as the track with smaller χ 2 p . With the above requirements on the d E/dx profile and optical filtering applied, this is the case for about 91% of the simulated CC1 p0π events. Consequently, we require l μ > l p , which results in a signal loss of about 9% and leads to a correct identification of the muon and the proton in about 99% of the remaining pairs. The application of this requirement rejects about 20% of the CR background remaining after the previous selection criteria are applied.

Broken track removal
As mentioned in Sect. 1, some of the two-track events induced by cosmic rays originate from a broken track. This may be caused by particles traveling across a region with inactive wires in the TPC, a soft scattering of CR muons off the argon nuclei, or an inefficiency in hit reconstruction. The mitigation of these effects is possible using the collinearity of the two tracks. Events originating from broken-tracks can be characterized by having a three-dimensional angle θ 12 (a) (b) (c) Fig. 6 The distribution of θ 12 , in MicroBooNE, for a CR background data, and b CC1 p0π simulated signal. A cut of δθ 12 = |θ 12 − 90 • | ≤ 55 • , to suppress background contributions is depicted in the figure. c The impact of a cut on δθ 12 on the number of simulated signal CC1 p0π events lost and the number of CR background events rejected between the two tracks close to 0 • or 180 • . This angle is useful in identifying background caused by such tracks. Figure  6 shows the individual distributions of the different samples before and after the application of the d E/dx, flash, and track length selections. We require that |θ 12 − 90 • | ≤ 55 • , as depicted in the figure.

Event selection based on quasielastic kinematics
Any cut on the detector response is interconnected with the kinematics and the available phase-space of the reaction. However, the selection cuts described above do not rely heavily on the specific two-body kinematical signature of CCQE interactions. They lead to a CR suppression of approximately 99.5% and maintain a signal purity of about 50%. Using the two-body kinematics of CCQE interactions, such as vertex activity, coplanarity, and the imbalance of the transverse momentum, p T , allows further CR rejection and increase of the CC1 p0π purity. These variables, unlike the detector observables, are not modeled by the detector simulation, but rely primarily on the model-dependent neutrino interaction generator. This dependence is reduced by using relatively loose cuts and performing cut-sensitivity studies. On  N ticks ), and the ratio R ΔQ for each, in the two induction planes (U and V) and in the collection plane (Y). In a CC1 p0π event, R ΔQ is expected to be close to unity for a box with dimensions that do not encapsulate too much noise from non-neutrino induced background the other hand, the effect of the cuts on the CR background rejection is directly measured with data, and hence contains no such model dependence.
We introduce three relevant cuts that are based on CCQE two-body kinematics. The first removes events with a large energy deposit near the vertex, where the energy is not associated with the muon and proton tracks. The second is based on the expected coplanarity angle between the plane spanned by the neutrino and muon and that spanned by the neutrino and the proton. The third is based on the imbalance of transverse momentum of the final state particles relative to the incoming neutrino direction. These cuts can be applied independently or together for enhanced CR rejection by focusing on a specific part of the CCQE-like interaction phase-space.

Removing events with large additional energy deposits near the vertex
Events in which multiple particles are produced but only one muon and one proton track are reconstructed have large charge depositions in the vertex region, not associated with the reconstructed muon and proton tracks. Such events can be identified and removed to enhance the signal purity. The PandoraNu algorithm combines hits to form tracks, and associates charge deposition, ΔQ, with a particle trajectory, allowing a measurement of the fraction of charge deposition that is not associated with reconstructed tracks in the vertex region.
For each reconstructed vertex, we project the position of the vertex onto each of the three wire plane views, and define a sequence of boxes of increasing size centered on the vertex. We study the vertex activity R ΔQ as a function of the size of the box: Figure 7 shows a typical simulated CC1 p0π event in which R ΔQ is close to unity for a box with dimensions ranging from 20 wires × 40 time-ticks to 150 wires × 300 timeticks. The size of the box should be small enough not to capture too much irrelevant noise or background.
We optimize the dimensions of the box to maximize the difference between the distribution of CC1 p0π events and background. Since the R ΔQ objects are three-dimensional ((R U , R V , R Y ) where (R ΔQ ) U ≡ R U , and similarly for V and Y planes), we need to use information from all three views to compare the distributions. Using the "energy test" [17], we find that the optimal box dimensions are around 50 wires × 100 time-ticks, corresponding to about 15 × 5.5 cm 2 . Notice that the boxes are two dimensional (wires × time-ticks) in each plane.
For simplicity, we optimize a one-parameter selection in the space of the three ratios, the radius r R ΔQ , of a sphere To optimize the cut on r R ΔQ , we define a figure of merit equal to the product of the CC1 p0π purity ( p) and efficiency ( ), (4) Figure 8 shows the figure of merit defined in Eq. 4 for the CC1 p0π events as a function of r R ΔQ . As can be seen in Fig.  8, the product of purity and efficiency is maximal for r R ΔQ ≈ 0.43 and thus we choose r R ΔQ ≤ 0.43 for the selection of CC1 p0π events.
CC 1p 0 (MC) eff. = 38.6% pur. = 65.3% Fig. 8 The product of the purity times efficiency for CC1 p0π events as a function of r R ΔQ , given by Eq. 3. The dashed line indicates the chosen value of r R ΔQ . See text for details

Coplanarity requirement
For each track we define φ as the azimuthal angle in the x − y plane, and the azimuthal difference between the two tracks as Δφ = φ p − φ μ . Assuming two-body kinematics, the muon and proton tracks should lie on a mutual plane with the beam axis direction (z), i.e., the azimuthal angular difference between the outgoing tracks is expected to be 180 • , as illustrated in the insert of Fig. 9. In CCQE-like events the Fermi motion of the nucleon, nuclear re-scattering, and the resolution of the detector produce small deviations around 180 • . Hence, a requirement of a large Δφ between the tracks suppresses contributions from events with multi-hadron production. To utilize the coplanarity, we use the reconstructed start and end points of the candidate muon and proton tracks, and correct the directionality of the tracks to exit from the reconstructed vertex. Figure 9a shows the difference between the reconstructed and generated Δφ for simulated CC1 p0π events after the application of all the selection criteria discussed above. The standard deviation of the distribution is depicted in the figure, and serves as an estimate of the detector resolution for Δφ.
The prominence of coplanarity in ν μ CCQE-like events was used in an analysis of similar events by the MINERνA collaboration [18]. We find that the MicroBooNE coplanarity resolution is about 7 • (see Fig. 9a), about twice as large as the MINERνA reported value of 3.8 • but still sufficient for our purposes. This difference is primarily due to the difference in typical event kinematics, which results from the lower energies of MicroBooNE as compared to MINERνA. Figure 9b shows the distributions of the reconstructed Δφ between the two tracks, after we apply the detection-based selection criteria described in Sect. 3. Figure 9c

Transverse momentum imbalance
For ν μ n → μ − p CCQE scattering off a single neutron, with no nuclear correlations, the component of the total reconstructed momentum transverse to the incoming neutrino, should be small and mainly due to the Fermi motion of the knocked out neutron, final state interactions of the emerging proton, and momentum reconstruction resolution.
We estimate the momenta of the final state particles from the stopping range of the tracks in liquid argon. Figure 10a shows the difference between the reconstructed and the generated p T for CC1 p0π events, after application of all selec-(a) (b) (c) Fig. 10 a The difference between the reconstructed and generated (truth-level) p T , for simulated CC1 p0π signal events. A fit for a Gaussian distribution function around the peak is depicted in the figure, as well as the width σ of the best-fit result. b The distributions of the reconstructed p T between the two tracks, after the application of the detection selection criteria. c The effect of a selection criterion on the maximal reconstructed p T tion criteria. The standard deviation of the distribution is given in the figure. This σ is only used to verify that our chosen cuts are far from the measurement resolution. Figure  10b shows the distributions of the reconstructed p T , after applying all previous selection criteria. Figure 10c shows the effect of imposing a maximum p T on the different samples. To enhance the contribution from CC1 p0π and suppress background, we require that p T ≤ 0.35 GeV/c. Table 1 shows the sequential impact of each of the applied selection criteria discussed above. The original number of pairs in the simulated-signal and data-background samples, labeled as "preselection", includes all events that survive CR rejection in MicroBooNE as discussed in Sect. 2.2. Applying the detector-based requirements described above we retain 45.1% of the CC1 p0π simulated signal, while suppressing about 99.5% of the CR background. After further selection based on the kinematics of CCQE-like interactions, including Δφ and the reconstructed p T , the CR background is reduced to about 0.07% of the number of cosmic events in the original sample.

Cosmic rejection summary
The purity of the CC1 p0π simulated signal, listed in Table 1, is 78.4% for the 1:1 cosmic overlay assumption, and is computed by N CC1 p0π /(N cosmic + N beam ), where N beam = 12, 676 is the number of close-track pairs induced by the simulated neutrino interactions, of which 78.4% (N CC1 p0π = 10, 020, see last row of Table 1) are contributed by the CC1 p0π after application of the event-section requirements, and N cosmic = 104. As discussed in Sect.2.1, for the real MicroBooNE case the cosmic contribution is about 10 times larger after the application of the software trigger. Consequently, the final purity would change from about 78.4% to N CC1 p0π /(10 × N cosmic + N beam ) ≈ 73%.

Sensitivity to the selection criteria parameters
The results shown in Table 1 were obtained using specific cut values. As part of cut optimization, we vary each of the cut parameters based on an arbitrary 10% variation of the parameter, or its resolution estimated from simulation, whichever is larger. To study the variation of the final efficiency and purity, an ensemble of 1000 combinations of cut parameters was generated, each parameter chosen at random from Gaussian distributions with the following means (μ) and standard deviations (σ ): The CR background rejection factor and the corresponding simulated signal efficiency were computed for each randomly sampled parameters combination. Then, the standard deviation σ of the resulting distributions served as a measure for the sensitivity to the selection criteria parameters values.
The resulting sensitivity of the detector-based selection criteria combination, is σ = 0.1% for the CR rejection factor and σ = 4.0% for the signal efficiency. The sensitivity of the combination of all selection criteria, is σ = 0.03% for the CR rejection factor and σ = 3.5% for the signal efficiency. Table 1 Application of the detector observable and kinematic selection criteria on close-proximity tracks contained in the detector fiducial volume. The requirements in the first five rows are made on variables insensitive to the neutrino interaction model (see Sect. 3), and in the four last rows are the kinematic selection criteria (see Sect. 4). The numbers in parentheses are the fractions of the original studied samples retained after the application of each selection criterion sequentially. The purity given in the table is the relative purity for the sample under study in this analysis, in which every event is forced to include a neutrino interaction. This purity takes into account, in addition to CR background, additional beam-related background simulated by GENIE (which is not CC1 p0π , which is not discussed in the paper

Summary
The use of state-of-the-art LArTPC detectors allows measurement of the final state characteristics of neutrino-argon interactions with unprecedented detail. While their use in current and future neutrino oscillation experiments will enable a new view into neutrino physics, measurements from these detectors can be limited by significant cosmogenic backgrounds. Rejecting such backgrounds is particularly challenging for LArTPC detectors positioned on the Earth's surface, such as those to be used in the Fermilab SBN program. Using a sample of cosmic data collected by MicroBooNE overlaid with simulated neutrino interactions generated using GENIE, we present, for the first time, methods for CR background removal in exclusive CCQE-like neutrino interactions. The event selection criteria remove CR backgrounds based on detector observables and the kinematics of CCQE interactions. The net result is a suppression of about three orders of magnitude in the CR background, while retaining 50-25% of the simulated signal events with a signal purity of about 50-80%, depending on the application of detector level cuts or the addition of kinematical cuts. The choice of cuts depends on the analysis goals and efficiency-purity combination required to meet these goals. While our study uses cosmic data collected by MicroBooNE and simulated neutrino interactions generated using GENIE, the methods presented are generic and can be adapted to other experiments that use LArTPC detectors. For example, the cosmic rejection procedure presented here was developed with the aim of testing nuclear physics models of the most basic CCQE process in well-defined kinematics. For that, and similar purposes, given the available MicroBooNE statistics, the focus was put on achieving high purity of the selected events and the price in efficiency is tolerable. For other purposes, different purity-efficiency combinations can be obtained by adopting different combination of the cuts [19].
Implementation of the external cosmic ray tagger in MicroBooNE [20], and other hardware improvements, are expected to allow comparable CR rejection with looser cuts that should result in higher signal selection efficiencies and comparable purities.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .