Measurement of the $\nu_{\mu}$ charged-current cross sections on water, hydrocarbon, iron, and their ratios with the T2K on-axis detectors

TL;DR

This paper reports precise measurements of muon neutrino charged-current cross sections on water, hydrocarbon, and iron at 1.5 GeV, including their ratios, confirming current interaction models with high accuracy.

Contribution

First high-precision measurement of neutrino cross sections on water and their ratios with hydrocarbon and iron at 1.5 GeV in the T2K experiment.

Findings

Cross sections measured with unprecedented precision.

Ratios of cross sections are consistent with theoretical models.

Results support current neutrino interaction models used in oscillation analyses.

Abstract

We report a measurement of the flux-integrated $\nu_{\mu}$ charged-current cross sections on water, hydrocarbon, and iron in the T2K on-axis neutrino beam with a mean neutrino energy of 1.5 GeV. The measured cross sections on water, hydrocarbon, and iron are $\sigma^{\rm{H_{2}O}}_{\rm{CC}}$ = (0.840$\pm 0.010$(stat.)$^{+0.10}_{-0.08}$(syst.))$\times$10$^{-38}$cm$^2$/nucleon, $\sigma^{\rm{CH}}_{\rm{CC}}$ = (0.817$\pm 0.007$(stat.)$^{+0.11}_{-0.08}$(syst.))$\times$10$^{-38}$cm$^2$/nucleon, and $\sigma^{\rm{Fe}}_{\rm{CC}}$ = (0.859$\pm 0.003$(stat.) $^{+0.12}_{-0.10}$(syst.))$\times$10$^{-38}$cm$^2$/nucleon respectively, for a restricted phase space of induced muons: $\theta_{\mu}<45^{\circ}$ and $p_{\mu}>$0.4 GeV/$c$ in the laboratory frame. The measured cross section ratios are ${\sigma^{\rm{H_{2}O}}_{\rm{CC}}}/{\sigma^{\rm{CH}}_{\rm{CC}}}$ = 1.028$\pm 0.016$(stat.)$\pm 0.053$(syst.), ${\sigma^{\rm{Fe}}_{\rm{CC}}}/{\sigma^{\rm{H_{2}O}}_{\rm{CC}}}$ = 1.023$\pm 0.012$(stat.)$\pm 0.058$(syst.), and ${\sigma^{\rm{Fe}}_{\rm{CC}}}/{\sigma^{\rm{CH}}_{\rm{CC}}}$ = 1.049$\pm 0.010$(stat.)$\pm 0.043$(syst.). These results, with an unprecedented precision for the measurements of neutrino cross sections on water in the studied energy region, show good agreement with the current neutrino interaction models used in the T2K oscillation analyses.

Figures

Tables

Table 1. Table 1: Summary of detector specifications of the Water Module, Proton Module, and one of the INGRID modules. The target masses are calculated inside the fiducial volumes, which correspond to the effective target masses and are specially tuned for this analysis, as described in Sec.  5 .

Parameter	Water Module	Proton Module	INGRID module
Target mass in fiducial volume (ton)	0.10	0.16	2.1
Main target materials and fraction	$H_{2}O$ (80%), $\mathrm{CH}$ (19%)	CH (98%)	Fe (96%)
Dimension of a scintillator (${\mathrm{cm}}^3$)	100$\times$2.5$\times$0.3	120$\times$2.5$\times$1.3 (SciBar-type),	120$\times$5$\times$1
		120$\times$5$\times$1 (INGRID-type)
Dimension of an iron plane (${\mathrm{cm}}^3$)	-	-	124$\times$124$\times$6.5
The number of readout channels	1280	1204	616
MPPC serial number	S13660	S10362-13-050C	S10362-13-050C
MPPC gain stability	10%	10%	10%
MPPC dark noise rate (hits/module/bunch)	0.2	12	6
Mean scintillator light yield for MIP	16	56 (SciBar-type),	23
(p.e. per scintillator thickness)		23 (INGRID-type)
Angular acceptance respect to beam axis	0∘ to 90∘	0∘ to 75∘	0∘ to 60∘
Period located at on-axis position	July 2016-	November 2010-May 2016	2009-

Table 2. Table 2: Nominal models of the neutrino-nucleus interactions implemented in NEUT used in this analysis.

Mode	Nominal model	Parameter
CCQE-like	Dipole type axial form factor	$M_A^{QE}=1.15$ GeV/$c^2$.
	RFG model by Smith-Moniz cite:rfg	$E_b$= 25, 27, 33 MeV and
	with binding energy ($E_b$) and Fermi surface momentum ($p_F$)	$p_F$= 217, 225, 250 MeV/c for
		12C , 16O, and 56Fe, respectively.
	RPA model by Nieves et al. cite:rpa	RPA is applied for 12O and 16C.
		RPA is not applied for 56Fe.
	2p2h model by Nieves et al. cite:mecnieves	Normalization
1$\pi$	Model by Rein-Sehgal cite:reinsehgal	$C_5^A\left(0\right)=1.01$,
		$M_A^{RES}=0.95$ GeV/$c^2$,
		$\mathrm{Isospin}\frac{1}{2}\mathrm{bg}=1.30$.
DIS	PYTHIA cite:pytia , Parton distribution function by	Energy dependent normalization
	GRV98 with Bodek and Yang correction Gluck:1998xa  Bodek:2003wc  Bodek:2004pc
Coherent	Model by Berger-Sehgal cite:cccohbarbarseghal	Normalization

Table 3. Table 3: Flux-integrated CC cross sections per nucleon for ν μ subscript 𝜈 𝜇 \nu_{\mu} on Fe, CH, and H 2 ​ O subscript H 2 O \rm H_{2}O simulated by NEUT. Neutrino interaction parameters used for the simulation are listed in Table  2 . Because RPA for Fe is not implemented in NEUT at present, the expectation of σ F ​ e subscript 𝜎 𝐹 𝑒 \sigma_{Fe} with RPA is not listed.

Cross section	NEUT expectation with RPA	NEUT expectation without RPA
${\sigma }_{H_{2}O}$	$0.819\times {10}^{-38}{\mathrm{cm}}^2$	$0.860\times {10}^{-38}{\mathrm{cm}}^2$
${\sigma }_{\mathrm{CH}}$	$0.832\times {10}^{-38}{\mathrm{cm}}^2$	$0.875\times {10}^{-38}{\mathrm{cm}}^2$
${\sigma }_{\mathrm{Fe}}$	not available	$0.904\times {10}^{-38}{\mathrm{cm}}^2$
${\sigma }_{H_{2}O}/{\sigma }_{\mathrm{CH}}$	0.984	0.983
${\sigma }_{\mathrm{Fe}}/{\sigma }_{H_{2}O}$	not available	1.051
${\sigma }_{\mathrm{Fe}}/{\sigma }_{\mathrm{CH}}$	not available	1.033

Table 4. Table 4: Parameters used for the event selection criteria for the on-axis detectors.

	Water Module	Proton Module	INGRID module
Time clustering	$\pm$50 nsec	$\pm$50 nsec	$\pm$50 nsec
Track matching with INGRID	$\pm$35∘	$\pm$35∘	-
	$\pm$150 mm	$\pm$150 mm	-
3D track matching	$\le$1 parallel plane	$\le$1 plane	$\le$1 plane
Vertexing	3 planes	2 planes	2 planes
	150 mm	150 mm	150 mm
Beam timing	$\pm$100 nsec	$\pm$100 nsec	$\pm$100 nsec
Upstream veto	$\ge$second parallel plane	$\ge$second plane	$\ge$first plane
Fiducial	700 mm$\times$700 mm	700 mm$\times$700 mm	700 mm$\times$700 mm
Reconstructed angle	45∘	45∘	45∘

Table 5. Table 5: Summary of the event selection for the Water Module. The purities of CC interactions are shown in parentheses.

Selection	Data					MC
		CC	NC	${\overline{\nu }}_{\mu }$, ${\nu }_e$, ${\overline{\nu }}_e$	CH B.G.	Wall B.G.	INGRID B.G.	All
Vertexing cut	1175980	$4.39\times {10}^4$ (4%)	1.66$\times {10}^2$	$1.12\times {10}^3$	$1.08\times {10}^4$	$9.10\times {10}^5$	$2.77\times {10}^5$	$1.24\times {10}^6$
Front veto cut	100790	$2.77\times {10}^4$ (21%)	1.04$\times {10}^3$	$9.38\times {10}^2$	$6.66\times {10}^3$	$8.09\times {10}^4$	$1.46\times {10}^4$	$1.32\times {10}^5$
Fiducial cut	17992	$1.25\times {10}^4$ (69%)	4.68$\times {10}^2$	$4.42\times {10}^2$	$3.51\times {10}^3$	$3.49\times {10}^2$	$5.84\times {10}^2$	$1.78\times {10}^4$
Track angle cut	17528	$1.20\times {10}^4$ (69%)	4.53$\times {10}^2$	$4.39\times {10}^2$	$3.39\times {10}^3$	$3.47\times {10}^2$	$5.64\times {10}^2$	$1.73\times {10}^4$

Table 6. Table 6: Summary of the event selection for the Proton Module. The purities of CC interactions are shown in parentheses.

Selection	Data					MC
		CC	NC	${\overline{\nu }}_{\mu }$, ${\nu }_e$, ${\overline{\nu }}_e$	Wall B.G.	INGRID B.G.	All
Vertexing cut	1321290	$5.56\times {10}^4$ (4%)	2.66$\times {10}^3$	$2.00\times {10}^3$	$1.03\times {10}^6$	$2.77\times {10}^5$	$1.36\times {10}^6$
Front veto cut	264550	$4.69\times {10}^4$ (15%)	2.25$\times {10}^3$	$1.72\times {10}^3$	$2.17\times {10}^5$	$3.63\times {10}^4$	$3.04\times {10}^5$
Fiducial cut	22930	$1.98\times {10}^4$ (85%)	9.52$\times {10}^2$	$7.31\times {10}^2$	$5.54\times {10}^2$	$9.97\times {10}^2$	$2.32\times {10}^4$
Track angle cut	22165	$1.92\times {10}^4$ (85%)	9.14$\times {10}^2$	$7.26\times {10}^2$	$5.51\times {10}^2$	$9.50\times {10}^2$	$2.23\times {10}^4$

Table 7. Table 7: Summary of the event selection for the INGRID module. The INGRID B.G. in the table represents backgrounds from the other INGRID modules. The purities of CC interactions are shown in parentheses.

Selection	Data				MC
		CC	NC	${\overline{\nu }}_{\mu }$, ${\nu }_e$, ${\overline{\nu }}_e$	CH B.G.	Wall B.G.	INGRID B.G.	All
Vertexing cut	3019430	$1.11\times {10}^6$ (44%)	$6.98\times {10}^4$	$3.20\times {10}^4$	$4.49\times {10}^4$	$9.45\times {10}^5$	$3.36\times {10}^5$	$2.54\times {10}^6$
Front veto cut	1468490	$1.07\times {10}^6$ (74%)	$6.74\times {10}^4$	$3.07\times {10}^4$	$3.97\times {10}^4$	$1.98\times {10}^5$	$4.33\times {10}^4$	$1.45\times {10}^6$
Fiducial cut	431211	$4.10\times {10}^5$ (88%)	$2.58\times {10}^4$	$1.14\times {10}^4$	$1.49\times {10}^4$	$1.52\times {10}^3$	$1.06\times {10}^2$	$4.65\times {10}^5$
Acceptance cut	308971	$2.88\times {10}^5$ (88%)	$1.81\times {10}^4$	$9.56\times {10}^3$	$1.07\times {10}^4$	$9.26\times {10}^2$	$6.73\times {10}^2$	$3.28\times {10}^5$
Track angle cut	293418	$2.74\times {10}^5$ (88%)	$1.72\times {10}^4$	$9.31\times {10}^3$	$1.02\times {10}^4$	$8.70\times {10}^2$	$6.38\times {10}^2$	$3.12\times {10}^5$

Table 8. Table 8: The number of selected events for the INGRID module before and after the pileup correction.

reconstructed angle bin	$N_{\mathrm{sel}}$	$N_{\mathrm{corr}}$	$N_{\mathrm{corr}}/N_{\mathrm{sel}}$
0-5∘	13106	13582.0	1.036
5-10∘	32928	33765.3	1.025
10-15∘	52272	53671.3	1.027
15-20∘	54205	55500.6	1.024
20-25∘	38540	39119.4	1.015
25-30∘	44097	45002.4	1.021
30-35∘	26615	26984.1	1.014
35-40∘	19709	20036.4	1.017
40-45∘	11946	12094.0	1.012
Total	293418	299755.5	1.022

Table 9. Table 9: Summary of the fraction of the backgrounds after the event selection. Non-target element backgrounds are neutrino interactions on neither CH nor H 2 ​ O subscript H 2 O \rm H_{2}O for the Water Module, on O, N, and Ti for the Proton Module, and on scintillators for the INGRID module.

Detector	angle bin	CC out of	Non-target	NC	${\overline{\nu }}_{\mu }$, ${\nu }_e$, ${\overline{\nu }}_e$	Wall	INGRID	All BG
		phase space	element
Water	0-5∘	44.5	26.1	28.6	43.5	4.90	55.2	216
Module	5-10∘	98.2	55.1	67.4	99.2	36.7	96.6	477
	10-15∘	145	72.0	83.7	103	73.8	10.3	615
	15-20∘	171	76.3	86.4	75.6	113	90.0	654
	20-25∘	165	58.2	76.7	51.7	58.9	77.2	527
	25-30∘	113	43.6	54.8	30.9	32.7	72.9	377
	30-35∘	84.4	27.6	32.4	19.6	13.2	33.4	229
	35-40∘	35.0	15.6	16.3	10.9	12.4	25.8	126
	40-45∘	40.2	7.70	6.99	4.72	1.24	9.74	82.4
	Total	896	382	453	439	3.47	564	3300
Proton	0-5∘	99.0	12.9	60.3	79.7	38.2	57.2	346
Module	5-10∘	255	35.7	145	172	47.4	154	905
	10-15∘	338	48.0	174	162	75.2	183	975
	15-20∘	352	49.1	177	129	145	150	997
	20-25∘	313	43.3	144	78.8	104	124	803
	25-30∘	243	34.5	101	50.6	83.7	107	616
	30-35∘	148	25.3	63.4	30.1	23.9	90.4	379
	35-40∘	67.6	16.6	32.4	15.2	20.5	56.1	207
	40-45∘	83.5	9.69	17.3	8.96	12.4	28.3	159
	Total	1870	275	914	726	551	950	5290
INGRID	0-5∘	1370	507	769	766	95.7	7.96	3540
module	5-10∘	2910	1310	1690	1740	145	101	7900
	10-15∘	4990	1990	2680	2020	147	122	11900
	15-20∘	5630	2020	3280	1720	114	216	13000
	20-25∘	3990	1440	2100	1010	109	49.0	8690
	25-30∘	5520	1680	3070	993	126	88.8	11500
	30-35∘	3320	997	1660	588	58.0	19.7	6650
	35-40∘	3650	702	1170	338	34.9	19.2	5920
	40-45∘	3080	456	801	144	40.2	13.2	4530
	Total	34500	11100	17200	9310	870	638	73600

Table 10. Table 10: Integrated ν μ subscript 𝜈 𝜇 \nu_{\mu} flux in the fiducial volume of each detector.

	Water Module	Proton Module	INGRID module
Integrated ${\nu }_{\mu }$ flux per ${10}^{21}$ POT ($/{\mathrm{cm}}^2$)	$5.13\times {10}^{13}$	$5.13\times {10}^{13}$	$5.08\times {10}^{13}$
Used POT in this analysis	$7.25\times {10}^{20}$	$5.89\times {10}^{20}$	$5.89\times {10}^{20}$
Integrated ${\nu }_{\mu }$ flux per used POT  ($/{\mathrm{cm}}^2$)	$3.72\times {10}^{13}$	$3.02\times {10}^{13}$	$2.99\times {10}^{13}$

Table 11. Table 11: Summary of the number of target nucleons.

	Number of target nucleons
$T_{\mathrm{WM}}^{H_{2}O}$	$4.939\times {10}^{28}$
$T_{\mathrm{WM}}^{\mathrm{CH}}$	$1.090\times {10}^{28}$
$T_{\mathrm{PM}}^{\mathrm{CH}}$	$9.230\times {10}^{28}$
$T_{\mathrm{ING}}^{\mathrm{Fe}}$	$1.206\times {10}^{30}$

Table 12. Table 12: Calculated cross sections using numbers of events expected by MC and their true values with the nominal model.

Target	Calculated cross sections	Expected cross sections
$H_{2}O$	$0.821\times {10}^{-38}{\mathrm{cm}}^2$	$0.819\times {10}^{-38}{\mathrm{cm}}^2$
CH	$0.832\times {10}^{-38}{\mathrm{cm}}^2$	$0.832\times {10}^{-38}{\mathrm{cm}}^2$
Fe	$0.904\times {10}^{-38}{\mathrm{cm}}^2$	$0.904\times {10}^{-38}{\mathrm{cm}}^2$

Table 13. Table 13: Summary of the systematic uncertainties for the cross section measurements (%).

Systematic uncertainty	${\sigma }_{H_{2}O}$	${\sigma }_{\mathrm{CH}}$	${\sigma }_{\mathrm{Fe}}$	${\sigma }_{H_{2}O}/{\sigma }_{\mathrm{CH}}$	${\sigma }_{\mathrm{Fe}}/{\sigma }_{H_{2}O}$	${\sigma }_{\mathrm{Fe}}/{\sigma }_{\mathrm{CH}}$
Flux-related	$+10.8$	$+11.5$	$+13$	$+0.6$	$+1.8$	$+1.1$
(hadron production and beam line)	$-8.9$	$-9.6$	$-11$	$-0.6$	$-1.8$	$-1.2$
Flux-related	$+0.3$	–	–	$+1.3$	$+1.1$	$+0.3$
(difference of running periods and location)	$-0.3$	–	–	$-1.3$	$-1.1$	$-0.3$
Interaction model-related	$+2.6$	$+3.1$	$+5.2$	$+2.3$	$+4.0$	$+2.7$
	$-2.6$	$-3.1$	$-5.2$	$-2.3$	$-4.0$	$-2.7$
Detector response-related	$+2.9$	$+2.5$	$+1.5$	$+4.5$	$+3.4$	$+2.8$
	$-2.9$	$-2.5$	$-1.5$	$-4.5$	$-3.4$	$-2.8$
Total	$+11.5$	$+13$	$+14$	$+5.2$	$+5.7$	$+4.1$
	$-9.7$	$-10$	$-12$	$-5.2$	$-5.7$	$-4.1$

Table 14. Table 14: List of the interaction model parameters and uncertainties used in the analysis.

Parameter	Nominal value	Uncertainties (1$\sigma$)
CCQE-like
$M_A^{\mathrm{QE}}$	1.15 GeV/$c^2$	0.18 GeV/$c^2$
$p_F{}_{}{}^{12}C$	217 MeV/c	31 MeV/c
$p_F{}_{}{}^{16}O$	225 MeV/c	31 MeV/c
$p_F{}_{}{}^{56}\mathrm{Fe}$	250 MeV/c	35 MeV/c
$E_b{}_{}{}^{12}C$	25 MeV/c	9 MeV/c
$E_b{}_{}{}^{16}O$	27 MeV/c	9 MeV/c
$E_b{}_{}{}^{56}\mathrm{Fe}$	33 MeV/c	11 MeV/c
$2p2h{\mathrm{normalization}}^{12}C$	100%	100%
$2p2h{\mathrm{normalization}}^{16}O$	100%	100%
$2p2h{\mathrm{normalization}}^{56}\mathrm{Fe}$	100%	100%
1$\pi$
$C_{\mathrm{A5}}$	1.01	0.12
$M_A^{\mathrm{Res}}$	0.95 GeV/$c^2$	0.15 GeV/$c^2$
$\mathrm{Isospin}\frac{1}{2}\mathrm{bg}$	1.30	0.20
CC multi-pion and DIS production
Normalization uncertainty is applied depending on neutrino energy by 0.4/$E_{\nu }$ (GeV)
CC coherent
$\mathrm{CC}\mathrm{coherent}\mathrm{normalization}{}_{}{}^{12}C$	100%	30%
$\mathrm{CC}\mathrm{coherent}\mathrm{normalization}{}_{}{}^{16}O$	100%	30%
Normalization of NC interactions
$\mathrm{NC}\mathrm{coherent}\mathrm{normalization}$	100%	30%
NC multi-pion and DIS production normalization	100%	30%
Secondary interaction of pions
Pion absorption normalization	100%	50%
Pion charge exchange normalization (500 MeV/$c$)	100%	50%
Pion charge exchange normalization ($p_{\pi }>$500 MeV/$c$)	100%	30%
Pion quasi elastic normalization (500 MeV/$c$)	100%	50%
Pion quasi elastic normalization ($p_{\pi }>$500 MeV/$c$)	100%	30%
Pion inelastic normalization	100%	50%

Table 15. Table 15: Summary of the neutrino interaction model-related uncertainties for each cross section measurement (%). Only the dominant systematic parameters are shown.

Parameter	${\sigma }_{H_{2}O}$	${\sigma }_{\mathrm{CH}}$	${\sigma }_{\mathrm{Fe}}$	${\sigma }_{H_{2}O}/{\sigma }_{\mathrm{CH}}$	${\sigma }_{\mathrm{Fe}}/{\sigma }_{H_{2}O}$	${\sigma }_{\mathrm{Fe}}/{\sigma }_{\mathrm{CH}}$
$M_A^{\mathrm{QE}}$	1.1	1.2	1.6	0.1	0.5	0.5
$M_A^{\mathrm{Res}}$	0.6	1.2	2.0	0.5	1.5	1.0
$C_{\mathrm{A5}}$	0.1	0.3	1.0	0.3	1.0	0.7
$\mathrm{Isospin}\frac{1}{2}\mathrm{bg}$	0.1	0.2	0.6	0.2	0.5	0.3
CC multi-pion and DIS production	0.2	0.7	1.3	0.5	1.1	0.6
NC multi-pion and DIS production	0.8	1.1	1.7	0.3	0.9	0.6
$2p2h{\mathrm{normalization}}^{12}C$	0.1	0.4	0.3	0.4	0.3	0.4
$2p2h{\mathrm{normalization}}^{16}O$	0.2	0	0	0.2	0.2	0
$2p2h{\mathrm{normalization}}^{56}\mathrm{Fe}$	0	0	0.4	0	0.4	0.4
Pion absorption normalization	0.2	0.1	0.2	0.2	0.4	0.2
Pion quasi elastic normalization (500 MeV/$c$)	0.2	0.1	0.1	0.2	0.4	0.4
Pion quasi elastic normalization ($p_{\pi }>$500 MeV/$c$)	0.2	0.1	0.3	0.1	0.1	0.1
Pion inelastic normalization	0.2	0.2	0.3	0.1	0.1	0.2
Backscattered protons and pions	1.5	0.8	1.0	1.8	1.8	1.3
$\mathrm{Total}$	2.6	3.1	5.2	2.3	4.0	2.7

Table 16. Table 16: Summary of the detector systematic uncertainties for the absolute cross section measurements (%).

cross section	${\sigma }_{H_{2}O}$		${\sigma }_{\mathrm{CH}}$	${\sigma }_{\mathrm{Fe}}$
Detector	Water Module	Proton Module	Proton Module	INGRID
Target mass	0.68	0.05	0.27	0.14
MPPC noise	0.01	0.09	0.39	0.09
Scintillator crosstalk	0.30	–	–	–
Hit efficiency	0.27	0.02	0.50	0.94
Event pileup	0.72	0.15	0.64	0.09
Beam-related background	1.09	0.31	1.31	0.38
Non-beam-related background	0.04	0.02	0.02	0.01
2D Track reconstruction	0.60	0.28	1.18	0.43
Track matching with INGRID	1.42	0.20	0.84	–
3D track matching	0.89	0.13	0.56	0.35
Vertexing	0.44	0.05	0.20	0.28
Beam timing cut	0.06	0.01	0.01	0.01
VETO and FV cut	1.19	0.18	0.72	0.52
Acceptance cut	–	–	–	0.61
Total	2.88		2.52	1.54

Table 17. Table 17: Summary of the cross sections predicted by NEUT with the various parameter values listed in Table  14 ( × 10 − 38 ​ cm 2 / nucleon absent superscript 10 38 superscript cm 2 nucleon \times 10^{-38}~{}\rm{cm^{2}/nucleon} for the absolute cross sections).

Parameter	${\sigma }_{H_{2}O}$	${\sigma }_{\mathrm{CH}}$	${\sigma }_{\mathrm{Fe}}$	${\sigma }_{H_{2}O}/{\sigma }_{\mathrm{CH}}$	${\sigma }_{\mathrm{Fe}}/{\sigma }_{H_{2}O}$	${\sigma }_{\mathrm{Fe}}/{\sigma }_{\mathrm{CH}}$
Nominal for $H_{2}O,\mathrm{CH},H_{2}O/\mathrm{CH}$	0.819	0.832	-	0.985	-	-
(RPA+RFG+MEC)
Nominal for $\mathrm{Fe},H_{2}O/\mathrm{Fe},\mathrm{Fe}/\mathrm{CH}$	0.860	0.875	0.904	0.982	1.052	1.034
(RFG+MEC)
$M_A^{\mathrm{QE}}-1\sigma$	0.781	0.792	0.857	0.986	1.046	1.031
$M_A^{\mathrm{QE}}+1\sigma$	0.855	0.869	0.948	0.984	1.058	1.041
$M_A^{\mathrm{Res}}-1\sigma$	0.779	0.791	0.861	0.985	1.050	1.034
$M_A^{\mathrm{Res}}+1\sigma$	0.859	0.873	0.948	0.984	1.054	1.037
$C_{\mathrm{A5}}-1\sigma$	0.789	0.800	0.871	0.986	1.061	1.046
$C_{\mathrm{A5}}+1\sigma$	0.853	0.866	0.941	0.984	1.054	1.037
$\mathrm{Isospin}\frac{1}{2}\mathrm{bg}-1\sigma$	0.806	0.818	0.889	0.985	1.051	1.035
$\mathrm{Isospin}\frac{1}{2}\mathrm{bg}+1\sigma$	0.835	0.848	0.923	0.985	1.053	1.037
$\mathrm{CCother}\mathrm{shape}-1\sigma$	0.797	0.809	0.880	0.985	1.052	1.036
$\mathrm{CCother}\mathrm{shape}+1\sigma$	0.841	0.854	0.928	0.985	1.053	1.037
$p_{F}C-1\sigma$	0.819	0.835	0.904	0.981	1.052	1.028
$p_{F}C+1\sigma$	0.819	0.825	0.904	0.993	1.052	1.045
$p_{F}O-1\sigma$	0.824	0.832	0.904	0.991	1.043	1.034
$p_{F}O+1\sigma$	0.812	0.832	0.904	0.976	1.065	1.034
$\mathrm{Eb}C-1\sigma$	0.819	0.831	0.904	0.986	1.052	1.035
$\mathrm{Eb}C+1\sigma$	0.819	0.833	0.904	0.984	1.052	1.032
$\mathrm{Eb}O-1\sigma$	0.818	0.832	0.904	0.984	1.054	1.034
$\mathrm{Eb}O+1\sigma$	0.820	0.832	0.904	0.986	1.051	1.034
$\mathrm{MEC}\mathrm{norm}C-1\sigma$	0.819	0.764	0.904	1.072	1.052	1.121
$\mathrm{MEC}\mathrm{norm}C+1\sigma$	0.819	0.900	0.904	0.910	1.052	0.959
$\mathrm{MEC}\mathrm{norm}O-1\sigma$	0.754	0.832	0.904	0.906	1.139	1.034
$\mathrm{MEC}\mathrm{norm}O+1\sigma$	0.884	0.832	0.904	1.063	0.978	1.034
$\mathrm{CCcoh}\mathrm{norm}C-1\sigma$	0.819	0.809	0.904	1.013	1.052	1.061
$\mathrm{CCcoh}\mathrm{norm}C+1\sigma$	0.819	0.855	0.904	0.958	1.052	1.007
$\mathrm{CCcoh}\mathrm{norm}O-1\sigma$	0.800	0.832	0.904	0.962	1.076	1.034
$\mathrm{CCcoh}\mathrm{norm}O+1\sigma$	0.838	0.832	0.904	1.007	1.030	1.034
$\mathrm{CCcoh}\mathrm{norm}\mathrm{Fe}-1\sigma$	0.819	0.832	0.896	0.958	1.043	1.023
$\mathrm{CCcoh}\mathrm{norm}\mathrm{Fe}+1\sigma$	0.819	0.832	0.913	0.958	1.062	1.044

Equations

$${\sigma }_A$$

$$N_{\mathrm{WM}}^{\mathrm{CH}\mathrm{BG}}$$

$$U_{ij}$$

$${\sigma }_{\mathrm{CC}}^{{\mathrm{H}}_2\mathrm{O}}=(0.840\pm 0.010{(stat.)}_{-0.08}^{+0.10}(syst.))\times 10^{-38}\mathrm{c}{\mathrm{m}}^2/nucleon,$$

$${\sigma }_{\mathrm{CC}}^{\mathrm{CH}}=(0.817\pm 0.007{(stat.)}_{-0.08}^{+0.11}(syst.))\times 10^{-38}\mathrm{c}{\mathrm{m}}^2/nucleon,$$

$${\sigma }_{\mathrm{CC}}^{\mathrm{Fe}}=(0.859\pm 0.003{(stat.)}_{-0.10}^{+0.12}(syst.))\times 10^{-38}\mathrm{c}{\mathrm{m}}^2/nucleon.$$

$$\frac{{\sigma }_{\mathrm{CC}}^{{\mathrm{H}}_2\mathrm{O}}}{{\sigma }_{\mathrm{CC}}^{\mathrm{CH}}}=1.028\pm 0.016(stat.)\pm 0.053(syst.),$$

$$\frac{{\sigma }_{\mathrm{CC}}^{\mathrm{Fe}}}{{\sigma }_{\mathrm{CC}}^{{\mathrm{H}}_2\mathrm{O}}}=1.023\pm 0.012(stat.)\pm 0.058(syst.),$$

$$\frac{{\sigma }_{\mathrm{CC}}^{\mathrm{Fe}}}{{\sigma }_{\mathrm{CC}}^{\mathrm{CH}}}=1.049\pm 0.010(stat.)\pm 0.043(syst.).$$

Full text

\subjectindex

C04, C32

1 Introduction

The Tokai-to-Kamioka (T2K) experiment Abe:2011ks is a long-baseline neutrino oscillation experiment that started taking physics data in 2010. The T2K experiment studies properties of neutrino oscillations via disappearance of muon (anti-)neutrinos and appearance of electron (anti-)neutrinos from a nearly pure muon (anti-)neutrino beam, which is produced by the J-PARC accelerator complex. The neutrino beam characteristics and neutrino-nucleus interactions are measured with a suite of near detectors, which are situated 280 m from the production target, consisting of the so-called INGRID Abe:2011xv and ND280 pod ; fgd ; Abgrall:2010hi ; Allan:2013ofa ; Aoki:2012mf . The INGRID is placed at the center of the neutrino beam (on-axis), while the ND280 is at an off-axis of 2.5*∘. The neutrino oscillation patterns are observed with the 2.5∘* off-axis far detector, Super-Kamiokande Fukuda:2002uc , which is located 295 km away from the production target. In order to precisely measure neutrino oscillations, understanding of the neutrino interactions with nuclei is essential. In the current T2K oscillation analysis Abe:2018wpn , data samples of charged-current candidates in which the interaction vertex is found in one of two Fine-Grained Detectors, FGD1 or FGD2 fgd , are used to constrain the neutrino flux prediction and cross section models. The former detector consists of 100% plastic scintillators (hydrocarbon) and the latter consists of a mixture of plastic scintillators and water, while the far detector consists of 100% water. The neutrino interaction model is used to extrapolate the near detector spectra to the (oscillated) far detector spectra in a few significant ways. First, the T2K off-axis near detector angular acceptance is more limited than the far detector. Second, the near detector event rate also includes significant interactions on materials other than the far detector (water) target. Finally, the interaction model is tuned at the near detector to predict the far detector energy spectra and this parameterization can be incomplete. Therefore testing the interaction model with different target materials and at various range of the neutrino energies are valuable to the T2K oscillation analysis. However, there are only a few publications of the neutrino cross sections on water so far cc1pifgd ; cc0pipod ; scifi . Two exclusive channels of charged current interactions are measured by the ND280 cc1pifgd ; cc0pipod with approximately 15% uncertainties with a mean neutrino energy of 0.6 GeV. There is only one measurement of axial vector mass scifi with 10% uncertainty with a mean neutrino energy above 1 GeV.

A new water-target neutrino detector, named the Water Module Koga:2015iqa , has been constructed for the precise measurements of neutrino interactions on water with a mean neutrino energy of 1.5 GeV. In this article, by using the Water Module and the other T2K detectors including the Proton Module ccincpm and INGRID Abe:2011xv , we measure the $\nu_{\mu}$ charged-current (CC) cross sections on water, hydrocarbon, iron, and their ratios. Dominant errors of the absolute cross section measurements come from the uncertainty of the T2K neutrino beam prediction, which largely cancels out when performing measurements on their cross section ratios. This method was established in the previous measurement of a cross section ratio between hydrocarbon and iron by using the Proton Module and INGRID ccincpm . In this article, measurements of neutrino interaction on water with the Water Module are conducted for the first time. In addition, in order to reduce the dependence on the Monte Carlo implemented model of neutrino-nucleus interactions in extracting the cross section values, a method for unfolding the total cross section as a function of muon scattering angles is implemented. Hereafter, we will describe the detector configuration, the Monte Carlo simulation, the used data sample, the event selection, the method to extract the cross sections, systematic uncertainties, and the results.

2 Detector configuration

We use the three detectors, INGRID, Proton Module, and Water Module as iron (Fe), hydrocarbon (CH), and water ($\rm H_{2}O$) interaction targets, respectively. Table 1 shows specifications of the three detectors. INGRID consists of fourteen identical modules arranged in a cross shape; each module has a sandwich structure comprised of nine iron planes and eleven tracking planes as shown in Fig. 1. INGRID has been operating since 2009 to monitor the neutrino beam rate, its direction, and stability in real-time. The tracking planes are formed from two layers of scintillator, each of which is composed of twenty-four bars that are oriented either horizontally or vertically. The thickness of the iron planes is 6.5 cm and the thickness of the scintillator is 1.0 cm. The iron planes, which play a role as the neutrino interaction target in this analysis, make up 96% of the total fiducial mass of the module. There are veto planes surrounding the module designed for tracking the charged particles entering into the detector. More detailed information about the INGRID can be found in Abe:2011xv . In this analysis, the central horizontal INGRID module is used as the iron target. The three horizontal INGRID modules surrounding the beam center are used for muon identification for the Proton Module and Water Module.

The Proton Module is a plastic scintillator target detector located between the horizontal and vertical INGRID modules, as shown in Fig. 2. It was built for the measurement of the neutrino cross section on hydrocarbon and it had been located at the on-axis position from November 2010 to May 2016. It consists of thirty-four tracking planes with each plane being an array of thirty-two scintillator bars that are oriented either horizontally or vertically. Two types (SciBar-type and INGRID-type) of scintillator bars, which have different sizes, are used in the inner and outer sections of each tracking plane. Hydrocarbon in the scintillators of the tracking planes serves as the neutrino interaction target and composes 98% of the total fiducial mass of the Proton Module. Similar to the INGRID modules, the Proton Module is composed of veto planes surrounding the tracking planes of the detector. More detailed information about the Proton Module can be found in ccincpm .

The Water Module is a neutrino detector with the interaction target region composed of 80% water and 20% plastic scintillators. The high fraction of water in the detector, in fact higher than the previous water-target neutrino detectors pod ; fgd , is essential to reduce the backgrounds induced by the neutrino interactions on non-water materials. The Water Module has been located at the on-axis position between the INGRID horizontal modules and vertical modules since June 2016, replacing the Proton Module. The Water Module consists of a stainless steel tank filled with water and sixteen scintillator tracking planes immersed in the water, as shown in Fig. 3. The eight tracking planes are placed alternately in the x-direction and y-direction along the z-direction so that three-dimensional tracks may be reconstructed. Each tracking plane is an array of eighty scintillator bars. The forty bars, called parallel scintillators, are placed along the xy-direction. The other forty bars, called grid scintillators, are placed along the z-direction with a grid-like structure in order to achieve a large angular acceptance. The plastic scintillators of dimension 100 cm (length) $\times$ 2.5 cm (width) $\times$ 0.3 cm (thickness) were produced in the Fermi National Accelerator Laboratory Beznosko:2005ba . The scintillators are made of polystyrene, infused with PPO (1%) and POPOP (0.03%). The manufactured scintillator, co-extruded with a white reflective coating of $\rm{TiO_{2}}$ infused in polystyrene, has a rectangular cross section with a groove to house a wavelength shifting (WLS) fiber (Kuraray Y-11 fiber ). The WLS fiber is glued onto the scintillator with optical cement (ELJEN TECHNOLOGY EJ-500 cement ). The surface of the scintillator is painted by a black cement of acrylic silicon to prevent optical crosstalk between the scintillators. Each layer of scintillator bars is affixed to a mechanical frame which sits inside a water tank. Spaces between scintillators are filled with water. Scintillation light from the scintillator is collected by the WLS fiber and detected by a multi-pixel photon counter (MPPC) mppc , similar to that for the INGRID and Proton Module. While the Hamamatsu S10362-13-050C MPPC was used in the INGRID and Proton Module, a newer type of MPPC, S13660 with higher gain, lower noise rate, crosstalk rate, and after-pulse rate, is used in the Water Module. The same Trip-t electronics Tript2007 are used for all the three detectors. To record data from the neutrino beam, delivered typically in eight bunches with a cycle of 581 ns for each 2.48 sec, a trigger from the J-PARC accelerator is provided to each detector. The integrated charge and hit timing of all channels are digitized and recorded with 2.5 photoelectrons (p.e.) threshold for each beam bunch.

3 Monte Carlo simulation

A Monte Carlo (MC) simulation is used for the estimation of background contamination and signal detection efficiency. Three pieces of software are used for the chain of simulation: JNUBEAM jnubeam for the neutrino flux prediction, NEUT NEUT for the neutrino interactions with nuclei, and a GEANT4 geant4 -based detector simulation. JNUBEAM simulates the interaction of 30 GeV primary protons on a graphite target, the propagation of the secondary and tertiary produced mesons in the magnetic fields induced by the magnetic horns and their decays in the decay volume. The simulation uses the proton beam profiles measured by the J-PARC neutrino beam line and is tuned with external hadron production measurements, mainly from the NA61/SHINE experiment na61det ; na612009 . We can select either a muon neutrino beam or a muon anti-neutrino beam by changing the current polarity of the focusing magnetic horns. In this analysis, data collected in the former beam configuration is used. The simulated on-axis neutrino beam has a mean energy of 1.5 GeV and a 1$\sigma$ standard deviation between $-$0.75 GeV and +0.85 GeV, as shown in Fig. 4.

For a given flux of incoming neutrinos, NEUT simulates the neutrino interactions with nuclei, including initial and final state interactions inside the nuclei, in order to provide the four-momenta of all induced particles. In this analysis, the version 5.3.3 of NEUT is used. CC quasi elastic (CCQE)-like, neutral-current (NC) elastic, CC and NC single pion production (1$\pi$), deep inelastic scattering (DIS), multi-pion production, and coherent interactions are simulated. The CCQE-like interactions, characterized by the inclusion of a single charged lepton and no mesons in the final state, are simulated with a relativistic Fermi Gas model (RFG) cite:rfg , random phase approximation (RPA) cite:rpa , and multi-nucleon (2p2h) interactions cite:mecnieves . In addition to the nominal NEUT model, we test the sensitivity of the analysis to determining alternate, available models cite:sf . Table 2 shows the nominal settings for each of the interaction models and tunable parameters in NEUT in this analysis. More details about the underlying neutrino interaction models implemented in NEUT that are used in the analysis can be found in Abe:2017vif . Figure 4 shows the energy of neutrinos that interacted with the target nuclei of the Water Module simulated by NEUT. The main modes of the CC interactions are CCQE111Here, 2p2h interactions are not included., CC1$\pi$, CC multi-pion and DIS production. The fraction of NC interaction is 30% of all interactions. Figure 5 shows the momentum and scattering angle distributions in the laboratory frame for muons produced by $\nu_{\mu}$ CC interactions. In this analysis, due to the limited acceptance of the horizontal INGRID modules to be used for muon identification for the Water Module and Proton Module as described in Sec. 5.5, we define the signal with a restricted phase space of muon kinematics, particularly CC interactions with $\theta_{\mu}<$45*∘* and $p_{\mu}>$0.4 GeV/$c$ in the laboratory frame. The cross section of the signal per nucleon is predicted by NEUT to be slightly different amongst $\rm H_{2}O$, CH, and Fe, as shown in Table 3. This is due to the target dependence of the total cross section of the CC coherent interaction, which is proportional to the square of atomic number, and the difference in the fraction of neutrons and protons per nucleus for the targets considered.

GEANT4 simulates the behavior of the secondary particles induced by the neutrino-nucleus interactions in the detector. The version v9r2p01n00 of GEANT4 and the physics list of QGSP BERT are used for the simulation. The geometry of the three detectors and the walls of the detector hall are modeled in GEANT4 based on the measurements performed during the detector construction. Responses of scintillator, MPPC, and electronics are modeled based on the measurements, as shown in Table 1. The energy deposited in the scintillators estimated by GEANT4 is converted to the observed number of p.e. by multiplying it by a constant determined from measurements with minimum ionization particles (MIP), performed during the detector operation. The following effects are taken into account: the quenching effect of the scintillator; position-dependent light collection efficiency of WLS fibers; attenuation and propagation time of the light in the WLS fiber; crosstalk between grid scintillators; MPPC noise; MPPC crosstalk and after-pulses; MPPC saturation; noise from electronics; gate width of the electronics; and statistical fluctuation of photon counting. For the physics analysis, the neutrino flux and interactions on detector targets, plastic scintillators, and main mechanical structures of the detector and the walls of the detector hall are simulated for the three detectors. Backgrounds from cosmic rays are negligible, as described in Sec. 7.3, and are not simulated for the physics analysis.

4 Data samples

In this article, the data samples recorded by both the INGRID and Proton Module were taken from November 2010 to May 2013. The total number of protons on target (POT) is $5.89\times 10^{20}$ with the neutrino-mode beam. In July 2016, after the Water Module construction and its commissioning were completed, the Water Module replaced the Proton Module for physics data taking. A total POT of 7.25$\times 10^{20}$ POT were collected with the neutrino-mode beam by the Water Module and INGRID during a period between October 2016 and April 2017.

5 Event selections

In this analysis, we define the signal with a restricted phase space of muon kinematics, particularly CC interactions with $\theta_{\mu}<$45*∘* and $p_{\mu}>$0.4 GeV/$c$. The main signature of the CC interactions is the presence of a muon-like track produced inside the detector. Neutrino interactions originating from outside the detectors, CC interactions with non-target materials inside the detectors (mainly scintillators for the studied case with the Water Module), $\overline{\nu}_{\mu},\nu_{e},\overline{\nu}_{e}$ CC interactions, and NC interactions are the main sources of backgrounds in this analysis. The backgrounds from the NC interactions do not produce muons. In order to identify the muons originating from the Water Module and Proton Module, events on the Water Module or Proton Module are required to have a track which penetrates at least two iron planes in one of the three horizontal INGRID modules near the beam center. This method for muon identification limits the phase space of the induced muon, because we reject the CC interactions with low momentum muons, which do not penetrate the iron planes, and high angle muons, which do not enter the three INGRID modules. The event selections applied to the three detectors are similar to that from a previous analysis ccincpm , achieving a similar selection performance for the cross section measurements in the three targets. Figure 6 shows an event display of a typical signal event passing the event selection criteria for the Water Module.

5.1 Event selections for the Water Module

5.1.1 Time clustering

Scintillator channels having charges larger than 2.5 p.e. are defined as a $``$hit”. Hits are clustered with the following criteria: if there are more than three hits within 100 nsec in the Water Module, all the hits within 50 nsec from the average time are grouped into a single cluster.

5.1.2 Two-dimensional track reconstruction

The two-dimensional tracks in the x-z and y-z views are reconstructed independently by using a cellular automaton algorithm scibar to cluster the hits. More details about the algorithm can be found in the reference. The hits in the neighbor scintillator planes are defined as a $``$cell”. Based on $\chi^{2}$ values given by the linear fitting of the relevant hits, it is judged if the pair of two cells having a common hit are merged into a new cell. This is repeated for all cells until no new cell is found and the long cells which have more than three hits are defined as tracks.

5.1.3 Two dimensional track matching with the horizontal INGRID modules

When two-dimensional tracks are reconstructed in the same beam bunch for both the Water Module and the three horizontal INGRID modules near the beam center, an attempt is made to match one to the other. The tracks are matched if they meet the following requirements:

• •

The upstream edge of the reconstructed track in the three INGRID modules is in the most upstream two layers of the INGRID modules.

• •

The difference between the reconstructed angle of the three INGRID modules and Water Module tracks with respect to the z-axis must be less than 35*∘*.

• •

At the halfway point between the three INGRID modules and Water Module, the distance between the three INGRID modules and Water Module track is less than 150 mm.

5.1.4 Three-dimensional track matching

Three-dimensional tracks are formed amongst the pairs of two-dimensional INGRID matched tracks in the x-z plane and in the y-z plane imposing that the difference between the two measurements of the z coordinates of the most upstream hits to be less than or equal to one plane of the parallel scintillators. If there are multiple candidates, we select a pair with the smallest difference of the most upstream hit point z. If there are still multiple candidates after the selection, we select a pair with the smallest difference of the most downstream hit point z.

Only events which have at least one INGRID-matched track are used for the analysis. Because the horizontal INGRID modules are located downstream of the Water Module, the angular acceptance is limited. In addition, the momentum acceptance is limited because the track is required to penetrate at least two iron planes of the INGRID modules for the matching.

5.1.5 Vertexing

After the three-dimensional track reconstruction, the most upstream z coordinate of each INGRID matched three-dimensional track is identified as a reconstructed vertex. If a pair of INGRID-matched three-dimensional tracks meet the following conditions they are identified as tracks coming from a common vertex:

• •

The difference between the most upstream z coordinate of the two tracks in the x-z view, added to the same difference in the y-z view, has to be less than three planes of the parallel scintillators.

• •

The distance between the upstream z coordinate of the two tracks in the x-y plane is less than 150 mm.

These cuts are applied to every vertex since each one is expected to correspond to a single neutrino interaction. The vertex position is re-defined as that of the longest INGRID-matched track amongst those that belong to the common vertex. The longest INGRID-matched track is defined as a muon-like track.

5.1.6 Beam timing cut

To reduce non-beam backgrounds, such as cosmic rays, only events within 100 nsec of the expected beam bunch timing are selected, as shown in Fig. 7. The individual event timing is defined as the time recorded by the MPPC channel with the largest light yield.

5.1.7 Upstream veto cut and fiducial cut

Two cuts are applied based on the position of the vertex to reduce beam-induced backgrounds from neutrino interactions outside the Water Module, mainly from the walls of the detector hall and the INGRID vertical modules. If the upstream point of a track is in the first or second plane of the parallel scintillators, then that event is rejected. The fiducial volume is defined as the central part of the Water module with dimensions of 70 cm (in x-coordinate) $\times$70 cm (in y-coordinate) $\times$ 21 cm (in z-coordinate).

The vertex is required to be within the fiducial volume for the neutrino event to be selected. Figure 8 shows distributions of the vertex used for these two cuts.

5.1.8 Reconstructed angle cut

The three-dimensional angle of the longest reconstructed track from a vertex is required to be smaller than 45*∘* to reduce large-angle muons since the detection efficiency for such kind of events is less than 10%, as described in Sec. 5.5.

5.1.9 Event selections summary

Table 4 shows a summary of the parameters used for the event selection. The numbers of selected events and the backgrounds in the Water Module at each selection step are summarized in Table 5. There are $1.73\times 10^{4}$ events expected in the MC after the event selection. The purity of the $\nu_{\mu}$ CC interactions on $\rm H_{2}O$ is 69.0% and the main background is from neutrino interactions on the scintillators (19.8%). The remaining background sources are NC interactions (2.9%) due to misidentification of pions, neutrino interactions of $\overline{\nu}_{\mu}$, $\nu_{e}$ and $\overline{\nu}_{e}$ (2.0%), photons from $\pi^{0}$ produced by neutrino interactions on the walls of the detector hall (2.4%), and backscattered production of neutrino interactions in the INGRID (3.1%). The muon-like tracks, identified as the longest INGRID-matched track, have 87% probability to be the true muons. Figure 9 shows the neutrino energy, muon momentum, and angle distributions of the selected events predicted by MC. The main interaction modes are CCQE, CC1$\pi$, CC multi-pion and DIS production. Figure 10 (upper left) shows the angle distribution of the reconstructed muon-like tracks for events which passed all event selection in the Water Module.

5.2 Event selections for the Proton Module

The event selections for the Proton Module and INGRID module are very similar to those for the Water Module. However, due to the difference in the scintillator layout, a few parameters for the cellular automaton algorithm and event selection have been optimized as listed in Table 4.

The numbers of selected events and the backgrounds in the Proton Module at each selection step are summarized in Table 6. After the event selection, a total of $2.23\times 10^{4}$ events are expected by MC. The purity of the CC interactions on CH is 85.4%. Background sources are NC interactions (4.2%), neutrino interactions of $\overline{\nu}_{\mu}$, $\nu_{e}$ and $\overline{\nu}_{e}$ (2.4%), photons from $\pi^{0}$ produced by neutrino interactions on the walls of the detector hall (2.1%) and backscattered events from neutrino interactions in the INGRID (5.2%). Figure 10 (upper right) shows the angle distribution of the reconstructed muon-like tracks for events which passed all event selection in the Proton Module.

5.3 Event selections for the INGRID module

The event selections are applied for the horizontal INGRID module located at the beam center with the parameters listed in Table 4. In addition, an $``$acceptance cut” is applied only for the INGRID module in order to achieve a similar angular acceptance with the Water Module and Proton Module. An imaginary module located directly behind the INGRID module is defined, as shown in Fig. 11. The distance between the INGRID module and the imaginary module is the same as between the Water Module and the INGRID horizontal modules. The reconstructed tracks are then projected further downstream, even if the track has stopped in the INGRID module. If at least one reconstructed track from the vertex reaches the imaginary module, that event is selected.

The numbers of selected events and the backgrounds in the INGRID module at each selection step are summarized in Table 7. After the event selection, a total of $3.12\times 10^{5}$ events are expected by MC. The purity of the $\nu_{\mu}$ CC interactions on Fe is 88.1%. Background sources are NC interactions (5.2%), neutrino interactions of $\overline{\nu}_{\mu}$, $\nu_{e}$ and $\overline{\nu}_{e}$ (2.9%), neutrino interactions on scintillator (3.3%), photons from $\pi^{0}$ produced by neutrino interactions on walls of the detector hall (0.3%) and the other INGRID modules (0.2%). Figure 10 (lower) shows the angle distribution of the reconstructed muon-like tracks for events which passed all event selection in the INGRID module.

5.4 Pileup correction for the INGRID module

If more than one neutrino event occurs in the detector at the same bunch timing, we sometimes fail to count them. Therefore, a correction must be applied to account for this event pileup effect. For the INGRID module, this effect is estimated in each bin of the reconstructed track angle by merging multiple bunches to enrich the pileup rate artificially. Table 8 shows the number of selected events before and after the pileup correction. For the Water Module and Proton Module, the effect of the pileup is small due to the small target mass therefore no correction is applied.

5.5 Selection efficiencies

Figure 12 shows the selection efficiency of CC interactions for the Water Module, Proton Module, and the INGRID module as a function of true muon scattering angle and momentum. Because the selection efficiencies for the CC interactions with $\theta_{\mu}>$45*∘* or $p_{\mu}<$400 MeV, are less than 10%, these events are excluded from the signal sample defined in this analysis. Figure 13 shows the efficiency of the signal for the three detectors and their ratios as a function of the muon scattering angle. The signal efficiency is almost constant as a function of muon momentum, while it depends on the muon scattering angle. In this analysis, the cross section is calculated by a sum of the differential cross sections as a function of the muon scattering angle, as described in Sec. 6. In this method, the efficiency is calculated for each bin of the scattering angle and the dependence of the signal efficiency on the MC models used in this analysis is reduced.

6 Cross section analysis

6.1 Analysis method

The flux-integrated $\nu_{\mu}$ cross sections of CC interactions on water ($\sigma_{\rm{H_{2}O}}$), hydrocarbon ($\sigma_{\rm{CH}}$), and iron ($\sigma_{\rm{Fe}}$) defined in a restricted phase space of the induced muon, $\theta_{\mu}<45^{\circ}$ and $p_{\mu}>$0.4 GeV/$c$, are measured as a sum of the differential cross sections as a function of the muon scattering angle. They are calculated as follows:

[TABLE]

where ${A}$ represents the type of the target material ($\rm{H_{2}O}$, CH, and Fe) and $D$ is the corresponding detector (Water Module, Proton Module, and INGRID). $N^{\rm{sel}}$ is the number of selected events, $N^{\rm{BG}}$ is the number of expected backgrounds, $\Phi$ is the integrated $\nu_{\mu}$ flux, $T$ is the number of target nucleons, and $\varepsilon$ is the detection efficiency of the signal. Subscript $i$ is a bin index of the true muon scattering angle and subscript $j$ is a bin index of the reconstructed angle of the muon-like track. The true and reconstructed muon scattering angle bins are defined as 9 bins from 0*∘* to 45*∘* with a bin width of 5*∘*, which are optimized based on the detector resolution. $U_{ij}$ is a probability that events in the reconstructed angle bin $j$ are in the true muon scattering angle bin $i$. The CC cross section ratios are estimated by taking the ratios of the $\sigma_{\rm{H_{2}O}}$, $\sigma_{\rm{CH}}$ and $\sigma_{\rm{Fe}}$.

The $N^{\rm sel}$ is estimated based on data as shown in Fig. 10 for the Water Module and Proton Module, and Table 8 for the INGRID module with the pileup correction. Except for the $\sigma_{\rm H_{2}{O}}$ measurement with the Water Module, in which the backgrounds from CC interactions on plastic scintillator ($N^{\rm{CH~{}BG}}_{\rm WM}$) are estimated with data from the Proton Module, other backgrounds $N^{\rm BG}$ are estimated by MC simulation. The $N^{\rm{CH~{}BG}}_{\rm WM}$ is estimated as follows:

[TABLE]

where $\sigma_{i~{}\rm{CH}}$ is the differential cross section on the CH target with $i$-th muon scattering angle bin. The other backgrounds are estimated by MC as summarized in Table 9 in detail. The integrated $\nu_{\mu}$ fluxes $\Phi$ are estimated to be $\Phi^{\rm H_{2}O}_{\rm WM}=3.72\times 10^{13}$ /cm2 with 7.25$\times$1021 POT, $\Phi^{\rm CH}_{\rm PM}=3.02\times 10^{13}$ /cm2 with 5.89$\times$1021 POT, and $\Phi^{\rm Fe}_{\rm INGRID}=2.99\times 10^{13}$ /cm2 with 5.89$\times$1021 POT by MC, as shown in Table 10. Although data samples used for the Proton Module and INGRID module are at the same delivered POT, the fact that the Proton Module is 1.2 m closer to the production target than the INGRID module, leads to a small difference in the integrated flux between them. The number of target nucleons, $T$, is calculated based on measurements performed during the detector construction as shown in Table 11. The detection efficiency of the signal, $\varepsilon$, is estimated by MC as shown in Fig. 13 in each true muon scattering angle bin.

The $U_{ij}$, probability that events in the reconstructed angle bin $j$ are in the true muon scattering angle bin $i$, is calculated as follows based on Bayes’s theorem:

[TABLE]

where $P(\theta^{\rm{recon}}_{j}|\theta^{\rm{true}}_{i})$ is calculated by MC as shown in Fig. 14. The $P$($\theta^{\rm{true}}_{i}$) is calculated by an iterative unfolding method DAgostini:1994fjx , which is briefly described as follows:

$P(\theta^{\rm{true}}_{i})$ is set to a flat prior, 2. 2.

calculate $U_{ij}$, 3. 3.

$P(\theta^{\rm{true}}_{i})$ is set to $\sum_{j}U_{ij}(N^{\rm{sel}}_{j}-N^{\rm{BG}}_{j})/\sum_{ij}U_{ij}(N^{\rm{sel}}_{j}-N^{\rm{BG}}_{j})$, 4. 4.

repeat (2)–(3).

The number of required iterations is set to 10 as described in Sec. 6.2.

6.2 Consistency test

From the number of selected events and the quantities described earlier in this section, the flux-integrated CC cross sections on $\rm H_{2}O$, CH, Fe, and their ratios are calculated based on Eq. 1. In this section, a consistency test is performed by replacing the number of selected events of data with that of the MC expectation, in order to check the consistency between the calculated cross section and MC expectation. Figure 15 shows the relation between the number of iterations and deviations of the calculated cross sections from MC expectation and the number of iterations when it is set to 10. Table 12 shows the calculated cross sections and their consistency with the MC expectation. The consistency test is performed with not only the nominal cross section model but also a few alternative models.

7 Systematic uncertainties

There are three main sources of systematic uncertainties for the cross section measurements: neutrino flux, neutrino interaction models, and detector response. The uncertainty evaluation for each source is detailed in this section.

7.1 Systematic uncertainties from the neutrino flux

The T2K neutrino flux simulation, based on JNUBEAM mentioned in Sec. 3, relies on several measurements as inputs, including the hadron production measurements and information from the J-PARC beam line monitors. The uncertainty on the flux prediction takes into account the uncertainties in the measurements of the hadron scattering experiments, the hadronic interaction models and the uncertainties in the beam profile measurements with the beam line monitors. Details about the sources of the flux uncertainty can be found in jnubeam . Figure 16 shows the calculated total on-axis flux uncertainty as a function of neutrino energy.

The uncertainty of the neutrino flux is related to the systematic uncertainties on the number of expected backgrounds ($N^{\rm{BG}}$), integrated flux ($\Phi$), detection efficiency ($\varepsilon$), and reconstructed-to-true transfer matrix ($U$). To evaluate the systematic effects on the cross section measurement, the number of produced and selected neutrino events in each bin of the reconstructed track angle and true muon scattering angle is varied by using the calculated flux uncertainty, including correlations between the true neutrino energy bins. Therefore the variations of $N^{\rm{BG}}$, $\Phi$, $\varepsilon$, and $U$ are calculated and the variation of the cross section result is determined. This is repeated for many toy data sets and the 68% range of the distribution of the cross section variation around the central value is taken as the size of the flux-related systematic uncertainty. The first row in Table 13 shows the calculated flux uncertainties. They are approximately 10% for the absolute cross section measurement and 1–2% for the cross section ratios.

In addition, uncertainties due to difference in position of the INGRID module compared with the Water Module and Proton Module, and difference of the running periods between the Water Module, Proton Module, and INGRID module are estimated separately. The former is estimated to be 0.31% based on measurement of the detector location. The latter is estimated to be 1.03% based on the beam stability measurements of the INGRID module between the different running periods. Their quadratic sums are summarized in the second row of Table 13.

7.2 Systematic uncertainties from the neutrino interaction models

The NEUT neutrino interaction model has a number of uncertainties that can affect the detection efficiency ($\varepsilon$), background contamination ($N^{\rm{BG}}$), and reconstructed-to-true transfer matrix ($U$). To evaluate the model-related effect on the cross section measurement, for each $\pm$ 1$\sigma$ variation of a given interaction model parameter, a deviation of the cross section from the nominal value calculated based on the induced variation of $\varepsilon$, $N^{\rm{BG}}$, and $U$ are set as a systematic uncertainty. Table 14 shows the nominal values and the uncertainties of the neutrino interaction parameters. More details about the simulation models used can be found in Abe:2017vif . In addition, uncertainties from pion final state interactions inside nuclei are taken into account: for each type of interaction, the uncertainties are assigned as normalization, as shown in Table 14.

When the uncertainty is calculated, no correlation amongst the different target nuclei for the Fermi momentum ($p_{F}$), binding energy ($E_{b}$), 2p2h, and CC coherent normalizations is assumed. Full correlation amongst the different nuclei is assumed for the other parameters. Table 15 shows the calculated uncertainties and they are in the range between 2.6% and 5.2%. The dominant ones come from the uncertainties of the axial vector mass of the CCQE, CC1$\pi$, and the energy-dependent normalization of the CC multi-pion and DIS production. The uncertainty of the beam-induced backgrounds coming from outside of the detector is not included here, although it affects the $N^{BG}$. It is calculated as one of the detector systematics, as described in Sec. 7.3.

In addition to the systematic effects estimated by NEUT, the uncertainties of backscattered protons and pions produced by neutrino interactions with nuclei, which mainly affects the position of the reconstructed vertex, are estimated independently. A fraction of the events generated inside the fiducial volume have reconstructed vertices outside the fiducial volume due to backscattered secondary protons or pions. The fraction of such events is 3.0% for the Water Module, 1.6% for the Proton Module, and 2.0% for the INGRID module with respect to the total number of selected events. The number and the uncertainty of such backscattered secondary particles may not be simulated well by NEUT, so a 50% conservative uncertainty is assumed, which leads to 1.5%, 0.8%, and 1.0% uncertainties for the Water Module, Proton Module, and INGRID module respectively in the total number of selected events. This is taken to be the 1$\sigma$ uncertainty for all reconstructed angle bins. In addition, no correlations between the target materials are assumed for this error.

7.3 Systematic uncertainties from the detector responses

Uncertainties of the detector response are estimated based on the difference between data and MC for the cosmic rays and beam-induced muons coming from outside of the detectors. We take into account the following errors: target mass, MPPC noise, scintillator crosstalk, hit efficiency of the scintillator, event pileup, beam-induced backgrounds from outside of the detector, two dimensional tracking efficiency, and three-dimensional tracking efficiency. In addition, the uncertainties of the reconstructed variables used for the event selections are taken into account as follows: two-dimensional track matching with the INGRID modules, three-dimensional track matching, vertexing, beam timing cut, veto and fiducial volume cut, and reconstructed angle cut. The effect from non beam-induced backgrounds is estimated to be less than 0.1% with beam-off data and is not included in the systematic uncertainties.

In order to evaluate these uncertainties on the cross section measurement, MC simulations are produced by varying detector parameters independently within their uncertainties by 1$\sigma$. The difference in the number of selected events in each bin of reconstructed track angle with respect to varying their uncertainty by 1$\sigma$ defines the 1$\sigma$ standard deviation systematic uncertainty in the number of events. Table 16 shows a summary of the uncertainties from the detector response for the absolute cross sections. For the measurements of the cross section ratios, no correlation is assumed between the three detectors except for the beam-induced backgrounds from the outside the detector, which is treated as a common uncertainty. The fourth row in Table 13 shows the total uncertainty from the detector response. They are approximately 2% for the absolute cross section measurement and 4% for the cross section ratios because most of the systematics do not cancel between the detectors.

The total systematic uncertainties of the cross section measurements are estimated as a quadratic sum of the uncertainties of the neutrino flux, neutrino interaction, and detector response. Table 13 shows the total systematic uncertainties and they are between 10% and 14% for the absolute cross section measurements and approximately 5% for the cross section ratios.

8 Results

The measured flux-integrated cross sections of $\nu_{\mu}$ CC interactions per nucleon at a mean neutrino energy of 1.5 GeV, defined in a restricted phase space of induced muon, $\theta_{\mu}<45^{\circ}$ and $p_{\mu}>$0.4 GeV/$c$, on $\rm H_{2}O$, CH, and Fe are

[TABLE]

[TABLE]

[TABLE]

The cross section ratios are

[TABLE]

[TABLE]

[TABLE]

The errors of both the measured absolute cross section and cross section ratios are dominated by the systematic uncertainties. This is the most precise measurement to date of neutrino cross sections on water in this energy region and the first measurement of neutrino cross section ratios of water-to-hydrocarbon and water-to-iron. Figure 17 shows the measured cross sections and their predictions by NEUT with nominal and varied parameters of the axial vector mass $M_{A}^{QE}$, normalization of 2p2h interaction, and Fermi momentum with 1$\sigma$, which have a relatively large effect on the cross section ratios in the parameters listed in Table 14. This is due to the fact that the variations of 2p2h normalization and Fermi momentum are only applied for $\rm H_{2}O$ but not CH and Fe as a conservative way to deal with our poor understanding of the target dependency of the neutrino interaction. The predictions agree with data within the uncertainties. This measurement validates the neutrino interaction models on the water-target and difference between water, plastic, and iron and confirms the reliability of the T2K oscillation analysis. Additional comparison of data and predictions with other parameters of the neutrino interaction listed in Table 14 are summarized in Fig. 18 and Table 17. All of the predictions agree with data within the estimated uncertainties. These results of the measurements and the neutrino flux at on-axis location are provided in text and ROOT format under the following link cite:link .

9 Conclusion

For the precise measurement of neutrino oscillation parameters, understanding of neutrino interactions with nuclei is essential. We reported measurements of the flux-integrated $\nu_{\mu}$ charged-current cross sections on water, hydrocarbon, iron, and their ratios in the T2K on-axis neutrino beam with a mean neutrino energy of 1.5 GeV in a restricted phase space for the kinematics of the induced muon with $\theta_{\mu}<$45*∘* and $p_{\mu}>$0.4 GeV/$c$ in the laboratory frame. This is the most precise measurement to date of neutrino cross sections on water in this energy region and the first measurement of neutrino cross section ratios of water-to-hydrocarbon and water-to-iron. The results agree with current neutrino interaction models used in the T2K oscillation analysis within their uncertainties.