MAcNLOPS for ZZ Pair Production at the LHC

This paper presents and validates a MAcNLOPS implementation for $pp \to ZZ$ production in MadGraph5_aMC@NLO + Pythia8 that eliminates negative H weights through a veto on the first shower emission, offering a promising alternative to MC@NLO with reduced negative weights and negligible computational overhead.

The Gist

Imagine you are trying to predict exactly how two heavy, invisible balloons (called Z bosons) will bounce around after a massive collision inside a giant, high-speed pinball machine (the Large Hadron Collider, or LHC).

To do this, physicists use a computer program that acts like a two-step recipe:

1. The Hard Calculation: A precise math formula that predicts the main collision.
2. The Shower: A simulation that adds the "messy" details, like tiny sparks and debris flying off, which the simple math formula misses.

The problem is that when you try to mix these two steps together perfectly, the computer sometimes gets confused. It tries to subtract the "messy" part from the "hard" part to avoid double-counting. But sometimes, the math gets so tricky that the computer assigns a "negative weight" to certain predicted events.

The Problem: The "Negative Weight" Glitch

Think of a "negative weight" like a debt in a bank account. If you have 100 events with a positive weight (money in the bank) and 10 events with a negative weight (debt), your total is 90.

While this is mathematically correct, it's a nightmare for the computer. To get a clear picture of the 90 events, the computer has to generate thousands of extra events just to cancel out the noise from the debts. It's like trying to hear a whisper in a room full of people shouting; you have to shout louder (generate more data) just to find the signal. This wastes time and computing power.

The Solution: The "MAcNLOPS" Fix

The authors of this paper, Yuxiao Che and Rikkert Frederix, tested a new method called MAcNLOPS to fix this debt problem for the specific case of creating two Z bosons.

Here is how their method works, using a simple analogy:

The Old Way (MC@NLO):
Imagine you are sorting a pile of mail. You have "Good Letters" (positive events) and "Bad Letters" (negative events). The old method says, "Keep both, but remember the Bad Letters cancel out some of the Good ones." You end up with a messy pile where you have to do extra math to figure out the final count.

The New Way (MAcNLOPS):
The new method says, "Let's just throw away the Bad Letters immediately."

• Step 1: They identify the "Bad Letters" (the negative events) and delete them from the pile.
• Step 2: But wait! If you just delete them, you lose information. So, they apply a "Veto" (a strict rule) to the "Good Letters" that are about to get messy.
• The Trick: They say, "If a Good Letter looks like it might have been a Bad Letter, we will randomly reject it with a specific probability." This rejection acts as a perfect mathematical compensation for the Bad Letters they threw away.

It's like a bouncer at a club. Instead of letting in a group of people who will cause trouble (negative events) and then asking them to leave later, the bouncer simply doesn't let them in. To make sure the crowd size stays the same, the bouncer occasionally turns away a few nice people at the door, but only in a way that perfectly balances the math.

What They Found

The authors tested this new "Bouncer" method against the old "Debt" method using the collision of two Z bosons.

1. The Results Match: When they looked at the final results—how the particles moved, how much energy they had, and where they went—the new method gave almost the exact same answer as the old method.
2. The "Low Energy" Zone: There was a tiny, tiny difference in the very lowest energy areas (the "soft" region). The paper explains this is like a rounding error that happens because of how the computer cuts off the simulation at very low speeds. It's a negligible difference that doesn't affect the big picture.
3. Efficiency: The best part? The new method removed all the "negative weight" events without making the computer work any harder. In fact, it made the process cleaner.

The Bottom Line

This paper proves that for creating pairs of Z bosons, you can use this new "Veto" trick to get rid of the confusing negative numbers that slow down simulations. The results are just as accurate as the old way, but the computer doesn't have to do the extra work of managing the "debts."

Note on Limitations: The paper specifically notes that this method only removes the "Bad Letters" (negative H events). It doesn't fix the "Bad Bank Accounts" (negative S events) that can happen if the starting math is negative. However, for the specific job of simulating Z boson pairs, this new method is a promising, cleaner alternative to the old way.

Technical Summary

Technical Summary: MAcNLOPS for ZZ Pair Production at the LHC

Problem Statement
Precision predictions for LHC observables require Next-to-Leading Order (NLO) calculations matched with Parton Showers (PS). The standard MC@NLO method achieves this by subtracting the shower approximation from the real-emission matrix element and adding the integrated contribution to the Born-like sample. While this preserves NLO accuracy, it frequently generates events with negative weights. These negative weights arise when the shower approximation locally exceeds the real matrix element ($K > R$). Although mathematically valid due to local cancellation, negative weights are practically undesirable: they reduce the statistical power of event samples, necessitate the generation of significantly more events to achieve a given precision, and complicate their use in experimental analyses.

Existing alternatives like POWHEG largely avoid negative weights by generating the hardest emission via a Sudakov formulation, but this alters the division of labor between the NLO calculation and the shower, potentially affecting the shower's logarithmic structure and recoil strategy. Other methods, such as MC@NLO-$\Delta$, reduce negative weights but do not eliminate them entirely.

Methodology: MAcNLOPS
This paper presents an implementation of the MAcNLOPS (MC@NLO with Positive Weights) prescription, a minimal modification of the MC@NLO framework designed to remove negative weights while preserving the shower's role in generating the first emission.

The core mechanism involves a two-step process applied to a standard MC@NLO event sample (generated via MadGraph5_aMC@NLO and interfaced with Pythia8):

1. Removal of Negative H Events: The "H" events (real-emission kinematics) with negative weights (occurring where $K > R$) are explicitly removed from the sample. The "H" events with positive weights ($R > K$) are retained.
2. Veto on S Events: To compensate for the removed negative H contribution, a veto is applied to the "S" events (Born-like kinematics) after the first shower emission.
   • The S events are showered up to the first QCD emission.
   • A veto probability $P_{veto}$ is calculated based on the ratio of the real matrix element ($R$) to the shower counterterm ($K$).
   • In the region where $K > R$ (the source of negative weights), S events are accepted with probability $R/K$.
   • In the region where $R \ge K$, S events are always accepted.

This procedure ensures that the negative H contribution is cancelled by a corresponding reduction in the S sample, resulting in a final sample with no negative H weights. The method retains the additive structure of MC@NLO for positive regions and the multiplicative treatment for negative regions.

Key Contributions and Implementation Details

• First Implementation: This work provides the first implementation of MAcNLOPS for $pp \to ZZ$ production.
• Process Selection: The $ZZ$ process was chosen because its Born-level final state is a color singlet. This isolates the initial-state radiation (ISR) complications, avoiding the added complexity of final-state radiation from colored particles, while remaining more complex than a simple $2 \to 1$ process.
• Workflow: The implementation splits the MC@NLO sample into S and H events. Negative H events are discarded. S events are passed to Pythia8 for the first emission; a veto is applied based on the $R/K$ ratio; the accepted S events are then combined with the positive H events and passed back to Pythia8 for the remainder of the shower evolution.
• Theoretical Accuracy: The authors demonstrate that the difference between MAcNLOPS and standard MC@NLO is proportional to terms that vanish at NLO accuracy. The discrepancy is of NNLO size for infrared-safe observables, ensuring that the NLO accuracy of the inclusive cross-section is preserved.

Results
The implementation was validated against standard MC@NLO predictions for $pp \to ZZ$ at $\sqrt{s} = 13$ TeV using Pythia8 for showering.

• Radiation-Sensitive Observables: Distributions sensitive to the first emission, such as the transverse momentum of the ZZ pair ($p_T^{ZZ}$) and the azimuthal separation of the Z bosons ($\Delta\phi_{ZZ}$), show agreement between MAcNLOPS and MC@NLO within statistical uncertainties.
• Inclusive Observables: Born-dominated observables, specifically the invariant mass ($M_{ZZ}$) and rapidity separation ($\Delta y_{ZZ}$), remain unchanged, confirming that the veto procedure does not bias the hard process kinematics.
• Low-$p_T$ Region: In the very low-$p_T$ region (near the shower cut-off), a small power-suppressed difference is observed. This is attributed to the mismatch between the removal of negative H events (which can occur below the cut-off via the $\Theta$ function) and the application of the veto (which only applies to S events with emissions above the cut-off). This effect is formally negligible but leads to a slight rate difference in the first few bins.

Significance and Claims
The paper claims that MAcNLOPS is a promising alternative to standard MC@NLO that successfully removes all negative H weights with negligible additional computational cost.

• It preserves the key practical advantage of MC@NLO: the first emission is still generated by the parton shower, maintaining the shower's native logarithmic structure and recoil strategy.
• The method does not eliminate negative S weights (which arise if the Born cross-section $eB$ is negative), but it isolates the negative H problem, which is often the dominant source of inefficiency.
• The authors conclude that while the current implementation is tailored to color-singlet final states, it offers a controlled environment to study the method's behavior, suggesting it could be extended to more complex processes in future work.