Event shapes and Inclusive Hadron Spectra at FCC-ee energies

This paper analyzes hadronic final states at Future Circular Electron-Positron Collider (FCC-ee) energies using PYTHIA simulations to investigate event shapes and inclusive spectra, extract the strong coupling constant $\alpha_s$ with NNLO accuracy, and assess systematic uncertainties and soft gluon dynamics for future QCD studies.

The Gist

Imagine you are a detective trying to figure out how strong a specific invisible glue is. This glue, called the Strong Force, is what holds the tiny building blocks of the universe (quarks) together to form particles like protons and neutrons. The strength of this glue is measured by a number called $\alpha_s$ (alpha-s).

This paper is a "practice run" for a future super-microscope called the FCC-ee (Future Circular Electron-Positron Collider). The scientists are simulating what will happen when they smash electrons and positrons together at incredibly high speeds to measure this glue strength with extreme precision.

Here is the breakdown of their study using everyday analogies:

1. The Experiment: Smashing Particles Like Billiard Balls

Think of the collider as a giant billiard table. They shoot two balls (an electron and a positron) at each other.

• The Crash: When they hit, they vanish and turn into pure energy, which instantly sprouts new particles (quarks and gluons).
• The Shower: These new particles don't stay alone. They are like a group of rowdy kids at a party who start throwing confetti (gluons) at each other. This creates a "shower" of particles.
• The Result: Eventually, the confetti settles into stable groups called hadrons (like pions and protons), which fly out in all directions. The scientists catch these with a giant detector.

2. The Clues: "Event Shapes"

How do you measure the strength of the glue? You look at the shape of the party.

• The "Thrust" (T): Imagine the particles flying out. If they all fly in two opposite straight lines (like a dumbbell), that's a "2-jet" event. This means the glue held them tight and they didn't scatter much. If they fly everywhere in a ball, that's a "spherical" event.
  • The scientists use a number called Thrust to measure this. A score of 1.0 is a perfect line; 0.5 is a perfect ball.
  • The Twist: They actually look at 1 minus Thrust. Why? Because if the glue is strong and causes a lot of "confetti" (gluons) to fly off, the shape changes. By measuring how much the shape deviates from a perfect line, they can calculate exactly how strong the glue is.
• The "C-Parameter": This is like measuring how "spread out" the party is. It helps them see if the particles are clumping together or spreading thin.

3. The Problem: The "Static" and the "Guests"

The scientists realized that at the high energies of the future collider, the data will be messy. They identified two main sources of "noise":

• The "Flashback" (Initial State Radiation):
  Imagine the two billiard balls are about to crash, but one of them sneezes a photon (a particle of light) before the crash. This sneeze steals some energy. Now, the crash happens at a lower speed than intended.

  • The Issue: This makes the "party" look different than it should. It creates fake shapes that look like the glue is weaker or stronger than it really is. The paper shows that at high energies, this "sneezing" happens a lot, and they have to filter out these events, which throws away a lot of their data.

• The "Imposter Guests" (Backgrounds):
  The collider isn't just making the specific particles they want. It's also making other things, like pairs of W bosons, Z bosons, or even Higgs bosons.

  • The Issue: These are like guests at the party wearing different costumes. A "W-pair" event looks like a 4-jet mess, while a "Z-pair" looks like a 6-jet mess. If the scientists don't separate these "imposters" from the real signal, their calculation of the glue strength will be wrong.

4. The Solution: Cleaning Up the Data

The paper proposes a strategy to clean the data:

• Radiative Cuts: They will mathematically "cut out" any event where the particles lost too much energy to that initial "sneeze."
• Subtractive Corrections: They will use computer simulations to guess how many "imposter guests" are in the room and subtract them from the total count.
• The Trade-off: The paper warns that being too strict with these cuts will throw away too much data (like throwing away 95% of your photos to get the perfect one). They need to find the perfect balance between a clean picture and having enough pictures to analyze.

5. The Results: What They Found

• Measuring the Glue: They successfully simulated how to extract the strength of the strong force ($\alpha_s$) at four different energy levels. They found that as the energy gets higher, the glue gets slightly weaker (a known property of the universe called "asymptotic freedom").
• The "Party Size" (Multiplicity): They also counted how many particles are created in each crash. They found that as you smash the balls harder, the "party" gets bigger (more particles), but it follows a predictable logarithmic pattern, just like a tree growing branches.
• The "Shape of the Crowd" (Momentum): They looked at how fast the particles are moving. They found that the distribution of speeds forms a specific "hump-backed" shape, which matches what theory predicts, though at the very highest energies, the computer simulation started to drift slightly from the theory.

The Big Picture

This paper is a rehearsal. The scientists are saying: "We know the FCC-ee will be amazing, but it will be messy. We have to be very careful about filtering out the 'sneezes' and 'imposters' to get a perfect measurement of the Strong Force."

If they succeed, they will know the strength of the universe's glue with 0.1% accuracy. This is crucial because if we don't know exactly how strong this glue is, our predictions for other things—like how the Higgs boson works or how heavy the top quark is—will be slightly off. It's about making sure our map of the universe is perfectly accurate.

Technical Summary

Here is a detailed technical summary of the paper "Event shapes and Inclusive Hadron Spectra at FCC-ee energies" by Philip Mathew, R. Aggarwal, and M. Kaur.

1. Problem Statement

The Future Circular Electron-Positron Collider (FCC-ee) aims to operate at center-of-mass (c.m.) energies ranging from the $Z$-pole (91.2 GeV) up to the top-quark pair threshold (365 GeV). While the primary goal is to determine the strong coupling constant ($\alpha_s$) with unprecedented precision (0.1%), higher collision energies introduce significant challenges not present at previous colliders like LEP:

• Initial State Radiation (ISR): Enhanced QED radiation from initial leptons reduces the effective c.m. energy ($\sqrt{s'}$) of the hadronic system. This leads to "radiative returns" where events mimic lower-energy physics (e.g., $Z$-pole events at 365 GeV), distorting event shape distributions.
• Electroweak Backgrounds: At energies above the $Z$-pole, significant backgrounds arise from $W^+W^-$, $ZZ$, $t\bar{t}$, and Higgs ($ZH$, $H\nu\bar{\nu}$) production. These processes create complex multi-jet topologies that mimic or distort the signal from the primary $e^+e^- \to \gamma/Z \to \text{hadrons}$ channel.
• Systematic Uncertainties: As statistical errors become negligible due to high luminosity, systematic uncertainties (ISR corrections, background subtraction, hadronization modeling, and detector effects) become the limiting factor for $\alpha_s$ extraction.

2. Methodology

The authors conducted a phenomenological study using PYTHIA 8.313 to simulate hadronic final states at four specific FCC-ee c.m. energies: 91.2 GeV, 160 GeV, 240 GeV, and 365 GeV.

• Event Generation: Samples of $5 \times 10^6$events were generated for each energy point. The study isolated the$e^+e^- \to \gamma/Z$channel and separately simulated electroweak backgrounds ($WW, ZZ, t\bar{t}, H$).
• Observables:
  • Event Shapes: The study focused on Thrust ($T$) and the C-parameter. These are infrared- and collinear-safe observables sensitive to hard gluon radiation and the geometry of the final state.
  • Inclusive Spectra: The study analyzed charged hadron multiplicity ($\langle N_{ch} \rangle$) and the scaled momentum distribution (via the variable $\xi = \ln(1/x_p)$).
• Analysis Techniques:
  • Radiative Cuts: Applied cuts on the effective energy $\sqrt{s'}$ to mitigate ISR distortions.
  • Background Subtraction: Modeled the contribution of electroweak backgrounds to determine the necessary subtraction corrections.
  • $\alpha_s$ Extraction: Fitted the simulated event shape distributions to Next-to-Next-to-Leading Order (NNLO) perturbative QCD predictions (using coefficients from Weinzierl) to extract $\alpha_s$.
  • Hadronization Modeling: Compared parton-level predictions with hadron-level PYTHIA results to quantify hadronization corrections.
  • Inclusive Evolution: Compared the energy evolution of multiplicity and $\xi$ peaks against 3NLO (for multiplicity) and Modified Leading Log Approximation (MLLA) (for $\xi$) theoretical predictions, invoking Local Parton-Hadron Duality (LPHD).

3. Key Contributions

• Quantification of High-Energy Distortions: The paper provides a comprehensive map of how ISR and specific electroweak backgrounds distort Thrust and C-parameter distributions at FCC-ee energies. It highlights that radiative returns create distinct peaks in distributions that must be filtered.
• Optimization of Correction Strategies: It demonstrates that while tight radiative cuts (e.g., $\sqrt{s'} \ge 0.95\sqrt{s}$) effectively remove ISR distortions, they result in a massive loss of statistics (e.g., only ~4.7% survival at 365 GeV). This underscores the need for sophisticated correction strategies rather than simple cuts.
• Background Characterization: The study details how different backgrounds ($WW, ZZ, t\bar{t}$) shift event shapes. For instance, $t\bar{t}$ events produce highly spherical topologies (high $1-T$and$C$), while$WW$ events are more collimated but still distinct from the signal.
• Validation of Fitting Procedures: The authors validated their NNLO fitting procedure at 91.2 GeV against LEP data, showing consistency within ~1%, before extrapolating the method to higher energies.

4. Key Results

• $\alpha_s$ Determination:
  • Extracted values of $\alpha_s$ decrease with energy as expected due to asymptotic freedom:
    • 91.2 GeV: $\alpha_s \approx 0.1284$
    • 160 GeV: $\alpha_s \approx 0.1186$
    • 240 GeV: $\alpha_s \approx 0.1123$
    • 365 GeV: $\alpha_s \approx 0.1061$
  • The deviation between $\alpha_s$ values extracted from Thrust and C-parameter is $\le 0.6\%$ across all energies.
  • Statistical errors are negligible; the dominant uncertainty is systematic.
• Impact of ISR and Backgrounds:
  • At 365 GeV, without cuts, the ISR "radiative return" creates a massive peak at $\sqrt{s'} \approx 91.2$ GeV, completely obscuring the 365 GeV signal in event shape distributions.
  • Electroweak backgrounds reduce the survival fraction of the pure $\gamma/Z$ signal to 49% at 365 GeV after necessary corrections.
• Inclusive Hadron Spectra:
  • Multiplicity: The mean charged multiplicity $\langle N_{ch} \rangle$ follows the 3NLO QCD prediction logarithmically. Values range from 20.37 (91.2 GeV) to 27.80 (365 GeV). Slight deviations appear above 240 GeV, potentially hinting at limits of fixed-order expansions.
  • Momentum Distribution ($\xi$): The peak position $\xi^*$ evolves with energy. While PYTHIA results match MLLA predictions at 91.2 GeV, they begin to deviate significantly above 160 GeV. This suggests that current fragmentation parameters (tuned at LEP) may not extrapolate perfectly to FCC-ee energies, or that higher-order logarithmic terms are required.

5. Significance

This study serves as a critical reference for the QCD program at the FCC-ee.

• Feasibility of Precision Physics: It confirms that while $\alpha_s$ can be extracted with high precision at FCC-ee energies, the path to 0.1% accuracy is blocked by systematic uncertainties related to ISR and backgrounds, not statistics.
• Strategic Guidance: The results argue against relying solely on hard kinematic cuts to remove backgrounds, as the loss of statistics is prohibitive. Instead, they advocate for advanced subtractive corrections and improved theoretical modeling of ISR and hadronization.
• Model Validation: The observed deviations in inclusive spectra at high energies provide a roadmap for refining Monte Carlo generators (like PYTHIA) and theoretical frameworks (MLLA/3NLO) to better describe the transition between perturbative and non-perturbative QCD regimes at multi-hundred GeV scales.

In conclusion, the paper establishes that the FCC-ee has the potential to revolutionize QCD precision measurements, but realizing this potential requires a re-evaluation of correction strategies and a deeper understanding of systematic effects at high energies.