[論文レビュー] Decay time integrals in neutral meson mixing and their efficient evaluation

T. M. Karbach, G. Raven+1 more|arXiv (Cornell University)|Jul 3, 2014

Particle physics theoretical and experimental studies参考文献 1被引用数 26

ひとこと要約

本稿では、複素誤差関数、特にFaddeeva関数 w(z) を含む中性メソン混合における崩壊時間積分の解析的解を提示する。正規化、受容率補正、モーメントのための式を導出し、従来のCERNLIB実装よりも高速かつ高精度な最適化されたC++ルーチンを新たに開発した。これらのルーチンは、ROOT 5.34/08以降、RooFitに統合されている。

ABSTRACT

In neutral meson mixing, a certain class of convolution integrals is required whose solution involves the error function $\mathrm{erf}(z)$ of a complex argument $z$. We show the the general shape of the analytic solution of these integrals, and give expressions which allow the normalisation of these expressions for use in probability density functions. Furthermore, we derive expressions which allow a (decay time) acceptance to be included in these integrals, or allow the calculation of moments. We also describe the implementation of numerical routines which allow the numerical evaluation of $w(z)=e^{-z^2}(1-\mathrm{erf}(-iz))$, sometimes also called Faddeeva function, in C++. These new routines improve over the old CERNLIB routine(s) WWERF/CWERF in terms of both speed and accuracy. These new routines are part of the RooFit package, and have been distributed with it since ROOT version 5.34/08.

研究の動機と目的

• 有限の崩壊時間分解能を有する時間に依存する中性メソン混合において生じる畳み込み積分の解析的解を導出すること。
• 受容率補正および分解能効果を含む実験的崩壊時間分布の正規化された確率密度関数を提供すること。
• これらの積分を評価するために不可欠なFaddeeva関数 w(z) = e^(-z²) erfc(-iz) の高性能な数値ルーチンを開発すること。
• 従来のCERNLIB WWERF/CWERFルーチンに代わる、より高速かつ高精度なC++実装をRooFitフレームワークに統合すること。
• 高エネルギー物理学実験におけるモーメントおよび受容率補正付き崩壊率の正確な計算を可能にすること。

提案手法

• 理論的崩壊率とガウス分解能関数の畳み込みに対して解析的表現を導出する。時間に依存する混合と崩壊を組み込む。
• 複素誤差関数を含む積分の解を表現するためにFaddeeva関数 w(z) を用い、計算を効率化する。
• 複素平面の異なる領域において級数展開と漸近近似を適用し、数値的安定性を確保する。
• |z| が小さい領域ではテイラー級数、|z| が大きい領域では漸近展開を用いたピecewiseアルゴリズムを実装し、最適化されたルックアップテーブルを併用する。
• RooFitにおけるC++テンプレートベースの実装を導入し、doubleおよびlong double精度を両方サポート。高速かつ高精度な評価を実現。
• さまざまな複素引数に対してベンチマークテストを実施し、CERNLIBルーチンと比較して速度と精度を検証する。

実験結果

リサーチクエスチョン

• RQ1 ガウス分解能関数と畳み込まれた場合の、中性メソン混合における崩壊時間積分は、どのように解析的に解けるか？
• RQ2 時間に依存する混合および分解能効果を考慮した場合、正規化、受容率補正、モーメントの閉形式表現は何か？
• RQ3 高エネルギー物理学の応用において、C++でFaddeeva関数 w(z) を効率的かつ高精度に評価するにはどうすればよいか？
• RQ4 レガシCERNLIBのWWERF/CWERFルーチンと比較して、速度と精度の面でどの程度の性能向上が達成できるか？
• RQ5 確率密度関数への応用を想定した複素誤差関数の数値評価を、複素平面全体にわたり安定化させるにはどうすればよいか？

主な発見

• 本稿では、分解能および受容率効果を含む中性メソン混合における崩壊率分布の正規化について、正確な解析的表現を導出した。
• Faddeeva関数 w(z) が、複素平面上でのこれらの積分評価を可能にする中心的な数学的道具であることが示された。
• 新規C++実装は顕著な性能向上を達成しており、一般的な用途ではCERNLIBのWWERF/CWERFと比較して2–3倍高速である。
• 特に動的範囲が広い領域や、漸近展開と級数展開の遷移領域付近において、数値的精度が向上している。
• 複素平面全体にわたり安定した実装であり、適切な引数還元により対称性および分岐カットを正しく扱っている。
• これらのルーチンは、ROOTバージョン5.34/08以降、RooFitパッケージに正常に統合され、高エネルギー物理学のソフトウェアスタックにおける標準的なツールとなっている。