Phenomenology of Matching Exponentiated Photonic Radiation to a Parton Shower in KKMChh

This paper introduces the KKMChh program, which extends the soft photon exponentiation of KKMC to hadron collisions by developing the NISR algorithm to consistently interface with parton showers and QED-inclusive parton distribution functions, along with presenting results on its phenomenological impact.

The Gist

Imagine you are trying to predict exactly how a crowd of people will behave at a massive concert. You have a very detailed map of the crowd (the Parton Distribution Functions, or PDFs) that tells you where everyone is standing before the music starts. You also have a super-accurate weather app (the KKMChh program) that predicts how the wind (photons) will blow and how people will react to it once the show begins.

The problem? Your crowd map might have already guessed how the wind would affect the people. If you use the map and the weather app together without checking, you might count the wind's effect twice. This is called "double-counting," and it ruins your prediction.

This paper is about building a clever "undo button" to fix that mistake. Here is the breakdown of their solution:

1. The Setup: The Crowd and the Wind

In the world of particle physics, when protons smash together (like at the Large Hadron Collider), they are actually made of smaller particles called quarks.

• The PDFs (The Map): These are mathematical descriptions of what the quarks are doing inside the proton. Some modern maps already include a little bit of "wind" (QED radiation) in their calculations.
• KKMChh (The Weather App): This is a sophisticated computer program that simulates the collision. It is famous for being incredibly precise about how the "wind" (photons) blows, especially the soft, gentle breezes that happen constantly.

2. The Problem: The Double-Counting Trap

If you use a map that already includes the wind, and then you run the weather app which also calculates the wind, you end up with a storm that is twice as strong as reality.

• The Analogy: Imagine you are baking a cake. Your recipe (the PDF) already says, "Add 1 cup of sugar." But your cooking assistant (KKMChh) also says, "Add 1 cup of sugar!" If you listen to both, your cake is way too sweet. You need a way to tell the assistant, "Hey, the sugar is already in the mix; don't add more."

3. The Solution: The "Negative ISR" (NISR)

The authors developed a new algorithm called NISR (Negative Initial State Radiation). Think of this as a "sugar subtractor."

• How it works: Before the weather app starts blowing its wind, NISR looks at the crowd map and says, "Okay, this part of the map already accounts for the wind. Let's mathematically subtract that wind from the map first."
• The "Undo" Button: It's like hitting Ctrl+Z on the wind effects in the map. Once the wind is removed from the map, the weather app (KKMChh) can safely add its own, highly accurate wind simulation without worrying about double-counting.

4. Why This Matters: The "Forward-Backward" Dance

The paper tests this new "sugar subtractor" by looking at a specific dance move called Forward-Backward Asymmetry.

• The Dance: When particles collide, they sometimes fly off to the left (forward) or the right (backward). Physicists want to know the exact ratio of left vs. right to understand the fundamental laws of the universe.
• The Test: They ran simulations with the new NISR tool. They found that for most of the "dance" (low energy), the difference was tiny—like a slight shift in weight that you wouldn't notice. However, at higher energies or when looking at extreme angles, the correction became important.
• The Result: The NISR tool successfully removed the extra "sugar" (double-counted wind) and allowed the weather app to give a perfect prediction, matching the physics of the real world.

5. The Catch: It Takes Time

The authors admit that running this "sugar subtractor" takes extra computer time.

• The Metaphor: It's like hiring a professional editor to check your recipe before you bake. If you are making a simple cupcake (a low-precision experiment), you might not need the editor; it's faster to just bake it. But if you are making a wedding cake (a high-precision experiment like the Z boson studies), the extra time is worth it to ensure the cake isn't ruined.

Summary

In short, this paper introduces a smart "undo" feature for particle physics simulations. It allows scientists to use modern, complex maps of protons (which already include some wind effects) while still using a super-accurate weather app to simulate the wind. By mathematically removing the wind from the map first, they ensure the final prediction is perfectly accurate, avoiding the mistake of counting the wind twice. This helps physicists understand the universe with greater precision than ever before.

Technical Summary

Here is a detailed technical summary of the paper "Phenomenology of Matching Exponentiated Photonic Radiation to a Parton Shower in KKMChh."

1. Problem Statement

The paper addresses a critical challenge in high-energy physics phenomenology: the consistent matching of the KKMChh Monte Carlo generator with Parton Distribution Functions (PDFs) that already include Quantum Electrodynamics (QED) effects.

• The Context: KKMChh is an extension of the precision $e^+e^-$ generator KKMC, adapted for hadronic collisions ($pp \to Z/\gamma^* \to \ell\ell + n\gamma$). It utilizes the Coherent Exclusive Exponentiation (CEEX) formalism to implement exact $\mathcal{O}(\alpha)$ and next-to-leading logarithm (NLL) $\mathcal{O}(\alpha^2 L)$ QED Initial State Radiation (ISR) with soft photon resummation at the amplitude level.
• The Conflict: Modern PDF sets (e.g., NNPDF3.1-LUXqed) often include QED evolution, meaning the parton densities themselves account for photon radiation.
  • If KKMChh generates ISR on top of these PDFs without correction, it leads to double-counting of QED effects.
  • Conversely, standard PDFs use a collinear approximation for photons, whereas KKMChh calculates ISR with full kinematic precision (including transverse momentum). Simply using a PDF with QED but ignoring KKMChh's superior ISR treatment loses precision.
• The Goal: Develop a method to "back out" the QED component from existing PDFs at a specific scale so that KKMChh can generate the ISR with its high-precision algorithm, ensuring no double-counting while maintaining theoretical accuracy.

2. Methodology: The NISR Algorithm

The authors introduce a novel algorithm termed Negative Initial State Radiation (NISR) to resolve the matching problem.

• Concept: Instead of modifying the PDFs offline, NISR is applied "on-the-fly" during event generation. It effectively convolves the quark PDFs with "reversed" radiator functions to subtract the QED evolution present in the PDF up to a specific scale $Q_0$.
• Mathematical Implementation:
  • The standard ISR radiator function $\rho_{ISR}$ is modified to include a "negative" component.
  • The PDF for a quark $q$ is convoluted with a reversed half-radiator function:
    $$\rho_{ISR}\left(-\frac{1}{2}\gamma_{ISR}(x_s), u\right)$$
  • Two new kinematic variables, $u_1$ and $u_2$, are introduced for the two incoming partons. The constraint $x = \hat{x}(1-u_1)(1-u_2)$ ensures momentum conservation.
  • The effective ISR energy fraction $v$ is modified to $v'$ such that:
    $$1 - v' = (1 - v)(1 - u_1)(1 - u_2)$$
• Scale Selection ($Q_0$):
  • For PDFs that include QED evolution (e.g., LUXqed), $Q_0$ is set to the hard process scale (e.g., the mass of the $Z$ boson). This removes the QED evolution accumulated from the starting scale up to the hard interaction.
  • For PDFs without explicit QED evolution (but containing QED contamination in input data), $Q_0$ is set to a low scale (e.g., a few GeV) to remove the contamination accumulated during QCD evolution.

3. Key Contributions

1. Development of NISR: The formulation of a practical algorithm to decouple QED effects from PDFs dynamically, allowing the use of current state-of-the-art PDF sets with high-precision QED generators.
2. Validation of Mass Independence: The authors demonstrate that the logarithmic dependence on quark masses in the ISR calculation cancels out with the NISR subtraction through second order. This confirms that the method is robust against variations in quark mass definitions (current vs. constituent masses).
3. Phenomenological Application: The integration of NISR into studies of Forward-Backward Asymmetry ($A_{FB}$), a precision observable sensitive to Initial-Final Interference (IFI) and QED effects.

4. Results

The paper presents several quantitative results validating the NISR approach:

• Cross-Section Consistency:
  • Table 1 and Table 2 compare the Born cross-section (no photons) with the QED-corrected cross-section (ISR + NISR, no FSR).
  • The fractional difference between the "Born" and "ISR+NISR" cases is small (ranging from ~0.07% to 0.27% depending on quark flavor), consistent with expected Next-to-Leading Order (NLO) non-logarithmic corrections ($\delta_Q \approx Q^2 \frac{\alpha}{\pi}(-\frac{1}{2} + \frac{\pi^2}{3})$).
  • Crucially, when quark masses were artificially set to 500 MeV (Table 2), the results remained nearly identical to the physical mass case. This proves the cancellation of logarithmic mass dependence, validating the theoretical consistency of the NISR subtraction.
• Parton Luminosity Analysis:
  • The authors analyzed the QED-corrected joint parton luminosity functions. They found that the up-quark contribution dominates the QED correction due to its electric charge ($2e/3$), rather than its small mass.
  • Figures 2 and 3 show that QED corrections modify the luminosity distribution, particularly at higher invariant masses ($M_{q\bar{q}}$).
• Forward-Backward Asymmetry ($A_{FB}$):
  • Applied to $pp \to \mu^+\mu^-$ at 8 TeV.
  • Findings: The effect of NISR on $A_{FB}$ is negligible for invariant masses up to 130 GeV.
  • At higher masses (>130 GeV) and high rapidities, the effect increases, but statistical limitations in that region make the precise magnitude difficult to pin down in this study.
  • The results are largely independent of whether NISR is applied at the hard scale (with LUXqed) or a low scale (2 GeV) with standard PDFs.

5. Significance and Conclusion

• Precision Physics: The NISR algorithm enables the use of KKMChh's superior exponentiated QED calculations with modern PDF sets without double-counting. This is essential for precision measurements at the LHC (e.g., $W$-mass, $Z$-pole observables).
• Efficiency Trade-off: The authors note that while NISR is theoretically necessary for consistency with QED-inclusive PDFs, it adds computational overhead to Monte Carlo runs.
  • Recommendation: If the required precision of a specific calculation is lower than the magnitude of the NISR effect (which is often small, <0.2%), it is more efficient to omit NISR.
  • Future Outlook: The authors suggest that generating "pruned" PDF sets (QED-removed) in advance would be the most efficient solution for future high-precision campaigns.
• Tribute: The paper acknowledges the foundational work of Stanisław Jadach (1947–2023), a pioneer in the field of radiative corrections and the creator of the KKMC framework.

In summary, this paper provides a robust technical solution for interfacing high-precision QED event generators with modern QED-inclusive PDFs, ensuring the integrity of precision electroweak measurements in hadron collisions.