[论文解读] A simple alternative to the Crystal Ball function

S. Das|arXiv (Cornell University)|Mar 6, 2016

Particle physics theoretical and experimental studies参考文献 4被引用 7

一句话总结

本文提出了GaussExp和ExpGaussExp函数——作为高能物理中拟合具有指数尾部的峰形时，对Crystal Ball函数的更简单替代方案。通过用指数尾部替代幂律尾部，新函数降低了拟合不稳定性，减少了参数数量，并在模拟和真实数据中均表现出更优的收敛性和稳定性，尤其在大规模数据集或高尾参数情况下表现更佳。

ABSTRACT

We present a simple alternative to the Crystal Ball function that has an exponential tail stitched to a Gaussian core. It has one parameter less than the Crystal Ball function and, where appropriate, offers more stable fits to peaks that continue into exponential tails. The function may also be extended with two exponential tails on each side of the Gaussian, and this has two parameters less than the corresponding double-shouldered Crystal Ball function. This function has been used to model background and signal processes in a recent Higgs pair production search and may be of versatile use in experimental physics and other fields.

研究动机与目标

• 为解决Crystal Ball函数在拟合具有长指数尾部的分布时的不稳定性及高参数敏感性问题，特别是在幂律指数n取极端值时的不稳定性。
• 开发一种参数更少、连续可微（C¹）的函数，以保持高斯核心与指数尾部之间的平滑拼接。
• 在高统计量数据集中提升拟合收敛性和参数可靠性，因为Crystal Ball函数在n的幂律指数取极端值时常因数值不稳定性而失效。
• 为实验粒子物理中不变质量分布的建模提供一种稳健替代方案，例如在希格斯玻色子搜索中的应用。
• 在实际应用中展示新函数的实用性，包括在CMS希格斯粒子对产生搜索中对背景和信号的建模。

提出的方法

• GaussExp函数在高斯核心的一侧采用指数尾部，并在过渡点处强制保证函数值和一阶导数的连续性。
• 过渡点由单一衰减常数k定义，替代了Crystal Ball函数中的两个参数（α和n），使参数数量减少一个。
• 双峰型的ExpGaussExp函数在高斯核心的两侧均采用指数尾部，使用kL和kH作为衰减常数，与双峰型Crystal Ball函数相比，参数数量减少两个。
• 这些函数的构造方式确保在过渡点处，指数尾部与高斯部分在函数值和斜率上完全匹配，从而实现C¹连续性。
• 使用ROOT/MINUIT框架中的最大似然法进行拟合，并在模拟和真实数据上与Crystal Ball函数进行比较。
• 方法在指数衰减分布和模拟的H→γγ峰上进行了验证，评估了在不同事件数量下的收敛性和参数不确定性。

实验结果

研究问题

• RQ1 能否使用具有指数尾部的C¹函数替代Crystal Ball函数中的幂律尾部，同时降低拟合不稳定性？
• RQ2 所提出的GaussExp函数是否在参数数量更少的前提下，相比Crystal Ball函数实现更稳定的拟合，特别是在高统计量数据集中？
• RQ3 ExpGaussExp函数是否能比双峰型Crystal Ball函数更可靠地建模信号峰的非对称尾部？
• RQ4 在具有长指数尾部的大数据集上拟合时，新函数在收敛性和参数不确定性方面表现如何？
• RQ5 GaussExp和ExpGaussExp函数是否能在真实实验数据中成功应用，例如在希格斯玻色子质量重建中？

主要发现

• GaussExp函数在所有测试案例中（1,000至1,000,000事件）均成功收敛，而Crystal Ball函数在100,000和1,000,000事件的数据集中因指数n的不稳定性而无法收敛。
• 在100,000事件时，Crystal Ball函数返回n = 127 ± 40，表明参数约束极差；而GaussExp函数返回与真实值一致且不确定性小的参数。
• 在H→γγ信号峰拟合中，Crystal Ball函数未能收敛，返回nL = 100 ± 120和nH = 63 ± 131，参数极不稳定且物理上不合理。
• ExpGaussExp函数实现了接近1的χ²/ndf，并报告了较小的参数不确定性，而Crystal Ball函数的χ²/ndf更差且未能收敛。
• GaussExp和ExpGaussExp函数在所有测试数据集（包括大规模样本和高尾参数区域）中均表现出更优的收敛性和稳定性。
• 新函数已在近期CMS希格斯粒子对产生搜索中成功应用，分别用于背景（GaussExp）和信号（ExpGaussExp）不变质量分布的建模。