[論文レビュー] PHOTOS Interface in C++; Technical and Physics Documentation

N. Davidson, T. Przedziński+1 more|arXiv (Cornell University)|Nov 3, 2010

Nuclear Physics and Applications被引用数 28

ひとこと要約

この論文は、PHOTOSモンテカルロジェネレータをFORTRANからC++に移行したものを提示しており、HepMCイベントレコードへのC++ネイティブインターフェースを新たに実装することで、保守性・拡張性・物理的精度の向上を実現している。再実装はHEPEVTとの完全な互換性を維持しており、プロセス依存の行列要素やイベントレコード構造のより良い制御といった高度な機能もサポートしており、行列要素補正を有効にしたLHC $Z$および$W$崩壊解析では、系統的不確実性が0.3%未満にまで低下している。

ABSTRACT

For five years now, PHOTOS Monte Carlo for bremsstrahlung in the decay of particles and resonances has been available with an interface to the C++ HepMC event record. The main purpose of the present paper is to document the technical aspects of the PHOTOS Monte Carlo installation and present version use. A multitude of test results and examples are distributed together with the program code. The PHOTOS C++ physics precision is better than its FORTRAN predecessor and more convenient steering options are also available. An algorithm for the event record interface necessary for process dependent photon emission kernel is implemented. It is used in Z and W decays for kernels of complete first order matrix elements of the decays. Additional emission of final state lepton pairs is also available. Physics assumptions used in the program and properties of the solution are reviewed. In particular, it is explained how the second order matrix elements were used in design and validation of the program iteration procedure. Also, it is explained that the phase space parametrization used in the program is exact.

研究の動機と目的

• FORTRANからC++へのイベントレコードインターフェースの再実装を通じて、PHOTOSモンテカルロジェネレータの近代化を図り、現代のヒッグズ・エレクトロウィーク物理ソフトウェアとの統合性を向上させること。
• レガシーフォーマットのHEPEVTに代わる、C++におけるHepMCイベントレコードの完全なサポートを実現し、現在のヒッグズ・エレクトロウィーク物理ソフトウェアの標準に合わせること。
• 崩壊過程、特に多光子およびレプトン対の放射において、スピンおよび行列要素情報への直接アクセスを提供することで、物理的精度を向上させること。
• 異常カップリングや検出器の受容率効果を含む、$W$および$Z$ボソン崩壊における系統的誤差の高精度な推定を可能にすること。
• 従来のHEPEVTとの後方互換性を保ちつつ、将来の物理的拡張を想定した拡張可能なC++クラスを導入すること。

提案手法

• オブジェクト指向設計の原則を用いて、コードの明確性と保守性を向上させるために、PHOTOSのコアアルゴリズムをC++に再実装した。
• HepMCイベントレコードを解析・操作する新しいインターフェース層を設計し、母粒子-娘粒子の崩壊チェーンおよび運動量情報の抽出を実現した。
• 崩壊を崩壊粒子の静止系に回転させるフレーム変換システムを実装し、量子化軸を母粒子の運動量に一致させることで、正確な行列要素計算を可能にした。
• プロセス依存の一次元QED行列要素を内部アルゴリズムに統合し、最終状態放射の正確な取り扱いを可能にした。
• 残りのFORTRANコアへのC++ラッパーを介してHEPEVTとの後方互換性を確保し、既存のワークフローに影響を与えないようにした。
• 強制処理、ログ出力、およびユーザー定義メソッド（例：processParticleおよびprocessBranch）を用いたカスタム解析を可能にする構成オプションを提供した。

実験結果

リサーチクエスチョン

• RQ1 HepMCのような現代的なC++イベントレコード標準をサポートするには、PHOTOSモンテカルロジェネレータをどのように近代化できるか？
• RQ2 FORTRANからC++への移行において、物理的精度を完全に維持するにはどのような変更が必要か？
• RQ3 イベントレコードインターフェースをどのように拡張すれば、プロセス依存の行列要素およびスピン情報へのアクセスを可能にし、精度を向上させられるか？
• RQ4 高エネルギー物理学の応用において、新しいC++ベースのアーキテクチャがもたらす性能および精度の向上はどのようなものか？
• RQ5 新しいインターフェースは、$W$および$Z$ボソン崩壊における系統的誤差推定といった高度な物理的研究をどのように支援できるか？

主な発見

• C++ベースのPHOTOS v3.60はHepMCイベントレコードと完全に互換性があり、現代のヒッグズ・エレクトロウィーク物理ソフトウェアスタックへのシームレスな統合を可能にした。
• 移行によりスピンおよび行列要素情報への完全なアクセスが可能となり、多光子およびレプトン対放射におけるプロセス依存QED補正の正確な取り扱いが実現した。
• 強制処理、ログ出力、およびprocessParticleおよびprocessBranchメソッドを用いたカスタム解析を含む、高度なユーザー設定が可能になった。
• ZおよびW崩壊シミュレーションにおける系統的不確実性は0.3%にまで低下し、行列要素補正を有効にした場合には0.2%まで低下した。これはLHC用途における高い精度を裏付けるものである。
• FORTRANコアへのC++ラッパーによりHEPEVTインターフェースが維持されており、パフォーマンスや機能を損なうことなく後方互換性が保証された。
• 新しい設計により、将来の物理的拡張（例：TAUOLAのようなインターフェースを通じた弱い補正との統合可能性を含む）が可能になった。