HadroTOPS: A Monte Carlo Event Generator For Hadron Production In Two-Photon Scattering In Electron Positron Collisions

The paper introduces HadroTOPS, a flexible Monte Carlo event generator designed to simulate hadronic two-photon fusion events in electron-positron collisions with exact leading-order QED coupling and customizable phase space decays, specifically supporting partial wave analyses for various final states relevant to collaborations like BESIII.

The Gist

Imagine a high-energy particle collider as a massive, ultra-fast billiard table where electrons and positrons (the antimatter twins of electrons) zoom around and crash into each other. Usually, when they collide, they smash together to create new particles directly. But sometimes, they don't hit each other head-on. Instead, they pass close by, and like two people waving at each other from passing trains, they exchange "invisible waves" of energy called photons.

These two photons then crash into each other and, surprisingly, turn into a burst of new matter—usually groups of particles called hadrons (like pions or kaons). This is called "two-photon scattering."

The Problem:
Physicists want to study these rare, messy bursts of new particles to understand the fundamental building blocks of the universe. To do this, they need a way to predict what should happen if their theories are correct. They need a "simulator" or a "digital twin" of the experiment.

The problem is that existing simulators were like old, broken maps. They were either too simple (ignoring the complex rules of how the particles interact) or too rigid (only able to simulate specific, simple crashes). They couldn't handle the messy, multi-particle explosions that happen in these photon collisions, making it hard to compare real data with theory.

The Solution: HadroTOPS
The authors of this paper have built a new, highly flexible computer program called HadroTOPS. Think of it as a state-of-the-art "virtual crash test" for these specific types of particle collisions.

Here is how it works, using simple analogies:

1. The Perfect Blueprint (QED):
   The program starts with the "perfect blueprint" of the collision. It calculates the exact rules of how the electrons and positrons emit their photon waves using the most precise math available (Quantum Electrodynamics). It's like knowing exactly how hard the billiard balls are moving before they even touch.

2. The "Flat" Canvas (Phase Space):
   Once the photons collide, they create a cloud of new particles. HadroTOPS generates these particles in a way that is "flat" or "uniform." Imagine a painter throwing paint onto a canvas. Instead of painting a specific picture (like a flower or a face), the program throws the paint evenly everywhere.
   Why is this useful? Because if the computer paints a specific picture, it might accidentally hide the real physics. By painting everything evenly, physicists can later look at the "paint splatters" (the real experimental data) and see exactly where the physics deviates from the flat background. This is crucial for "Partial Wave Analysis," a technique used to figure out the shape and spin of the particles created.

3. The Recipe Book (Theoretical Inputs):
   The program isn't just a random generator; it has a recipe book. For some common particle combinations (like two pions), it uses advanced mathematical recipes called "dispersive formalism" that account for how these particles bounce off each other and form resonances (temporary, unstable particles). For other combinations (like kaons), it uses real-world data from previous experiments as a guide.

4. The Flexible Chameleon:
   One of the best features of HadroTOPS is its flexibility.

   • Customizable: If a physicist wants to study a new type of particle collision, they can easily add it to the program.
   • Adjustable: It can simulate collisions where the photons are almost real (quasi-real) or highly virtual (carrying a lot of energy), just like real experiments at facilities like BESIII or Belle II.
   • Multi-Particle: It can handle simple crashes (two particles) and complex explosions (three or more particles) with the same ease.

What It Can Do Right Now:
The paper states that the program currently includes specific "recipes" for creating:

• Pion pairs ($\pi^+\pi^-$, $\pi^0\pi^0$)
• Mixed pairs like $\pi^0\eta$
• Kaon pairs ($K^+K^-$, $K^0_S K^0_S$)
• Eta pairs ($\eta\eta$)
• A specific complex decay chain involving the $f_1(1285)$ particle turning into an eta and two pions.

Why It Matters:
This tool allows physicists to run millions of "virtual experiments" on their computers. They can compare the results of these virtual crashes with the actual data coming from real particle colliders. If the real data looks different from the HadroTOPS simulation, it tells the physicists that there is something new or interesting happening in nature that their current theories haven't explained yet.

In short, HadroTOPS is a versatile, high-precision digital simulator that helps physicists decode the messy, beautiful explosions of matter created when light collides with light.

Technical Summary

Technical Summary of HadroTOPS: A Monte Carlo Event Generator for Hadron Production in Two-Photon Scattering

Problem Statement

The analysis of hadronic two-photon interactions at electron-positron ($e^+e^-$) colliders requires Monte Carlo (MC) generators capable of modeling complex hadronic processes across a broad kinematic spectrum. A specific limitation exists in the study of multi-hadron production, where suitable simulation tools are lacking. Partial wave analyses (PWAs) of various hadronic final states necessitate MC simulations that precisely predict the quantum electrodynamics (QED) components of the process while simultaneously generating a uniform phase space distribution for two-photon to multi-hadron transitions. Existing generators, such as Ekhara3.2, have limitations: while Ekhara3.2 includes full dispersive amplitudes for $\pi\pi$ final states, it is restricted to specific channels and does not support the flexible, flat phase-space generation required for general PWA tools across diverse final states. Furthermore, many available generators possess different kinematic limitations that prevent their use in all experimental configurations, particularly those involving single-tag or double-tag modes where photon virtualities ($Q^2$) are substantial.

Methodology

The authors present HadroTOPS, a C++ Monte Carlo event generator designed to address these gaps. The methodology integrates exact leading-order QED coupling with a flexible phase-space generation algorithm.

Theoretical Framework

The generator models the process $e^+e^- \to e^+e^-X$, where $X$ is a hadronic final state formed via $\gamma^*\gamma^*$ fusion.

• QED Calculation: The code extends the leading-order QED calculation for inclusive and single-meson production to fully exclusive two-meson ($M_1M_2$) production. It utilizes a helicity amplitude decomposition of the hadronic tensor $H_{\mu\nu}$ into helicity amplitudes $H_{\lambda_1\lambda_2}$ for the fusion of two virtual photons.
• Cross Section Formulation: The unpolarized cross section is expressed in terms of 25 helicity-dependent cross sections and interference response functions (e.g., $\sigma_0, \sigma_2, \sigma_{TL}, \tau_{TT}$, etc.). The formulation explicitly retains the full kinematic information, including polar and azimuthal angular distributions of the final state particles, which is crucial for differential analysis.
• Phase Space Generation: The generator employs a highly efficient algorithm based on Schuler's method (used in GALUGA2.0 and Ekhara). It generates final-state vectors using five Lorentz invariants ($W^2, t_1, t_2, s_1, s_2$) mapped from a six-dimensional unit hypercube.
  • Logarithmic Mapping: To handle the $(t_1 t_2)^{-1}$ dependence in the cross section and ensure numerical stability in the region of small photon virtualities, the code implements a logarithmic mapping for momentum transfers. This cancels the pole at $t_1 t_2 = 0$ and mimics the cross-section behavior, significantly improving generation efficiency.
  • Flat Phase Space: The hadronic state decays into an arbitrary number of final-state particles with a flat phase space distribution, allowing the QED couplings to be applied as weights. This design facilitates the use of PWA tools where the MC sample must cover the phase space uniformly.

Implementation of Specific Channels

The generator incorporates both theoretical inputs and experimental data for specific channels:

1. $\pi\pi$ and $\pi^0\eta$ Channels: These are modeled using a dispersive formalism (Muskhelishvili–Omnès representation) that enforces analyticity and unitarity. This approach accounts for strong final-state interactions (FSI) and coupled-channel effects (e.g., $\pi\pi$ and $K\bar{K}$ for the scalar channel). The input covers the kinematic range $2m_\pi < W < 2$ GeV and $0 \le Q^2_{1,2} \le 4$ GeV$^2$.
2. $K\bar{K}$ and $\eta\eta$ Channels: These are modeled using available experimental data from Belle and preliminary BESIII data. Since these measurements are primarily untagged ($Q^2 \approx 0$), the generator uses experimental cross sections for the transverse-transverse ($\sigma_{TT}$) component and employs user-configurable form-factor models (e.g., vector-pole or VMD) to extrapolate to finite virtualities.
3. $f_1(1285) \to \eta\pi^+\pi^-$: This three-body final state is simulated via intermediate states $a_0(980)\pi$ and $f_0(500)\eta$. The model uses an effective Lagrangian framework, incorporating Breit-Wigner propagators and Omnès functions for rescattering, extended to include a second virtuality $Q^2_2$.
4. Custom Single-Resonance Mode: Users can simulate arbitrary hadronic final states through a single resonance using phenomenological models for mass dependence and transition form factors (TFFs), allowing for the estimation of reconstruction efficiencies.

Software Architecture

HadroTOPS is implemented in C++ with support for both double and quadruple precision floating-point arithmetic (using __float128) to ensure numerical stability. It is designed as a class-based structure for easy integration into detector software frameworks. It supports weighted and unweighted event generation, with the unweighting efficiency determined by a pre-run weighted loop to establish the maximum event weight.

Key Results

The paper presents validation results and performance benchmarks:

• Luminosity Functions: The calculated two-photon luminosity functions agree with analytical results from Ref. [1] within statistical fluctuations across a wide range of energies and momentum transfers.
• Comparison with Ekhara3.2: For the $e^+e^- \to e^+e^-\pi^0\pi^0$ process, HadroTOPS shows excellent agreement with Ekhara3.2 in terms of invariant mass ($W$), momentum transfer ($Q^2$), and polar angle ($\cos\theta_\pi$) distributions. Minor discrepancies in the modulation angle $\tilde{\phi}_1$ are attributed to differences in how the differential cross section is evaluated (analytic interpolation vs. direct matrix element evaluation).
• Exclusive vs. Inclusive: The generator successfully demonstrates the difference between inclusive and exclusive formulations, particularly in the modulation angle dependence, which is only captured in the exclusive formula.
• Performance: Weighted event generation times are consistently below 0.1 ms. Unweighting efficiencies vary significantly by channel:
  • High efficiency ($>30:1$) for two-body states like $\pi^+\pi^-$, $K^+K^-$, and $\eta\eta$.
  • Lower efficiency ($\sim 1000:1$ to $2600:1$) for the $f_1(1285) \to \eta\pi^+\pi^-$ channel, reflecting the complexity of the three-body phase space and resonance structure.
• Physical Observables: The generator reproduces expected physical behaviors, such as the dominance of the Born contribution in charged pion production at low masses and the suppression of spin-1 production (Landau-Yang theorem) for real photons, which is recovered only when virtualities are included.

Significance and Claims

The authors claim that HadroTOPS provides a flexible and precise tool for the study of two-photon physics, specifically addressing the need for MC generators that:

1. Support Partial Wave Analysis: By providing flat phase-space samples with exact QED couplings, it allows for the correct description of leading-order QED parts while enabling the replacement of the two-photon cross section with relevant observables in PWA fits.
2. Handle Multi-Hadron Final States: It fills a gap in simulation tools for multi-hadron production, which is currently hindered by a lack of suitable generators.
3. Integrate Theory and Data: It uniquely combines state-of-the-art dispersive theoretical inputs (for $\pi\pi$ and $\pi\eta$) with experimental data (for $K\bar{K}$ and $\eta\eta$), allowing for a unified simulation framework.
4. Enable Differential Studies: The exclusive formulation retains full azimuthal-angle modulations, enabling the study of interference response functions that are typically washed out in inclusive analyses. This is particularly relevant for single-tag and double-tag experimental configurations at facilities like BESIII and Belle II.

The paper positions HadroTOPS as a baseline implementation based on the best currently available theoretical and experimental inputs, with a modular design intended for future extensions to additional hadronic final states and the inclusion of irreducible background contributions (e.g., Bhabha scattering interference) in future updates.