[论文解读] SUSY Les Houches Accord 2 I/O made easy

Thomas Hahn|ArXiv.org|May 4, 2006

Mathematics, Computing, and Information Processing被引用 7

一句话总结

本文提出一个 Fortran 77 库，简化了在超对称（SUSY）粒子物理计算中读写 SLHA2 格式文件的操作。该库使用单个双精度复数数组，通过预处理器宏访问，简化了数据处理，使与谱计算器和事件生成器的集成更加便捷，同时保持与 SLHALib 1 的向后兼容性，并支持 SLHA2 中的复数参数。

ABSTRACT

A library for reading and writing data in the SUSY Les Houches Accord 2 format is presented. The implementation is in native Fortran 77. The data are contained in a single array conveniently indexed by preprocessor statements.

研究动机与目标

• 简化高能物理应用中超对称 Les Houches 协议 2（SLHA2）格式数据的输入/输出。
• 提供一个轻量级、可移植的 Fortran 77 库，用于处理 SLHA2 文件的读写，而无需修改物理参数。
• 在支持 SLHA2 中引入的复数参数的同时，保持与 SLHALib 1 的向后兼容性。
• 通过使用预处理器定义的符号名称替代手动数组索引，降低编程复杂度，提升开发效率。
• 通过标准化、人类可读的 SLHA2 I/O，实现谱计算器与蒙特卡洛事件生成器之间的无缝集成。

提出的方法

• 该库使用一个名为 `slhadata` 的双精度复数数组作为所有 SLHA2 输入和输出的中心数据容器。
• 通过预处理器宏（在 `SLHA.h` 中定义）将符号名称映射到数组索引，避免用户记忆或计算偏移量。
• SLHA 块（如 SMInputs、MinPar）通过类似 `BlockSMInputs(i)` 的宏访问，其索引根据 SLHA2 规范自动计算。
• 衰变信息通过专用子程序管理：`SLHANewDecay`、`SLHAAddDecay`、`SLHAGetDecay` 和 `SLHADecayTable`，支持动态构建衰变表。
• 通过双精度复数数据类型支持复数参数，与仅支持实数的 `SLHALib 1` 形成区分。
• `SLHAInfo` 子程序用于添加元数据（如程序名称和版本），替代 `SLHAWrite` 中先前的硬编码参数，提升可扩展性。

实验结果

研究问题

• RQ1 如何使 Fortran 基础的 SUSY 软件开发者更容易使用且更少出错地进行 SLHA2 I/O 操作？
• RQ2 在可移植、向后兼容的 Fortran 77 库中，表示和访问 SLHA2 数据结构的最高效方式是什么？
• RQ3 如何在遗留的 Fortran 77 环境中本机支持 SLHA2 中的复数参数？
• RQ4 什么样的接口设计可以最小化将 SLHA2 集成到现有谱计算器和事件生成器工作流中的学习曲线和编码错误？
• RQ5 通过预处理器宏实现的单一数组抽象能否有效替代复杂数据结构，同时保持性能和正确性？

主要发现

• 该库通过使用预处理器宏，成功将 SLHA2 数据访问与底层数组索引解耦，显著减少编码错误并提高可维护性。
• 通过 `double complex slhadata` 实现从实数到复数数据的过渡，全面支持 SLHA2 中的复数参数，包括 MSSM 及其扩展中的参数。
• `SLHAExist` 函数现在返回一个整数，指示条目为实数、复数或无效，增强了类型安全性和错误检测能力。
• `SLHAInfo` 子程序允许灵活注入元数据（如程序名称、版本），而无需修改 `SLHAWrite` 接口，提升了可扩展性。
• 该库保持与 SLHALib 1 的向后兼容性，从 SLHA1 升级到 SLHA2 工作流时，仅需极少代码修改。
• 构建过程简单明了，支持常见 Fortran 编译器，使用 `./configure && make`，源代码可自由获取，遵循 GNU LGPL 许可，地址为 http://www.feynarts.de/slha。