PHOKHARA

PHOKHARA
Name	PHOKHARA
Title	PHOKHARA
Developer	H. Czyż, J.H. Kühn, G. Rodrigo, M. Szopa, T. Czyż, A. B. Arbuzov
Released	2000
Latest release version	10.0
Programming language	Fortran
Operating system	Unix-like, Linux
Genre	Monte Carlo event generator
License	Proprietary / academic

PHOKHARA is a precision Monte Carlo event generator developed for radiative return studies and precision measurements in high-energy particle physics. It provides next-to-leading order and higher-order corrections for processes in electron-positron annihilation relevant to precision determinations of the anomalous magnetic moment of the muon, the running of the electromagnetic coupling constant, and hadronic cross sections measured at facilities such as DAΦNE, KEK, SLAC, and LEP. PHOKHARA interfaces with experimental analyses performed by collaborations including KLOE, BaBar, Belle, CMD-2, and SND.

Overview

PHOKHARA is designed to simulate radiative processes like e+e- → hadrons + nγ with high theoretical control, targeting measurements connected to the hadronic vacuum polarization contribution to the muon g−2 and determinations of the hadronic R-ratio used in global fits including those by the Particle Data Group. It supports channels with two-pion, three-pion, four-pion, kaon, and baryonic final states relevant for experiments at VEPP-2M, VEPP-2000, BEPCII, and high-luminosity B-factories. The code has been benchmarked against generators such as KKMC, BabaYaga, and MCGPJ to ensure compatibility with analyses from CMD-3 and BaBar.

Theoretical Foundations

PHOKHARA implements quantum electrodynamics radiative corrections at next-to-leading order and beyond, combining exact matrix elements with factorization approaches used in precision calculations by groups associated with CERN, DESY, and INFN. The generator incorporates form factor parametrizations informed by dispersive treatments developed in studies related to the Vector Meson Dominance approach and constraints from measurements by NA7, CELLO, TASSO, and ALEPH. Treatment of initial-state radiation and final-state radiation follows techniques validated against predictions from QED calculations used in analyses from LEP collaborations such as ALEPH, DELPHI, L3, and OPAL to control infrared singularities and soft-photon resummation used in comparisons with ResBos-type formalisms.

Implementation and Features

The codebase is written primarily in Fortran and structured to provide modular implementations of channels, enabling deployment at computing centers including CERN and university clusters associated with MIT, University of Rome La Sapienza, and University of Warsaw. Features include event generation with exact NLO matrix elements, optional structure-function approximations, multi-photon emission with exponentiation, and user-definable acceptance cuts matching detectors like KLOE-2, Belle II, BaBar, and CLEO. PHOKHARA provides differential distributions such as invariant mass spectra and angular observables used in unfolding procedures employed by collaborations like BESIII and NA62. Interfaces allow comparison with fitting frameworks used by collaborations including PDG global fits and phenomenology tools developed at IHEP, JINR, and SLAC National Accelerator Laboratory.

Validation and Comparisons

PHOKHARA's predictions have been cross-checked against experimental data from KLOE, BaBar, Belle, CMD-2, and SND, and validated through tuned comparisons with Monte Carlo generators KKMC and BabaYaga@NLO. Benchmark studies performed in workshops at La Thuile and Radcor meetings compared PHOKHARA outputs to analytic calculations used by theory groups at MPI for Physics and Max Planck Institute for Physics. Validation campaigns addressed radiative return observables critical to the evaluations by Muon g−2 Theory Initiative and global electroweak fits led by groups at CERN and SLAC.

Applications in Experimental Physics

PHOKHARA has been used to extract the e+e− → π+π− cross section employed by KLOE to constrain the hadronic contribution to muon g−2 and to provide templates for ISR-based measurements by BaBar and Belle. It supports studies of resonance parameters for states such as the ρ(770), ω(782), φ(1020), and higher vector resonances analyzed by CMD-3 and SND. Experiments at DAΦNE and VEPP-2000 utilize PHOKHARA for luminosity-normalized measurements and systematic uncertainty estimates, while B-factory analyses at SLAC and KEK exploit its multi-photon emissions modeling to study radiative returns to low-energy hadronic states and backgrounds relevant for searches performed by Belle II and BaBar.

Development History and Contributors

PHOKHARA development began around 2000 and evolved through contributions from theorists and experimentalists affiliated with institutions such as University of Silesia, University of Kraków, University of Zurich, CERN, INFN, and DESY. Key contributors include authors who published foundational papers and updates presented at conferences organized by EPS, ICHEP, and the European Physical Society meetings. Ongoing maintenance and channel additions have been coordinated with experimental collaborations KLOE-2, BaBar, and Belle to ensure relevance for current precision programs at Fermilab and J-PARC initiatives. The project continues to influence precision studies feeding into global efforts by the Muon g−2 Theory Initiative and phenomenology networks across Europe and North America.

Category:Monte Carlo event generators