宇宙射线和强子相互作用蒙特卡洛前端

此软件包提供了一个简单通用的用户界面，用于宇宙射线和高能粒子物理学中使用的流行事件生成器。通过去除复杂的Fortran风格接口、ASCII输入卡和C++依赖性，该软件包简化了粒子相互作用的模拟，使得更广泛的用户能够更容易、更快地访问。

用法

Python用户界面

要使用支持的事件生成器之一进行交互模拟，请导入该软件包并定义碰撞的参数。然后，创建事件生成器的一个实例，并生成事件。

import numpy as np
import chromo

# Define the parameters of the collisions
kinematics = chromo.kinematics.CenterOfMass(
    13 * chromo.constants.TeV,
    "proton", "proton")
# Create an instance of an event generator
generator = chromo.models.Sibyll23d(kinematics)

nevents = 0
average_pt = 0

# Generate 10000 events
for event in generator(10000):
    # Filter event
    event = event.final_state_charged()
    # do something with event.pid, event.eta, event.en, event.pt, etc.
    # these variables are numpy arrays, that can be histogrammed or counted like
    pt = event.pt[np.abs(event.pid) == 211]
    # The list could be empty
    if len(pt) > 0:
        nevents += 1
        average_pt += np.mean(pt)

average_pt = average_pt / nevents
print("Average pT for charged pions {0:4.3f}".format(average_pt))

更多示例，例如这个可以在示例文件夹中找到。

命令行界面（CLI）

安装chromo也会提供命令行界面。如果您的Python运行环境设置正确，您可以进行以下操作：

 chromo --help

查看命令行界面的帮助。如果不起作用，只需将 chromo 替换为 python -m chromo。命令行界面是为 CRMC 的用户设计的，感觉非常熟悉。CLI 可以以多种输出格式写入事件，具体如下。HepMC 输出可以管道输入到 RIVET 和许多支持该格式的其他工具。

输出格式

• HepMC（通过 pyhepmc，可选 gzip 压缩）
• ROOT（通过 uproot）
• 事件 SVG 图像（通过 pyhepmc 包）

支持的模式和引用方法

请注意，chromo 仅提供以下模式的用户界面，本身不包含任何粒子物理模型。在公开工作中使用这些模型时，请务必按照以下链接正确引用原始模型参考文献。另外，如果您在工作中发现 chromo 有用，我们非常感谢您的认可、脚注或链接到 chromo。

h = hadron，N = nucleon（p 或 n），A = nucleus，γ = photon，e = electron/positron

通过 PyPI 安装

支持的平台

• Python 3.8+
• Linux，Mac OS X（x86 和 M1/M2），Windows

推荐安装 chromo 的方法是使用预编译的二进制 wheel，这对于大多数常见架构和 Python 版本都可用

pip install chromo

高级和开发者安装说明可在此处找到：这里。

作者

• Anatoli Fedynitch
• Hans Dembinski
• Anton Prosekin
• Sonia El Hadri
• Keito Watanabe

授权

chromo 的源代码根据 BSD 3-clause 许可证授权（请参阅 LICENSE 以获取详细信息）。事件生成器的源代码根据不同的条件分别授权（请参阅子目录中位于 COPYING 文件中的条件）。