## First Measurement of Neutrino and Antineutrino Coherent Charged Pion Production on Argon

R. Acciarri, 1 C. Adams, 2 J. Asaadi, 3 B. Baller, 1 T.Bolton, 4 C. Bromberg, 5 F. Cavanna, 2, 6 E. Church, 2 D. Edmunds, 5 A. Ereditato, 7 S. Farooq, 4 B. Fleming, 2 H. Greenlee, 1 R. Hatcher, 1 G. Horton-Smith, 4 C. James, 1 E. Klein, 2 K. Lang, 8 P. Laurens, 5 R. Mehdiyev, 8 B. Page, 5 O. Palamara, 2, 9 K. Partyka, 2 G. Rameika, 1 B. Rebel, 1 E. Santos, 10 A. Schukraft, 1 M. Soderberg, 3, 1 J. Spitz, 2 A.M. Szelc, 2 M. Weber, 7 T. Yang, 1 and G.P. Zeller 1

(ArgoNeuT Collaboration)

1 Fermi National Accelerator Laboratory, Batavia, IL 60510 USA 2 Yale University, New Haven, CT 06520 USA 3 Syracuse University, Syracuse, NY 13244 USA 4 Kansas State University, Manhattan, KS 66506 USA 5 Michigan State University, East Lansing, MI 48824 USA 6 Universit` a dell'Aquila e INFN, 67100 L'Aquila, Italy 7 University of Bern, 3012 Bern, Switzerland 8 The University of Texas at Austin, Austin, TX 78712 USA 9 INFN - Laboratori Nazionali del Gran Sasso, 67100 Assergi, Italy 10 Imperial College London, London, SW7 2AZ, UK (Dated: April 6, 2015)

We report on the first cross section measurements for charged current coherent pion production by neutrinos and antineutrinos on argon. These measurements are performed using the ArgoNeuT detector exposed to the NuMI beam at Fermilab. The cross sections are measured to be 2 6 . +1 2 . -1 0 . ( stat ) +0 3 . -0 4 . ( syst ) × 10 -38 cm 2 / Ar for neutrinos at a mean energy of 9 6 GeV . and 5 5 . +2 6 . -2 1 . ( stat ) +0 6 . -0 7 . ( syst ) × 10 -39 cm 2 / Ar for antineutrinos at a mean energy of 3 6 GeV. .

Neutrinos can produce single pion final states by coherently scattering from the entire nucleus. Both neutral current (NC) and charged current (CC) processes are possible. In these interactions, the squared fourmomentum transfer to the target nucleus, | t | , is small so the nucleus remains unchanged. In this Letter, we focus on the CC coherent pion production from muon neutrinos and antineutrinos on argon:

$$\nu _ { \mu } + \text{Ar} \to \mu ^ { - } + \pi ^ { + } + \text{Ar} ; \text{ \quad \ \ } ( 1 ) \text{ \ \ } \text{is}$$

$$\bar { \nu } _ { \mu } + \text{Ar} \rightarrow \mu ^ { + } + \pi ^ { - } + \text{Ar} ; \text{ \quad \ \ } ( 2 ) \text{ \ \ } \text{cur}$$

where the low | t | condition entails that the pions and muons are forward going with respect to the incoming neutrino direction.

There are several models from which one can extract cross sections and kinematical predictions for this interaction. The Rein-Seghal [1] model has been used to successfully describe high energy data within experimental uncertainties since the first observation of coherent pion production at the Aachen-Padova spark chamber [2] in 1983. This approach is based on Adler's Partially Conserved Axial Current (PCAC) theorem [3], which relates the pion production cross section to the cross section for the pion-nucleus scattering. This model is still the standard for neutrino generators today, such as genie [4], nuwro [5], and neut [6], with continued updates to the formalism and the pion-nucleus scattering data that is used. With recent interest in coherent pion production in the theoretical community, other PCAC models have been proposed [7, 8]. Microscopic models [9-11]

have also been suggested, which employ a full quantum mechanical treatment that explores the excitation and decay of the ∆ resonance. While the PCAC based models are a simple approach, tailored for the description of high energy data, their extension to the few GeV regime is not straightforward. Notably, the K2K [12] and SciBooNE [13] collaborations found cross section upper limits for the CC coherent pion production well below Rein-Seghal's estimation. The microscopic models are better motivated at lower neutrino energies but currently cannot be used to describe high energy data. Given the differences in these models, more experimental measurements are necessary to validate and tune the models and, in particular, better understand the transition region between microscopic and PCAC validity at E ν ∼ 3 -5 GeV.

In this Letter, a measurement of CC coherent pion production from the ArgoNeuT (Argon Neutrino Test) experiment is presented. ArgoNeuT [14] is a 170 L liquid argon time projection chamber (LArTPC), with dimensions 47 × 40 × 90 cm . The electric field inside the 3 TPC is 481V cm, and the drifted charge from parti-/ cle interactions is read out in two planes of 240 wires with 4 mm pitch (the induction and collection planes). The angle between the induction and collection plane wires is 60 degrees. ArgoNeuT is exposed to the NuMI beam [15] set in an antineutrino-enhanced mode, which provides a flux that is mostly muon antineutrino but still rich in muon neutrinos. The total number of protons on target (POT) accumulated during a 5-month

run is 1 2 . × 10 20 and the estimated integrated fluxes are 6 6 . × 10 11 muon neutrinos per cm 2 and 3 0 . × 10 12 muon antineutrinos per cm . 2 The differential flux can be found in reference [16]. Neutrino interactions comprise almost 60% of all the neutrino/antineutrino-induced events in the detector [16]. During this run, the MINOS near detector [17] placed downstream of ArgoNeuT is also operational. The muons that exit ArgoNeuT's TPC volume are matched to MINOS, in which the momentum and charge are reconstructed.

Using the LArSoft software [18], (anti)neutrino interactions are reconstructed, rendering a full characterization of the charged particles emerging in the ArgoNeuT detector. The software also provides the framework for a Monte Carlo (MC) simulation of the experiment. This is achieved by employing genie [4] as the neutrino event generator and geant4 [19] for the simulation of the propagation of products in the detector. The complete ArgoNeuT geometry is simulated along with the signal formation processes and taking into account electronic noise. The simulated events are fully reconstructed in the same way as for data. The propagation of particles in the MINOS near detector is simulated with geant3 [20]. A standalone version of MINOS simulation and reconstruction is used to characterize the matching of tracks passing from ArgoNeuT into MINOS.

The search for CC coherent pion production starts with an event selection which is used to find the two track topology of Eqs. (1) and (2). Each of the selection criteria described below is chosen in order to maximise the significance, defined as s/ √ s + , where b s and b are the numbers of signal and background events which pass the selection in the MC simulation. The MC used assumes the signal as modeled by Rein-Seghal. We start by requiring that two tracks are reconstructed in the event, originating from the same vertex. One track, identified as the muon, must be reconstructed in both ArgoNeuT and MINOS and matched between the two detectors. The unmatched track is the pion candidate. ArgoNeuT's precise calorimetry is used to discriminate pions from protons by defining an acceptance window for the mean dE/dx of the unmatched track. While the dE/dx of a pion will correspond to a Minimum Ionizing Particle (1 MIP), a proton track will leave an energy deposition several times higher ( &gt; 2 MIP). By applying a selection criteria on the dE/dx of the pion-candidate track, the CC quasi-elastic background is almost fully removed. The calorimetry capabilities of the detector are further exploited by investigating the Analog-toDigital (ADC) readout at the wires at the vertex. Low energy protons emerging at the vertex induce high ADC readouts at the first wire hits which are used to exclude the event. This selection reduces the background of interactions where multiple low-energy protons are produced, added either by nuclear effects or a result of deep inelastic scattering.

The lack of any particles other than the muon and the pion emerging from the vertex is further reinforced by another selection criteria. For each event, the charge readout inside a ∼ 20 cm × 20 cm box defined in the wire number versus drift time view of the collection plane is counted; the fraction of this charge that is associated with the two outgoing tracks must amount to at least 86%(84%) for antineutrino(neutrino) events. The collection plane is used because its response is better calibrated compared with the induction plane. This verification is crucial since it removes background events with activities around the interaction vertex that are not originated from the muon and the pion.

The event selection defined makes the most of the precise calorimetry and the high imaging resolution the ArgoNeuT detector is capable of, which are a characteristic of LArTPCs. We estimate the selection efficiencies to be (18 4 . ± 1 8)% for neutrino and (21 8 . . ± 0 8)% . for antineutrino events. The inefficiency is dominated by the track reconstruction inefficiency for overlapping tracks or complex topologies when the pion interacts with the argon nucleus. The systematic uncertainties associated to the assumptions on the kinematics of the signal events are accessed by estimating the efficiency using a different generator ( nuwro ). The difference between the efficiencies obtained with the two generators is kept as the systematic uncertainty.

A total of 167 antineutrino and 150 neutrino events have the two-track topology in the TPC with one track matched to a reconstructed track in MINOS. After applying the event selection described, 30 antineutrino and 24 neutrino candidate events remain. This event sample contains a background fraction, predominantly resonant and deep inelastic interactions, that ideally would be reduced by selecting events with low | t | = ∣ ∣ ( q -p π ) 2 ∣ ∣ , where q represents the momentum transfer from the neutrino and p π is the momentum carried by the pion. This approach is not feasible because most pions are not contained in the ArgoNeuT TPC so their momentum can't be estimated. Instead, we achieve signal from background separation by applying a multivariate method which exploits the topological and calorimetric information reconstructed in each event. The ROOT Toolkit for Multivariate Analysis [21] was used to create a Boosted Decision Tree (BDT) which is trained using genie signal and background samples. The classification is based on the angles of the pion and muon tracks, the visible energy loss of the pion from the TPC's calorimetry, the reconstructed muon momentum from MINOS and the mean stopping power of the first third of the muon track. The last of these parameters was added to help distinguish events where the start of the muon and pion tracks is overlapping. The angular parameters have the highest discrimination power. An example of a neutrino

FIG. 1: An example of CC Coherent pion production from a neutrino in ArgoNeuT. The neutrino's incoming direction is along the horizontal coordinate; the muon track corresponds to the most forward going one, making an angle of 1 2 . ◦ with the incoming neutrino direction. The opening angle between the muon and the pion track is 10 6 . ◦ . A kink in the pion trajectory is visible.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0598v3-figure-b5479ff3-6774-4352-aa0e-408781967d32.jpg)

interaction classified as signal by the BDT is show in Figure 1.

To estimate the rate of signal events, the BDT distribution in data is fitted to a linear combination of templates for signal and background obtained from simulation. The fit preserves the shape of the signal and background BDT distributions and finds the scale of these which best agrees with the data by minimising the effective χ 2 = -2 ln L , where L represents the likelihood assuming Poisson-distributed counts in each bin. The statistical error is found by evaluating the 1 σ interval, determined by ∆ χ 2 = χ 2 -χ 2 min = 1. Figure 2 shows the data and the best-fit signal and background distributions. The antineutrino signal is estimated to be 7 9 . +3 7 . -3 0 . events and the neutrino signal is 7 0 . +3 3 . -2 6 . events. The background contamination in the signal region (BDT Classification &gt; 0) is small.

The systematic uncertainties affecting the measurement are listed in Table I. These are dominated by the flux-scale uncertainty (10 -12%). Reconstruction effects have their impact estimated by adjusting the reconstructed values by ± 1 σ , where σ is the uncertainty on the reconstructed parameter. The absolute muon momentum estimated from the track curvature in the MINOS detector has a 4% systematic uncertainty [22] and the angular uncertainty assigned to tracks recon-

FIG. 2: Best-fit of the signal and background templates to the data. A BDT classification value of -1 means background-like and 1 is signal-like. The background and signal shapes are scaled to minimize an effective χ 2 function from which the statistical error is also extracted.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0598v3-figure-e3d86677-b29b-4254-ab73-2e476b99fe9c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0598v3-figure-ee6b4df5-9b5f-4284-bb50-53da4a326477.jpg)

structed in ArgoNeuT is 1 ◦ [23]. The contribution of background uncertainties is found by adjusting the contribution from each individual background process by ± 20% [24]. The rate at which the charge of the muon is mis-identified is also estimated and treated like the other backgrounds, though its contribution was found to be negligible. The effect of nuclear interactions affecting the production of background events is also considered. This is done by evaluating the fraction of background events added by final state interactions and reweighting this sample by a conservative factor ( ± 20%). Finally, the systematic error associated with the signal modeling is investigated by generating a signal template using nuwro . The difference in the number of signal events found after repeating the fit is our estimation of the systematic uncertainty.

The flux-averaged cross section is found by dividing the number of signal events by the efficiency of the selection, the number of target nuclei in the fiducial volume and the integrated (anti)neutrino flux. The measurements we report are

$$\xi _ { - } \quad \left \langle \sigma _ { \bar { \nu } _ { \mu } } \right \rangle = 5. 5 ^ { + 2. 6 } _ { - 2. 1 } ( s t a t ) ^ { + 0. 6 } _ { - 0. 7 } ( s y s t ) \times 1 0 ^ { - 3 9 } \text{cm} ^ { 2 } \quad ( 3 )$$

$$\frac { \psi _ { i = 1 } } { i = 1 } \Big < \sigma _ { \nu _ { \mu } } \Big > = 2. 6 ^ { + 1. 2 } _ { - 1. 0 } ( s t a t ) ^ { + 0. 3 } _ { - 0. 4 } ( s y s t ) ) \times 1 0 ^ { - 3 8 } \text{cm} ^ { 2 } \quad ( 4 )$$

per argon nuclei at 〈 E ¯ ν µ 〉 = 3 6 . ± 1 5 GeV and . 〈 E ν µ 〉 = 9 6 . ± 6 5 GeV, where the . ± 1 5(6 5) GeV represents the . . range that contains 68% of the flux. A comparison between these measurements, existing data, and the ReinSeghal model are shown in Figure 3. The antineutrino measurement agrees well with the Rein-Seghal model while the neutrino one deviates by ∼ 1 2 . σ .

In this Letter, we have presented the first cross section measurement of CC coherent pion production on argon. This is also the first time that machine learning techniques have been applied to LArTPC data analysis. The large uncertainties on the final cross section

TABLE I: Contributions to the total systematic uncertainty on the flux-averaged cross sections. The dominant backgrounds in this analysis are the CC quasi-elastic (QE), resonant (RES), and deep inelastic scattering (DIS) interactions.

|                         | Cross Section Uncertainty [%]   | Cross Section Uncertainty [%]   |
|-------------------------|---------------------------------|---------------------------------|
| Systematic Effect       | ¯ ν µ                           | ν µ                             |
| Flux normalization      | +10 . 0 - 12 . 0                | +10 . 0 - 12 . 0                |
| MINOS momentum res.     | ± 4 . 1                         | ± 4 . 3                         |
| ArgoNeuT angle res.     | ± 1 . 6                         | ± 2 . 7                         |
| CC QE background        | +0 . 3 - 0 . 4                  | +1 . 2 - 0 . 6                  |
| CC RES background       | +0 . 2 - 0 . 5                  | +0 . 4 - 0 . 3                  |
| CC DIS background       | ± 0 . 1                         | ± 0 . 3                         |
| Nuclear Effects         | ± 0 . 3                         | ± 0 . 7                         |
| POT                     | ± 0 . 1                         | ± 0 . 1                         |
| Number of Argon Targets | ± 2 . 2                         | ± 2 . 2                         |
| Efficiency              | ± 0 . 8                         | ± 1 . 8                         |
| Signal modeling         | ± 0 . 8                         | ± 5 . 7                         |
| Total systematics       | +11 . 3 - 13 . 1                | +12 . 9 - 14 . 5                |

values are dominated by the statistical errors. Using the precise calorimetry and the high resolution of the interaction vertex which are fundamental for this analysis, future LArTPC experiments will be able to provide decisive measurements for the understanding of neutrino induced coherent pion production.

We gratefully acknowledge the cooperation of the MINOS collaboration in providing their data for use in this analysis. We wish to acknowledge the support of Fermilab, the Department of Energy, and the National Science Foundation in ArgoNeuT's construction, operation, and data analysis. We also wish to acknowledge the support of STFC, FNAL and FCT's grant SFRH/BD/69814/2010.

- [1] D. Rein, and L. M. Seghal, Phys. Lett. B 657 , 207-209 (2007).
- [2] H. Faissner, et al. , Phys. Lett. B 125 , 230-236 (1983).
- [3] S. Adler, Phys. Rev. 135 , B963-B966 (1964).
- [4] C. Andreopoulos, et al. , Nucl. Instr. &amp; Meth. A 614 , 87 (2010).
- [5] T. Golan, C. Juszczak and J. T. Sobczyk, Phys. Rev. C 86 , 015505 (2012).
- [6] Y. Hayato, Acta Phys. Polon. B 40 , 2477-2489 (2009).
- [7] C. Berger, L. M. Seghal, Phys. Rev. D 79 , 053003 (2009).
- [8] E. A. Paschos and D. Schalla, Phys. Rev. D 80 , 033005, (2009).
- centeVacas, Phys. Rev. C 75 , 055501 (2007).
- [10] E. Hernandez, J. Nieves and M. Valverde, Phys. Rev. D 76 , 033005 (2007).
- [11] S. X. Nakamura, et al. , arXiv:0901.2366.
- [12] M. Hasegawa, et al. , Phys. Rev. Lett. 95 , 252301 (2005).
- [13] K. Hiraide, et al. , Phys. Rev. D 78 , 112004 (2008).
- [14] C. Anderson, et al. (ArgoNeuT Collaboration) JINST 7 , P10019 (2012).
- [15] K. Anderson, et al. , FERMILAB-DESIGN-1998-01 (1998).
- [16] R. Acciarri, et al. (ArgoNeuT Collaboration), Phys. Rev. D 89 , 112003 (2014).
- [17] D. G. Michael, et al. (MINOS Collaboration), Nucl. Instr. &amp; Meth. A 596 , 190 (2008).
- [18] E. Church, arXiv:1311.6774.
- [19] S. Agostinelli, et al. , Nucl. Instr. &amp; Meth. A 506 , 250 (2003).
- [20] Application Software Group, CERN Program Library Long Writeup W5013, CERN, 1994.
- [21] A. Hocker, et al. , PoS ACAT , 040 (2007).
- [22] P. Adamson, et al. , Phys. Rev. D 81 , 072002 (2010).
- [23] J. Spitz, Ph.D. thesis, Yale University, [FERMILABTHESIS-2011-36] (2011).
- [24] J. G. Morfin, and J. Nieves, and J. T. Sobczyk, Adv. High Energy Phys. 2012 , 934597 (2012).
- [25] H.-J. Grabosch, et al. , Zeitschrift fur Physik C Particles and Fields 31 , 203-211 (1986).
- [26] P. Vilain et al. , Phys. Lett. B 313 , 267-275 (1993).
- [27] A. Higuera, et al. , Phys. Rev. Lett. 113 , 261802 (2014).

FIG. 3: ArgoNeuT's CC coherent pion cross section measurements (open circle and filled circle) compared to ReinSeghal's model as implemented in genie and nuwro [5]. The statistical error is dominant (the systematic uncertainty is shown alone for comparison). Data from other experiments in the same energy range is also shown. These consist in measurements made by SKAT (filled square, open square), CHARM II (diamond) and MINERvA (triangle, inverted triangle) [25-27]. These measurements are scaled to Argon assuming the A 1 / 3 dependance from the Rein-Seghal model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0598v3-figure-8313f69d-4c4e-47cb-be8f-9af9d415fe7b.jpg)
R. Acciarri C. Adams J. Asaadi B. Baller T. Bolton C. Bromberg F. Cavanna E. Church D. Edmunds A. Ereditato S. Farooq B. Fleming H. Greenlee R. Hatcher G. Horton-Smith C. James E. Klein K. Lang P. Laurens R. Mehdiyev B. Page O. Palamara K. Partyka G. Rameika B. Rebel E. Santos M. Soderberg J. Spitz A. M. Szelc M. Weber T. Yang G. P. Zeller
First Measurement of Neutrino and Antineutrino Coherent Charged Pion Production on Argon
We report on the first cross section measurements for charged current
coherent pion production by neutrinos and antineutrinos on argon. These
measurements are performed using the ArgoNeuT detector exposed to the NuMI beam
at Fermilab. The cross sections are measured to be
$2.6^{+1.2}_{-1.0}(stat)^{+0.3}_{-0.4}(syst) \times 10^{-38}
\textrm{cm}^{2}/\textrm{Ar}$ for neutrinos at a mean energy of $9.6$ GeV and
$5.5^{+2.6}_{-2.1}(stat)^{+0.6}_{-0.7}(syst) \times 10^{-39}
\textrm{cm}^{2}/\textrm{Ar}$ for antineutrinos at a mean energy of $3.6$ GeV.