[论文解读] MATCHIG: A program for matching charged Higgs boson production at hadron colliders

Johan Alwall|ArXiv.org|Mar 14, 2005

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

一句话总结

本文提出 MATCHIG，一种蒙特卡洛事件生成器，用于在强子对撞机中实现带电希格斯玻色子产生过程中 $2\to2$ $gb\to tH^{\pm}$ 与 $2\to3$ $gg\to tbH^{\pm}$ 过程之间的匹配方案。通过减去源自胶子分裂为 $b\bar{b}$ 的共线近似所导出的双重计数项，该方法消除了伴随 $b$ quark 在低横动量区域的过计数，从而在所有 $p_{\perp,b}$ 区域实现了精确的微分截面。

ABSTRACT

This manual describes how to use the MatCHig code for matching the charged Higgs boson production processes gg->tbH+/- and gb->tH+/-. A negative term, correcting for the double-counting between these processes, is implemented as an external process to PYTHIA, allowing the Monte Carlo generation of matched events. Results from the matching were published in JHEP 0412, 050. The code can be downloaded from http://www.isv.uu.se/thep/MC/matchig/ .

研究动机与目标

• 解决强子对撞机中带电希格斯玻色子产生过程中 $2\to2$ $gb\to tH^{\pm}$ 与 $2\to3$ $gg\to tbH^{\pm}$ 过程之间的双重计数问题。
• 提高伴随 $b$ quark 横动量较低区域的带电希格斯玻色子微分截面的计算精度。
• 通过 Pythia 的外部过程接口，实现部分子矩阵元与部分子喷注之间的稳定匹配方案。
• 为如 MSSM 等模型提供高精度模拟工具，支持更准确的实验信号预测。
• 通过引入负的双重计数项，校正共线近似中的过计数，确保事件生成的可靠性。

提出的方法

• 双重计数项 $\sigma_{\mathrm{DC}}$ 由胶子分裂 $g\to b\bar{b}$ 的 leading-order $b$-夸克密度导出，与 $2\to2$ 过程的共线贡献相匹配。
• 该双重计数项通过矩阵元 $\mathcal{M}_{2\to2}$、部分子分布函数以及分裂函数 $P_{qg}(z)$ 计算得出，同时考虑有限 $m_b$ 和 $\hat{s}$ 的运动学约束。
• 该方法通过 Les Houches 接口在 Pythia 中作为外部过程实现，过程 ID 为 LPRUP(1)=10000，事件通过 KFPR(MSTI(1),2) 进行标识。
• 采用减法方案 $\sigma = \sigma_{2\to2} + \sigma_{2\to3} - \sigma_{\mathrm{DC}}$，确保在低 $p_{\perp,b}$ 区域无双重计数。
• 该实现支持壳模型部分子喷注（MSTP(63)=0），用于部分子层级的比较，推荐设置 MSTP(32)=12 和 MSTP(39)=8 以实现一致的尺度选择。
• 因子化与重整化尺度通过 PARP(193) 和 PARP(194) 设置，推荐 $\mu_F \approx 0.5(m_{H^\pm} + m_t)/2$ 以获得最佳精度。

实验结果

研究问题

• RQ1 如何在伴随 $b$ quark 横动量较低时，一致地匹配 $2\to2$ 与 $2\to3$ 带电希格斯玻色子产生过程，以避免双重计数？
• RQ2 在 $2\to2$ 过程中，胶子分裂为 $b\bar{b}$ 的共线近似下，所产生双重计数项的正确解析形式为何？
• RQ3 如何通过外部过程与 Les Houches 接口，在 Pythia 蒙特卡洛框架中实现该匹配方案？
• RQ4 哪些尺度选择（因子化与重整化）能产生最精确的带电希格斯玻色子微分截面？
• RQ5 在低 $p_{\perp,b}$ 区域，与仅使用 $2\to2$ 过程相比，包含 $2\to3$ 过程并进行正确匹配是否能改善微分截面？

主要发现

• 双重计数项 $\sigma_{\mathrm{DC}}$ 被解析导出为 $g\to b\bar{b}$ 分裂产生的 $b$-夸克密度对 $2\to2$ 过程的 leading-order 贡献，确保与 $2\to3$ 矩阵元的一致性。
• 仅使用 $2\to3$ 过程在高 $p_{\perp,b}$ 区域会低估 $p_{\perp,b} \lesssim 100\,\mathrm{GeV}$ 区域的微分截面，因此必须与 $2\to2$ 过程匹配。
• 通过 Pythia 的外部过程接口实现双重计数项，可实现匹配事件生成，并正确实现干涉项抵消。
• 推荐设置 MSTP(32)=12 与 MSTP(39)=8，配合 PARP(193) 与 PARP(194)，可实现一致的因子化与重整化尺度选择，且 $\mu_F \approx 0.5(m_{H^\pm} + m_t)/2$ 为最优选择。
• 壳模型部分子喷注（MSTP(63)=0）对于准确的部子层级比较是必要的，尽管在喷注层级分析中预期无可观测差异。
• 该代码通过组合 $gb\to tH^{\pm}$（ISUB=161）与 $gg\to tbH^{\pm}$（ISUB=401）并加入负的双重计数项，成功生成了匹配事件，确保了总截面的正确性。