Estimating detector systematic uncertainties for the T2K far detector

Basic Information

• Paper ID: 2501.00147
• Title: Estimating detector systematic uncertainties for the T2K far detector
• Author: Michael Reh (on behalf of the T2K Collaboration)
• Classification: hep-ex (High Energy Physics - Experiment)
• Publication Time/Conference: 25th International Workshop on Neutrinos from Accelerators
• Paper Link: https://arxiv.org/abs/2501.00147

Abstract

T2K is a long-baseline neutrino oscillation experiment that measures oscillation parameters related to νμ(ν̄μ) disappearance and νe(ν̄e) appearance. T2K uses Super-Kamiokande (SK) as its far detector, and the systematic uncertainties of the SK detector constitute one of the primary sources of systematic uncertainty in the T2K oscillation analysis. The detector error estimation procedure quantifies uncertainties by fitting SK atmospheric neutrino Monte Carlo simulations to data within a Markov Chain Monte Carlo (MCMC) framework. The procedure has been upgraded to a 540-dimensional MCMC fit, incorporating systematic error treatment for newly added νeCC1π± and NCπ0 enhanced analysis samples.

Research Background and Motivation

Core Problem

The T2K experiment requires precise measurement of neutrino oscillation parameters, but systematic uncertainties from the Super-Kamiokande far detector have become one of the primary limiting factors for measurement precision. Detector modeling errors and event reconstruction uncertainties directly impact the extraction accuracy of neutrino oscillation parameters.

Significance

1. Substantial Physical Importance: Precise measurement of neutrino oscillation parameters is crucial for understanding fundamental particle physics and cosmic evolution
2. Systematic Error Dominance: SK detector systematic uncertainties represent the primary source of uncertainty in current T2K analyses
3. Precision Requirements: As statistical precision improves, control of systematic uncertainties becomes increasingly critical

Limitations of Existing Methods

Traditional detector systematic error estimation methods often rely on specific physical model assumptions, lacking sufficient universality and robustness to adapt to detector condition changes and improvements in modeling choices.

Research Motivation

To develop a model-independent, robust detector systematic error estimation method capable of accurately quantifying SK detector modeling uncertainties and effectively propagating them to the T2K beam oscillation analysis.

Core Contributions

1. Established an MCMC-based detector systematic error estimation framework using atmospheric neutrino data as a control sample
2. Proposed a universal parameterization method using shift (β) and broadening (α) parameters to describe systematic offsets in particle identification variables
3. Developed a 540-dimensional MCMC fitting program supporting newly added νeCC1π± and NCπ0 analysis samples
4. Achieved significant computational efficiency improvements, reducing analysis convergence time from months to weeks
5. Constructed a complete error propagation chain from atmospheric neutrino fitting to covariance matrices in T2K beam analysis

Detailed Methodology

Task Definition

Input: SK atmospheric neutrino data and Monte Carlo simulated events Output: Detector systematic error covariance matrices for various energy intervals and sample types in T2K beam analysis Constraints: Maintain physical reasonableness and ensure statistical significance of error estimates

Core Methodological Architecture

1. Systematic Error Parameterization

T2K employs a universal parameterization method to describe detector modeling errors:

L → α · L + β

Where:

• L: Nominal value of particle identification (PID) variable
• α: Broadening parameter controlling distribution width changes
• β: Shift parameter controlling distribution center position offset

2. Atmospheric Neutrino Fitting

Parameter estimation using MCMC framework:

Steps:

1. MCMC proposes {α, β} parameter combinations
2. Modify MC event PID values according to the formula
3. Construct histograms from modified MC events binned by topology
4. Calculate shape likelihood function between data and modified MC
5. Repeat to obtain posterior distribution of {α, β}

3. Error Matrix Construction

Propagate atmospheric neutrino fitting results to T2K beam analysis via toy MC program:

1. Randomly sample {α, β} parameters from MCMC posterior distribution
2. Modify T2K beam MC event PID values
3. Observe event migration between different analysis samples
4. Quantify neutrino energy spectrum variations
5. Construct covariance matrix

Technical Innovations

1. Model Independence

The method does not depend on specific physical model assumptions, capable of capturing various potential detector modeling errors, including:

• Inaccurate light scattering intensity modeling
• Photomultiplier tube response characteristic deviations
• Water optical property variations

2. Event Topology Separation

Uncertainty separation based on true event topologies above Cherenkov threshold, ensuring systematic errors correctly affect different signal and background samples.

3. High-Dimensional Parameter Space

The upgraded 540-dimensional MCMC fit enables finer description of systematic error characteristics across different event types and kinematic regions.

Experimental Setup

Datasets

1. SK Atmospheric Neutrino Data: Nearly 30 years of atmospheric neutrino observations covering T2K-relevant energy ranges and event types
2. T2K Beam Data: Run 1-4 datasets including νμ and ν̄μ beam modes
3. Monte Carlo Simulations: New MC production employing updated neutrino interaction models

Evaluation Metrics

1. Shape Likelihood Function: Assesses consistency between data and MC distributions
2. Posterior Distribution Convergence: MCMC diagnostics ensure reliable parameter estimation
3. Covariance Matrix Physical Reasonableness: Validates error correlation rationality

Analysis Samples

• Single ring electron-like sample (Single ring e-like)
• Single ring muon-like sample (Single ring μ-like)
• νeCC1π+ sample (newly added)
• NCπ0 enhanced sample (newly added)

Experimental Results

Main Results

1. PID Variable Systematic Offset Quantification

Successfully quantified systematic offsets in critical PID variables such as electron/muon identification through {α, β} parameters, as shown in Figure 2 depicting the impact of different parameter combinations on distribution shapes.

2. Energy Spectrum Uncertainty Estimation

Figure 3 presents posterior predicted energy spectra for different T2K analysis samples, with bin errors reflecting 1σ uncertainties from the SK detector error matrix. Significant differences in bin errors appear between 2022 and 2023 analyses through removal of certain "temporary" systematic errors.

3. Computational Efficiency Improvements

Achieved significant performance enhancements:

• Convergence Time: Reduced from O(months) to O(weeks)
• Parameter Dimensionality: Extended to 540-dimensional MCMC fitting
• Sample Coverage: Supports systematic error estimation for newly added analysis samples

Method Validation

Validated the rationality of applying atmospheric neutrino fitting results to T2K beam analysis by comparing overlaps in kinematic ranges, event topologies, and detector conditions between atmospheric and beam neutrinos.

Related Work

Neutrino Oscillation Experiment Systematic Uncertainties

Long-baseline neutrino oscillation experiments universally face detector systematic error challenges, with experiments such as NOνA and DUNE developing corresponding error estimation methods.

Detector Modeling Techniques

Precise modeling of water Cherenkov detectors remains an important topic in the field, involving multiple aspects including optical properties, PMT response, and event reconstruction algorithms.

Advantages of This Work

Compared to other methods, the proposed framework offers:

1. Greater Universality: Model-independent parameterization method
2. Higher Precision: Fine description through 540-dimensional parameter space
3. Better Robustness: Adaptability to detector condition changes

Conclusions and Discussion

Main Conclusions

1. Successfully established an MCMC-based SK detector systematic error estimation framework
2. Achieved effective error propagation from atmospheric neutrino data to T2K beam analysis
3. Provided important systematic error control means for forthcoming T2K analyses

Limitations

1. Assumption Constraints: Assumes similar systematic error characteristics between atmospheric and beam neutrinos
2. Model Dependence: Still relies on basic accuracy of MC simulations
3. Computational Complexity: 540-dimensional MCMC fitting requires substantial computational resources

Future Directions

1. Further optimize MCMC algorithms to improve convergence efficiency
2. Explore finer event classification and systematic error modeling
3. Combine other control samples to validate error estimation accuracy

In-Depth Evaluation

Strengths

1. Strong Methodological Innovation: Proposes model-independent systematic error parameterization method
2. Reasonable Experimental Design: Using atmospheric neutrinos as control sample has strong physical motivation
3. Sophisticated Technical Implementation: 540-dimensional MCMC fitting demonstrates high technical capability
4. High Practical Value: Directly serves precision enhancement in T2K physics analysis

Weaknesses

1. Limited Theoretical Analysis: Lacks in-depth discussion of theoretical foundations for parameterization method
2. Single Verification Approach: Primarily relies on atmospheric neutrino data, lacking independent verification
3. Insufficient Uncertainty Quantification: Inadequate discussion of uncertainties introduced by the method itself

Impact

1. Academic Contribution: Provides important systematic error control method for neutrino physics experiments
2. Practical Value: Directly influences precision of T2K experimental physics results
3. Technical Generalization: Method generalizable to other similar particle physics experiments

Applicable Scenarios

This method is particularly suitable for:

1. Long-baseline neutrino oscillation experiments
2. Large-scale water Cherenkov detectors
3. Particle physics measurements requiring high-precision systematic error control

References

1 K. A. et al., The T2K experiment, NIM-A 659, 106 (2011) 2 S. F. et al., The Super-Kamiokande detector, NIM-A 501, 418 (2003) 3 K. A. et al., Measurements of neutrino oscillation parameters from the T2K experiment using 3.6×10²¹ protons on target, Eur. Phys. J. C 83, 782 (2023)

Overall Assessment: This is a technically rigorous experimental physics paper that provides important methodological contributions to systematic error control in the T2K experiment. While somewhat limited in theoretical analysis, its practical value and technical innovation make it an important work in the neutrino physics field.