Analysis of quarkonium polarization in proton-proton (p-p) collisions at LHC using PYTHIA model

This paper utilizes the PYTHIA8 Monte Carlo model to simulate and analyze the polarization of $J/\psi$ and $\Upsilon(1S)$ quarkonia in proton-proton collisions at 7 and 13 TeV, incorporating detector effects to compare theoretical predictions with recent ALICE experimental measurements.

The Gist

The Big Picture: The "Spinning Top" Mystery

Imagine you have a collection of tiny, heavy spinning tops (called quarkonia, specifically J/psi and Upsilon). These tops are made of a heavy particle and its anti-particle stuck together. When they break apart, they shoot out two smaller particles (muons) in opposite directions.

Physicists want to know: Are these tops spinning in a specific direction, or are they spinning randomly?

• Random spin: The tops are like a bag of marbles rolling around; they point everywhere.
• Ordered spin: The tops are like a synchronized dance troupe, all pointing their "heads" in the same direction.

This "direction" is called polarization. Knowing how they spin helps scientists understand the fundamental rules of how matter is built (Quantum Chromodynamics).

The Problem: The "Blindfolded Photographer"

The scientists in this paper used a computer program called PYTHIA to simulate these collisions. Think of PYTHIA as a video game engine that creates a perfect, theoretical world.

In this perfect world, the tops are completely random (unpolarized). They spin in every direction equally. If you took a photo of them in this perfect world, the angles of the broken pieces would look like a perfect circle.

However, real life isn't perfect.

In a real experiment (like at the LHC), the "cameras" (detectors) aren't perfect. They have two main flaws:

1. Blurry Vision (Momentum Smearing): The camera can't measure the speed of the broken pieces perfectly. It's a little blurry.
2. Bad Angles (Inefficiency): The camera has blind spots. It misses some pieces depending on which way they are flying.

The Discovery: The "Fake Spin" Illusion

The researchers asked a crucial question: "If we take our perfect, random simulation and run it through a 'blurry camera' with 'blind spots,' does it start to look like the tops are spinning in a specific direction?"

The answer was a resounding YES.

Here is the analogy:
Imagine you have a bag of marbles rolling randomly in all directions (unpolarized). Now, imagine you put a filter over your eyes that only lets you see marbles rolling to the right if they are moving slowly, but lets you see marbles rolling to the left if they are moving fast.

Suddenly, the marbles you see don't look random anymore. They look like they have a pattern! You might think, "Wow, the marbles are all spinning to the right!" But that's an illusion created by your filter.

The paper found that:

• When they added "blur" and "blind spots" to their computer simulation, the data started showing a "fake" polarization.
• For the lighter tops (J/psi), this fake effect was very strong. It looked like they were spinning one way at low speeds and the other way at high speeds.
• For the heavier tops (Upsilon), the effect was smaller because they are heavier and easier to track, but it was still there.

The Solution: Cleaning Up the Mess

The researchers then tried to fix the "blurry photo." They applied mathematical corrections to account for the camera's bad angles and blurry vision.

The Result:
Once they cleaned up the data, the "fake spin" disappeared. The tops went back to looking perfectly random, just as the computer simulation originally intended.

Why This Matters

This paper is a warning label for scientists. It says:

"If you see a pattern in how these particles spin, don't assume it's a new law of physics. It might just be your camera is blurry or has a blind spot."

They showed that if you don't correct for these detector errors, you can invent a "puzzle" that doesn't actually exist. By fixing the errors, the data matches the theory (which says they are random), and the "puzzle" vanishes.

Summary in One Sentence

This paper proves that the strange "spinning patterns" seen in particle physics data can sometimes be an optical illusion caused by imperfect detectors, and that once you correct for those imperfections, the particles are actually spinning randomly as expected.

Technical Summary

1. Problem Statement

The production mechanism of quarkonia (bound states of heavy quark-antiquark pairs, specifically $J/\psi$ and $\Upsilon(1S)$) in hadronic collisions remains a significant challenge in Quantum Chromodynamics (QCD). While theoretical models like the Color-Singlet Model (CSM) and Non-Relativistic QCD (NRQCD) attempt to describe production cross-sections, they often fail to consistently predict the polarization (spin alignment) of these vector mesons. This discrepancy is known as the "$J/\psi$ polarization puzzle."

A critical issue in experimental measurements is the potential for artificial polarization signals induced by detector effects. Previous studies using the PYTHIA8 Monte Carlo generator reported non-zero polarization parameters, which contradicts the theoretical expectation that quarkonia produced in PYTHIA8 are inherently unpolarized and decay isotropically. The authors investigate whether detector inefficiencies, momentum smearing, and kinematic selection cuts can create spurious anisotropies in the angular distribution of decay muons, leading to false polarization measurements.

2. Methodology

The study utilizes the PYTHIA8 event generator (Tune 4C) to simulate $p-p$ collisions at center-of-mass energies of $\sqrt{s} = 7$ and $13$ TeV. The analysis focuses on the forward rapidity region ($2.5 < y < 4.0$), mimicking the acceptance of the ALICE Muon Spectrometer.

• Physics Processes: Simulations include gluon fusion ($gg \to Q\bar{Q}$) and quark-antiquark annihilation ($q\bar{q} \to Q\bar{Q}$) within the NRQCD framework, accounting for both color-singlet and color-octet states. The study systematically varies Multi-Parton Interactions (MPI) and Color Reconnection (CR) settings to isolate their effects.
• Reference Frames: Polarization is analyzed in two standard frames:
  • Helicity (HE) Frame: The $z$-axis aligns with the quarkonium momentum in the center-of-mass frame.
  • Collins-Soper (CS) Frame: The $z$-axis aligns with the bisector of the angle between the two incoming proton beams.
• Polarization Extraction: The angular distribution of decay muons ($W(\theta, \phi)$) is parameterized to extract polarization parameters ($\lambda_\theta, \lambda_\phi, \lambda_{\theta\phi}$). Three independent extraction methods are employed:
  1. 1D Angular Fit: Fitting integrated angular distributions.
  2. Population Counting: Bin-by-bin ratio analysis.
  3. Average Method: Computing mean values of trigonometric functions of decay angles.
• Detector Simulation: To mimic realistic experimental conditions, the authors introduce:
  • Momentum Smearing: Gaussian smearing of single muon momentum ($p_T$) with resolutions ranging from 0.5% to 3%.
  • Efficiency Maps: Probabilistic acceptance/rejection based on a constructed efficiency matrix dependent on $p_T$ and pseudorapidity ($\eta$).
  • Mass Window Cuts: Selection of dimuon candidates within $\pm 3\sigma$ of the fitted invariant mass peak.

3. Key Contributions

• Resolution of the PYTHIA Polarization Discrepancy: The paper clarifies why previous PYTHIA-based studies reported non-zero polarization. It demonstrates that the generator itself produces unpolarized quarkonia; the observed anisotropies are artifacts of detector-level distortions.
• Quantification of Detector-Induced Bias: The study provides a quantitative analysis of how momentum resolution and mass-window cuts distort angular distributions. It specifically highlights that the $J/\psi$ meson is more susceptible to these effects than $\Upsilon(1S)$ due to its lower mass and larger decay opening angle.
• Validation of Correction Strategies: The authors demonstrate that applying rigorous acceptance and resolution corrections can successfully recover the original isotropic (unpolarized) distribution, validating the necessity of such corrections in experimental analyses.
• Frame Invariance Check: The study confirms that while individual polarization parameters ($\lambda_\theta, \lambda_\phi$) are frame-dependent and sensitive to smearing, the frame-invariant parameter $\lambda_{inv}$ remains robust and consistent across different frames.

4. Key Results

• Generator Level: In the absence of detector effects, all extracted polarization parameters ($\lambda_\theta, \lambda_\phi, \lambda_{\theta\phi}$) are statistically consistent with zero across the entire $p_T$ range (0–14 GeV/c) for both $J/\psi$ and $\Upsilon(1S)$. Variations in MPI and CR settings did not induce significant polarization.
• Detector Effects (Smearing + Mass Cut):
  • Applying momentum smearing and mass-window cuts to the unpolarized sample induced significant artificial polarization.
  • For $J/\psi$, the $\lambda_\theta$ parameter showed a transition from apparent longitudinal polarization ($\lambda_\theta < 0$) at low $p_T$ to transverse polarization ($\lambda_\theta > 0$) at high $p_T$ in the Helicity frame.
  • The effect was more pronounced for $J/\psi$ than $\Upsilon(1S)$ because the lower mass of $J/\psi$ leads to a broader reconstructed mass distribution, making the angular distribution more sensitive to the $\pm 3\sigma$ mass cut.
• Correction Recovery: After applying efficiency corrections derived from a 40% subset of the data to the remaining 60%, the extracted polarization parameters returned to values consistent with zero, effectively restoring the isotropic nature of the decay.
• Comparison with Experiment: The corrected simulation results for $\Upsilon(1S)$ align well with preliminary ALICE measurements at $\sqrt{s}=13$ TeV. No direct comparison was made for $J/\psi$ as corresponding ALICE measurements were not yet available at the time of writing.

5. Significance

This work is crucial for the interpretation of quarkonium polarization data at the LHC. It establishes that detector resolution and selection biases are major sources of systematic error that can mimic physical polarization signals.

• For Experimentalists: It underscores the absolute necessity of detailed detector modeling and rigorous correction procedures to avoid misinterpreting instrumental artifacts as new physics or deviations from QCD predictions.
• For Theorists: It clarifies that the "puzzle" of non-zero polarization in some Monte Carlo studies may stem from uncorrected detector effects rather than flaws in the underlying production models (like NRQCD).
• Methodological Impact: The study validates the use of frame-invariant quantities ($\lambda_{inv}$) as a robust observable when comparing results across different reference frames or experimental setups with varying detector characteristics.

In conclusion, the paper confirms that quarkonia produced in PYTHIA8 are unpolarized and that the non-zero polarization observed in some analyses is an artifact of experimental conditions that can be corrected, thereby reinforcing the need for high-precision detector modeling in future LHC analyses.