[论文解读] First Results with HIJING++ on High-energy Heavy Ion Collisions

Gábor Papp, G. G. Barnaföldi+7 more|Repository of the Academy's Library (Library of the Hungarian Academy of Sciences)|May 7, 2018

High-Energy Particle Collisions Research被引用 3

一句话总结

本论文介绍了 HIJING++ 3.1.1，这是一个基于 C++ 的全模块化、并行化高能重离子碰撞事件生成器，基于 PYTHIA8 构建，并集成了更新的核效应。通过线程级并行化，其性能提升了约 24 倍（非并行化部分占 4.2%），在 √sNN = 5.02 TeV 的 pp 和 pPb 碰撞中，与 ALICE 实验数据在带电强子谱和核子修改因子方面表现出良好一致性，但喷射淬火建模仍有改进空间。

ABSTRACT

We present preliminary results with HIJING++ (3.1.1) for identified hadron production in high-energy heavy ion collisions at LHC energies. The recently developed HIJING++ version is based on the latest version of PYTHIA8 and contains all the nuclear effects that have been included in the HIJING2.552, which will be improved by a new version of the shadowing parametrization and jet quenching module. Here, we summarize the structure and the speed gain due to parallelization of the new program code, also presenting some comparison between experimental data.

研究动机与目标

• 将原有的基于 FORTRAN 的 HIJING 事件生成器现代化，转化为一种可移植、模块化且支持并行计算的 C++ 框架，兼容现代 HEP 软件（如 AliRoot）。
• 利用最新的 PYTHIA8 和 LHAPDF6 库，实现更新的核效应（如 $Q^2$ 依赖的遮蔽效应和喷射淬火）
• 通过多线程并行化实现显著的计算加速，同时保持与原始 HIJING 2.552 物理精度相当
• 支持未来集成到大规模探测器模拟框架中，并实现 GPU 加速的随机数生成
• 为未来对撞机的下一代模拟提供基础，具备改进的物理建模能力和可扩展性

提出的方法

• 使用 PYTHIA8 的类结构作为基础，将原始 FORTRAN 的 HIJING 2.552 代码重写为面向对象的 C++ 代码
• 通过模块化类集成了核效应，包括遮蔽（带 $Q^2$ 依赖性）、Cronin 效应、喷射淬火和多重部分子相互作用
• 利用 HijQueue 和 HijManager 类实现基于线程的并行化，将事件分发到 CPU 核心上
• 通过 HijAnalysis 类支持用户自定义分析，配备可配置的直方图收集模块
• 通过将默认随机数生成器替换为 GPU 兼容版本，在初步测试中优化了性能
• 应用 Amdahl 定律量化加速效果，并识别出 4.2% 的非并行化部分（主要为队列管理）

实验结果

研究问题

• RQ1 能否在保留原始 HIJING 核心物理特性的前提下，从遗留的 FORTRAN 版本中构建出一个现代化、基于 C++ 且支持并行计算的事件生成器？
• RQ2 新版本 HIJING++ 在多大程度上能再现 √sNN = 5.02 TeV 的 pp 和 pPb 碰撞中带电强子谱的实验数据？
• RQ3 线程级并行化在减少计算时间方面的效率如何？其理论加速极限是多少？
• RQ4 为何当前的喷射淬火模型导致核子修改因子 R_{pPb} 过度抑制？如何改进？
• RQ5 模块化的 C++ 设计能否实现与 AliRoot 等大规模 HEP 模拟框架的无缝集成？”

主要发现

• 根据 Amdahl 定律拟合结果，HIJING++ 3.1.1 在 20 核系统上实现了最大理论加速约 24 倍，仅 4.2% 的代码为非并行化部分
• 在 √s = 5.02 TeV 的 pp 碰撞中，带电强子产额与 ALICE 数据相差在 50% 以内，且在中等和高 pT 区域，HIJING++ 的表现优于 PYTHIA8（取决于参数设置）
• 在 √sNN = 5.02 TeV 的 pPb 碰撞中，HIJING++ 的结果与 ALICE 数据在强子谱和核子修改因子方面相差在 30% 以内
• 核子修改因子 R_{pPb} 在 HIJING++ 中普遍被抑制，表明能量损失被高估，可能由于 N_bin 缩放相对于 N_part 的处理不准确
• 代码结构支持未来集成更先进的喷射淬火模型（如 GLV）以及 GPU 加速的随机数生成，初步测试显示性能有轻微提升
• 模块化的 C++ 设计使得代码更易于维护、可扩展，并与现代 HEP 软件栈（如 AliRoot）兼容，减少了代码膨胀和交叉依赖