T2K neutrino flux prediction

The Tokai-to-Kamioka (T2K) experiment studies neutrino oscillations using an off-axis muon neutrino beam with a peak energy of about 0.6 GeV that originates at the J-PARC accelerator facility. Interactions of the neutrinos are observed at near detectors placed at 280 m from the production target and at the far detector -- Super-Kamiokande (SK) -- located 295 km away. The flux prediction is an essential part of the successful prediction of neutrino interaction rates at the T2K detectors and is an important input to T2K neutrino oscillation and cross section measurements. A FLUKA and GEANT3 based simulation models the physical processes involved in the neutrino production, from the interaction of primary beam protons in the T2K target, to the decay of hadrons and muons that produce neutrinos. The simulation uses proton beam monitor measurements as inputs. The modeling of hadronic interactions is re-weighted using thin target hadron production data, including recent charged pion and kaon measurements from the NA61/SHINE experiment. For the first T2K analyses the uncertainties on the flux prediction are evaluated to be below 15% near the flux peak. The uncertainty on the ratio of the flux predictions at the far and near detectors is less than 2% near the flux peak

Evidence of electron neutrino appearance in a muon neutrino beam

The T2K Collaboration reports evidence for electron neutrino appearance at the atmospheric mass splitting, |Δm322|≈2.4×10-3 eV2. An excess of electron neutrino interactions over background is observed from a muon neutrino beam with a peak energy of 0.6 GeV at the Super-Kamiokande (SK) detector 295 km from the beam’s origin. Signal and background predictions are constrained by data from near detectors located 280 m from the neutrino production target. We observe 11 electron neutrino candidate events at the SK detector when a background of 3.3±0.4(syst) events is expected. The background-only hypothesis is rejected with a p value of 0.0009 (3.1σ), and a fit assuming νμ→νe oscillations with sin⁡22θ23=1, δCP=0 and |Δm322|=2.4×10-3 eV2 yields sin⁡22θ13=0.088-0.039+0.049(stat+syst)