[論文レビュー] PYTHIA 6.3 Physics and Manual

Torbjörn Sjöstrand, Leif Lönnblad+2 more|ArXiv.org|Aug 14, 2003

Particle physics theoretical and experimental studies被引用数 168

ひとこと要約

PYTHIA 6.3 は、量子色力学（QCD）プロセス、パートンシャワー、ハドロン化、およびボーズ・アインシュタイン相関を含む量子干渉効果を詳細にモデル化する、高エネルギー物理学のための包括的なモンテカルロイベントジェネレータである。Bose-Einstein相関を処理するための3つの異なるアルゴリズム（BE₀、BE_λ、BE_m）を実装しており、それぞれ運動量保存と対相関の精度を最適化している。BE_m は、同一粒子対の近接を抑制し、4粒子補正を最適化することで、入力の f₂(Q) 形状への一致が向上している。

ABSTRACT

The PYTHIA program can be used to generate high-energy-physics `events', i.e. sets of outgoing particles produced in the interactions between two incoming particles. The objective is to provide as accurate as possible a representation of event properties in a wide range of reactions, with emphasis on those where strong interactions play a role, directly or indirectly, and therefore multihadronic final states are produced. The physics is then not understood well enough to give an exact description; instead the program has to be based on a combination of analytical results and various QCD-based models. This physics input is summarized here, for areas such as hard subprocesses, initial- and final-state parton showers, beam remnants and underlying events, fragmentation and decays, and much more. Furthermore, extensive information is provided on all program elements: subroutines and functions, switches and parameters, and particle and process data. This should allow the user to tailor the generation task to the topics of interest.

研究の動機と目的

• 高エネルギー粒子衝突を、イベントレベルの詳細を完全に含む統合的で正確かつ拡張可能なフレームワークとしてシミュレートすること。
• マルチハドロン最終状態における、パートンシャワーリング、ビーム残渣、多重相互作用、ハドロン化などの非摂動的QCD効果をモデル化すること。
• 特に同一粒子対におけるボーズ・アインシュタイン相関を扱うための複数のアルゴリズムを実装・比較すること。
• λ や R といったパラメータのチューニングと、BE補正におけるエネルギー運動量保存の維持を通じて、実験的観測量との整合性を確保すること。
• 特に e⁺e⁻、pp、ep 衝突において、完全なドキュメンテーションと設定可能성을備えた、現実的なイベントサンプルの生成を研究者に支援すること。

提案手法

• ハード散乱過程、パートンシャワー、ハドロン化をシミュレートするために、解析的QCD計算と現象論的モデルの組み合わせを用いる。
• ハドロン化には、カラーテンソル状態のストリング形成と、ハドロンへの逐次崩壊を特徴とする Lund ストリング断片化モデルを採用する。
• ボーズ・アインシュタイン相関による運動量シフトを補正するために、3つの異なるアルゴリズム（BE₀、BE_λ、BE_m）を適用し、エネルギー運動量保存を確保する。
• BE_m では、相関のある粒子対 (i,j) の近くに配置された補正粒子対 (k,l) を、近接度をペナルティとする重み付き測度 W_ijkl を用いて選択する。このペナルティは、(1−exp(−Q²R²)) を用いて同一粒子に近い場合の寄与を抑制する。
• 確率的微分断面積に従ってイベントをサンプリングするためのボイドアルゴリズムを乱数生成とともに用い、統計的正確性を保証する。
• HEPEVT標準と互換性のある階層的イベント記録構造を実装し、モジュラリティと拡張性を実現する。

実験結果

リサーチクエスチョン

• RQ1 エネルギー運動量保存を保ちながら、イベントジェネレータにおけるボーズ・アインシュタイン相関を正確にモデル化する方法は何か？
• RQ2 ハドロン最終状態における同一粒子相関をシミュレートする際、BE₀、BE_λ、BE_m といった異なるアルゴリズムの相対的利点と制限は何か？
• RQ3 同一粒子の近接が、BE相関モデルの精度に与える影響は何か？ また、その影響をどのように軽減できるか？
• RQ4 複数の粒子対による異なる運動量シフトが、最終的な相関関数に与える影響は何か？ そして、その補正はどのように行われるか？
• RQ5 補償粒子対を用いて運動量保存を導入する際、BE相関形状 f₂(Q) をどのようにして維持できるか？

主な発見

• BE_m アルゴリズムは、同一パートナーに近い粒子の寄与を抑制する重み付き測度 W_ijkl を最大化することで、入力の f₂(Q) 形状への一致が最良となる。
• BE_m アルゴリズムにおける平均スケーリング係数 α は 0.73 であり、複数の運動量シフトが有効相関強度を低下させるため、補正が必要であることを示している。
• BE₀ アルゴリズムは、入力の f₂(Q) よりも低い観測二粒子相関を生じさせるため、データとの一致を得るためには λ パラメータのチューニングが必要である。
• BE_λ アルゴリズムは、Lund モデルにおけるストリング長にインspiredされており、ストリング断片化に内在するカラーフラグメント接続性のおかげで、相関構造が改善されている。
• 抑制因子 (1−exp(−Q²R²)) は、同一粒子に近い領域での物理的でない運動量シフトを効果的に防止し、安定性と正確性を向上させている。
• 特に密度の高い最終状態において、BE_m アルゴリズムは BE₀ や BE_λ よりも、BE相関関数のピーク高さと幅をよりよく保持している。