[论文解读] RADGEN 1.0. Monte Carlo Generator for Radiative Events in DIS on Polarized and Unpolarized Targets

Igor Akushevich, H. Boettcher+1 more|ArXiv.org|Jun 17, 1999

Radiation Detection and Scintillator Technologies被引用 14

一句话总结

RADGEN 1.0 是一种用于极化与非极化靶的深度非弹性散射（DIS）事件的蒙特卡罗生成器，通过两种形式（Mo-Tsai 和 Bardin-Shumeiko）实现 QED 辐射修正。它能够精确模拟弹性、准弹性及非弹性过程中的辐射尾部，经验证的结果显示对能量截断 Δ 的依赖性极小，约为束流能量的 0.1%，这对 HERMES 等精密实验至关重要。

ABSTRACT

A new Monte-Carlo generator including real radiated photons in DIS on polarized and unpolarized targets is presented. Analytical and numerical tests are performed and discussed in details.

研究动机与目标

• 开发一种蒙特卡罗生成器，以准确模拟在极化与非极化核靶上发生的深度非弹性散射（DIS）中的真实光子辐射。
• 实现并比较两种 QED 辐射修正形式——Mo-Tsai 与 Bardin-Shumeiko——用于包含式 DIS 过程。
• 通过引入运动学修正和光子能量阈值，确保对辐射尾部（弹性、准弹性、非弹性）的可靠模拟。
• 通过解析与数值比较验证生成器的性能，特别是在极化不对称性修正方面的应用。
• 为 HERMES 等精密实验提供工具，其中辐射修正显著影响截面测量结果。

提出的方法

• 生成器使用两个独立模块：FERRAD（非极化情况）与 POLRAD 2.0（极化情况），两者均基于 FORTRAN 代码。
• 通过光子发射后的真正四动量转移与能量转移重构运动学，进行如下替换：$ q \to q - k $，其中 $ k $ 为辐射光子的动量。
• 截面被划分为非辐射部分（包含圈修正与软光子修正）与辐射贡献部分：$ \sigma_{\text{obs}} = \sigma_{\text{non-rad}}(\Delta) + \sigma_{\text{in}} + \sigma_{\text{q}} + \sigma_{\text{el}} $。
• 截断参数 $ \Delta $ 定义软光子的阈值；若 $ E_\gamma < \Delta $，则事件被视为非辐射。
• 生成器分别计算非极化与极化情况下的辐射修正（RC）因子 $ \delta_{u,p} = \sigma^{\text{obs}}_{u,p} / \sigma^{\text{Born}}_{u,p} $。
• 采用查表法、样条插值或 FINT 插值，并对光子能量与方位角进行分箱，以高效采样微分截面。

实验结果

研究问题

• RQ1 真实光子发射的辐射修正如何影响极化与非极化靶上 DIS 的观测截面？
• RQ2 为使辐射修正对软光子阈值的依赖性最小化，能量截断 $ \Delta $ 的最优值是多少？
• RQ3 Mo-Tsai 与 Bardin-Shumeiko 形式在预测包含式 DIS 中辐射尾部贡献方面有何异同？
• RQ4 考虑到总截面修正量较大，辐射修正对极化不对称性的影响程度如何？
• RQ5 该生成器能否可靠模拟辐射光子的运动学分布及其对 $ x $、$ y $、$ Q^2 $ 与 $ W^2 $ 的影响？

主要发现

• 截断参数 $ \Delta $ 的最优值约为束流能量的 0.1%，此时结果对 $ \Delta $ 的依赖性可忽略，确保了结果的稳健性。
• 极化不对称性修正因子 $ \delta_p $ 显著小于非极化情况下的 $ \delta_u $，这是由于不对称性比值中因子化修正项发生了抵消。
• 生成器成功再现了 $ \delta_{u,p} $ 随 $ x $ 与 $ y $ 的变化行为，且在 $^3$He 中当 $ x \sim 0.4 $ 时，两者几乎相等。
• 极化不对称性修正因子主要由来自辐射尾部的非因子化修正主导，而非因子化部分。
• 对极化（POLRAD）与非极化（FERRAD）模块的数值比较显示良好一致性，验证了两种情况下的实现可靠性。
• 该生成器已成功集成至 HERMES 蒙特卡罗程序 HMC 中，经请求可提供使用。