Data-driven method to estimate contamination from light ion beam transmutation at colliders

This paper proposes a data-driven method that utilizes the distinct time dependence and smaller size of beam contaminants to define control regions for quantifying and mitigating the impact of light-ion beam transmutation on physics analyses at colliders like the LHC and RHIC.

The Gist

Imagine you are a chef trying to bake the perfect, giant chocolate cake (representing a Quark-Gluon Plasma, the hot soup of the early universe). To do this, you plan to smash two massive, heavy cakes together in a giant kitchen (a particle collider).

But recently, scientists decided to try something new: smashing smaller, lighter cakes together (like Oxygen or Neon ions). They want to see if the "soup" still forms, just in a smaller pot. This helps them understand if the size of the pot matters for how the soup behaves.

The Problem: The "Rotten Fruit" in the Bowl

Here is the catch: When you spin these light cakes around a track at near-light speed, they are fragile. As they zoom around, they bump into the invisible magnetic fields of other cakes. Sometimes, these bumps cause the light cakes to break apart.

• The Original Cake: A whole Oxygen atom.
• The Broken Pieces: Smaller fragments like Helium, Carbon, or Nitrogen.

In a heavy cake (like Lead), if it breaks, the pieces fly off the track and disappear. But with light cakes, the broken pieces are still the right shape to stay on the track! They keep spinning around, mixing with the fresh, whole cakes.

Over time, your bowl of "pure" Oxygen cakes gets contaminated with a growing pile of broken Helium and Carbon pieces. When you smash them together, you aren't just smashing Oxygen vs. Oxygen; you're smashing Oxygen vs. Helium, or Helium vs. Helium. This ruins your experiment because you can't tell if the results are from the "pure" collision or the "rotten" contamination.

The Solution: A Data-Driven Detective Trick

The authors of this paper propose a clever, data-driven method to figure out exactly how much "rotten fruit" is in the bowl, without needing to simulate every single crash with a supercomputer (which is too hard to do perfectly).

They use a trick similar to sorting laundry by time and size.

1. The Two Clues

They look at two things for every collision:

• Time: How long has the machine been running? (At the start, there is no contamination. Later, there is more.)
• Size (The "Track Count"): How many particles came out of the crash? A head-on collision of two big cakes makes a huge mess (many tracks). A collision involving a tiny broken piece makes a smaller mess (fewer tracks).

2. The "Safe Zone" (Control Regions)

Imagine a graph where the X-axis is Time and the Y-axis is Mess Size.

• The Reference Zone (Start of the run): At the very beginning (Time = 0), the bowl is pure. The "Mess Size" distribution here is your Gold Standard. It shows what a clean collision looks like.
• The High-Purity Zone (Big Messes): Even later in the run, if you see a huge mess (lots of tracks), it must be a collision between two big, whole cakes. Broken pieces are too small to make a huge mess. This area is still "pure" and tells you how the number of collisions is dropping over time (because the beam gets weaker).

3. The Magic Math

Here is the detective work:

1. Look at the High-Purity Zone later in the run. Use it to calculate a "scaling factor." This tells you: "Okay, the beam is weaker now, so we expect 50% fewer big collisions than at the start."
2. Take your Gold Standard (from the start) and shrink it by that 50%.
3. Now, look at the Messy Zone (smaller messes) later in the run.
4. Subtract the shrunk Gold Standard from the actual data.
5. Whatever is left over? That is the Contamination!

It's like saying: "I know I started with 100 apples. I know I lost 20 to rot. If I see 90 apples on the table, but I expected 80 based on my math, then the extra 10 must be the rotten ones I didn't account for."

Why This Matters

This method is like having a self-cleaning filter for your data.

• It's Robust: It doesn't need to know the exact physics of how the atoms break; it just looks at the patterns in the data.
• It's Flexible: It can tell you how the contamination grows minute-by-minute.
• It's Practical: It helps scientists at the Large Hadron Collider (LHC) and RHIC clean up their data so they can finally answer the big question: Does a tiny pot of soup behave like a giant one?

The "Gotchas" (Complicating Factors)

The paper also warns about a few things that could mess up the math, like:

• Pile-up: Sometimes two collisions happen at the exact same time, making the "mess" look bigger than it is. (Like two people dropping their laundry baskets at once).
• Multiple Rotten Fruits: What if there are different types of broken pieces (Helium, Carbon, etc.)? The method can still handle this, but you have to be careful about where you draw the line between "clean" and "dirty."
• Late Start: If the scientists turn on the detectors a few minutes late, they miss the "pure" start. But they can use math to guess what the start looked like.

The Bottom Line

This paper gives scientists a simple, smart way to separate the signal from the noise. By using the time the machine runs and the size of the collision, they can mathematically peel away the "broken pieces" to see the true physics underneath. It turns a messy, confusing problem into a clean, solvable puzzle.

Technical Summary

1. Problem Statement

In high-energy nuclear physics, collisions of relativistic light ions (e.g., Oxygen-16, Neon-20) are increasingly used to study the system-size dependence of Quark-Gluon Plasma (QGP) dynamics. However, a significant experimental challenge arises from beam transmutation.

• Mechanism: As light ions circulate in a collider, they undergo electromagnetic dissociation due to interactions with the strong electromagnetic fields of other ions. Unlike heavy ions (where daughter products often have different charge-to-mass ratios and are lost from the beam), light ion dissociation products (e.g., $^{16}\text{O} \to {}^{12}\text{C} + {}^{4}\text{He}$) often retain the same charge-to-mass ratio ($Z/A = 0.5$).
• Consequence: These daughter ions (contaminants) remain trapped in the beam and accumulate over time. This leads to a growing population of contaminant ions (e.g., Helium-4 in an Oxygen beam) that collide with the primary beam ions.
• Impact: These "contaminant collisions" (e.g., He-O, C-O) have different collision sizes and impact parameters compared to the intended primary collisions (O-O). They introduce asymmetric backgrounds that are difficult to distinguish from glancing primary collisions, potentially skewing measurements of QGP-like effects and system-size dependence.
• Limitation of Current Methods: Simulating these effects from first principles is difficult due to complex variables involving beam optics, decay kinematics, and cross-sections. There is a lack of robust, data-driven methods to quantify this contamination in real-time.

2. Methodology

The authors propose a data-driven estimation method analogous to the ABCD method used in particle physics, utilizing two uncorrelated variables to define control regions:

1. Time ($t$): The transmutation effect is time-dependent. Contaminants are negligible at the start of a beam "fill" but accumulate as the fill progresses.
2. Collision System Size Proxy ($N_{\text{trk}}$): The total number of charged tracks in an event. Primary head-on collisions produce the highest $N_{\text{trk}}$, while contaminant collisions (involving lighter nuclei) produce lower $N_{\text{trk}}$.

The Algorithm:

• Reference Control Region ($t_0 < t < t_1$): Defined at the beginning of the fill. Here, contamination is assumed negligible. This region establishes the "pure" shape of the $N_{\text{trk}}$ distribution for the primary beam.
• High Purity Control Region ($N_{\text{trk}} > N^{\text{cut}}_{\text{trk}}$): Defined by a track count threshold. Since contaminant collisions involve fewer nucleons, they cannot produce $N_{\text{trk}}$ values above this cut. This region tracks the overall normalization (intensity decay) of the primary beam over time.
• Region of Interest (ROI): Defined as $t > t_1$ and $N_{\text{trk}} < N^{\text{cut}}_{\text{trk}}$. This region contains the mixture of primary and contaminant collisions.

Calculation:

1. Calculate a scaling factor $f$ by comparing the integral of the $N_{\text{trk}}$ distribution in the High Purity region at time $t$ versus the reference time. This accounts for beam intensity decay.
2. Scale the Reference Control Region distribution by $f$.
3. Subtract the scaled Reference distribution from the total distribution in the ROI.
4. The residual represents the time-dependent distribution of contaminant collisions.

3. Simulation Study & Results

The authors validated the method using a "toy model" simulation:

• Setup: Simulated $^{16}\text{O}+^{16}\text{O}$ collisions with a single contaminant species ($^{4}\text{He}$).
• Physics Models: Used HG-Pythia (Glauber model for initial state + Pythia for final state) to generate $N_{\text{trk}}$ distributions.
• Time Evolution: Modeled primary beam decay (exponential) and contaminant growth (saturation curve) over a 14-hour fill, assuming a 5% contamination level at the end.
• Key Findings:
  • The method successfully extracted the $N_{\text{trk}}$ shape and rate of the contaminant (He-O) collisions.
  • Closure Test: Comparing extracted rates to input rates showed sub-percent level agreement across the entire fill duration.
  • The method correctly identified that contaminant contributions are concentrated in the low-$N_{\text{trk}}$ region and grow visibly over time.

4. Complicating Factors & Mitigation Strategies

The paper addresses three major experimental challenges and proposes solutions:

• A. Pileup: Multiple interactions per bunch crossing can distort the $N_{\text{trk}}$ distribution and violate the assumption of a constant shape.
  • Mitigation: Use vertex reconstruction to tag/veto pileup events, or employ luminosity leveling (separating beams to maintain constant instantaneous luminosity) to stabilize the pileup profile during the fill.
• B. Multiple Contaminants: If multiple daughter species exist (e.g., C, N, B), the contaminant $N_{\text{trk}}$ distribution may have multiple "knees."
  • Mitigation: Perform an upward scan of the $N^{\text{cut}}_{\text{trk}}$ threshold. If new structures appear in the extracted contaminant distribution, raise the cut to ensure all contaminants are excluded from the High Purity region.
• C. Delayed Data Collection: Detectors often require ramp-up time, meaning the "Reference Control Region" ($t_0$) may already contain some contamination.
  • Mitigation: Perform the extraction with varying start times (delays) and extrapolate the results back to $t=0$ using a linear or simulated fit to estimate the true initial contamination (which is effectively zero).

5. Key Contributions

1. Novel Methodology: Introduces a robust, data-driven technique to quantify beam contamination without relying on complex, uncertain first-principles simulations.
2. Time-Dependent Resolution: Unlike static background estimates, this method provides a dynamic view of how contamination evolves throughout a beam fill.
3. Validation: Demonstrates high precision (sub-percent) in a realistic simulation environment.
4. Practical Framework: Provides specific strategies (leveling, vertex tagging, threshold scanning) to mitigate experimental systematic errors.

6. Significance

• Immediate Application: This method is critical for analyzing recent and upcoming datasets from the LHC (Run 3 and beyond) and RHIC, specifically for Oxygen-Oxygen (OO) and Neon-Neon (NeNe) collision programs.
• Physics Impact: Accurate background subtraction is essential to correctly interpret system-size dependence in QGP formation. Without this, observed effects could be misattributed to QGP dynamics rather than beam contamination.
• Future Outlook: The method enables the LHC to push for higher instantaneous luminosities in light-ion runs by providing a way to monitor and correct for the inevitable buildup of beam contaminants. Furthermore, the authors suggest that understanding these contaminants could eventually allow them to be used as a novel physics probe themselves.