[论文解读] ManeParse: Mathematica Toolbox for PDF Uncertainties and Application to New Physics Searches

E. Godat|arXiv (Cornell University)|Oct 20, 2015

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

一句话总结

ManeParse 是一个基于 Mathematica 的工具箱，可实现对多个合作组（包括 LHAPDF6 和 CTEQ PDS 格式）部分子分布函数（PDFs）的高精度、用户友好的分析。它支持插值、误差传播（Hessian 和蒙特卡洛方法）以及亮度计算，通过在本征底夸克分析和超高能中微子的低-x 外推等应用中展示了其在新物理搜索中的实用性。

ABSTRACT

As the LHC begins Run 2 at an even higher energy, one of the top priorities will be to search for new particles (possibly from SUSY) at the highest energy scales. In addition to direct production of new particles, they can mix with Standard Model (SM) particles to yield discrepancies from the usual predictions. To distinguish this new physics from old uncertainties, we need tools to precisely quantify the uncertainties of the SM predictions. Here we present a versatile set of utility functions built with Mathematica that is capable of working with a variety of PDF formats including the recent LHAPDF6 format. This software can perform PDF calculations within the Mathematica framework and compare results from different PDF collaborations. The package includes both the central PDF value as well as the full error sets needed for PDF uncertainty analysis; a variety of sample error definitions are implemented. We demonstrate this package for the case of a new heavy scalar particle production at the LHC.

研究动机与目标

• 开发一种轻量级、易用的工具，用于在高级计算环境中读取和分析来自多个合作组的 PDF，包括 LHAPDF6 和 CTEQ PDS 格式。
• 通过支持 Hessian 和蒙特卡洛误差传播技术，实现标准模型预测中不确定性的精确量化。
• 通过在 Mathematica 环境中实现亮度和截面的精确计算，促进新物理搜索。
• 提供一种快速、可靠的 PDF 网格插值方法，采用自定义的四点拉格朗日插值例程。
• 支持向低-x 区域的外推，适用于超高能中微子物理等应用，该区域实验数据稀疏。

提出的方法

• 该工具包使用自定义的四点拉格朗日插值例程，将离散的 PDF 网格视为连续函数，实现在任意 x 和 Q 值处的平滑评估。
• 它将 LHAPDF6 和 CTEQ PDS 格式的 PDF 数据文件解析为 x、Q 和部分子种类的三维网格，并将数据存储在内存中以实现快速访问。
• 插值函数集成于 Mathematica 的原生环境，可直接与内置的绘图、积分和符号计算工具协同使用。
• 误差传播通过 Hessian 和蒙特卡洛两种方法实现，完整的误差集被存储并可供访问，以支持不确定性量化。
• 亮度计算通过数值积分卷积公式实现：$ \frac{d\mathcal{L}_{ij}}{d\tau} = \frac{1}{\delta_{ij}+1} \int_{\tau}^{1} \frac{1}{x} \left[ f_i(x) f_j(\tau/x) + f_j(x) f_i(\tau/x) \right] dx $。
• 该工具包包含一个演示笔记本，提供用于绘图、插值和误差带可视化的示例程序，提升了新用户的易用性。

实验结果

研究问题

• RQ1 如何在高级计算环境中开发一种轻量级、用户友好的工具，用于读取和分析来自多个合作组的 PDF？
• RQ2 当应用于不同合作组的多样化 PDF 网格时，自定义插值例程在保持精度和速度方面能达到何种程度？
• RQ3 ManeParse 是否能有效重现已知结果（如亮度计算），并支持不同 PDF 集之间的直接比较？
• RQ4 该工具包在实验数据有限的低-x 区域对 PDF 进行外推的能力如何？这对高能物理应用有何影响？
• RQ5 ManeParse 是否能够以一种可为新物理搜索提供信息的方式，区分本征底夸克贡献与标准 PDF 不确定性？

主要发现

• ManeParse 在数秒内即可解析完整的 PDF 家族，单个文件加载仅需零点几秒，性能与基于 C++ 和 FORTRAN 的 LHAPDF6 及 CTEQ 代码相当。
• 插值例程每秒可执行约 2,000 次调用，支持在 PDF 网格上高效进行数值积分和绘图。
• PDF 的求和规则成功得到验证，确认在给定 Q 值下，插值后的 PDF 在所有部分子种类上总和为一，验证了插值的准确性。
• 该工具包成功重现了先前工作的亮度计算结果，并支持在不同 PDF 集之间进行直接比较，包括与自定义本征底夸克 PDF 的比较。
• ManeParse 通过可调的 x 幂律外推方法，可靠地实现了向低-x 区域的扩展，支持超高能中微子相互作用的研究。
• ManeParse 与 Mathematica 的绘图和计算工具集成，可生成复杂可视化效果，包括三维相空间图和带误差带的比值图。