PhyStat-v Workshop on Statistical Issues in Experimental Neutrino Physics
→
Asia/Tokyo
Yoshi Uchida
(Imperial College London)
Description
Overview
Neutrino experiments are reaching exciting new levels of precision: percent-level systematic uncertainties in oscillation experiments, the eV scale in direct mass measurements, 1025 year half-lives in double beta decay experiments, and neutrino-nucleus cross section measurements at the 5% level. With this new precision comes an important focus on the statistical understanding of data from neutrino experiments, and techniques to ensure the quality of the results produced.
The PhyStat series of workshops has been a successful institution in the world of collider physics over the past 15 years, providing a forum for discussing topics such as setting upper limits on or discovering rare processes, and combining results from different experiments. We now aim to turn this success to the field of neutrino physics.
PhyStat-ν will address statistical issues in the broad range of modern neutrino physics, including making measurements of model parameters, setting limits, discovery criteria, discrete choice of models, Bayesian vs. frequentist inference, and combining experiments. Two workshops are planned for 2016: this workshop at the Kavli Institute for the Physics and Mathematics of the Universe in Kashiwa, Japan in May, and a second workshop at Fermilab in the autumn. These workshops will consist of invited and contributed talks, and poster and discussion sessions, and will bring together experts in the analysis of neutrino data with experts in statistics to explore the statistical issues in neutrino physics.
PhyStat-ν will address statistical issues in the broad range of modern neutrino physics, including making measurements of model parameters, setting limits, discovery criteria, discrete choice of models, Bayesian vs. frequentist inference, and combining experiments. Two workshops are planned for 2016: this workshop at the Kavli Institute for the Physics and Mathematics of the Universe in Kashiwa, Japan in May, and a second workshop at Fermilab in the autumn. These workshops will consist of invited and contributed talks, and poster and discussion sessions, and will bring together experts in the analysis of neutrino data with experts in statistics to explore the statistical issues in neutrino physics.
We are looking forward to seeing you in Kashiwa.
Local Organizing Committee Members:
Chair: Mark HARTZ (Kavli IPMU/TRIUMF)
Christophe BRONNER (Kavli IPMU)
Richard CALLAND (Kavli IPMU)
Yoshinari HAYATO (ICRR)
Yasuhiro NISHIMURA (ICRR)
Kimihiro OKUMURA (ICRR)
Scientific Organizing Committee Members:
Chair: Yoshi UCHIDA (Imperial College London)
Jun CAO (IHEP, CAS)
Daniel CHERDACK (Colorado State University)
Robert COUSINS (UCLA)
David VAN DYK (Imperial College London)
Joe FORMAGGIO (MIT)
Mark HARTZ (Kavli IPMU (WPI), The University of Tokyo)
Pilar HERNANDEZ (University of Valencia)
Thomas JUNK (Fermilab)
Asher KABOTH (Royal Holloway, University of London)
Louis LYONS (Imperial College London)
Shun SAITO (Kavli IPMU (WPI), The University of Tokyo)
Subir SARKAR (University of Oxford, Neils Bohr Institute)
Elizabeth WORCESTER (BNL)
Kai ZUBER (TU Dresden)
Jun CAO (IHEP, CAS)
Daniel CHERDACK (Colorado State University)
Robert COUSINS (UCLA)
David VAN DYK (Imperial College London)
Joe FORMAGGIO (MIT)
Mark HARTZ (Kavli IPMU (WPI), The University of Tokyo)
Pilar HERNANDEZ (University of Valencia)
Thomas JUNK (Fermilab)
Asher KABOTH (Royal Holloway, University of London)
Louis LYONS (Imperial College London)
Shun SAITO (Kavli IPMU (WPI), The University of Tokyo)
Subir SARKAR (University of Oxford, Neils Bohr Institute)
Elizabeth WORCESTER (BNL)
Kai ZUBER (TU Dresden)
Contact
Participants
-
-
Registration
-
10:30
Registration and Tea and Coffee
-
Opening Session
Introductory Overview, with talks on general neutrino physics and the topics to be discussed at the workshop
- 2
- 3
- 4
-
12:10
Lunch
-
Statistical Methods at Current Experiments
Overview talks on how current and recent experiments are treating their data
- 5
- 6
-
14:50
Tea and Coffee
- 7
- 8
-
Short Poster Talks
Short introductory talks from poster presenters
-
9
Estimation of neutrino oscillation parameters with Markov Chain Monte CarloThe T2K experiment is a long-baseline accelerator neutrino experiment with the ability to measure Δm32, θ23, θ13 and δCP from (anti-)νμ disappearance as well as (anti-)νe appearance. In this poster we show a Bayesian framework using the Metropolis-Hastings algorithm to estimate the oscillation parameters. The advantage of Markov Chain Monte Carlo is its ability to scan a multidimensional parameter space, allowing to evaluate at the same time the oscillation parameters and the systematic uncertainties on the neutrino interaction model, accelerator flux model and detector reconstruction. We present the results obtained for a simultaneous fit of νμ disappearance and νe appearance, and a probability of anti-νe appearance.Speaker: Ms Leïla HAEGEL (Université de Genève)
-
10
Statistical techniques used for neutrino oscillation analyses with the VALOR frameworkVALOR is a well established neutrino fitting group that is leading several neutrino oscillation analyses in T2K, producing world leading results. The analysis framework has been validated and optimised for both a hybrid Bayesian-Frequentist approach to fitting parameters, where the systematic parameters are marginalised, as well as a frequentist analysis, that consists of the likelihood ratio maximisation as a function of all the fit parameters. The VALOR group is involved in many other neutrino oscillation experiments such as DUNE and Hyper-K, which goal is the measurement of δCP , and the Fermilab Short Baseline Neutrino program, a novel experiment aiming to search for sterile neutrinos. In this talk the statistical techniques used in the neutrino oscillation analyses performed by the VALOR group are shown.Speaker: Davide Sgalaberna (University of Geneva)
-
11
Anti-𝜈e Appearance at T2K using VALORVALOR is a well established neutrino fitting group that is leading several neutrino oscillation analyses in T2K, producing world leading results. The analysis framework has been validated and optimised for both a hybrid Bayesian-Frequentist approach to fitting parameters, where the systematic parameters are marginalised, as well as a frequentist analysis, that consists of the likelihood ratio maximisation as a function of all the fit parameters. The VALOR group is involved in many other neutrino oscillation experiments such as DUNE and Hyper-K, which goal is the measurement of δCP , and the Fermilab Short Baseline Neutrino program, a novel experiment aiming to search for sterile neutrinos. In this talk the statistical techniques used in the neutrino oscillation analyses performed by the VALOR group are shown.Speaker: Mr raj shah (STFC, UK)
-
12
Tuning generators to the global neutrino cross-section dataset.The T2K Neutrino Oscillation Experiment uses a prior estimation of a scattering cross-section model when generating nominal Monte-Carlo (MC) neutrino interaction events. A selection of the available models in the NEUT MC generator are fit to a range of historical neutrino interaction cross-section measurements made external to T2K. A nominal model is then chosen which is capable of finding the best agreement with the fitted data. One difficulty that must be overcome when performing these fits is the need to include datasets where only limited information on the observed event kinematics and associated systematic uncertainties are provided. In this poster we review the current efforts to produce tighter constraints on the neutrino interaction cross-section model. We also highlight recent attempts to use multiple kinematic distributions from a single experiment when performing likelihood fits, whilst using the MC model to approximate the true correlations between the different kinematic variables.Speaker: Mr John Stowell (University of Sheffield)
-
13
Empirical Parametrisations of Neutrino OscillationsThe discovery of neutrino oscillations driven by the U[e3] element of the PMNS mixing matrix has prompted a step change in the way neutrino oscillation experiments are analysed, commonly described by the shorthand of "switching from a 2-neutrino to 3-neutrino picture". The difference is less clear cut than this. On the one hand, some older measurements used parametrisations that implicitly assumed at least three neutrino flavours (and less commonly required 2 independent mass splittings). The more interesting change relates to the way oscillations are parametrised for fitting and inference. "Two flavour" analyses used parameters that -- regardless of how they were named -- directly described the shape of the observed data (the location and depth of a 'dip'), whereas newer analyses are explicitly concerned with mixing angles that are not direct observables and only related to the data through mathematical transformations. This has non-trivial consequences for inferences from the data. This poster highlights some of these unintended consequences and considers whether alternative formulations for analysing and presenting data are possible.Speakers: Dr Phillip Litchfield (Imperial College London), Dr Phillip Litchfield (Imperial College London)
-
14
Neutrino Oscillation on EpicyclesGlobal fit of neutrino oscillations is shown to be equivalent to Fourier transformation in a complex plane. Adding new sterile neutrinos or new oscillation modes is equivalent to adding Fourier components. Therefore new sterile neutrinos should be introduced with physics motivations beyond phenomenology.Speaker: Benda Xu (RCNS, Tohoku Univ., Japan)
-
15
Study of the Neutrino Oscillation at Very Short BaselineVerification of neutrino oscillations at a very short baseline of a nuclear reactor is one way to check the possible existence of the sterile neutrinos. Challenging features of this short baseline experiment come from: the high muon rates at shallow depth of overburden, the escaping gamma energy due to the small size of the detector in a limited space, and the inaccurate predictions of the reactor anti-neutrino flux. The oscillation parameters of the sterile neutrino have been studied using the energy spectrum from the NEOS (Neutrino Experiment for Oscillation at Short baseline) data collected for 6 months at 24 m baseline of a 3 GW_th commercial reactor.Speaker: Mr Youngju Ko (Chung-Ang University)
-
16
Optimization of detector design of NuPRISM, a new water-Cherenkov neutrino near-detectorNuPRISM is a proposed novel water Cherenkov detector in the J-PARC neutrino beam line that will make neutrino cross section measurements critical to both T2K and T2HK as well as searches for short baseline neutrino oscillations. By observing neutrino interactions over off-axis angles ranging from 1 to 4 degrees with NuPRISM detector, systematic uncertainties for T2K and T2HK oscillation measurements can be significantly reduced. Evaluation of the sensitivity to neutrino interaction and oscillation physics with the NuPRISM detector is ongoing. In this poster, the event reconstruction performance for different configurations of photosensor sizes and photocathode coverages are shown.Speaker: Tomoyo Yoshida (Tokyo Institute of Technology)
-
17
Boosted Decision Tree Approach to Track Finding in the COMET ExperimentThe Coherent Muon to Electron Transition (COMET) experiment is designed to search for muon to electron conversion, a process which has very good sensitivity to Beyond the Standard Model physics. The first phase of the experiment is currently under construction at J-PARC. This phase is designed to probe muon to electron conversion 100 times better than the current limit. The experiment will achieve this sensitivity by directing a high intensity muon beam at a stopping target. The detectors probe the resulting events for the signal 105 MeV electron from muon to electron conversion. A boosted decision tree (BDT) algorithm has been developed to find this signal track. This BDT is used to combine energy deposition and timing information with a reweighted inverse hough transform to filter out background hits. The resulting hits are fit using a RANdom SAmple Consensus (RANSAC) fit, which chooses the best fit parameters for an optimized selection of the filtered hits. These hits are then passed to the track fitting algorithm.Speaker: Mr Ewen Gillies (Imperial College London)
-
9
-
Pre-Reception Talk
-
18
Welcome from the DirectorSpeaker: Hitoshi Murayama (Kavli IPMU)
-
18
-
18:00
Poster Session and Drinks Reception
-
-
-
Statistical Methods at Current Experiments
Overview talks on how current and recent experiments are treating their data
- 19
- 20
-
10:20
Tea and Coffee
- 21
-
12:40
Lunch
-
Contributed and invited talks and discussions on statistical methods and issues
-
24
Searching for PeV neutrinos with IceCubeThe IceCube Neutrino Observatory at the geographic South Pole has been taking high-quality data in the deep-ice since 2010. Extremely high energy neutrinos that carry energies above 1 PeV have been observed. The most constraining limits on the flux of ultra-high energy neutrinos, which could come from cosmogenic origins have been obtained based on a binned Poisson-likelihood-ratio method. In this talk, we will discuss the implementation of this statistical tool in light of an extremely rare event scenario.Speaker: Dr Lu Lu (Chiba University)
- 25
- 26
-
15:45
Tea and Coffee
- 27
- 28
- 29
-
24
-
18:00
Travel to Dinner Venue
-
18:30
Workshop Dinner
-
-
-
Contributed and invited talks and discussions on statistical methods and issues
- 30
- 31
-
10:10
Tea and Coffee
-
32
Sensitivity to the Neutrino Mass HierarchyPrecisely one neutrino mass hierarchy is realized in nature. How can this crucial fact be incorporated into a statistical analysis of the mass hierarchy? It implies that the problem is not hypothesis testing, but rather model selection. The corresponding tools are standard, but different. We will discuss some statistical issues that will be relevant in this kind of experiment, present a simple formula for the corresponding sensitivity and confront it with simulations of the experiment JUNO.Speaker: Emilio Ciuffoli (IMP, CAS)
-
33
The KM3NeT/ORCA detectorThe ORCA detector is the low-energy branch of the new generation neutrino telescope KM3NeT which is currently being constructed in the Mediterranean Sea. Interactions of neutrinos that are created in cosmic ray interactions travel and through the Earth can be observed in the detector. By measuring the matter-induced oscillation effects on this flux, KM3NeT/ORCA aims at settling the question on the neutrino mass ordering and to constrain other oscillation parameters. In this contribution an introduction to the ORCA detector will be given, together with the Monte-Carlo simulation efforts and statistical methods used to estimate the sensitivity to determine whether the neutrino mass hierarchy is normal or inverted.Speaker: Dr Ronald Bruijn (University of Amsterdam/Nikhef)
- 34
-
35
Multiple Hypothesis Testing and Testing One Hypothesis Multiple Times: a unified (re)viewIn statistics, the problem of testing one hypothesis multiple times can be formulated in terms of hypothesis testing when a nuisance parameter is present only under the alternative, also known as “look-elsewhere effect”. Each possible value of the nuisance parameter specifies a different alternative hypothesis and a unique global p-value is provided to summarize the statistical evidence in support (or against) the null hypothesis. From a physics perspective, this scenario occurs quite often in the searches for new signals over an energy or mass spectrum, and in both nested and non-nested frameworks. An alternative way to search for new emissions is to refer to the classical and widely known multiple hypothesis testing approach. Separate tests of hypothesis are conducted at different locations producing an ensemble of local p-values, the smallest is reported as evidence for the new resonance, once adequately adjusted to control the false detection rate (type I error rate). The aim of this talk is to review both approaches, and propose a graphical tool to identify recurrent scenarios where a simple multiple hypothesis testing procedure can be used to provide valid inference with respect to stringent significance requirements, without encountering the usual problem of over-conservativeness.Speaker: Ms Sara Algeri (Imperial College London)
-
13:15
Lunch
-
Panel Discussion
Discussion on matters arising during the workshop
-
36
Panel DiscussionSpeakers: Dean Karlen (University of Victoria), Dr Elizabeth Worcester (BNL), Michael Betancourt (University of Warwick), Ms Sara Algeri (Imperial College London), Dr Yoshi Uchida (Imperial College London)
-
15:30
Tea and Coffee
-
36
-
Closing Session
Experimental Summary talkby Bob Cousins, Theoretical Summary
- 37
- 38
-
39
Closing Statements
-