An effective topology-transfer nuisance for Genie-based neutrino-argon semi-inclusive analyses: a methodological proposal with MicroBooNE-constrained diagnostic

技术摘要：一种针对基于 Genie 的中微子-氩核半包含分析的高效拓扑转移干扰参数

问题陈述
液氩时间投影室（LArTPC）振荡分析正越来越多地根据约 35 MeV 的质子重建阈值，将中微子相互作用样本拆分为质子标记型（$N_p$）和无质子型（$0p$）拓扑结构。在标准的 GENIE Monte Carlo 重加权库中存在一个关键空白：它缺乏一个能够将事件在 35 MeV 边界两侧进行转移的“拨盘”。虽然 GENIE 的内核级联（INC）模型通过恒定的移除能量（$W_s \approx 25$ MeV）来模拟核结合，但它缺乏传输代码（如 GiBUU）中存在的与位置相关的平均场势能。这种缺失导致 GEN_IE 对近阈值质子（$T_p \approx 35\text{--}100$ MeV）的能量损失估计不足，从而导致 GENIE 将本应属于 $0p$ 的事件错误分类为 $N_p$。这种“质子可见性迁移”扭曲了 $0p/N_p$ 的比例，直接影响能量重建策略和振荡系统误差。现有的 GENIE 重加权拨盘（如 MFP$_\pi$, MFP$_N$）在特定拓扑内起作用，但无法在阈值之间重新分配事件。

方法论
作者提出了一个单一的经验干扰参数 $f_{mig}$，其定义为对标称 GENIE 输出的逐事件重加权。该参数作为一个有效的拓扑转移代理，吸收了未指明的多种近阈值建模误差组合（包括质子侧末态相互作用、核子动量分布、RES/MEC 运动学形状以及残余平均场能量损失）。

该方法论包括：

1. 物理动机： 使用 $^{40}$Ar Woods-Saxon 密度分布和 Walecka 型相对论平均场势能（GiBUU NL3* 深度）进行的半解析估计表明，GENIE 缺失了约 7.8 MeV 的残余能量损失。这为 $f_{mig}$ 的量级提供了一个参考尺度。
2. 参数化： 迁移通过将事件从 $N_p$ 拓扑线性转移到 $0p$ 拓扑来实现：
   • $N_{0p} \to N_{0p} + f_{mig} N_{Np}$
   • $N_{Np} \to N_{Np} (1 - f_{mig})$
     这在改变拓扑划分的同时保持了总事件率不变。
3. 数据约束： 该参数使用两组公开的 MicroBooNE 数据集进行约束：
   • 主要约束： 直接针对同步 CC $0p/N_p$ 微分截面发布（PRD 110）进行扫描。此分析使用未经过 sector-weight 正则化的展开截面，仅对整体归一化进行剖析（profiling）。
   • 次要约束： 对 NC $\Delta$ 辐射侧带协方差（64 个 bin，4 个通道）进行拟合。该拟合包含五个 sector 归一化权重（CCQE+MEC, CCRES, CCDIS, NCRES, NCDIS）以及岭正则化（高斯先验）和两个迁移分数（$f^{CC}_{mig}, f^{NC}_{mig}$）。
4. 统计处理： 通过剖析 $\Delta\chi^2$ 扫描评估显著性。由于侧带协方差的有效自由度较低（$n_{eff} \approx 2.5$）且存在物理边界（$f_{mig} \ge 0$），作者采用了感知边界的检验统计量，并通过 5,000 次伪实验（toys）来校准经验 p 值。

核心贡献

• 识别缺失的干扰参数： 本文识别了 GENIE 重加权基准中的一个特定缺口：即无法模拟靠近质子检测阈值的拓扑转换末态相互作用。
• 提出 $f_{mig}$： 引入了一个作为有效拓扑转移代理的方法论诊断参数，它可以实现为平坦的逐事件权重而不涉及事件再生，这使得它在经过验证后适用于未来的振荡分析（如 DUNE）。
• 诊断副产品： 分析表明，引入 $f_{mig}$ 缓解了在迁移盲拟合中观察到的 CCRES sector 权重补偿性抑制现象。在缺乏迁移拨盘的情况下，拟合会人为地抑制 CCRES（$w_{CCRES} \approx 0.69$}）；而引入 $f_{mig}$ 后，CCRES 权重向单位值回归（$w_{CCRES} \approx 1.07$）。
• 鲁棒性检查： 文中包含了关于先验敏感性、协方差结构（Peelle's Pertinent Puzzle）以及与现有强子学拨盘（如 MFP$_N$）正交性的广泛检查。

结果

• 主要约束 (PRD 110)： 直接 CC 扫描得出最佳拟合迁移分数 $\hat{f}_{mig} = 0.093 \pm 0.024$。这对应于在没有 sector-weight 正则化的情况下，对 $f_{mig} > 0$ 的偏好达到了 3.83$\sigma$ (Wilks) 和 3.24$\sigma$ (经验)。
• 次要约束 (侧带)： NC $\Delta$ 侧带拟合得出 $\hat{f}^{CC}_{mig} = 0.050^{+0.029}_{-0.022}$。在 Wilks 检验下，显著性范围为 2.00$\sigma$（无信息先验）至 2.92$\sigma$（基于外部截面数据的物理信息先验）；在经验校准下，显著性为 2.36$\sigma$–3.16$\sigma$。
• 一致性： 两个结果在统计上是兼容的（分离度约为 $\approx 1.2\sigma$）。PRD 110 扫描中的值较高，归因于缺乏用于吸收归一化偏差的 sector-weight 自由度，从而迫使所有拓扑张力都进入了 $f_{mig}$。
• 物理阐释： 拟合的量级（侧带中约为 5%，直接扫描中约为 9%）与级联生成器（GENIE）与传输生成器（GiBUU）在近阈值质子方面的透明度差异（报告为 5%–10%）是一致的。

意义与主张
本文声称 $f_{mig}$ 是一个方法论诊断参数和有效的拓扑转移代理，而非对特定微观机制（如仅由平均场能量损失引起）的唯一测量。该参数吸收了包括 FSI、运动学建模误差和重建阈值在内的多种近阈值效应。

作者强调：

1. 结果取决于侧带拟合中所使用的先验模型以及 MicroBooNE 数据集的特定协方差结构。
2. 主要经验证据来自无先验的直接 CC 扫描，该扫描确立了对拓扑转移效应的强烈偏好（$>3\sigma$）。
3. 该参数尚未经过验证，不能直接用于 DUNE 振荡分析；它需要在重建层面通过独立的探测器级验证（例如使用 SBND 或 DUNE-ND）来确认迁移分数，并区分不同的物理起源。
4. 该参数作为一个候选的 Genie 拨盘，旨在解决 $0p$ 拓扑结构中存在的系统性缺失，这种缺失目前正迫使对 sector 归一化进行补偿性调整，进而可能导致振荡结果的偏差。

论文结论认为，尽管目前的数据尚无法分离出具体的物理起源（是残余平均场还是其他运动学效应），但 MicroBooNE 数据有力支持了存在 $\sim 5\text{--}10\%$ 拓扑转移效应的事实，这使得在未来的 LArTPC 分析中引入新的干扰参数成为必要。