Citation for this paper:
Abe, K., Amey, J., Andreopoulos, C., Anthony, L., Antonova, M., Karlen, D., … Zito, M. (2018). First measurement of the νμ charged-current cross section on a water target without pions in the
UVicSPACE: Research & Learning Repository
_____________________________________________________________
Faculty of Science
Faculty Publications
_____________________________________________________________
First measurement of the νμ charged-current cross section on a water targetwithout pions in the final state
K. Abe, J. Amey, C. Andreopoulos, L. Anthony, M. Antonova, D. Karlen, … & M. Zito January 2018
© 2018 K. Abe et al. This is an open access article distributed under the terms of the Creative Commons Attribution License. https://creativecommons.org/licenses/by/4.0/
This article was originally published at:
First measurement of the
ν
μcharged-current cross section on a water target
without pions in the final state
K. Abe,45J. Amey,15C. Andreopoulos,43,25L. Anthony,25M. Antonova,20S. Aoki,22A. Ariga,1 Y. Ashida,23S. Ban,23 M. Barbi,37G. J. Barker,53G. Barr,33C. Barry,25M. Batkiewicz,11V. Berardi,16S. Berkman,3,49S. Bhadra,57S. Bienstock,34 A. Blondel,10S. Bolognesi,5S. Bordoni,13,*B. Bourguille,13,* S. B. Boyd,53D. Brailsford,24A. Bravar,10C. Bronner,45
M. Buizza Avanzini,9 J. Calcutt,27R. G. Calland,21 T. Campbell,7 S. Cao,12S. L. Cartwright,41M. G. Catanesi,16 A. Cervera,14A. Chappell,53C. Checchia,18D. Cherdack,7 N. Chikuma,44 G. Christodoulou,25J. Coleman,25 G. Collazuol,18D. Coplowe,33A. Cudd,27 A. Dabrowska,11G. De Rosa,17T. Dealtry,24P. F. Denner,53S. R. Dennis,25 C. Densham,43F. Di Lodovico,36S. Dolan,33O. Drapier,9K. E. Duffy,33J. Dumarchez,34P. Dunne,15S. Emery-Schrenk,5 A. Ereditato,1 T. Feusels,3,49A. J. Finch,24G. A. Fiorentini,57G. Fiorillo,17C. Francois,1,† M. Friend,12,† Y. Fujii,12,† D. Fukuda,31Y. Fukuda,28A. Garcia,13C. Giganti,34F. Gizzarelli,5 T. Golan,55M. Gonin,9D. R. Hadley,53L. Haegel,10 J. T. Haigh,53D. Hansen,35J. Harada,32M. Hartz,21,49T. Hasegawa,12,†N. C. Hastings,37T. Hayashino,23Y. Hayato,45,21
A. Hillairet,50T. Hiraki,23 A. Hiramoto,23S. Hirota,23M. Hogan,7 J. Holeczek,42F. Hosomi,44 K. Huang,23 A. K. Ichikawa,23M. Ikeda,45J. Imber,9 J. Insler,26R. A. Intonti,16T. Ishida,12,† T. Ishii,12,† E. Iwai,12K. Iwamoto,44 A. Izmaylov,14,20B. Jamieson,54M. Jiang,23S. Johnson,6P. Jonsson,15C. K. Jung,30,‡M. Kabirnezhad,29A. C. Kaboth,39,43 T. Kajita,46,‡H. Kakuno,47J. Kameda,45D. Karlen,50,49T. Katori,36E. Kearns,2,21,‡ M. Khabibullin,20A. Khotjantsev,20 H. Kim,32J. Kim,3,49S. King,36J. Kisiel,42A. Knight,53A. Knox,24T. Kobayashi,12,†L. Koch,40T. Koga,44P. P. Koller,1
A. Konaka,49L. L. Kormos,24Y. Koshio,31,‡ K. Kowalik,29Y. Kudenko,20,§ R. Kurjata,52T. Kutter,26J. Lagoda,29 I. Lamont,24M. Lamoureux,5 P. Lasorak,36M. Laveder,18M. Lawe,24M. Licciardi,9 T. Lindner,49Z. J. Liptak,6 R. P. Litchfield,15X. Li,30 A. Longhin,18 J. P. Lopez,6 T. Lou,44L. Ludovici,19X. Lu,33L. Magaletti,16K. Mahn,27 M. Malek,41S. Manly,38L. Maret,10A. D. Marino,6J. F. Martin,48P. Martins,36S. Martynenko,30T. Maruyama,12,† V. Matveev,20K. Mavrokoridis,25W. Y. Ma,15E. Mazzucato,5 M. McCarthy,57N. McCauley,25K. S. McFarland,38
C. McGrew,30 A. Mefodiev,20C. Metelko,25 M. Mezzetto,18A. Minamino,56 O. Mineev,20S. Mine,4 A. Missert,6 M. Miura,45,‡ S. Moriyama,45,‡ J. Morrison,27Th. A. Mueller,9 T. Nakadaira,12,† M. Nakahata,45,21K. G. Nakamura,23 K. Nakamura,21,12,†K. D. Nakamura,23Y. Nakanishi,23S. Nakayama,45,‡ T. Nakaya,23,21K. Nakayoshi,12,†C. Nantais,48
C. Nielsen,3,49K. Nishikawa,12,†Y. Nishimura,46P. Novella,14J. Nowak,24H. M. O’Keeffe,24K. Okumura,46,21 T. Okusawa,32W. Oryszczak,51S. M. Oser,3,49T. Ovsyannikova,20R. A. Owen,36Y. Oyama,12,† V. Palladino,17 J. L. Palomino,30 V. Paolone,35N. D. Patel,23P. Paudyal,25M. Pavin,34D. Payne,25Y. Petrov,3,49L. Pickering,15 E. S. Pinzon Guerra,57C. Pistillo,1B. Popov,34,∥M. Posiadala-Zezula,51J.-M. Poutissou,49A. Pritchard,25P. Przewlocki,29
B. Quilain,23T. Radermacher,40E. Radicioni,16 P. N. Ratoff,24M. A. Rayner,10E. Reinherz-Aronis,7C. Riccio,17 E. Rondio,29B. Rossi,17S. Roth,40A. C. Ruggeri,17A. Rychter,52K. Sakashita,12,† F. Sánchez,13S. Sasaki,47,† E. Scantamburlo,10K. Scholberg,8,‡J. Schwehr,7M. Scott,49Y. Seiya,32T. Sekiguchi,12,†H. Sekiya,45,21,‡D. Sgalaberna,10
R. Shah,43,33A. Shaikhiev,20F. Shaker,54D. Shaw,24M. Shiozawa,45,21T. Shirahige,31M. Smy,4 J. T. Sobczyk,55 H. Sobel,4,21J. Steinmann,40T. Stewart,43P. Stowell,41Y. Suda,44S. Suvorov,20A. Suzuki,22S. Y. Suzuki,12,†Y. Suzuki,21 R. Tacik,37,49M. Tada,12,†A. Takeda,45Y. Takeuchi,22,21R. Tamura,44H. K. Tanaka,45,‡H. A. Tanaka,48,49,¶T. Thakore,26 L. F. Thompson,41S. Tobayama,3,49W. Toki,7 T. Tomura,45 T. Tsukamoto,12,† M. Tzanov,26W. Uno,23M. Vagins,21,4
Z. Vallari,30 G. Vasseur,5 C. Vilela,30T. Vladisavljevic,33,21T. Wachala,11 J. Walker,54,‡ C. W. Walter,8,‡ Y. Wang,30 D. Wark,43,33M. O. Wascko,15A. Weber,43,33R. Wendell,23,‡ M. J. Wilking,30C. Wilkinson,1 J. R. Wilson,36 R. J. Wilson,7 C. Wret,15Y. Yamada,12,† K. Yamamoto,32S. Yamasu,31 C. Yanagisawa,30,**T. Yano,22S. Yen,49 N. Yershov,20M. Yokoyama,44,‡ T. Yuan,6 M. Yu,57A. Zalewska,11J. Zalipska,29L. Zambelli,12,† K. Zaremba,52
M. Ziembicki,52E. D. Zimmerman,6 and M. Zito5 (T2K Collaboration)
1
University of Bern, Albert Einstein Center for Fundamental Physics, Laboratory for High Energy Physics (LHEP), Bern, Switzerland
2
Boston University, Department of Physics, Boston, Massachusetts, USA
3University of British Columbia, Department of Physics and Astronomy,
Vancouver, British Columbia, Canada
4University of California, Irvine, Department of Physics and Astronomy, Irvine, California, USA 5
IRFU, CEA Saclay, Gif-sur-Yvette, France
6University of Colorado at Boulder, Department of Physics, Boulder, Colorado, USA 7
8Duke University, Department of Physics, Durham, North Carolina, USA 9
Ecole Polytechnique, IN2P3-CNRS, Laboratoire Leprince-Ringuet, Palaiseau, France
10University of Geneva, Section de Physique, DPNC, Geneva, Switzerland 11
H. Niewodniczanski Institute of Nuclear Physics PAN, Cracow, Poland
12High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki, Japan 13
Institut de Fisica d’Altes Energies (IFAE), The Barcelona Institute of Science and Technology, Campus UAB, Bellaterra (Barcelona) Spain
14
IFIC (CSIC & University of Valencia), Valencia, Spain
15Imperial College London, Department of Physics, London, United Kingdom 16
INFN Sezione di Bari and Universit`a e Politecnico di Bari, Dipartimento Interuniversitario di Fisica, Bari, Italy
17
INFN Sezione di Napoli and Universit `a di Napoli, Dipartimento di Fisica, Napoli, Italy
18INFN Sezione di Padova and Universit `a di Padova, Dipartimento di Fisica, Padova, Italy 19
INFN Sezione di Roma and Universit `a di Roma“La Sapienza”, Roma, Italy
20Institute for Nuclear Research of the Russian Academy of Sciences, Moscow, Russia 21
Kavli Institute for the Physics and Mathematics of the Universe (WPI), The University of Tokyo Institutes for Advanced Study, University of Tokyo,
Kashiwa, Chiba, Japan
22Kobe University, Kobe, Japan 23
Kyoto University, Department of Physics, Kyoto, Japan
24Lancaster University, Physics Department, Lancaster, United Kingdom 25
University of Liverpool, Department of Physics, Liverpool, United Kingdom
26Louisiana State University, Department of Physics and Astronomy,
Baton Rouge, Louisiana, USA
27Michigan State University, Department of Physics and Astronomy, East Lansing, Michigan, USA 28
Miyagi University of Education, Department of Physics, Sendai, Japan
29National Centre for Nuclear Research, Warsaw, Poland 30
State University of New York at Stony Brook, Department of Physics and Astronomy, Stony Brook, New York, USA
31
Okayama University, Department of Physics, Okayama, Japan
32Osaka City University, Department of Physics, Osaka, Japan 33
Oxford University, Department of Physics, Oxford, United Kingdom
34UPMC, Universit´e Paris Diderot, CNRS/IN2P3, Laboratoire de Physique Nucl´eaire et de Hautes
Energies (LPNHE), Paris, France
35University of Pittsburgh, Department of Physics and Astronomy,
Pittsburgh, Pennsylvania, USA
36Queen Mary University of London, School of Physics and Astronomy,
London, United Kingdom
37University of Regina, Department of Physics, Regina, Saskatchewan, Canada 38
University of Rochester, Department of Physics and Astronomy, Rochester, New York, USA
39
Royal Holloway University of London, Department of Physics, Egham, Surrey, United Kingdom
40
RWTH Aachen University, III. Physikalisches Institut, Aachen, Germany
41University of Sheffield, Department of Physics and Astronomy, Sheffield, United Kingdom 42
University of Silesia, Institute of Physics, Katowice, Poland
43STFC, Rutherford Appleton Laboratory, Harwell Oxford, and Daresbury Laboratory,
Warrington, United Kingdom
44University of Tokyo, Department of Physics, Tokyo, Japan 45
University of Tokyo, Institute for Cosmic Ray Research, Kamioka Observatory, Kamioka, Japan
46
University of Tokyo, Institute for Cosmic Ray Research, Research Center for Cosmic Neutrinos, Kashiwa, Japan
47
Tokyo Metropolitan University, Department of Physics, Tokyo, Japan
48University of Toronto, Department of Physics, Toronto, Ontario, Canada 49
TRIUMF, Vancouver, British Columbia, Canada
50University of Victoria, Department of Physics and Astronomy,
Victoria, British Columbia, Canada
51University of Warsaw, Faculty of Physics, Warsaw, Poland 52
53University of Warwick, Department of Physics, Coventry, United Kingdom 54
University of Winnipeg, Department of Physics, Winnipeg, Manitoba, Canada
55Wroclaw University, Faculty of Physics and Astronomy, Wroclaw, Poland 56
Yokohama National University, Faculty of Engineering, Yokohama, Japan
57York University, Department of Physics and Astronomy, Toronto, Ontario, Canada
(Received 24 August 2017; published 8 January 2018; corrected 22 March 2018) This paper reports the first differential measurement of the charged-current interaction cross section ofνμ on water with no pions in the final state. This flux-averaged measurement has been made using the T2K experiment’s off-axis near detector, and is reported in doubly differential bins of muon momentum and angle. The flux-averaged total cross section in a restricted region of phase space was found to be σ ¼ ð0.95 0.08ðstatÞ 0.06ðdet systÞ 0.04ðmodel systÞ 0.08ðfluxÞÞ × 10−38cm2=n.
DOI:10.1103/PhysRevD.97.012001
I. INTRODUCTION
The Tokai-to-Kamioka (T2K) experiment[1]is a long-baseline neutrino oscillation experiment, with a beam originating at Japan Proton Accelerator Complex (J-PARC) which consists primarily of muon neutrinos. T2K has measured the disappearance of muon neutrinos[2]
and the appearance of electron neutrinos[3], using the off-axis ND280 near detector on the J-PARC site, and the Super-Kamiokande detector [4], located 295 km away.
At the energies of the T2K beam line, the main neutrino interaction process is charged-current quasielastic (CCQE) interactions (νμþ n → μ−þ p). Because these neutrino interactions occur within nuclear targets and not on free nuclei, additional nuclear effects and final state interactions can modify the composition and kinematics of the particles that are observed to be exiting the interaction. This paper focuses on a measurement of CCQE-like events (in bins of the muon angle and momentum) in which no pions (charged or neutral) are observed in the final state (CC0π). The active target regions of the ND280 near detector, as will be discussed in Sec. II, are primarily composed of plastic scintillator, but the far detector is water based. While the near detector measurements on hydrocarbon targets can help to greatly constrain the flux and cross section uncertainties for the oscillation analyses, one of the
dominant remaining uncertainties is due to potential differences between the oxygen and carbon cross sections that are not currently well constrained by the ND280 detector [5]. The ND280 detector also contains water targets, and this paper presents a measurement of the νμ CC0π interaction cross section on water. This process is very important for T2K’s neutrino oscillation measure-ments since this is the dominant reaction in the far detector. While there are differential measurements of the CCQE cross sections on carbon or hydrocarbon[6–10], there are none on water. The K2K experiment has published a measurement of the axial vector mass in neutrino-oxygen CCQE interactions[11]using the SciFi detector, but not a differential cross section measurement.
II. THE T2K EXPERIMENT
A high-intensity beam is produced at J-PARC at Tokai village, Ibaraki, Japan and directed to Super-Kamiokande. In order to enhance the sensitivity to T2K’s primary physics goals, the appearance of electron neutrinos and the dis-appearance of muon neutrinos, the beam energy peaks at the oscillation maximum of 0.6 GeV, and this is achieved by adopting the off-axis method[1].
A. The T2K beam line
High-intensity 30 GeV protons from the J-PARC accel-erator strike a graphite target every 2.48 seconds and produce charged pions and kaons. These pions and kaons are focused to the forward direction by three horn magnets
[12]and decay in the 96-m long decay volume. Neutrinos (antineutrinos) are produced from decays of positively (negatively) charged pions. Horn magnets select pions of either sign by flipping the current direction.
T2K started data taking in 2010 with a beam of primarily muon neutrinos until May 2014 and then started data taking with a primarily muon antineutrino beam until May 2016, and is continuing to take data. The analysis presented here uses data taken in the neutrino mode.
*Now at CERN.
†Also at J-PARC, Tokai, Japan.
‡Also at Kavli IPMU (WPI), the University of Tokyo, Japan. §Also at National Research Nuclear University“MEPhI” and
Moscow Institute of Physics and Technology, Moscow, Russia.
∥Also at JINR, Dubna, Russia.
¶Also at Institute of Particle Physics, Canada.
**Also at BMCC/CUNY, Science Department, New York,
New York, USA.
Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI. Funded by SCOAP3.
The neutrino flux is calculated based on the measurement of primary proton beam profiles and hadron production data, including the measurement of the pion and kaon yields[13]
by the NA61/SHINE experiment[14]at CERN. The total absolute flux uncertainty is about 10% at the peak energy, and details of the flux calculation are described in Ref.[15].
B. The T2K ND280 near detector
This analysis looks for neutrino interactions in the T2K off-axis near detector, shown in Fig. 1, specifically for events in the pi-zero (PØD) subdetector. This subdetector, described in more detail in [16], consists of alternating planes of scintillator bars, sandwiched between lead or brass radiator layers. A more detailed schematic of the PØD is shown in Fig. 2. There are a total of 40 scintillator modules in the PØD, each one composed of a vertical layer of triangular scintillator bars and a layer of horizontal bars. In the upstream and downstream portions of the PØD (the “upstream ECal” and the “central ECal”), the scintillator layers alternate with lead sheets. The central 25 layers alternate the scintillator layers with brass sheets and bladders which can be filled with water or air. This analysis will use data taken in both configurations. The total mass of water in the PØD in the fiducial volume is approximately 1902 kg. Downstream of the PØD is a tracker with three TPCs[17]and two fine-grained scintillator detectors[18]. The PØD and tracker are surrounded by electromagnetic calorimeters [19]inside a magnet with a 0.2 T field that is instrumented with muon detectors[20].
III. ANALYSIS A. Event selection
Hadrons produced within the nucleus are subject to final state interactions (FSIs) that can reabsorb the hadron or
alter their kinematics as they emerge from the nucleus. Pion absorption for example can make a charged-current reso-nant (CCRES) interaction appear as a CCQE final state. This forces the selection to be based on a topology classified by the number and type of outgoing particles from the nucleus in order to be less model dependent. A single muon and zero pion topology is called CCQE-like or CC0π. The event selection identifies CC0π-enhanced samples by selecting events where a single track was reconstructed in the PØD. Identical selections are applied to both the water-in and water-out samples.
The result in this paper relies on a subset of the neutrino mode data as broken down in Table I and totaling 5.52 × 1020 protons on the target (POT). The selection starts by identifying beam spills where the spill information and all ND280 subdetectors are known to have high-quality data. Each spill of the proton beam consists of several FIG. 1. An exploded view of T2K’s off-axis near detector.
FIG. 2. A schematic of the PØD, showing the structure of the layers. The beam travels from the left to the right. The red box outlines the approximate fiducial volume used for this analysis.
TABLE I. T2K runs and their associated POT, filtered for spills where all ND280 detectors were flagged with good data quality.
T2K Run POT 2 Water 4.29 × 1019 2 Air 3.55 × 1019 3c Air 1.35 × 1020 4 Water 1.63 × 1020 4 Air 1.76 × 1020 Total 5.52 × 1020
clearly separated bunches of protons. The selection then selects bunches containing tracks reconstructed in the TPCs associated with vertices reconstructed within the PØD fiducial volume, which extends from the middle of the scintillator of the first water layer to the middle of the scintillator in the last water layer and 25 cm in from the edges in the xy-plane (Fig. 2). In practice, these tracks are reconstructed from segments found in the PØD and TPC subdetector reconstruction processes. After identifying these bunches, the analysis identifies the highest-momentum negatively charged track as the muon candidate. If no negatively charged track is found, the bunch is removed from the selection. The selection then applies a TPC quality cut, removing tracks with no more than 18 nodes in the TPC reconstruction. The final cut for the selected CC0π sample requires that only a single track was reconstructed anywhere in the PØD in that bunch. The cross section reported here is restricted to bins in the region of muon kinematics where cosθμ≥ 0 and pμ≤ 5 GeV. The total number of selected events in the water-in data sample is 12 777 (and based on simulations it is expected that approximately 3860 of these events are true inter-actions on water) and the total number in the water-out data (which has a larger POT exposure) sample is 13 370. Due to the requirement of a negatively charged track, the ¯νμ contamination is this sample is very small and is estimated to be approximately 0.1% for both samples.
Figure 3 shows distributions of events that pass the selection for the water-in sample for data and a Monte Carlo (MC) detector simulation. The MC simulation chain primarily uses the NEUT neutrino event generator [21]
to provide the kinematics for particles emerging from neutrino interactions and a GEANT4-based[22,23]package to simulate these particles moving through the detector
(usingGEANT version 4.9.4). The QGSP_BERT model is used for hadronic interactions. The MC distribution is separated by interaction channel. A detailed review of neutrino interactions can be found in[24]for example. The two largest sources of background are charged-current interactions with a single outgoingπþ (CC1πþ) and any other charged-current interaction (CCOther) not catego-rized as CC0π or CC1πþ. The CC1πþ topology is due primarily to pion resonance production and CCOther to deep inelastic scattering[24]. Additionally, neutral current backgrounds are classified as “BKG” and interactions occurring outside the PØD fiducial volume are classified as“out of FV” or “OOFV.” Distributions for the water-out sample look similar.
To provide a data-driven constraint on the background, sideband selections are used for the CC1πþ and CCOther topologies. Their selection is identical to the signal selec-tion with the excepselec-tion of the final cut. For the CC1πþ sideband, events with exactly two PØD reconstructed tracks in a bunch along with at least one PØD reconstructed Michel electron are selected. For the CCOther sideband, events with greater than two PØD reconstructed tracks are selected. These cuts reduce overlap between the signal and sideband selections as shown in Figs.4and5for the water-in samples. The water-out sideband samples look very similar. The numbers of selected events in the signal and sidebands regions are summarized in TablesII andIII.
B. Cross section extraction
The event selections described above are binned in the reconstructed double differentialðpμ; cos θμÞ phase space. Reconstructed kinematics are approximations to the true initial state of the muon. To extract the true kinematics from the reconstructed, an unfolding technique is used based
0 0.5 1 1.5 2 2.5 3 3.5 4 5 ]] 0 0 0 0.2 0.4 0.6 0.8 1 1.2 3 10 × Data π CC-0 + π CC-1 CC-other BKG out FV Events [GeV] μ p 4.5
FIG. 3. Data and MC distributions of the water-in CC0π signal selection.NEUT5.3.2MC has been normalized to data POT and is sorted into various truth topologies.
0 0.5 1 1.5 2 2.5 3 3.5 4 [G [G 5]] 0 0 0 10 20 30 40 50 60 Data π CC-0 + π CC-1 CC-other BKG out FV Events [GeV] μ p 4.5
FIG. 4. Data and MC distributions of the water-in CC1πþ sideband selection.NEUT5.3.2MC has been normalized to data POT and is sorted into various truth topologies.
on D’Agostini’s method with the MC truth as the prior[25]. The purpose of unfolding is to remove detector reconstruction related imperfections to achieve a more accurate representation of how the muon emerged from the interaction. By correcting the data samples from PØD water-in and water-out configurations separately, a subtraction procedure ultimately gives a cross section on
water. Additionally, the MC is tuned to account for flux, interaction modeling, and detector corrections. The final tuned MC is then used to calculate efficiency corrections, purity corrections, and the unfolding matrix.
The exclusive CC0π signal is evaluated by correcting the data selection based on the CC0π signal purity in the MC. The signal purity is calculated using the MC truth information and the data from the sideband samples. As described above, CC1πþinteractions are the largest sources of background, followed by CCOther interactions. Therefore, sideband samples for these two backgrounds were selected as a data-driven background constraint. The ratio of the overall data sideband normalization to the overall MC sideband normalization is calculated and used to constrain the corresponding background in the CC0π MC selection as
B0j¼ Bj Sd
Sm; ð1Þ
where j denotes a bin index, B0 the sideband constrained background, B the original background in the MC signal selection, and Sdand Sm the sideband normalizations from data and MC respectively. This affects the signal purity used in the background correction,
p0j¼ 1 − B0j N0m;j
; ð2Þ
where N0m;j¼ Nm;j− Bjþ B0j and Nm is the original number of events in the CC0π MC selection.
For this analysis, the initial prior used in the D’Agostini unfolding technique was taken from the MC truth and a single iteration is used. Fake-data studies showed that multiple iterations did not improve results but increased uncertainties. This increase is due to the correlation introduced between the data and successive priors [26]
and can produce large fluctuations between neighboring bins. Additionally while in fake-data tests the totalχ2 for the unfolded result was slightly larger than expected (∼1.7 not including any systematics), it did not significantly decrease with additional iterations. Thus a single iteration was chosen. This must be regarded as a regularized result. With the MC truth as the prior, a single iteration Bayesian unfolding matrix is equivalent to directly constructing an unfolding matrix based on the MC. This means that the unfolding matrix is calculated based on the mapping between truth and reconstructed kinematics in the MC.
To extract the neutrino cross section on water an addi-tional subtraction step after the unfolding is required. The PØD fiducial volume contains plastic scintillator layers and thin brass sheets sandwiched between layers of water[16]. The water layers act as a passive target making it difficult to know whether an interaction occurred on water or on some other target nucleus. To work around this, the PØD was designed to be drained and filled during different run
0 0.5 1 1.5 2 2.5 3 3.5 4 [G [G 5]] 0 0 0 50 100 150 200 250 Data π CC-0 + π CC-1 CC-other BKG out FV [GeV] μ p Events 4.5
FIG. 5. Data and MC distributions of the water-in CCOther sideband selection.NEUT5.3.2MC has been normalized to data POT and is sorted into various truth topologies.
TABLE II. Number of water-out data events contained in the selected signal and sideband samples after each selection cut is applied. The cut on the number of Michel electrons is only applied to the CC1πþsample.
Water-out data Number of selected events
Cut CC0π CC1πþ CC Other
Data quality and PØDþ TPC
μ−candidate in fiducial volume 33 083 33 083 33 083
Number of PØD tracks 13 370 8107 11 606
Number of Michel e− NA 1710 NA
TABLE III. Number of water-in data events contained in the selected signal and sideband samples after each selection cut is applied. The cut on the number of Michel electrons is only applied to the CC1πþsample.
Water-in data Number of selected events
Cut CC0π CC1πþ CC Other
Data quality and PØDþ TPC
μ−candidate in fiducial volume 27 713 27 713 27 713
Number of PØD tracks 12 777 6297 8639 Number of Michel e− NA 1434 NA
periods. Everything being equal except for the inclusion or exclusion of water, a subtraction of the true water-in and water-out distributions should give the number of inter-actions on water. Since the PØD is in two different detector configurations with and without water, it is necessary to correct for detector smearing and inefficiencies before the subtraction. A direct subtraction of the reconstructed event rates would give an incorrect estimate for the actual event rate on water. Therefore, we first unfold the reconstructed distribution for water-in and water-out separately to get an approximation of their true distributions, then subtract the unfolded distributions to get the distribution of interactions on water. Specifically, this is given by
NH2O i ¼ Uw ijNwj ϵw i − RUaijNaj ϵa i ; ð3Þ
where the indexes i and j indicate true and reconstructed bins respectively, w and a indicate water-in and water-out periods respectively, N is the number of purity-corrected, signal events measured in the data signal selection, ϵ the selection efficiency, and R the flux normalization factor between water-in and water-out periods. Uijrepresents the unfolding matrix. From this, the differential cross section on water can be expressed as
dσ dxi¼ NH2O i FwNnΔxi ; ð4Þ
where Fwis the integrated flux over the water-in period, Nn the number of neutrons, andΔxithe area of bin i across the phase space x. Fw is the POT-weighted flux of all water-in periods and calculated using flux simulations and from the constraints described above. As the water-in and water-out periods have different beam exposures, total flux for the water-out periods, Fa, is used to scale the flux normali-zation ratio, R ¼ Fw=Fa. Additionally, correlated errors between the water-in and water-out periods are taken into account in the error propagation and discussed in Sec.III C.
C. Systematic uncertainties
As in previous T2K analyses, the systematic uncertain-ties can be separated into three principal categories: flux uncertainties, cross section model uncertainties, and detec-tor uncertainties[9,27].
The flux and cross section errors are calculated by reweighting the individual simulated events. Once the events are reweighted, the unfolding matrix is regenerated and the cross section is recalculated using the new weights. The uncertainties for flux and cross section parameters are taken from the same covariance matrices used as inputs to the near detector fit in the T2K oscillation analysis [27]. There are 46 total parameters in the covariance matrix: 25 for flux, 6 for final-state interactions, and 15 other cross section parameters. By generating many throws for the
model parameters and calculating the resulting event weights, we can calculate a covariance matrix for the final result. If dσO
dxi;k gives the cross section result for bin i and
parameter throw k which has n throws, then the covariance between bins i and j is
covði; jÞ ¼1 n Xn−1 k¼0 dσO dxi;k −dσO¯ dxi dσO dxj;k −dσO¯ dxj ð5Þ where ¯ dσO dxi ¼ 1 n Xn−1 k¼0 dσO dxi;k: ð6Þ
The flux uncertainty is dominated by the uncertainty in the total flux normalization and is found to be 8.76%. Flux uncertainties are due in large part to uncertainties in the hadron production model but are affected by beam line uncertainties as well [15]. A parametrization in neutrino energy and flavor is used to propagate flux uncertainties. Cross section model uncertainties include both uncertain-ties in basic neutrino-nucleus scattering and uncertainuncertain-ties from FSI [27]. Parameters that govern the neutrino inter-action and nuclear description models used in the NEUT generator[21]are used to propagate cross section uncertainties.
Neutrino-nucleus scattering depends on nucleon form factors and the nuclear medium model. Neutrino generators typically assume a dipole form factor governed by an axial mass parameter for QE and resonant charged-current interactions. Additionally, many of the interaction param-eters tune the normalization on certain interaction channels that are poorly understood. These include 2p2h contribu-tions that arise from multinucleon correlacontribu-tions inside the nucleus. The current understanding of these multinucleon interactions is based on the meson exchange model where two nucleons exchange a meson current such that the charged-current interaction involves both nucleons. This is expected to enhance the CCQE-like cross section. Thus, instead of a single final state nucleon, two or more nucleons may be ejected out of the nucleus. Other parameters, such as the nuclear binding energy and the Fermi momentum describe the nuclear medium.
An intranuclear cascade model describes secondary pion propagation within the nucleus[28]. Uncertainties on FSI are calculated by tuning the pion production, absorption, rescattering and charge-exchange probabilities.
The primary neutrino generator used for this analysis is a tuned version of NEUT 5.3.2 [21,29]. The tuning is applied based on best-fit values from a fit of three models implemented in NEUT to external MINERνA and MiniBooNE CCQE measurements[30]. The NEUT model uses the Smith-Moniz relativistic Fermi gas model [31]
and a multinucleon exchange (2p2h) model [32,33]. An axial mass of MQEA ¼ ð1.15 0.03Þ GeV=c2, the amplitude of the multinucleon effects 2p2h norm¼ ð27 12Þ%, and Fermi momentum pF ¼ ð223 5Þ MeV=c were used. The uncertainties in these parameters are propagated to the final cross section measurement. A secondary neutrino generator used for testing isGENIE 2.8.0 [34].
Table IV shows the individual contributions of several interaction model parameters to the total cross section. Note that the 2p2h normalization, binding energy Eb, and Fermi momentum pF are nuclei dependent. These are calculated with 100 throws without taking correlations into account. This illustrates and identifies the effects of a single parameter on the cross section result. For the final result, 800 throws were used along with covariances across the full parameter-space to generate the interaction and other systematic uncertainties, and 2000 throws were used for statistical errors. The lower number of throws used in Table IV means there is a larger statistical error on the uncertainty itself.
Detector systematics describe uncertainties in recon-structed properties of events rather than the underlying physics. These are separated into weight systematics and variation systematics. Weight systematics modify the weight to be applied to simulated events. These include effects such as differences in tracking efficiencies between data and simulation. Variation systematics modify observ-ables such as the reconstructed particle momentum of a simulated event. These allow for migration of events between different kinematic bins. Systematics related to the uncertainty on the fiducial mass must be treated in a different manner, as these affect both the expected number of events and the total mass used to normalize the cross section measurement.
The event reconstruction in the PØD uses a Kalman filter algorithm to calculate the properties of tracks found in the detector. The reconstructed momentum will primarily be influenced by two main effects: the curvature of the track in the tracker and the energy loss within the PØD. Because relatively little information about the track curvature can be obtained by the PØD and because the small amount of energy loss in the tracker will have little effect on the
reconstructed initial momentum, the systematics for these two effects are treated separately. Uncertainties in the energy loss in the PØD are estimated from a sample of cosmic ray muons that pass through the tracker and stop in the PØD. To propagate these uncertainties, the recon-structed momentum in the simulation is altered by
ðΔpÞPØD
scale¼ ðxsσsÞðpreco− precoTPC1Þ ð7Þ for the scale uncertainty and
ðΔpÞPØD
res ¼ ðxrσrÞðpreco− precoTPC1− ptruthþ ptruthTPC1Þ ð8Þ for the resolution uncertainty. In these equationsσ repre-sents the uncertainty and x is a normally distributed random number that changes for each systematics throw. Similar numbers are obtained for both the water-out and water-in geometries, and a 1.4% uncertainty is applied to the energy loss scale and a 7% uncertainty to the resolution. Tracker curvature systematics follow the same procedure as for previous T2K analyses [27]. Scale and resolution uncer-tainties are applied to the momentum extrapolated from the track curvature. Additional uncertainties are applied to account for distortions in the magnetic field, charge misidentification, TPC cluster reconstruction, and TPC track reconstruction. Studies have shown that the PØD track reconstruction and matching PØD tracks to TPC tracks have an efficiency greater than 99.8%, so no correction or uncertainty is applied due to this.
In some cases, an interaction that occurred outside the fiducial volume will be reconstructed in the fiducial volume. Using the nominal NEUT simulation, such out of fiducial volume events represent 2.89% of events in the water-in sample and 3.95% of events in the water-out sample. An analysis of the spatial distribution of events in the layers of the PØD central electromagnetic calorimeter, immediately downstream of the water target, suggests that there may be as much as a factor of 2 discrepancy in the amount of migration between data and simulation (but this is a small fraction of the total number of events). To account for this, we increase the weights of events originating in the central electromagnetic calorimeter by a factor of 2 and add an uncertainty of 50% of the initial weight. The distribu-tions of these events in the two PØD geometries are identical, as there is no water in these layers. These events will cancel out up to small differences in statistics and the detector response, so a large uncertainty corresponds to a small effect on the final result. The detector systematics are summarized in TableV
The one-track selection is sensitive to pileup where particles from several interactions are reconstructed in the PØD in the same bunch. The main pileup sources are cosmic ray muons, multiple neutrino interactions in the same beam bunch, and interactions occurring in the material outside ND280 (called sand events) that produce TABLE IV. The fractional uncertainties on the restricted phase
space cross section due to several of the individual cross section parameters.
Source Uncertainty [%] on cross section
MQEA 0.8
MRESA 0.3
CA5 single-pion interaction 0.5 2p2h normalization (16O) 2.8
Eb (16O) 0.05
tracks in the PØD. The cosmic ray muon pileup rate is estimated to be approximately2 × 10−4per bunch, which is negligible. The beam-related pileup rate is dominated by sand events, and is calculated separately for each ND280 run period. The sand muon corrections are given in Table VI.
Finally, there is some uncertainty in the total target mass. Mass uncertainties change the total number of events expected as well as the total number of target nuclei used to obtain the correct cross section normalization. Because the PØD is composed of several materials, the mass uncertainty also has a small effect on the smearing matrix, since the cross section shape will not necessarily have the exact same shape for different materials. In the case of the water in the fiducial volume, the total mass of the water is known very well, but there is still uncertainty in how the water is distributed throughout the different layers. Measurements have shown that there is approximately 15% more water in the most downstream water layer compared to more centrally located layers. There is also a smaller uncertainty on the amount of water in individual layers away from the edges. To estimate the uncertainty due to the mass, a 0.8% uncertainty on the total water mass is applied with additional 15% uncertainties in the layers at the ends of the water target and 5% uncertainties on the
remaining layers. The layer-by-layer variations are per-formed keeping the total mass fixed. These give us an uncertainty of 1.5% due to the target mass.
For all systematic errors, correlations between the water-in and water-out periods were taken water-into consideration water-in Eq. (3). This was done by ensuring identical seeding for throws in calculating Eq. (5). Since flux uncertainties do not distinguish between water-in and water-out samples the flux uncertainties should be fully correlated. The oxygen binding energy parameter affects only the water-in sample, which means that the overall cross section uncertainties are not fully correlated. Some PØD-specific detector systematics, such as the PØD target mass system-atic, have dependence on the water state. Thus detector systematics are not fully correlated either. Further, any differences between water-in and water-out detectors are taken into account with different standard deviations on the underlying variations. Statistical errors are treated as uncorrelated.
IV. RESULTS
The result shown in Fig.6uses data from T2K Runs 2–4. It is reported as a double differential cross section in the outgoing muon kinematicsðpμ; cos θμÞ. Black data points taken from Eq.(4)show the double differential result with full error bars. The colored error bars show the cumulative uncertainties from various sources, starting with the data statistics and ending with the detector systematics. Errors from each source are added in quadrature. MC predictions fromNEUT5.3.2(tuned) andGENIE2.8.0are shown as solid and dashed blue lines.
The fractional contribution from each source of uncer-tainty is shown in Fig. 7, binned and plotted in the same scheme as Fig.6. In most regions, the statistical error from the data is the single most dominant source of bin-by-bin uncertainty, but overall the statistical errors are comparable to the total systematic uncertainty. Aside from the low sensitivity bins, the fractional bin-by-bin errors lie on the order of 10%–20%.
This result is compared to the T2K CC0π cross section on C8H8 [9] in Fig. 8. Since both results were obtained using the T2K flux, the error bars shown do not include the flux uncertainty. The binnings do not exactly match between the two analyses, but when overlaid, both results appear to be consistent within uncertainties for most of the phase space. The largest areas of discrepancy are in the high-angle regions, where the CC0π water cross section is higher than the hydrocarbon cross section.
Finally, the result is compared to the Martini et al.[35]
model predictions on carbon in Fig. 9 and the SuSAv2 model predictions on water[36,37]in Fig.10. Both models are shown with and without 2p2h. Overall, our result agrees better with an inclusion of 2p2h in most regions of phase space. This is consistent with a similar comparison per-formed for the measurement on C8H8[9].
TABLE V. Detector systematic uncertainties considered in this analysis. The uncertainty in each category gives the overall uncertainty on the restricted phase space cross section, calculated for 100 throws. The final result in Eq.(9)throws all uncertainties simultaneously and uses a larger number of throws.
Systematic Uncertainty [%]
PØD energy loss scale 1.3
PØD energy loss resolution 6.7
Tracker momentum scale 1.5
Tracker momentum resolution 0.2 Magnetic field distortion 0.04
Charge ID efficiency 0.1
TPC cluster efficiency 0.3
TPC track efficiency 0.4
Out of fiducial volume 0.8
Detector mass 1.5
TABLE VI. Sand event pileup rates for each ND280 run period. The weights of simulated events are reduced by a factor of 1 − psand.
Run psand [bunch−1]
2 water 0.0144
2 air 0.0173
3c 0.0188
4 water 0.0217
The unfolding procedure allows events to migrate into and out of regions in the phase space outside the binned area. Specifically, this“out-of-range” region refers to where cosθμ< 0 or pμ> 5 GeV. As the selection efficiency here
is very low, we calculated a reduced, total cross section by integrating over only to the explicitly binned areas of the ðpμ; cos θμÞ phase space. The reduced, total cross section was found to be 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 NEUT 5.3.2 (tuned) T2K result GENIE 2.8.0 < 0.600 μ θ cos ≤ 0.000 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.5 1 1.5 2 2.5 3 3.5 < 0.700 μ θ cos ≤ 0.600 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.5 1 1.5 2 2.5 < 0.800 μ θ cos ≤ 0.700 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 < 0.850 μ θ cos ≤ 0.800 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 < 0.900 μ θ cos ≤ 0.850 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 < 0.925 μ θ cos ≤ 0.900 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2d 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 < 0.975 μ θ cos ≤ 0.925 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2d 0 0.2 0.4 0.6 0.8 1 < 1.000 μ θ cos ≤ 0.975 [GeV] μ p [GeV] μ p [GeV] μ p [GeV] μ p [GeV] μ p [GeV] μ p [GeV] μ p [GeV] μ p
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 -2 10 -1 10 detector mass flux cross-section fsi statistical (mc) statistical (data) < 0.600 μ θ cos ≤ 0.000 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 -2 10 -1 10 < 0.700 μ θ cos ≤ 0.600 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 -2 10 -1 10 < 0.800 μ θ cos ≤ 0.700 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 -2 10 -1 10 < 0.850 μ θ cos ≤ 0.800 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 -2 10 -1 10 < 0.900 μ θ cos ≤ 0.850 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 -2 10 -1 10 < 0.925 μ θ cos ≤ 0.900 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 -2 10 -1 10 < 0.975 μ θ cos ≤ 0.925 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 -2 10 -1 10 < 1.000 μ θ cos ≤ 0.975 [GeV] μ p [GeV] μ p [GeV] μ p [GeV] μ p [GeV] μ p [GeV] μ p [GeV] μ p [GeV] μ p
Fractional Error Fractional Error
Fractional Error Fractional Error
Fractional Error Fractional Error
Fractional Error Fractional Error
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2d 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 on water π CC0 on carbon (FGD1) π CC0 < 0.600 μ θ cos ≤ 0.000 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2d 0 0.5 1 1.5 2 2.5 3 3.5 < 0.700 μ θ cos ≤ 0.600 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.5 1 1.5 2 2.5 < 0.800 μ θ cos ≤ 0.700 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 < 0.850 μ θ cos ≤ 0.800 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.5 1 1.5 2 2.5 < 0.900 μ θ cos ≤ 0.850 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 < 0.94 μ θ cos ≤ 0.90 < 0.925 μ θ cos ≤ 0.900 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 < 0.98 μ θ cos ≤ 0.94 < 0.975 μ θ cos ≤ 0.925 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 < 1.00 μ θ cos ≤ 0.98 < 1.000 μ θ cos ≤ 0.975 [GeV] μ p pμ [GeV] [GeV] μ p [GeV] μ p [GeV] μ p [GeV] μ p [GeV] μ p [GeV] μ p
FIG. 8. Here, the CC0π water cross section from the PØD is overlaid with the CC0π C8H8 cross section previously measured by T2K[9]. The error bars shown include all sources of uncertainty except the flux. In regions where cosθμslices are different between the water and hydrocarbon result, the bin edges are noted in the legend.
0 0.5 1 1.5 2 2.5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 on water π CC0
Martini CCQE w. RPA on C
Martini CCQE w. RPA+2p2h on C < 0.600 μ θ cos ≤ 0.000 0 0.5 1 1.5 2 2.5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.5 1 1.5 2 2.5 3 3.5 < 0.700 μ θ cos ≤ 0.600 0 0.5 1 1.5 2 2.5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.5 1 1.5 2 2.5 < 0.800 μ θ cos ≤ 0.700 0 0.5 1 1.5 2 2.5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 < 0.850 μ θ cos ≤ 0.800 0 0.5 1 1.5 2 2.5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 < 0.900 μ θ cos ≤ 0.850 0 0.5 1 1.5 2 2.5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2 d 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 < 0.925 μ θ cos ≤ 0.900 0 0.5 1 1.5 2 2.5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2d 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 < 0.975 μ θ cos ≤ 0.925 0 0.5 1 1.5 2 2.5 /GeV/neutron) 2 cm -38 (10 dpθ dcos σ 2d 0 0.2 0.4 0.6 0.8 1 < 1.000 μ θ cos ≤ 0.975 [GeV] μ p [GeV] μ p [GeV] μ p [GeV] μ p [GeV] μ p pμ [GeV] [GeV] μ p [GeV] μ p
FIG. 9. A comparison of the CC0π water cross section against two Martini model predictions on carbon, one with 2p2h contributions and one without.
FIG. 10. A comparison of the CC0π water cross section against the SuSAv2 model. Predictions for CCQE and 2p2h on oxygen are shown separately along with their sums and average over the bins.
σCC0π
νμH2O¼ 0.95 0.08ðstatÞ 0.06ðdet systÞ
0.04ðmodel systÞ 0.08ðfluxÞ × 10−38cm2=n: ð9Þ This is significantly higher than the NEUT (GENIE) prediction of 0.66ð0.68Þ × 10−38cm2=n primarily due to the disagreement between data and MC in the high-angle regions, which cover a large portion of the reduced phase space. The breakdown of the fractional uncertainties on Eq. (9)are given in Table VII.
V. CONCLUSION
Using the T2K near detector, ND280, a CCQE-like, flux-integrated cross section on water was reported here in the double-differential phase space of the outgoing muon. Since the ¯νμ contribution to this sample is insignificant, this is effectively a measurement on oxygen. The result complements previous T2K cross section measurements. A comparison to a double-differential cross section result on carbon[9]shows good agreement with the exception of a few low momentum bins in the high-angle region. Comparisons were also performed to T2K simulations using NEUT[21]and GENIE[34]. Overall it appears that the tuned NEUT prediction is favored over GENIE, but in
the angular regions of0.7 < cos θμ< 0.85 it shows better agreement with GENIE. Another comparison to the Martini et al.[35]CCQE prediction with and without 2p2h prefers the 2p2h contribution. A comparison to the SuSAv2 model predictions [36,37] also generally agrees with the data within errors, though the data points tend to be lower than the model in the more forward regions.
New techniques for CCQE measurements, such as exploring the transverse kinematics space, or incorporating proton kinematics are ongoing in T2K. Future analyses from T2K will have the benefit of higher statistics and can also include antineutrino cross sections.
The fluxes used for this result, the extracted cross section values and covariance matrices for the errors are available at[38].
ACKNOWLEDGMENTS
We thank the J-PARC staff for superb accelerator performance. We thank the CERN NA61/SHINE Collaboration for providing valuable particle production data. We acknowledge the support of MEXT, Japan; NSERC (Grant No. SAPPJ-2014-00031), NRC and CFI, Canada; CEA and CNRS/IN2P3, France; DFG, Germany; INFN, Italy; National Science Centre (NCN) and Ministry of Science and Higher Education, Poland; RSF, RFBR, and MES, Russia; MINECO and ERDF funds, Spain; SNSF and SERI, Switzerland; STFC, UK; and DOE, USA. We also thank CERN for the UA1/NOMAD magnet, DESY for the HERA-B magnet mover system, NII for SINET4, the WestGrid and SciNet consortia in Compute Canada, and GridPP in the United Kingdom. In addition, participation of individual researchers and institutions has been further supported by funds from ERC (FP7), H2020 Grant No. RISE-GA644294-JENNIFER, EU; JSPS, Japan; Royal Society, UK; the Alfred P. Sloan Foundation and the DOE Early Career program, USA.
[1] K. Abe et al. (T2K Collaboration),Nucl. Instrum. Methods Phys. Res., Sect. A 659, 106 (2011).
[2] K. Abe et al. (T2K Collaboration),Phys. Rev. D 85, 031103 (2012).
[3] K. Abe et al. (T2K Collaboration), Phys. Rev. Lett. 112, 061802 (2014).
[4] Y. Fukuda et al. (Super-Kamiokande Collaboration),Nucl. Instrum. Methods Phys. Res., Sect. A 501, 418 (2003). [5] K. Abe et al. (T2K Collaboration), Phys. Rev. Lett. 112,
061802 (2014).
[6] A. A. Aguilar-Arevalo et al. (MiniBooNE Collaboration), Phys. Rev. D 81, 092005 (2010).
[7] J. L. Alcaraz-Aunion and J. Walding (SciBooNE Collabo-ration),AIP Conf. Proc. 1189, 145 (2009).
[8] G. A. Fiorentini et al. (MINERvA Collaboration), Phys. Rev. Lett. 111, 022502 (2013).
[9] K. Abe et al. (T2K Collaboration),Phys. Rev. D 93, 112012 (2016).
[10] K. Abe et al. (T2K Collaboration),Phys. Rev. D 92, 112003 (2015).
[11] R. Gran et al. (K2K Collaboration), Phys. Rev. D 74, 052002 (2006).
[12] T. Sekiguchi et al. (T2K Collaboration), Nucl. Instrum. Methods Phys. Res., Sect. A 789, 57 (2015).
TABLE VII. Fractional uncertainties on the total cross section.
Source Uncertainty [%] Statistics (data) 7.7 Statistics (MC) 1.8 Detector mass 1.5 Detector other 6.3 FSI 3.4 Cross section 3.2 Flux 8.7
[13] N. Abgrall et al. (NA61/SHINE Collaboration),Eur. Phys. J. C 76, 84 (2016).
[14] N. Abgrall et al. (NA61/SHINE Collaboration),J. Instrum. 9, P06005 (2014).
[15] K. Abe, N. Abgrall, H. Aihara et al. (T2K Collaboration), Phys. Rev. D 87, 012001 (2013).
[16] S. Assylbekov et al. (T2K ND280 PØD Collaboration), Nucl. Instrum. Methods Phys. Res., Sect. A 686, 48 (2012). [17] N. Abgrall et al. (T2K ND280 TPC Collaboration),Nucl.
Instrum. Methods Phys. Res., Sect. A 637, 25 (2011). [18] P. A. Amaudruz et al. (T2K ND280 FGD Collaboration),
Nucl. Instrum. Methods Phys. Res., Sect. A 696, 1 (2012). [19] D. Allan et al. (T2K ND280 ECAL Collaboration),J.
Ins-trum. 8, P10019 (2013).
[20] S. Aoki et al.,Nucl. Instrum. Methods Phys. Res., Sect. A 698, 135 (2013).
[21] Y. Hayato,Nucl. Phys. B, Proc. Suppl. 112, 171 (2002). [22] S. Agostinelli et al. (GEANT4 Collaborations), Nucl.
Instrum. Methods Phys. Res., Sect. A 506, 250 (2003). [23] J. Allison et al.,IEEE Trans. Nucl. Sci. 53, 270 (2006). [24] J. A. Formaggio and G. P. Zeller,Rev. Mod. Phys. 84, 1307
(2012).
[25] G. D’Agostini,Nucl. Instrum. Methods Phys. Res., Sect. A 362, 487 (1995).
[26] T. Adye,arXiv:1105.1160.
[27] K. Abe et al. (T2K Collaboration),Phys. Rev. D 91, 072010 (2015).
[28] L. Salcedo, E. Oset, M. Vicente-Vacas, and C. Garcia-Recio, Nucl. Phys. A484, 557 (1988).
[29] Y. Hayato, Acta Phys. Pol. B 40, 2477 (2009). [30] C. Wilkinson et al.,Phys. Rev. D 93, 072010 (2016). [31] R. A. Smith and E. J. Moniz,Nucl. Phys. B43, 605 (1972);
B101, 547(E) (1975).
[32] J. Nieves, I. R. Simo, and M. J. Vicente Vacas,Phys. Rev. C 83, 045501 (2011).
[33] R. Gran, J. Nieves, F. Sanchez, and M. J. Vicente Vacas, Phys. Rev. D 88, 113007 (2013).
[34] C. Andreopoulos, C. Barry, S. Dytman, H. Gallagher, T. Golan, R. Hatcher, G. Perdue, and J. Yarba, arXiv: 1510.05494.
[35] M. Martini, M. Ericson, G. Chanfray, and J. Marteau,Phys. Rev. C 80, 065501 (2009).
[36] J. E. Amaro, M. B. Barbaro, J. A. Caballero, T. W. Donnelly, A. Molinari, and I. Sick, Phys. Rev. C 71, 015501 (2005).
[37] G. D. Megias, J. E. Amaro, M. B. Barbaro, J. A. Caballero, and T. W. Donnelly,Phys. Rev. D 94, 013012 (2016). [38]
http://t2k-experiment.org/wp-content/uploads/nd280data-numu-cc0pi-xs-on-h2o-2017.tar.gz.
Correction: The names of ten authors were missing in the original publication and have been inserted.
