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

技術要約：Genieベースのニュートリノ・アルゴン半包含解析における効果的なトポロジー転送によるスニセンス（不確実性）

問題提起
液体アルゴン時報投影チャンバー（LArTPC）の振動解析では、$\sim 35$ MeVの陽子再構成閾値に基づき、ニュートリノ相互作用サンプルを陽子タグ付き（$N_p$）および陽子なし（$0p$）のトポロジーに分割することが増えている。標準的なGENIEモンテカルロ・リウェイト・ライブラリには、この$35$MeVの境界を越えてイベントを転送する「ダイヤル」が欠けているという、決定的なギャップが存在する。GENIEの核内カスケード（INC）モデルは、一定の除去エネルギー（$W_s \approx 25$MeV）によって原子核の結合をモデル化しているが、輸送コード（例：GiBUU）に見られるような位置依存の平均場ポテンシャルを欠いている。この欠落により、閾値付近の陽子（$T_p \approx 35\text{--}100$MeV）に対するエネルギー損失が過小評価され、本来$0p$であるべきイベントを$N_p$と誤分類してしまう。この「陽子の可視性マイグレーション（移動）」は、$0p/N_p$比を歪ませ、エネルギー再構成戦略や振動系の系統誤差に直接的な影響を与える。既存のGENIEリウェイト・ダイヤル（例：MFP$_\pi$, MFP$_N$）は、あるトポロジー内では機能するが、閾値を越えてイベントを再分配することはできない。

手法
著者は、NominalなGENIE出力に対するイベントごとのリウェイトとして定義される、単一の経験的なスニセンス・パラメータ $f_{mig}$ を提案している。これは、実効的なトポロジー転送・代理変数として機能し、記述不明な近閾値の誤モデリング効果（陽子側の最終状態相互作用、核子運動量分布、RES/MECの運動学的形状、および残留平均場エネルギー損失の組み合わせ）を吸収するものである。

手法は以下の通りである：

1. 物理的動機： $^{40}$ArのWoods-Saxon密度プロファイルとWalecka型相対論的平均場ポテンシャル（GiBUU NL3*の深さ）を用いた半解析的な推定により、GENIEに欠落している残留エネルギー損失が $\approx 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. データ制約： パラメータは、2つの公開されたMicroBooNEデータセットを用いて制約される：
   • 主要な制約： 同時CC $0p/N_p$ 微分断面積リリース（PRD 110）に対する直接スキャン。この解析では、セクター重みによる正則化を行わないアンフォールドされた断面積を使用し、全体の正規化のみをプロファイリングしている。
   • 二次的な制約： NC $\Delta$ ラジエティブ・サイドバンド共分散（4チャネル、64ビン）へのフィット。これには、5つのセクター正規化重み（CCQE+MEC, CCRES, CCDIS, NCRES, NCDIS）とリッジ正則化（ガウス事前分布）、および2つのマイグレーション・フラクション（$f^{CC}_{mig}, f^{NC}_{mig}$）が含まれる。
4. 統計的処理： 有意性はプロファイル $\Delta\chi^2$ スキャンによって評価される。サイドバンド共分散の有効自由度が低い（$n_{eff} \approx 2.5$）こと、および物理的な境界（$f_{mig} \ge 0$）があるため、著者は境界を考慮したテスト統計量を用い、経験的なp値を較正するために5,000回の擬似実験（トイ・モンテカルロ）を用いて検証を行っている。

主な貢献

• 欠落しているスニセンスの特定： 本論文は、GENIEのリウェイト基底における特定のギャップ、すなわち陽子検出閾値付近でのトポロジー横断的な最終状態相互作用をモデル化できない点を特定している。
• $f_{mig}$ の提案： これは、イベントの再生成を伴わずにフラットなイベントごとの重みとして実装可能な、実効的なトポロジー転送・代理変数としてのメソッド論的診断パラメータである。これは、検証を条件として、将来の振動解析（例：DUNE）への適用が可能である。
• 診断的副産物： 本解析は、$f_{mig}$ を導入することで、マイグレーションを考慮しないフィットで見られたCCRESセクターの補償的な抑制が緩和されることを示している。$f_{mig}$ がない場合、フィットは$0p$の過剰を説明するためにCCRESを人工的に抑制する（$w_{CCRES} \approx 0.69$）。しかし、$f_{mig}$ を用いると、CCRESの重みは1へと緩和される（$w_{CCRES} \approx 1.07$）。
• 頑健性のチェック： 本論文には、事前分布への感度、共分散構造（Peelle's Pertinent Puzzle）、および既存のハドロン・ダイヤル（例：MFP$_N$）との直交性に関する広範なチェックが含まれている。

結果

• 主要な制約 (PRD 110)： 直接的なCCスキャンは、$\hat{f}_{mig} = 0.093 \pm 0.024$ という最良フィット・マイグレーション・フラクションを与えた。これは、$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$ となった。
• 一貫性： これら2つの結果は統計的に互換性がある（約 $1.2\sigma$ の分離）。PRD 110スキャンにおけるより高い値は、正規化の不一致を吸収するためのセクター重みの自由度が欠如しているため、トポロジーのテンションがすべて $f_{mig}$ に押し付けられたことに起因する。
• 物理的解釈： フィットされた大きさ（サイドバンドで $\sim 5\%$、直接スキャンで $\sim 9\%$）は、カスケード（GENIE）と輸送（GiBUU）のジェネレーター間で報告されている、閾値付近の陽子の透明性の違い（$\sim 5\text{--}10\%$）と一致している。

意義と主張
本論文は、$f_{mig}$ は特定の微視的なメカニズム（例：平均場エネルギー損失のみ）の固有の測定値ではなく、メソッド論的な診断パラメータであり、実効的なトポロジー転送・代理変数であると主張している。このパラメータは、FSI、運動学的モデリングの誤り、および再構成閾値を含む、近閾値の効果の組み合わせを吸収するものである。

著者は以下の点を強調している：

1. 結果は、サイドバンド・フィットで使用される事前モデルおよびMicroBooNEデータセットの特定の共分散構造に条件付けられている。
2. 主要な経験的証拠は、事前分布に依存しない直接CCスキャンから得られており、これはトポロジー転送効果に対する強い選好（$>3\sigma$）を確立している。
3. このパラメータは、まだDUNEの振動解析に直接使用するための検証は完了していない。再構成レベルでのマイグレーション・フラクションを確認し、競合する物理的起源を区別するために、独立した検出器レベルの検証（例：SBNDまたはDUNE-ND）が必要である。
4. このパラメータは、現在セクター正規化の補償的な調整を強いており、振動解析にバイアスを与える可能性がある $0p$ トポロジーの系統的な欠損に対処するための、候補となるGenieダイヤルとして機能する。

結論として、現在のデータでは特定の物理的起源（残留平均場か他の運動学的効果か）を孤立させることはできないが、$\sim 5\text{--}10\%$ のトポロジー転送効果が存在するという事実はMicroBooNEのデータによって強く支持されており、将来のLArTPC解析において新しいスニセンス・パラメータが必要であることを示唆している。