Abstract Measuring observables to constrain models using maximum-likelihood estimation is fundamental to many physics experiments. Wilks' theorem provides a simple way to construct confidence intervals on model parameters, but it only applies under certain conditions. These conditions, such as nested hypotheses and unbounded parameters, are often violated in neutrino oscillation measurements and other experimental scenarios. Monte Carlo methods can address these issues, albeit at increased computational cost. In the presence of nuisance parameters, however, the best way to implement a Monte Carlo method is ambiguous. This paper documents the method selected by the NOvA experiment, the profile construction. It presents the toy studies that informed the choice of method, details of its implementation, and tests performed to validate it. It also includes some practical considerations which may be of use to others choosing to use the profile construction.
more »
« less
Grid-based minimization at scale: Feldman-Cousins corrections for light sterile neutrino search
High Energy Physics (HEP) experiments generally employ sophisticated statistical methods to present results in searches of new physics. In the problem of searching for sterile neutrinos, likelihood ratio tests are applied to short-baseline neutrino oscillation experiments to construct confidence intervals for the parameters of interest. The test statistics of the form Δχ 2 is often used to form the confidence intervals, however, this approach can lead to statistical inaccuracies due to the small signal rate in the region-of-interest. In this paper, we present a computational model for the computationally expensive Feldman-Cousins corrections to construct a statistically accurate confidence interval for neutrino oscillation analysis. The program performs a grid-based minimization over oscillation parameters and is written in C++. Our algorithms make use of vectorization through Eigen3, yielding a single-core speed-up of 350 compared to the original implementation, and achieve MPI data parallelism by employing DIY. We demonstrate the strong scaling of the application at High-Performance Computing (HPC) sites. We utilize HDF5 along with HighFive to write the results of the calculation to file.
more »
« less
- Award ID(s):
- 2013070
- PAR ID:
- 10338866
- Editor(s):
- Biscarat, C.; Campana, S.; Hegner, B.; Roiser, S.; Rovelli, C.I.; Stewart, G.A.
- Date Published:
- Journal Name:
- EPJ Web of Conferences
- Volume:
- 251
- ISSN:
- 2100-014X
- Page Range / eLocation ID:
- 02065
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
A bstract We evaluate the statistical significance of the 3+1 sterile-neutrino hypothesis using ν e and $$ \overline{\nu} $$ ν ¯ e disappearance data from reactor, solar and gallium radioactive source experiments. Concerning the latter, we investigate the implications of the recent BEST results. For reactor data we focus on relative measurements independent of flux predictions. For the problem at hand, the usual χ 2 -approximation to hypothesis testing based on Wilks’ theorem has been shown in the literature to be inaccurate. We therefore present results based on Monte Carlo simulations, and find that this typically reduces the significance by roughly 1 σ with respect to the naïve expectation. We find no significant indication in favor of sterile-neutrino oscillations from reactor data. On the other hand, gallium data (dominated by the BEST result) show more than 5 σ of evidence supporting the sterile-neutrino hypothesis, favoring oscillation parameters in agreement with constraints from reactor data. This explanation is, however, in significant tension (∼ 3 σ ) with solar neutrino experiments. In order to assess the robustness of the signal for gallium experiments we present a discussion of the impact of cross-section uncertainties on the results.more » « less
-
In the presence of extra neutrino states at high scales, the low-energy effective leptonic mixing matrix (LMM) is in general nonunitary. We revisit the question of what is our current knowledge of individual LMM matrix elements without assuming unitarity. We first demonstrate that a minimal set of experimental constraints suffices in bounding LMM nonunitarity parameters to the level of , without the use of neutrino oscillation data. We then revisit oscillation results as a complementary cross-check, using different physics and different experimental techniques to probe a similar parameter space. We correct some common misconceptions in the neutrino nonunitarity literature resulting from an incautious treatment of input parameters. We find that neutrino oscillation experiments can constrain LMM nonclosure, but, contrary to claims in the literature, are completely insensitive to the overall normalization of the LMM. Thus, we conclude that oscillation experiments, including the future DUNE experiment, have no power in excluding nonunitarity altogether. Published by the American Physical Society2024more » « less
-
A<sc>bstract</sc> We derive new constraints on effective four-fermion neutrino non-standard interactions with both quarks and electrons. This is done through the global analysis of neutrino oscillation data and measurements of coherent elastic neutrino-nucleus scattering (CEνNS) obtained with different nuclei. In doing so, we include not only the effects of new physics on neutrino propagation but also on the detection cross section in neutrino experiments which are sensitive to the new physics. We consider both vector and axial-vector neutral-current neutrino interactions and, for each case, we include simultaneously all allowed effective operators in flavour space. To this end, we use the most general parametrization for their Wilson coefficients under the assumption that their neutrino flavour structure is independent of the charged fermion participating in the interaction. The status of the LMA-D solution is assessed for the first time in the case of new interactions taking place simultaneously with up quarks, down quarks, and electrons. One of the main results of our work are the presently allowed regions for the effective combinations of non-standard neutrino couplings, relevant for long-baseline and atmospheric neutrino oscillation experiments.more » « less
-
null (Ed.)A bstract Our herein described combined analysis of the latest neutrino oscillation data presented at the Neutrino2020 conference shows that previous hints for the neutrino mass ordering have significantly decreased, and normal ordering (NO) is favored only at the 1 . 6 σ level. Combined with the χ 2 map provided by Super-Kamiokande for their atmospheric neutrino data analysis the hint for NO is at 2 . 7 σ . The CP conserving value δ CP = 180° is within 0 . 6 σ of the global best fit point. Only if we restrict to inverted mass ordering, CP violation is favored at the ∼ 3 σ level. We discuss the origin of these results — which are driven by the new data from the T2K and NOvA long-baseline experiments —, and the relevance of the LBL-reactor oscillation frequency complementarity. The previous 2 . 2 σ tension in ∆ m 2 21 preferred by KamLAND and solar experiments is also reduced to the 1 . 1 σ level after the inclusion of the latest Super-Kamiokande solar neutrino results. Finally we present updated allowed ranges for the oscillation parameters and for the leptonic Jarlskog determinant from the global analysis.more » « less