One-, two-, and three-dimensional photon femtoscopy

This paper advocates for the use of the $C(\Delta E, Q_{\rm inv})$ kinematic variable over $C(Q_{\rm inv})$ for two-photon correlation functions in femtoscopy, highlighting its growing experimental feasibility due to advancements in detector technology.

The Gist

Imagine you are trying to take a photograph of a crowded dance floor to figure out how big the dancers are and how close they stand to each other. In the world of high-energy physics, scientists do something similar, but instead of dancers, they are looking at photons (particles of light) created when heavy atomic nuclei smash into each other at nearly the speed of light.

This technique is called Femtoscopy. It's like using the "dance moves" of particles to measure the size of the explosion that created them.

Here is the story of this paper, broken down into simple concepts:

1. The Problem: The "Ghost" Dancers

When two heavy nuclei collide, they create a fireball. Most of the light (photons) we see comes from the very beginning of this explosion. These are the "direct photons" we want to study.

However, there is a massive problem: Background Noise.
Most of the photons detected actually come from the decay of a particle called a pion (specifically $\pi^0$). Think of these as "ghost dancers" that show up later and mess up the photo. They dilute the signal, making it hard to see the real pattern of the direct photons.

Because there are so many of these "ghosts," the signal is very weak. Scientists have been struggling to get a clear picture because they don't have enough data (statistics) to separate the real dancers from the ghosts.

2. The Old Way: Looking at the Wrong Angle

For years, scientists tried to analyze this data using a one-dimensional map. They measured a single number called $Q_{inv}$ (Invariant Momentum Difference).

The Analogy:
Imagine you are trying to find two twins who are standing very close together in a crowd.

• The old method ($Q_{inv}$) was like measuring only the distance between them.
• The problem is that in a crowd, many unrelated people might happen to be the same distance apart by pure chance.
• For photons, this method mixes up the "real twins" (direct photons that are correlated) with "strangers" (uncorrelated photons).
• The Result: The peak of the signal (the evidence that the twins are together) gets flattened and washed out. It's like trying to hear a whisper in a noisy room; the old method made the whisper even quieter.

The authors of this paper argue that this old method is not optimal. It's like trying to solve a 3D puzzle by looking at a flat 2D shadow.

3. The New Idea: A Two-Dimensional Map

The authors propose a new way to look at the data. Instead of just measuring the distance, they suggest measuring two things at once:

1. The distance between the photons (still using $Q_{inv}$).
2. The difference in their energy ($\Delta E$).

The Analogy:
Now, imagine you are looking at the twins again, but this time you have a special camera that takes a 3D photo.

• You can see how close they are standing (distance).
• You can also see if they are wearing the exact same outfit (energy).
• Real "twin" photons (direct photons) tend to be close and have very similar energies.
• The "ghost" photons (background noise) might be close in distance, but their energies are all over the place.

By using this two-dimensional map ($C(\Delta E, Q_{inv})$), the scientists can separate the signal from the noise much better.

• The "ghosts" get pushed to the side of the map where they don't interfere.
• The "real twins" stay right in the center, forming a sharp, tall peak.

4. Why This Matters

• Sharper Vision: This new method preserves the full height of the signal peak. In the old method, the peak was squashed down, making it hard to measure.
• Better Geometry: By seeing the peak clearly, scientists can accurately measure the size of the "fireball" created in the collision. This tells us about the geometry of the universe at the moment of the Big Bang (or at least, the tiny version created in the lab).
• Counting the Ghosts: Interestingly, the amount the peak is squashed in the old method can actually tell us how many "ghost" photons there are. So, the new method helps us measure the signal and the background simultaneously.

The Conclusion

The paper is a call to action for experimental physicists. It says: "Stop using the old, blurry one-dimensional map. Switch to the new, sharp two-dimensional map."

With modern detectors getting faster and more sensitive, we finally have the technology to use this better method. It's like upgrading from a grainy black-and-white photo to a high-definition 3D video. We can finally see the structure of the photon-emitting system clearly, without the background noise blurring the picture.

In short: The authors found a better way to organize the data so that the "signal" stands out clearly against the "noise," allowing us to measure the size of the subatomic explosion with much greater precision.

Technical Summary

1. Problem Statement

Photon Femtoscopy and Statistical Limitations:
Femtoscopy (HBT interferometry) is a standard technique for measuring the space-time geometry of particle-emitting sources in high-energy nuclear collisions. While well-established for pions, its application to direct photons is challenging. Direct photons are produced early in the collision and do not undergo final-state interactions, offering a unique probe of the early stages. However, they are heavily contaminated by a background of decay photons (primarily from $\pi^0$).

The Core Issue:
The presence of decay photons dilutes the Bose-Einstein (BE) correlation signal. The amplitude of the BE peak is reduced by a factor of $f^2$ (where $f$ is the fraction of direct photons). Historically, experimental analyses have used a one-dimensional correlation function, $C(Q_{inv})$, where $Q_{inv}$ is the invariant momentum difference.

• The Flaw: The authors argue that $C(Q_{inv})$ is suboptimal for photons. Unlike massive particles, photons travel at the speed of light regardless of energy. Therefore, two photons emitted close in space-time but with different energies "fly together" in the lab frame.
• The Consequence: $Q_{inv}$ is primarily sensitive to the transverse ($q_{side}$) and longitudinal ($q_{long}$) momentum differences but is only weakly related to the "out" component ($q_{out}$), which is the direction of the pair momentum. Since the BE enhancement for photons is restricted to small $q_{out}$, projecting the data onto $Q_{inv}$ mixes correlated pairs (small $q_{out}$) with uncorrelated pairs (large $q_{out}$). This dilutes the BE peak amplitude, leading to systematic underestimations of the correlation strength ($\lambda$) and requiring large, model-dependent corrections.

2. Methodology

The paper proposes a shift in the kinematic variables used to construct the correlation function, moving from a 1D to a 2D representation.

Proposed Variables:
Instead of $C(Q_{inv})$, the authors advocate for a two-dimensional correlation function:
$$C(\Delta E, Q_{inv})$$
Where:

• $\Delta E = E_2 - E_1$ is the energy difference between the photon pair.
• $Q_{inv} = \sqrt{-(p_2 - p_1)^2}$ is the invariant momentum difference (equivalent to the invariant mass for photons).

Theoretical Justification:

• Vector Mapping: The authors demonstrate that the vector $(\Delta E, Q_{inv})$ is mathematically equivalent to the momentum difference vector in the Longitudinally Co-Moving System (LCMS). Specifically, $\Delta E$ maps directly to the "out" component ($q_{out}$), while $Q_{inv}$ maps to the transverse/longitudinal components ($\sqrt{q_{side}^2 + q_{long}^2}$).
• Preservation of Signal: By separating $\Delta E$ (which isolates the BE peak at $q_{out} \approx 0$) from $Q_{inv}$, the analysis avoids mixing uncorrelated pairs with large energy differences into the peak region.
• Background Separation: This 2D approach allows for better discrimination against backgrounds. For instance, $\pi^0$ decays produce photons with small opening angles (small $Q_{inv}$) but specific energy correlations, whereas anisotropic flow or detector effects (like cluster splitting) may manifest differently in the $\Delta E$ vs. $Q_{inv}$ plane.

Statistical Handling:
The authors address the concern that adding a dimension increases statistical noise. They argue that modern fitting techniques (Maximum Likelihood methods assuming Poissonian statistics) can handle low-count bins and empty bins effectively. Therefore, the 2D analysis does not incur a penalty in fit precision compared to 1D projections.

3. Key Contributions

1. Critique of $C(Q_{inv})$: The paper provides a rigorous argument that the standard 1D invariant mass difference representation dilutes the BE signal for photons because it fails to isolate the $q_{out}$ component where the correlation exists.
2. Proposal of $C(\Delta E, Q_{inv})$: It introduces a specific 2D kinematic representation that preserves the full height of the BE peak while maintaining the narrowness of the $\pi^0$ background peak (crucial for normalization).
3. Mathematical Equivalence: The authors derive the transformation equations showing that $C(\Delta E, Q_{inv})$ in the LCMS is practically equivalent to the standard 3D pion femtoscopy coordinates $C(q_{out}, \sqrt{q_{side}^2 + q_{long}^2})$ in the region of interest ($q < 0.1$ GeV/c).
4. Quantitative Impact: Through Monte Carlo simulations (Appendix), they demonstrate that the reduction of the peak strength $\lambda$ in $C(Q_{inv})$ is significant, particularly at high collision energies and high pair transverse momenta ($K_T$).

4. Results

• Simulation Findings: The Monte Carlo study (Fig. 2) shows that fitting a Gaussian to $C(Q_{inv})$ yields a significantly reduced correlation strength $\lambda$ compared to the true input value. This reduction worsens as the pair transverse momentum ($K_T$) and temperature ($T$) increase.
• Peak Preservation: In the proposed $C(\Delta E, Q_{inv})$ representation, the BE peak retains its full theoretical height because the integration over uncorrelated $q_{out}$ values is eliminated.
• Background Management: The 2D representation keeps the $\pi^0$ peak narrow (as $Q_{inv}$ is still the primary variable for the invariant mass), facilitating the normalization of the correlation function without the need for complex corrections.

5. Significance

• Revitalizing Photon Femtoscopy: With advancements in detector technology and electronics, the statistical limitations that previously hindered photon femtoscopy are being overcome. This paper provides the necessary theoretical framework to utilize this new data effectively.
• Accurate Geometry Extraction: By removing the dilution of the BE peak, this method allows for a more precise extraction of the source size ($R_{out}$) and the direct photon fraction ($f$) without relying on large, uncertain corrections.
• Standardization: The authors argue that the community should move away from the 1D $C(Q_{inv})$ standard (borrowed from pion physics) and adopt the 2D $C(\Delta E, Q_{inv})$ approach for photon pairs to ensure accurate physical interpretation of high-energy nuclear collisions.

In summary, the paper serves as a critical methodological guide, asserting that the choice of kinematic variables is paramount in photon femtoscopy and that a two-dimensional analysis is essential to unlock the full potential of direct photon correlations in studying the early stages of heavy-ion collisions.