2025-11-10T03:09:47.574386

Estimating detector systematic uncertainties for the T2K far detector

Reh
Tokai to Kamioka (T2K) is a long-baseline neutrino oscillation experiment that measures oscillation parameters related to both $ν_μ(\barν_μ)$ disappearance and $ν_e(\barν_e)$ appearance in a $ν_μ(\barν_μ)$ beam. T2K uses Super-Kamiokande (SK) as its far detector, and SK detector systematic errors are currently among the leading sources of systematic uncertainty in the T2K oscillation analysis. Therefore, accurate understanding of detector mis-modelling and event reconstruction uncertainties in SK is crucial to the extraction of the neutrino oscillation parameters. The detector error estimation procedure quantifies uncertainties by fitting SK atmospheric MC to data in a Markov Chain Monte Carlo (MCMC) framework. Uncertainties are separated based on true event topology above Cherenkov threshold in the SK atmospheric MC in order to capture how the uncertainties should affect the different signal and background samples when propagated to the T2K beam oscillation analysis. The procedure has been upgraded for the upcoming analysis cycle to include systematics targeting newly added $ν_eCC1π^\pm$ and $NCπ^0$-enhanced analysis samples, among other changes, culminating in a 540 dimensional MCMC fit between atmospheric MC and data in SK. The far detector analysis chain and its integration into the broader T2K oscillation analysis will be discussed in these proceedings.
academic

Estimating detector systematic uncertainties for the T2K far detector

基本信息

  • 论文ID: 2501.00147
  • 标题: Estimating detector systematic uncertainties for the T2K far detector
  • 作者: Michael Reh (代表T2K合作组)
  • 分类: hep-ex (高能物理实验)
  • 发表时间/会议: 第25届国际加速器中微子研讨会 (25th International Workshop on Neutrinos from Accelerators)
  • 论文链接: https://arxiv.org/abs/2501.00147

摘要

T2K是一个长基线中微子振荡实验，测量与νμ(ν̄μ)消失和νe(ν̄e)出现相关的振荡参数。T2K使用Super-Kamiokande (SK)作为远探测器，SK探测器的系统误差是T2K振荡分析中系统不确定性的主要来源之一。探测器误差估计程序通过在马尔可夫链蒙特卡罗(MCMC)框架中拟合SK大气中微子蒙特卡罗模拟与数据来量化不确定性。该程序已升级为540维MCMC拟合，包含新增的νeCC1π±和NCπ0增强分析样本的系统误差处理。

研究背景与动机

核心问题

T2K实验需要精确测量中微子振荡参数，但Super-Kamiokande远探测器的系统误差已成为限制测量精度的主要因素之一。探测器建模误差和事件重建不确定性直接影响中微子振荡参数的提取精度。

重要性

  1. 物理意义重大: 中微子振荡参数的精确测量对理解基本粒子物理和宇宙演化具有重要意义
  2. 系统误差主导: SK探测器系统误差是当前T2K分析中的主要不确定性来源
  3. 精度要求: 随着统计精度的提高，系统不确定性的控制变得更加关键

现有方法局限性

传统的探测器系统误差估计方法往往依赖于特定的物理模型假设，缺乏足够的普适性和鲁棒性，难以适应探测器条件的变化和建模选择的改进。

研究动机

开发一种模型无关的、鲁棒的探测器系统误差估计方法，能够准确量化SK探测器的建模不确定性，并将其有效传递到T2K束流振荡分析中。

核心贡献

  1. 建立了基于MCMC的探测器系统误差估计框架，使用大气中微子数据作为控制样本
  2. 提出了通用的参数化方法，使用平移(β)和展宽(α)参数来描述粒子识别变量的系统偏移
  3. 开发了540维MCMC拟合程序，支持新增的νeCC1π±和NCπ0分析样本
  4. 实现了显著的计算效率提升，将分析收敛时间从月级别缩短到周级别
  5. 构建了完整的误差传递链，从大气中微子拟合到T2K束流分析的协方差矩阵

方法详解

任务定义

输入: SK大气中微子数据和蒙特卡罗模拟事件 输出: T2K束流分析中各能量区间和样本类型的探测器系统误差协方差矩阵 约束: 保持物理合理性，确保误差估计的统计显著性

核心方法架构

1. 系统误差参数化

T2K采用通用的参数化方法来描述探测器建模误差：

L → α · L + β

其中：

  • L: 粒子识别(PID)变量的标称值
  • α: 展宽参数，控制分布的宽度变化
  • β: 平移参数，控制分布的中心位置偏移

2. 大气中微子拟合

使用MCMC框架进行参数估计：

步骤:

  1. MCMC提出{α, β}参数组合
  2. 根据公式修改MC事件的PID值
  3. 将修改后的MC事件分bin构建直方图
  4. 计算数据与修改MC之间的形状似然函数
  5. 重复过程获得{α, β}的后验分布

3. 误差矩阵构建

通过玩具MC程序将大气中微子拟合结果传递到T2K束流分析：

  1. 从MCMC后验分布中随机抽取{α, β}参数
  2. 修改T2K束流MC事件的PID值
  3. 观察事件在不同分析样本间的迁移
  4. 统计中微子能谱的变化
  5. 构建协方差矩阵

技术创新点

1. 模型无关性

该方法不依赖于特定的物理模型假设，能够捕获各种潜在的探测器建模误差，包括：

  • 光散射强度建模不准确
  • 光电倍增管响应特性偏差
  • 水的光学性质变化

2. 事件拓扑分离

基于Cherenkov阈值以上的真实事件拓扑进行不确定性分离，确保系统误差能够正确影响不同的信号和本底样本。

3. 高维参数空间

升级后的540维MCMC拟合能够更精细地描述不同事件类型和运动学区域的系统误差特征。

实验设置

数据集

  1. SK大气中微子数据: 近30年的大气中微子观测数据，覆盖T2K相关的能量范围和事件类型
  2. T2K束流数据: Run 1-4数据集，包含νμ和ν̄μ束流模式
  3. 蒙特卡罗模拟: 新的MC产品，采用更新的中微子相互作用模型

评价指标

  1. 形状似然函数: 评估数据与MC分布的一致性
  2. 后验分布收敛性: 通过MCMC诊断确保参数估计的可靠性
  3. 协方差矩阵的物理合理性: 检验误差相关性的合理性

分析样本

  • 单环电子样本 (Single ring e-like)
  • 单环缪子样本 (Single ring μ-like)
  • νeCC1π+样本 (新增)
  • NCπ0增强样本 (新增)

实验结果

主要结果

1. PID变量系统偏移量化

通过{α, β}参数成功量化了电子/缪子识别等关键PID变量的系统偏移，如图2所示的不同参数组合对分布形状的影响。

2. 能谱不确定性估计

图3展示了不同T2K分析样本的后验预测能谱，bin误差反映了来自SK探测器误差矩阵的1σ不确定性。从2022年到2023年分析，通过移除若干"临时"系统误差，bin误差有显著差异。

3. 计算效率提升

实现了显著的性能改进：

  • 收敛时间: 从O(月)缩短到O(周)
  • 参数维度: 扩展到540维MCMC拟合
  • 样本覆盖: 支持新增分析样本的系统误差估计

方法验证

通过比较大气中微子和束流中微子在运动学范围、事件拓扑和探测器条件方面的重叠，验证了将大气中微子拟合结果应用于T2K束流分析的合理性。

相关工作

中微子振荡实验系统误差

长基线中微子振荡实验普遍面临探测器系统误差的挑战，包括NOνA、DUNE等实验都在开发相应的误差估计方法。

探测器建模技术

水Cherenkov探测器的精确建模一直是该领域的重要课题，涉及光学性质、PMT响应、事件重建算法等多个方面。

本文优势

相比其他方法，本文提出的框架具有：

  1. 更强的普适性: 模型无关的参数化方法
  2. 更高的精度: 540维参数空间的精细描述
  3. 更好的鲁棒性: 能够适应探测器条件变化

结论与讨论

主要结论

  1. 成功建立了基于MCMC的SK探测器系统误差估计框架
  2. 实现了从大气中微子数据到T2K束流分析的有效误差传递
  3. 为即将到来的T2K分析提供了重要的系统误差控制手段

局限性

  1. 假设限制: 假设大气中微子和束流中微子的系统误差特征相似
  2. 模型依赖: 仍然依赖于MC模拟的基本准确性
  3. 计算复杂度: 540维MCMC拟合仍需要相当的计算资源

未来方向

  1. 进一步优化MCMC算法，提高收敛效率
  2. 探索更精细的事件分类和系统误差建模
  3. 结合其他控制样本验证误差估计的准确性

深度评价

优点

  1. 方法创新性强: 提出了模型无关的系统误差参数化方法
  2. 实验设计合理: 使用大气中微子作为控制样本具有很强的物理动机
  3. 技术实现完善: 540维MCMC拟合展现了高水平的技术能力
  4. 实用价值高: 直接服务于T2K物理分析的精度提升

不足

  1. 理论分析有限: 缺乏对参数化方法理论基础的深入讨论
  2. 验证手段单一: 主要依赖大气中微子数据，缺乏独立验证
  3. 不确定性量化: 对方法本身引入的不确定性讨论不足

影响力

  1. 学术贡献: 为中微子物理实验提供了重要的系统误差控制方法
  2. 实用价值: 直接影响T2K实验的物理结果精度
  3. 技术推广: 方法可推广到其他类似的粒子物理实验

适用场景

该方法特别适用于：

  1. 长基线中微子振荡实验
  2. 大型水Cherenkov探测器
  3. 需要高精度系统误差控制的粒子物理测量

参考文献

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)

总体评价: 这是一篇技术性很强的实验物理论文，为T2K实验的系统误差控制提供了重要的方法学贡献。虽然在理论分析方面有所不足，但其实用价值和技术创新性使其成为中微子物理领域的重要工作。