[论文解读] 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)）的中性介子混合衰变时间积分，提出了分析解。推导了归一化、接受度修正和矩的表达式，并引入了高度优化的 C++ 程序，其速度和精度均优于以往的 CERNLIB 实现，现已集成至 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) 的高性能数值程序，该函数是计算这些积分的关键。
• 用更快、更精确的 C++ 实现替代旧版 CERNLIB 的 WWERF/CWERF 程序，集成至 RooFit 框架中。
• 实现实验中矩和接受度修正后衰变速率的精确计算，适用于高能物理实验。

提出的方法

• 推导理论衰变速率与高斯分辨率函数卷积的解析表达式，同时考虑时间依赖的混合与衰变效应。
• 利用 Faddeeva 函数 w(z) 表达涉及复误差函数的积分解，实现高效计算。
• 在复平面的不同区域采用级数展开与渐近逼近，确保数值稳定性。
• 采用分段算法：小 |z| 时使用泰勒级数，大 |z| 时使用渐近展开，并结合优化的查表机制。
• 在 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 框架中，成为高能物理软件栈中的标准工具。