Background-Free Rapidity Asymmetry

• Background-free rapidity asymmetry is a method that isolates antisymmetric rapidity structures by canceling symmetric contributions using techniques like factorial moments and rapidity gaps.
• It provides robust observables in heavy-ion collision studies, including measures from long-range correlations, saturation-scale fluctuations, and jet-hadron analyses.
• The approach refines hydrodynamic models and signal extractions by enabling unbiased comparisons of forward-backward asymmetries while mitigating Poisson noise and detector effects.

Background-free rapidity asymmetry denotes a family of observables that isolate an antisymmetric structure in rapidity while removing, canceling, or strongly suppressing contaminating contributions such as Poisson counting noise, short-range dynamical correlations, uncorrelated bulk background, or explicit signal-background admixture. In the literature, the concept appears in several technically distinct settings: event-by-event rapidity-density asymmetry extracted from long-range correlations in heavy-ion collisions, odd rapidity modes generated by saturation-scale fluctuations, rapidity-asymmetry signals of jet-induced diffusion wake, and forward-backward asymmetry measurements in which signal and background are separated statistically rather than geometrically (Bialas et al., 2011, Bzdak et al., 2016, Yang et al., 6 Jan 2025, Yang et al., 6 Aug 2025, Pretz et al., 2008).

1. Observable classes and common logic

The central idea is always to construct an odd-in-rapidity quantity whose symmetric backgrounds either vanish identically or are removed by the measurement procedure. In symmetric collisions, this often means projecting onto the first odd rapidity mode. In jet-hadron observables, it means comparing samples that differ only by the rapidity of the associated jet while holding the trigger rapidity fixed or symmetric. In signal-extraction problems, it means weighting or subtracting background contributions so that the remaining asymmetry is attributable to the signal component.

Context	Observable	Background treatment
Long-range rapidity correlations	$D_-^2=\tfrac12[F_{20}-F_{11}]$	factorial moments remove Poisson noise; forward-backward gap suppresses short-range correlations
Saturation-scale fluctuations	$c_n\{n\}$ from $a_1$	in symmetric collisions, background-symmetric pieces drop out
Jet-hadron correlations	$A^{near,away}(\Delta\eta;Y)$ or $A_y(y)$	subtraction between rapidity classes cancels uncorrelated background
Longitudinal asymmetry	ratio of $dN_+/dy$ and $dN_-/dy$	ratio of two classes cancels many detector inefficiencies and most backgrounds

A recurring misconception is that “background-free” has a single operational meaning. The cited work uses the phrase in several non-equivalent senses: exact cancellation of an identical background, suppression of short-range sources by a rapidity gap, removal of Poisson noise through factorial moments, or simultaneous determination of signal and background asymmetries through optimal weighting. This suggests that the phrase is methodological rather than universal.

2. Long-range rapidity correlations and hidden asymmetry

A particularly direct construction starts from the single-particle rapidity density $\rho(y)$, defined so that in a small rapidity interval $\Delta y$ around $y$ the average number of particles is $c_n\{n\}$0. The full two-particle density is written as $c_n\{n\}$1, and the connected correlation is

$c_n\{n\}$2

After dividing the rapidity axis into bins $c_n\{n\}$3, one assumes that for a fixed underlying density profile $c_n\{n\}$4 the occupancies $c_n\{n\}$5 fluctuate independently and Poisson-distributed about their means. Under this “Poissonian” assumption, the factorial moments

$c_n\{n\}$6

directly equal moments of the dynamical densities in each bin, $c_n\{n\}$7 (Bialas et al., 2011).

For a symmetric collision at $c_n\{n\}$8, the event-by-event density can be expanded on $c_n\{n\}$9 in Legendre polynomials:

$a_1$0

with $a_1$1. Even $a_1$2 correspond to symmetric fluctuations, whereas odd $a_1$3 correspond to antisymmetric fluctuations and therefore measure violation of boost invariance event by event. The two-particle cumulant decomposes as

$a_1$4

The first odd mode is isolated by the dynamical variance of the forward-backward difference. Using two narrow bins centered at $a_1$5 and separated by a rapidity gap, with occupancies $a_1$6 and $a_1$7, one defines

$a_1$8

where $a_1$9 and $A^{near,away}(\Delta\eta;Y)$0. Because the bins are separated by $A^{near,away}(\Delta\eta;Y)$1, short-range correlations such as resonance decays or jet fragmentation are strongly suppressed, while factorial moments remove pure statistical fluctuations. In this sense, $A^{near,away}(\Delta\eta;Y)$2 is a direct, background-free measurement of rapidity-density asymmetry (Bialas et al., 2011).

Applied to STAR Au+Au data at $A^{near,away}(\Delta\eta;Y)$3 with bins at $A^{near,away}(\Delta\eta;Y)$4 and $A^{near,away}(\Delta\eta;Y)$5, the reported values for top-central collisions are

$A^{near,away}(\Delta\eta;Y)$6

which give

$A^{near,away}(\Delta\eta;Y)$7

With $A^{near,away}(\Delta\eta;Y)$8, the asymmetric component is clearly nonzero. The interpretation is that the rapidity profile in single events is not an even function of $A^{near,away}(\Delta\eta;Y)$9, so boost invariance is violated on the event-by-event basis even at central rapidity. The paper further notes that this effect may seriously influence the hydrodynamic expansion of the system (Bialas et al., 2011).

3. Odd rapidity modes from saturation-scale fluctuations

A closely related but more explicitly microscopic formalism writes the event-by-event single-particle distribution as

$A_y(y)$0

with $A_y(y)$1 by construction. Dropping $A_y(y)$2 and higher even terms isolates the antisymmetric fluctuation in the coefficient $A_y(y)$3. In a KLN-type model with proton saturation scales $A_y(y)$4 and $A_y(y)$5, one finds in the limit of small $A_y(y)$6-fluctuations

$A_y(y)$7

where $A_y(y)$8 (Bzdak et al., 2016).

The connected two-particle rapidity correlation satisfies

$A_y(y)$9

and in symmetric collisions all background-symmetric pieces such as $dN_+/dy$0 and $dN_+/dy$1 drop out. More generally, the genuine $dN_+/dy$2-particle cumulant obeys

$dN_+/dy$3

The first cumulants are

$dN_+/dy$4

with analogous expressions given for $dN_+/dy$5 and $dN_+/dy$6 (Bzdak et al., 2016).

The underlying fluctuation model takes

$dN_+/dy$7

so that $dN_+/dy$8. For small $dN_+/dy$9,

$dN_-/dy$0

In the color-domain picture, the variance of the mean $dN_-/dy$1 decreases with the number of domains, giving

$dN_-/dy$2

and therefore

$dN_-/dy$3

This was reported to be in excellent accord with the preliminary ATLAS result $dN_-/dy$4 (Bzdak et al., 2016).

For higher-order cumulants, the small-$dN_-/dy$5 behavior is

$dN_-/dy$6

up to overall factors of order one. Since $dN_-/dy$7 decreases with multiplicity, the paper predicts that $dN_-/dy$8 and $dN_-/dy$9 change sign at an intermediate multiplicity where $\rho(y)$0. This sign change is presented as a sharp, background-free signature of genuine saturation-scale-fluctuation-induced rapidity asymmetry (Bzdak et al., 2016).

4. Jet-hadron rapidity asymmetry and diffusion wake

In jet-medium studies, the relevant object is the per-trigger jet-hadron correlation in pseudorapidity and azimuth, defined relative to the leading jet by

$\rho(y)$1

with $\rho(y)$2 or $\rho(y)$3, $\rho(y)$4, and $\rho(y)$5. Projection onto near and away regions yields $\rho(y)$6 and $\rho(y)$7. The medium-induced modification is

$\rho(y)$8

and the rapidity asymmetry for a large gap $\rho(y)$9 relative to a reference small gap $\Delta y$0 is

$\Delta y$1

This observable was proposed as a robust signal of the diffusion wake induced by di-jets in Pb+Pb collisions (Yang et al., 6 Jan 2025).

Two background-removal schemes were presented. In the theoretical subtraction, CoLBT-hydro is run twice event by event, with and without the embedded di-jet, and the hydro-only charged-hadron spectrum is subtracted bin-by-bin from the jet+hydro spectrum. In the mixed-event subtraction, one forms the same-event pair distribution $\Delta y$2 and a mixed-event reference $\Delta y$3, with the trigger jet paired with charged hadrons from a different event of the same centrality class, and defines $\Delta y$4. After either subtraction, the projected asymmetries are very similar in the simulation (Yang et al., 6 Jan 2025).

For centrality 0–10% Pb+Pb at $\Delta y$5, with leading jet $\Delta y$6, sub-leading jet $\Delta y$7, $\Delta y$8, and soft hadrons $\Delta y$9, the predicted near-side asymmetry for $y$0 shows a negative dip of about $y$1 per trigger at $y$2 and a positive excess of $y$3 near $y$4. For $y$5, the dip moves to $y$6 with amplitude about $y$7. The away-side asymmetry has a similar shape but a larger peak-trough difference of about $y$8 per trigger. The physical interpretation is that a finite rapidity separation shifts the depletion region created by the diffusion wake of one jet relative to the enhancement associated with the other jet (Yang et al., 6 Jan 2025).

A later formulation generalizes the construction to dijets and $y$9-jets and removes the background by design rather than by subtraction. With trigger rapidity fixed or restricted to a symmetric interval around $c_n\{n\}$00, one defines the opposite-side hadron rapidity distribution

$c_n\{n\}$01

and the rapidity asymmetry

$c_n\{n\}$02

If

$c_n\{n\}$03

with the background $c_n\{n\}$04 independent of $c_n\{n\}$05, then

$c_n\{n\}$06

so the background cancels exactly and no mixed-event subtraction or event-plane correction is required (Yang et al., 6 Aug 2025).

In CoLBT-hydro simulations of Pb+Pb at $c_n\{n\}$07, the $c_n\{n\}$08-jet asymmetry for $c_n\{n\}$09 exhibits a valley at $c_n\{n\}$10 for forward jets and a peak at $c_n\{n\}$11 for central jets, while for $c_n\{n\}$12 the asymmetry is negligible. In dijets, selecting soft hadrons gives a negative dip at $c_n\{n\}$13 and a positive peak at $c_n\{n\}$14; again, the asymmetry tends to zero for $c_n\{n\}$15–$c_n\{n\}$16. Subtracting a pp baseline has a negligible effect on $c_n\{n\}$17 (Yang et al., 6 Aug 2025).

5. Longitudinal asymmetry and participant-zone rapidity shift

In collisions of identical nuclei, the numbers of participants from the two nuclei can differ event by event. Writing these as $c_n\{n\}$18 and $c_n\{n\}$19, the participant asymmetry is

$c_n\{n\}$20

and the participant-zone rapidity shift is

$c_n\{n\}$21

For small $c_n\{n\}$22,

$c_n\{n\}$23

This longitudinal asymmetry causes a shift in the center-of-mass rapidity of the participant zone and therefore modifies the final-state rapidity distribution (Raniwala et al., 2016).

The analysis strategy is to bin events by $c_n\{n\}$24, or by the experimentally accessible spectator asymmetry, compute the average $c_n\{n\}$25 in each bin, and fit

$c_n\{n\}$26

Analytically one expects $c_n\{n\}$27, $c_n\{n\}$28, and $c_n\{n\}$29, while finite binning, fluctuations in $c_n\{n\}$30, and details of the rapidity profile generate effective deviations. Toy simulations with shifted Gaussian rapidity distributions show that the ratio of two distributions shifted by $c_n\{n\}$31 expands into a polynomial whose linear coefficient directly measures $c_n\{n\}$32. Full event generators—HIJING, AMPT default, AMPT String-Melting, and DPMJET—were then used to fit the ratio of positive- and negative-asymmetry samples with a third-order polynomial in $c_n\{n\}$33, with the linear coefficient $c_n\{n\}$34 remaining dominant (Raniwala et al., 2016).

Experimentally, one cannot count $c_n\{n\}$35 and $c_n\{n\}$36 directly, but one can measure spectator neutrons in the two Zero-Degree Calorimeters and define

$c_n\{n\}$37

A Monte Carlo Glauber response matrix is then used to unfold $c_n\{n\}$38 or directly $c_n\{n\}$39. The final-state rapidity or pseudorapidity distributions for $c_n\{n\}$40 and $c_n\{n\}$41 are ratioed and fit with a third-order polynomial. Because the measurement is a ratio of two event classes, many detector inefficiencies, acceptance effects, and most backgrounds cancel, yielding an effectively background-free measure of the participant-zone rapidity shift (Raniwala et al., 2016).

The construction is not entirely model-independent. The paper states that the mapping $c_n\{n\}$42 relies on Glauber-MC tuning, that ZDC resolution is poor in very central collisions, and that the small nonlinear coefficients are sensitive to the shape of $c_n\{n\}$43. This suggests that the observable is robust at the level of cancellation but still calibration-dependent.

6. Statistical extraction, collider applications, and interpretive limits

A general statistical framework for background-free rapidity asymmetry treats signal and background asymmetries simultaneously. Let $c_n\{n\}$44 be a discriminating variable and $c_n\{n\}$45 the rapidity. If $c_n\{n\}$46 and $c_n\{n\}$47 are the signal and background shapes, and $c_n\{n\}$48 for $c_n\{n\}$49 and $c_n\{n\}$50 for $c_n\{n\}$51, one defines per-event weights

$c_n\{n\}$52

In the small-asymmetry limit, the weighted moments satisfy a $c_n\{n\}$53 linear system for the signal and background asymmetries. When one assumes $c_n\{n\}$54, the optimal estimator reduces to

$c_n\{n\}$55

For vanishing asymmetries, this estimator reaches the minimal variance bound given by the Cramér-Rao inequality (Pretz et al., 2008).

A complementary treatment includes explicit background subtraction, acceptance-luminosity asymmetry, and unfolding. In a binned forward-backward analysis with counts $c_n\{n\}$56, background estimates $c_n\{n\}$57, relative acceptance$c_n\{n\}$58luminosity factor $c_n\{n\}$59, and normalization $c_n\{n\}$60, the asymmetry estimator is

$c_n\{n\}$61

An unbinned extended-likelihood form replaces discrete bin counts by a signal PDF $c_n\{n\}$62 and background PDF $c_n\{n\}$63, while detector smearing is handled through a response matrix $c_n\{n\}$64 or a folded reconstruction-level PDF. This formulation makes the separation between background subtraction and unfolding explicit (Pate et al., 2 Feb 2026).

A collider application with a different meaning of “background-free” is the Tevatron top-antitop forward-backward asymmetry. There, the rapidity difference is

$c_n\{n\}$65

or in reconstructed lepton+jets events,

$c_n\{n\}$66

and

$c_n\{n\}$67

Matrix-element-level studies showed that the $c_n\{n\}$68 jets background has a positive slope in $c_n\{n\}$69, with $c_n\{n\}$70 for $c_n\{n\}$71 and $c_n\{n\}$72 for $c_n\{n\}$73, while the leading-order $c_n\{n\}$74 signal sample has $c_n\{n\}$75 by construction. The proposed remedies were tighter kinematic cuts, template fits in $c_n\{n\}$76, and bin-by-bin subtraction using the predicted background fractions and asymmetries (Hagiwara et al., 2012).

A plausible implication is that “background-free rapidity asymmetry” spans a spectrum of constructions. At one end are observables such as $c_n\{n\}$77 in which the background cancels exactly by symmetry and sample selection. At the other are statistically optimal extractions in which a background model remains necessary but the asymmetry estimator is unbiased, variance-efficient, or explicitly corrected for background and detector response. Across these uses, the common target is the same: an odd rapidity structure that can be interpreted dynamically rather than as an artifact of counting noise, combinatorial contamination, or acceptance imbalance.