 MINOS/MINOS+ STERILE 3+1 NEUTRINO DISAPPEARANCE RESULT: TWODETECTOR ANALYSIS (February 2019)

MINOS/MINOS+ searches for sterile neutrinos using chargedcurrent and neutralcurrent events in 10.56×10^{20} POT ν_{μ} running data collected by MINOS, and 5.80×10^{20} POT ν_{μ} running data collected by MINOS+. The previous analyses were based on the ratio between the measured neutrino energy spectra in the two detectors (FaroverNear ratio), whereas this analysis employs a twodetector fit method, directly fitting the reconstructed neutrino energy spectra in the two detectors to significantly improve the sterile neutrino sensitivity for Δm^{2}_{41} > 5×10^{3} eV^{2}. This plot shows the combined MINOS and MINOS+ FeldmanCousins 90% C.L. data contour using ν_{μ} disapearance and searching for a deficit in NC events. The combined MINOS/MINOS+ coverage is compared to the published MINOS FaroverNear data coverage (green; see also below) and coverages of other experiments. Note the axis is sin^{2}θ_{24}. All experiments bar SuperK and MINOS(+) used a twoflavour approach and therefore have had their contours made symmetric when converting from sin^{2}(2θ_{24}) to sin^{2}θ_{24}.
The PRL paper describing the details of this analysis can be found here: Phys. Rev. Lett. 122, 091803.
The arXiv paper can be found here: arXiv:1710.06488 [hepex].
Under Ancillary files for the arXiv paper, you can find a data release root file and corresponding example macro as well as a document containing supplemental info on this analysis. The root file contains the MINOS/MINOS+ data and Monte Carlo histograms, as well as the systematic uncertainty matrices (it can also be downloaded here). The content of this root file is listed in this pdf. The example macro for a MINOS/MINOS+ fit (which can also be downloaded here) is executed as follows: root l b q dataRelease_chi2Calc_compile.C+, which prints the parameters and corresponding χ^{2} output and then saves the results in a root file.
Published Results
 a listing of all journal publications by the MINOS collaboration
 a listing of theses written by MINOS graduate students, available as a web page with links or a text list
The Public Materials Page provides downloadable PDF files of plots and diagrams from recent results, formatted as presentation slides. These files are released by MINOS for public use by nonMINOS scientists. All have been presented at recent conferences and represent published or soontobe published results.
Highlights from Recent Results
Plots and brief explanatory text for recent MINOS and MINOS+ results. Click on the plot image for a fullpage version.
 MINOS/MINOS+ LARGE EXTRA DIMENSIONS RESULT: TWODETECTOR ANALYSIS (May 2018)

MINOS/MINOS+ searches for large extra dimensions using chargedcurrent and neutralcurrent events in 10.56×10^{20} POT ν_{μ} running data collected by MINOS, and 5.80×10^{20} POT ν_{μ} running data collected by MINOS+. The previous analyses were based on the ratio between the measured neutrino energy spectra in the two detectors (FaroverNear ratio), whereas this analysis employs a twodetector fit method, directly fitting the reconstructed neutrino energy spectra in the two detectors to significantly improve the sensitivity to large extra dimensions. This plot shows the combined MINOS and MINOS+ FeldmanCousins 90% C.L. data contour using ν_{μ} disapearance and searching for a deficit in NC events. No evidence for large extra dimensions is found. In the limit of a vanishing lightest neutrino mass m_{0}, the large extra dimension size R is constrained to be smaller than 0.30 μm at 90% C.L. by the MINOS and MINOS+ data sets.
 MINOS NONSTANDARD INTERACTIONS ANALYSIS RESULT (January 2017)

The 90% C.L. in the ε_{eτ},(δ_{CP}+δ_{eτ}) parameter space for normal (top) and inverted (bottom) neutrino mass hierarchy using ν_{e} and ν̅_{e} appearance candidates in the MINOS Far Detector. The shaded areas to the left of the solid contours indicate the MINOS allowed regions where additional oscillation parameters, including NSI, were included in the fit. The dotted contours show the limits where the additional oscillation parameters were fixed.
The PRL paper describing the details of this analysis can be found here: Phys. Rev. D 95, 012005.
The arXiv paper can be found here: arXiv:1605.06169 [hepex].
 MINOS LARGE EXTRA DIMENSIONS RESULT: FAROVERNEAR ANALYSIS (December 2016)

MINOS searches for large extra dimensions using chargedcurrent and neutralcurrent events in 10.56×10^{20} POT ν_{μ} running data collected by MINOS. This plot shows the FeldmanCousins 90% C.L. data contour for MINOS obtained from a FaroverNear ratio analysis. No evidence for large extra dimensions is found. In the limit of a vanishing lightest neutrino mass m_{0}, the large extra dimension size R is constrained to be smaller than 0.45 μm at 90% C.L. by the MINOS data set.
The PRD paper describing the details of this analysis can be found here: Phys. Rev. D. 94, 111101(R).
The arXiv paper can be found here: arXiv:1608.06964 [hepex].
 MINOS/DAYA BAY/BUGEY3 STERILE 3+1 RESULT (October 2016)

In a joint analysis with the Daya Bay collaboration, the MINOS sterile 3+1 FaroverNear result (see also below) is combined with the Daya Bay and Bugey3 sterile 3+1 results to allow a direct comparison of with the LSND and MiniBooNE results. This plot shows the MINOS and Daya Bay/Bugey3 combined 99% CL_{S} limit on sin^{2}2θ_{μe} compared to the LSND and MiniBooNE antineutrino 99% C.L. allowed regions. Regions of parameter space to the right of the red contour are excluded. The MiniBooNE neutrino 99% exclusion region is also shown. We note that the excursion to small mixing in the exclusion contour at around Δm^{2}_{41}∼5×10^{3} eV^{2} is due to the degenerate island in the MINOS exclusion limit. Also shown are the results of two global fits. The Kopp et al. 2013 fit included ν_{e} appearance experiments, and the Gariazzo et al. 2016 fit included all shortbaseline experiments with the exception of MiniBooNE data below 475 MeV.
The PRL paper describing the details of this joint analysis can be found here: Phys. Rev. Lett. 117, 151801.
The arXiv paper can be found here: arXiv:1607.01177 [hepex].
The results are also shown on the Daya Bay Collaboration website.
The root file containing the CL_{S} contours from MINOS, Daya Bay/Bugey3, and the combination of MINOS and Daya Bay/Bugey3 can be downloaded here. The root file containing the MINOS data and Monte Carlo histograms, as well as the systematic uncertainty matrices, can be downloaded here (the content of this root file is listed in this pdf and an example macro to run a MINOS fit can be downloaded here; root l b q dataRelease_chi2Calc_compile.C+ prints the parameters and corresponding χ^{2} output and then saves the results in a root file). See also below for more information on the MINOS sterile analysis result.
 MINOS/MINOS+ STERILE 3+1 NEUTRINO DISAPPEARANCE RESULT: FAROVERNEAR ANALYSIS (October 2016)
 MINOS/MINOS+ searches for sterile neutrinos using chargedcurrent and neutralcurrent events in 10.56×10^{20} POT ν_{μ} running data collected by MINOS, and 5.80×10^{20} POT ν_{μ} running data collected by MINOS+. This plot shows the combined MINOS and MINOS+ FeldmanCousins 90 and 95% C.L. contours for data using ν_{μ} disapearance and searching for a deficit in NC events. The combined MINOS and MINOS+ contours are compared to other experiments. Note the axis is sin^{2}θ_{24}. All experiments bar SuperK and MINOS(+) used a twoflavour approach and therefore have had their contours made symmetric when converting from sin^{2}(2θ_{24}) to sin^{2}θ_{24}.
 MINOS STERILE 3+1 NEUTRINO DISAPPEARANCE RESULT: FAROVERNEAR ANALYSIS (October 2016)

This plot shows the MINOS FeldmanCousins 90% C.L. for data using ν_{μ} disapearance and searching for a deficit in NC events. This contour is compared to other experiments. Note the axis is sin^{2}θ_{24}. All experiments bar SuperK and MINOS reported results in terms of sin^{2}2θ_{24}, so their contours were made symmetric when converting to sin^{2}θ_{24}. Also shown are the 90% C.L. allowed regions of two global fits. The Kopp et al. 2013 fit included ν_{e} appearance experiments, and the Gariazzo et al. 2016 fit included all shortbaseline experiments with the exception of MiniBooNE data below 475 MeV. To compare these global fits to disappearance data, we convert sin^{2}2θ_{μe} to sin^{2}θ_{24} by assuming θ_{14} = 0.15, the best fit value from the Kopp et al. global fit to appearance data.
The PRL paper describing the details of this analysis can be found here: Phys. Rev. Lett. 117, 151803
The arXiv paper can be found here: arXiv:1607.01176 [hepex]
The root file containing the MINOS data and Monte Carlo histograms, as well as the systematic uncertainty matrices, can be downloaded here (the content of this root file is listed in this pdf and an example macro to run a MINOS fit can be downloaded here; root l b q dataRelease_chi2Calc_compile.C+ prints the parameters and corresponding χ^{2} output and then saves the results in a root file).
 MINOS STERILE 3+1 ANTINEUTRINO DISAPPEARANCE RESULT: FAROVERNEAR ANALYSIS (October 2016)
 This plot shows the MINOS FeldmanCousins 90% C.L. for data using ν̅_{μ} disapearance. Antineutrino events from both ν̅_{μ} and ν_{μ} running data were used. This contour is compared to other experiments.
 MINOS/MINOS+ STERILE 3+1 ν_{e} APPEARANCE RESULT (October 2016)
 The 90% C.L. using ν_{e} and ν̅_{e} appearance candidates is compared to the MINOS and Daya Bay/Bugey3 combined 90% CL_{S} limit on sin^{2}2θ_{μe}, and the LSND and MiniBooNE antineutrino 90% C.L. allowed regions. The MiniBooNE neutrino 90% exclusion region is also shown. Regions of the parameter space to the right of the dashed brown contour are excluded by the appearance analysis.
 MINOS/MINOS+ STANDARD THREE FLAVOR ANALYSIS RESULT (October 2016)
 The 68% and 90% confidence limits resulting from a fit to 48.67 ktyrs of atmospheric data combined with disappearance and appearance data from the MINOS beam and disappearance data from the MINOS+ beam are shown as solid lines. The dashed orange line is the NOvA Neutrino 2016 result. The dashed pink line is the T2K Neutrino 2016 result.
 MINOS LARGE EXTRA DIMENSION ANALYSIS RESULT (May 2015)
 MINOS search for large extra dimensions using chargedcurrent and neutralcurrent events in 10.56x10^{20} POT ν_{μ} running data. This figure shows the 90% C.L. excluded region in the parameter space m_{0} vs. a, where m_{0} is the smallest neutrino mass and a is the large extra dimension size. The shaded area indicates that the excluded region is to the right of the plot.
 FIRST MINOS+ RESULTS  FAR DETECTOR CHARGED CURRENT SPECTRUM (June 2014)
 The MINOS+ experiment has operated since September 2013 using the upgraded NuMI accelerator neutrino beam, which is now running in a mediumenergy configuration. This figure shows the first spectrum of ν_{μ} and ν_{μ} chargedcurrent interactions in the MINOS+ Far Detector, summed with the existing lowenergy data from MINOS. The hatched histograms show the predicted spectra for MINOS and MINOS+, calculated using the bestfit oscillation parameters from MINOS. The two predictions are summed to give a combined spectrum, which is shown by the blue histogram along with its 1σ systematic uncertainty band. In addition, the red histogram shows the combined spectrum for the case of no oscillations. The observed data, indicated by the points with errors, are welldescribed by the oscillation model. The MINOS+ data significantly increase the event yield in the mediumenergy region, enabling precision measurements of the ν_{μ} survival probability curve.
 FIRST MINOS+ RESULTS  ν_{μ} SURVIVAL PROBABILITY CURVE (June 2014)
 This figure shows the ratio of the observed spectrum of accelerator ν_{μ} CC and ν_{μ} CC interations from MINOS and MINOS+, calculated with respect to the predicted spectrum without oscillations. The blue curve shows the corresponding oscillated prediction and its 1σ systematic uncertainty band, calculated using the bestfit oscillation parameters from MINOS. The observed data are welldescribed by the oscillation model.
 UPDATED OSCILLATION RESULTS  COMBINED ANALYSIS OF ν_{μ} DISAPPEARANCE AND ν_{e} APPEARANCE (June 2014)
 This figure shows the latest oscillation results from MINOS, which have been updated for the Neutrino 2014 conference. The results are based on a combined analysis of ν_{μ} disappearance and ν_{e} appearance, and use both the accelerator and atmospheric data from MINOS. The atmospheric data set has been increased by 30% with respect to the previous results from MINOS. The red contours in this figure show the 2D confidence limits on the parameters Δm^{2}_{32} and sin^{2}θ_{23} obtained for the normal (top panel) and inverted (bottom panel) hierarchy. The singleparameter 68% confidence intervals are: Δm^{2}_{32}= (2.34^{+0.09}_{0.09} )x10^{3} eV^{2} and sin^{2}θ_{23}= 0.43^{+0.16}_{0.04} for the normal hierarchy; and Δm^{2}_{32}= (2.37^{+0.11}_{0.07} )x10^{3} eV^{2} and sin^{2}θ_{23}= 0.43^{+0.19}_{0.05} for the inverted hierarchy. For comparison, the blue contours show the published ν_{μ} disappearance results from T2K (PRL 112, 181801, 2014). The 2D contours from MINOS are available here as a ROOT file.
 STERILE ANALYSIS  COMBINATION WITH BUGEY (June 2014)
 Unitarity constraints generate relationships between ν_{μ}↔ν_{e} transitions and ν_{e} and ν_{μ} disappearance. In this figure, the ν_{μ} disappearance results from MINOS are combined with ν_{e} disappearance results from the Bugey reactor experiment to yield 90% C.L. limits on the sterile mixing parameter sin^{2}2θ_{μe} = 4U_{e4}^{2}U_{μ4}^{2}, relevant to ν_{μ}↔ν_{e} transition experiments also shown. Regions of parameter space to the right of the red contour are excluded at 90% C.L. The MINOS data correspond to a 10.56x10^{20} POT exposure in neutrino running mode. The Bugey results are obtained from a GLoBES 2012 fit provided by P. Huber. It accounts for the new calculation of reactor fluxes, as described in P. Huber, Phys. Rev. C 85 029901 (2011). The MiniBooNE contours are provided by the MiniBooNE Collaboration and were published in A. A. AguilarArevalo et al., Phys. Rev. Lett. 110, 161801 (2013).
 STERILE ANALYSIS  90% C.L. CONTOUR (February 2015)
 MINOS searches for light sterile neutrinos that mix with the active neutrino flavors using chargedcurrent and neutralcurrent neutrino interactions. No evidence for sterile neutrinos is found in the MINOS data. This figure shows the 90% C.L. excluded regions on the parameter space Δm^{2}_{43} vs. sin^{2}(2θ_{24}) obtained from the MINOS analysis compared to CDHS, CCFR, SciBooNE and MiniBooNE experiments. MINOS takes advantage of its long baseline and broad energy spectrum to probe previously unexplored regions of this parameter space.
 STERILE ANALYSIS  SPECTRUM OF CHARGED CURRENT EVENTS (February 2015)
 This figure shows the far detector prediction of chargedcurrent reconstructed neutrino energy spectrum including backgrounds and the total systematic uncertainty and compares it to the observed MINOS data. The prediction assumes a threeflavor oscillation scenario with parameters Δm^{2}_{21}=7.59x10^{5} eV^{2}, Δm^{2}_{32}=2.41x10^{3} eV^{2}, sin^{2}θ_{13}=0.024, sin^{2}θ_{12}=0.319, and sin^{2}θ_{23}=0.388. No evidence of deviations from the standard threeflavor oscillation picture is found. MINOS uses these results to set limits on the existence of light sterile neutrinos.
 STERILE ANALYSIS  SPECTRUM OF NEUTRAL CURRENT EVENTS (February 2015)
 This figure shows the far detector prediction of neutralcurrent reconstructed neutrino energy spectrum including backgrounds and the total systematic uncertainty and compares it to the observed MINOS data. The prediction assumes a threeflavor oscillation scenario with parameters Δm^{2}_{21}=7.59x10^{5} eV^{2}, Δm^{2}_{32}=2.41x10^{3} eV^{2}, sin^{2}θ_{13}=0.024, sin^{2}θ_{12}=0.319, and sin^{2}θ_{23}=0.388. No evidence of deviations from the standard threeflavor oscillation picture is found. MINOS uses these results to set limits on the existence of light sterile neutrinos.
 NONSTANDARD INTERACTION (NSI) LIMITS  ν_{e} APPEARANCE (June 2014)
 The appearance of ν_{e} events at the Far Detector is sensitive to the NSI parameters ε_{eτ} and δ_{eτ}. These plots show the 90% C.L. allowed values of ε_{eτ} and δ_{eτ}+δ_{CP} assuming a normal hierarchy (top plot) and an inverted hierarchy (bottom plot). The appearance rate also depends on standard oscillation parameters which are assumed to be the following: sin^{2}θ_{23}=0.5, sin^{2}θ_{13}=0.025, and Δ m^{2}_{32}=2.43×10^{3} eV^{2}. The contour is produced by marginalizing over the values of δ_{CP}.
Previous Results reported in 2013
 COMBINED ANALYSIS OF ν_{μ} DISAPPEARANCE AND ν_{e} APPEARANCE USING BEAM AND ATMOSPHERIC DATA (August 2013)
 MINOS has performed a combined analysis of ν_{μ} disappearance and ν_{e} appearance using its complete set of accelerator and atmospheric neutrino data. The analysis uses a full threeflavour framework, with the oscillation parameters Δm^{2}_{32}, θ_{23}, θ_{13} and δ_{CP} varied in the fit. The analysis also incorporates an external constraint of sin^{2}θ_{13}=0.0242±0.0025, based on an average of reactor neutrino measurements. The left panels of this figure show the 2D confidence levels on Δm^{2}_{32} and θ_{23} for the normal and inverted hierarchy. The contours are calculated relative to the overall best fit point of Δm^{2}_{32}=2.41x10^{3} eV^{2} and sin^{2}θ_{23}=0.41, which is indicated by the star. The right panels show the corresponding 1D loglikelihood profiles as a function of Δm^{2}_{32} and θ_{23}. The singleparameter confidence limits are: Δm^{2}_{32}=[2.282.46]x10^{3} eV^{2} (68% C.L.) and sin^{2}θ_{23}=[0.350.65] (90% C.L.) for the normal hierarchy; and Δm^{2}_{32}=[2.322.53]x10^{3} eV^{2} (68% C.L.) and sin^{2}θ_{23}=[0.340.67] (90% C.L.) for the inverted hierarchy. The 2D contours are available here as a ROOT file. The corresponding 2D likelihood surfaces are here. More details on the analysis can be found in arXiv:hepex/1403.0867  published in Phys. Rev. Lett. 112, 191801 (2014).
 COMBINED ANALYSIS OF ν_{μ} DISAPPEARANCE AND ν_{e} APPEARANCE USING BEAM AND ATMOSPHERIC DATA (August 2013)
 This figure shows the 1D likelihood profile for the δ_{CP} parameter resulting from the combined analysis of ν_{μ} disappearance and ν_{e} appearance by MINOS using its complete set of accelerator and atmospheric neutrino data. Separate profiles are plotted for each combination of mass hierarchy and θ_{23} octant. The best fit occurs in the inverted hierarchy and lower octant; the worst fit is the normal hierarchy and upper octant, which is disfavoured by 2ΔlogL=1.74. The dashed horizontal lines indicate the 68\% (90\%) singleparameter confidence limits, which disfavour 36\% (11\%) of the threeparameter space defined by the mass hierarchy, θ_{23} octant, and δ_{CP}. More details on the analysis can be found in arXiv:hepex/1403.0867  published in Phys. Rev. Lett. 112, 191801 (2014).
 COMBINED ANALYSIS  OBSERVED SPECTRA OF BEAM AND ATMOSPHERIC NEUTRINOS (August 2013)
 This figure shows the observed distributions of muon neutrino and muon antineutrino chargedcurrent interactions in the MINOS Far Detector. The top row shows the energy spectra for MINOS accelerator data, separated by event type and beam configuration. The bottom rows shows the zenith angle distributions of MINOS atmospheric data, separated by energy and event type. These spectra form the ν_{μ} CC and ν_{μ} CC inputs to the combined analysis of ν_{μ} disappearance and ν_{e} appearance described in Phys. Rev. Lett. 112, 191801 (2014). In each case, the gray histograms show the predictions without oscillations; the red histograms show the bestfit oscillations; and the black points show the observed data.
 ν_{μ} DISAPPEARANCE  OSCILLATION MEASUREMENT USING BEAM AND ATMOSPHERIC DATA (May 2013)
 This figure shows the 90% confidence level contours for the neutrino oscillation parameters obtained by MINOS and other experiments. The MINOS analysis uses all ν_{μ} and ν_{μ} charged current interactions, assuming identical oscillation parameters for neutrinos and antineutrinos. The black line shows the 90% CL contour obtained from a combined analysis of the MINOS beam and atmospheric neutrino data. The red line shows the analysis of MINOS beam data only. These results are compared to the the measurements from SuperKamiokande (blue) and T2K (green). The best fit point for the MINOS combined analysis is Δm^{2}=(2.41^{+0.09}_{0.10} )x10^{3} eV^{2}, sin^{2}(2Θ)=0.950^{+0.035}_{0.036}, or sin^{2}(2Θ) > 0.890 (90% CL). A tabular form of the likelihood surface from this analysis can be found here as a csv file, which lists sin^{2}(2Θ), Δm^{2} and Δln(L) in three columns. The likelihood surface is also available here as a ROOT file. More details on this analysis can be found in arXiv:hepex/1304.6335  published in Phys. Rev. Lett. 110, 2518011 (2013).
 ν_{μ} DISAPPEARANCE  ANTINEUTRINO OSCILLATION MEASUREMENT USING BEAM AND ATMOSPHERIC DATA (May 2013)
 This figure shows the 90% confidence level contours for the antineutrino oscillation parameters. The black line shows the 90% CL contour obtained from the combined analysis of all the MINOS beam and atmospheric neutrino data, allowing the neutrino and antineutrino oscillation parameters to float independently. Also shown are the MINOS antineutrino measurements using beam data only (red) and atmospheric data only (blue), as well as the SuperKamiokande antineutrino measurement (dashed). The best fit antineutrino oscillation parameters obtained by the MINOS combined analysis are Δm^{2}=(2.50^{+0.23}_{0.25} )x10^{3} eV^{2}, sin^{2}(2Θ)=0.97^{+0.03}_{0.08}, or sin^{2}(2Θ) > 0.83 (90% CL). More details on this analysis can be found in arXiv:hepex/1304.6335  published in Phys. Rev. Lett. 110, 2518011 (2013).
 ν_{μ} DISAPPEARANCE  OBSERVED SPECTRA OF BEAM AND ATMOSPHERIC NEUTRINOS (May 2013)
 This figure shows the observed distributions of muon neutrino and muon antineutrino chargedcurrent interactions in the MINOS Far Detector. The top row of graphs shows the energy spectra for MINOS beam data, the bottom row shows the L/E distributions for MINOS atmospheric data. In each case, the black line shows the prediction without oscillations; the red line shows the best fit to oscillations; and the black points show the observed data. This plot is available here as a ROOT file
 ν_{e} APPEARANCE ANALYSIS  LIMITS and BOUNDS
 ν_{e} events appearing at the Far Detector can be used as a measure of the θ_{13} mixing angle. The plots show the allowed ranges and best fits for 2sin^{2}(2θ_{13})sin^{2}(θ_{23}) as a function of the CPviolating phase δ. The upper (lower) panel assumes the normal (inverted) neutrino mass hierarchy. (The value of δ and the mass hierarchy are fixed in the fit). This fit combines data taken in the neutrinoenhanced and antineutrinoenhanced beam modes. The contour is available here as a ROOT file along with a description of the contents.