Drell-Yan Angular Analysis

• The paper extracts angular coefficients (A0, A2, etc.) using a reweighted maximum-likelihood fit to rigorously test QCD predictions and the Lam–Tung relation.
• Drell–Yan angular analysis is a quantitative study of lepton-pair distributions that reveals key aspects of partonic dynamics and gluon spin behavior.
• The measurement confirms the Lam–Tung relation, supporting the vector nature of gluons and constraining models of QCD spin-correlation and parton-level processes.

The Drell-Yan angular analysis is the quantitative study of angular distributions of lepton pairs produced via the annihilation of quark and antiquark in hadron–hadron collisions (specifically $p\bar{p}\to\gamma^*/Z\to e^+e^-+X$ at $\sqrt{s}=1.96$ TeV, as executed by the CDF II experiment). It provides a direct probe of the partonic-level dynamics, the vector nature of the gluon, QCD radiative effects, and spin–correlation phenomena, by extracting a set of angular coefficients $A_i$ from the observed lepton angular spectrum in the Collins–Soper frame. The complete measurement disentangles the polarization states and their modulations in transverse momentum $P_T$, enabling a rigorous test of the Lam–Tung relation and the underlying theory.

1. Collins–Soper Angular Decomposition and Coefficient Definitions

In the Collins–Soper rest frame, the polar angle $\theta$ and azimuthal angle $\phi$ of the $e^+$ (or $e^-$) are defined with respect to the bisector of the beam directions. The differential cross section is expressed as: $$\frac{d\sigma}{d\Omega} \propto (1+\cos^2\theta) + \frac{A_0}{2}(1-3\cos^2\theta) + A_1\sin2\theta\cos\phi + \frac{A_2}{2}\sin^2\theta\cos2\phi + A_3\sin\theta\cos\phi + A_4\cos\theta + A_5\sin^2\theta\sin2\phi + A_6\sin2\theta\sin\phi + A_7\sin\theta\sin\phi$$ Each coefficient $A_i$ parameterizes a specific physical correlation:

• $\sqrt{s}=1.96$0: Longitudinal versus transverse polarization admixture of $\sqrt{s}=1.96$1
• $\sqrt{s}=1.96$2: $\sqrt{s}=1.96$3 azimuthal modulation from spin–spin and orbital angular momentum
• $\sqrt{s}=1.96$4: $\sqrt{s}=1.96$5 term (spin–flip interference)
• $\sqrt{s}=1.96$6, $\sqrt{s}=1.96$7: Parity-violating terms from $\sqrt{s}=1.96$8–$\sqrt{s}=1.96$9 interference (forward–backward and $A_i$0 asymmetry)
• $A_i$1, $A_i$2, $A_i$3: T-odd/higher-twist terms (predicted to vanish at leading order in perturbative QCD)

The extraction of these coefficients provides a minimal set necessary for full angular characterization.

2. Experimental Realization: CDF II Detector and Event Selection

CDF II employed a multi-layer tracking system (SVX, COT) in a solenoidal magnetic field, projective-tower EM/HAD calorimetry, and high-$A_i$4 triggers. Drell–Yan events were selected in the $A_i$5 mass window $A_i$6 GeV/c$A_i$7 and divided into CC, CP, and PP topologies according to electron rapidities and tracker associations. Electrons were required to exceed $A_i$8 thresholds and be matched to SVX tracks; QCD and electroweak backgrounds were carefully subtracted (total background $A_i$9). The dataset comprised $P_T$0140,000 events over $P_T$1 fb$P_T$2 of integrated luminosity.

3. Angular Coefficient Extraction and Systematics

The analysis binned events in five $P_T$3 intervals ($P_T$4 from $P_T$5 to $P_T$6 GeV/c):

• $P_T$7, $P_T$8, $P_T$9, $\theta$0, $\theta$1 GeV A re-weighted maximum-likelihood fit was performed on the $\theta$2 distributions using PYTHIA+GEANT simulated events. Each $\theta$3 bin’s data was fitted for $\theta$4 and $\theta$5 (from $\theta$6), and $\theta$7, $\theta$8 (from $\theta$9), floating normalizations. Acceptance, efficiency, resolution, and remaining backgrounds were fully modeled. Systematic uncertainties (electron ID/tracking, background subtraction, material) were propagated through repeated fits, found to be subdominant to statistical errors.

Measured coefficients ($\phi$0, $\phi$1, $\phi$2, $\phi$3 in units of $\phi$4):

$\phi$5 bin (GeV)	$\phi$6	$\phi$7	$\phi$8	$\phi$9	$e^+$0
0–10	4.8	0.17±0.14±0.07	0.16±0.26±0.06	-0.04±0.12±0.01	1.10±0.10±0.01
10–20	14.1	0.42±0.25±0.07	-0.01±0.35±0.16	0.18±0.16±0.01	1.01±0.17±0.01
20–35	26.0	0.86±0.39±0.08	0.52±0.51±0.29	0.14±0.24±0.01	1.56±0.26±0.01
35–55	42.9	3.11±0.59±0.10	2.88±0.84±0.19	-0.19±0.41±0.04	0.52±0.42±0.03
>55	73.7	4.97±0.61±0.10	4.83±1.24±0.02	-0.47±0.56±0.02	0.85±0.50±0.05

4. Lam–Tung Relation and Validation

The Lam–Tung relation, $e^+$1, is a consequence of the spin-1 nature of gluon emission in both $e^+$2 and $e^+$3 subprocesses. Experimental data yield

$e^+$4 = $e^+$5 $e^+$6 $e^+$7 $e^+$8 $e^+$9 across $e^-$0 bins, with average $e^-$1. This result is robustly consistent with zero, confirming the Lam–Tung relation to within $e^-$2 at high $e^-$3, and directly verifying that the gluon behaves as a vector boson in hadronic collisions.

5. Comparison with Theory and QCD Mechanisms

The measured $e^-$4 and $e^-$5 at low $e^-$6 conform to the annihilation formula: $e^-$7, while at high $e^-$8 they exceed this, evidencing $e^-$9 Compton scattering contributions. Comparison with predictions from NLO generators (DYRAD, POWHEG, FEWZ, MG+PS), matched parton shower models (PYTHIA+1jet), and NNLO corrections indicate excellent agreement with the data, whereas LO-only implementations and resummation-only generators (VBP) underestimate $$\frac{d\sigma}{d\Omega} \propto (1+\cos^2\theta) + \frac{A_0}{2}(1-3\cos^2\theta) + A_1\sin2\theta\cos\phi + \frac{A_2}{2}\sin^2\theta\cos2\phi + A_3\sin\theta\cos\phi + A_4\cos\theta + A_5\sin^2\theta\sin2\phi + A_6\sin2\theta\sin\phi + A_7\sin\theta\sin\phi$$0 and $$\frac{d\sigma}{d\Omega} \propto (1+\cos^2\theta) + \frac{A_0}{2}(1-3\cos^2\theta) + A_1\sin2\theta\cos\phi + \frac{A_2}{2}\sin^2\theta\cos2\phi + A_3\sin\theta\cos\phi + A_4\cos\theta + A_5\sin^2\theta\sin2\phi + A_6\sin2\theta\sin\phi + A_7\sin\theta\sin\phi$$1 at large $$\frac{d\sigma}{d\Omega} \propto (1+\cos^2\theta) + \frac{A_0}{2}(1-3\cos^2\theta) + A_1\sin2\theta\cos\phi + \frac{A_2}{2}\sin^2\theta\cos2\phi + A_3\sin\theta\cos\phi + A_4\cos\theta + A_5\sin^2\theta\sin2\phi + A_6\sin2\theta\sin\phi + A_7\sin\theta\sin\phi$$2.

6. Implications for Gluon Spin, QCD Spin-Correlation, and Global Fits

The precise verification of $$\frac{d\sigma}{d\Omega} \propto (1+\cos^2\theta) + \frac{A_0}{2}(1-3\cos^2\theta) + A_1\sin2\theta\cos\phi + \frac{A_2}{2}\sin^2\theta\cos2\phi + A_3\sin\theta\cos\phi + A_4\cos\theta + A_5\sin^2\theta\sin2\phi + A_6\sin2\theta\sin\phi + A_7\sin\theta\sin\phi$$3 rules out a sizable scalar-gluon component, confirming the vector nature of the gluon. The $$\frac{d\sigma}{d\Omega} \propto (1+\cos^2\theta) + \frac{A_0}{2}(1-3\cos^2\theta) + A_1\sin2\theta\cos\phi + \frac{A_2}{2}\sin^2\theta\cos2\phi + A_3\sin\theta\cos\phi + A_4\cos\theta + A_5\sin^2\theta\sin2\phi + A_6\sin2\theta\sin\phi + A_7\sin\theta\sin\phi$$4 dependence of $$\frac{d\sigma}{d\Omega} \propto (1+\cos^2\theta) + \frac{A_0}{2}(1-3\cos^2\theta) + A_1\sin2\theta\cos\phi + \frac{A_2}{2}\sin^2\theta\cos2\phi + A_3\sin\theta\cos\phi + A_4\cos\theta + A_5\sin^2\theta\sin2\phi + A_6\sin2\theta\sin\phi + A_7\sin\theta\sin\phi$$5 and $$\frac{d\sigma}{d\Omega} \propto (1+\cos^2\theta) + \frac{A_0}{2}(1-3\cos^2\theta) + A_1\sin2\theta\cos\phi + \frac{A_2}{2}\sin^2\theta\cos2\phi + A_3\sin\theta\cos\phi + A_4\cos\theta + A_5\sin^2\theta\sin2\phi + A_6\sin2\theta\sin\phi + A_7\sin\theta\sin\phi$$6 quantitatively demonstrates the role and importance of quark–gluon Compton scattering ($$\frac{d\sigma}{d\Omega} \propto (1+\cos^2\theta) + \frac{A_0}{2}(1-3\cos^2\theta) + A_1\sin2\theta\cos\phi + \frac{A_2}{2}\sin^2\theta\cos2\phi + A_3\sin\theta\cos\phi + A_4\cos\theta + A_5\sin^2\theta\sin2\phi + A_6\sin2\theta\sin\phi + A_7\sin\theta\sin\phi$$7), as predicted by pQCD at Tevatron energies. The measurement constrains parton-level angular correlations employed in global QCD fits and bolsters the validation of Monte Carlo tools incorporating NLO+parton shower matching. Future analyses at the LHC and other higher-energy experiments will further elucidate the balance of $$\frac{d\sigma}{d\Omega} \propto (1+\cos^2\theta) + \frac{A_0}{2}(1-3\cos^2\theta) + A_1\sin2\theta\cos\phi + \frac{A_2}{2}\sin^2\theta\cos2\phi + A_3\sin\theta\cos\phi + A_4\cos\theta + A_5\sin^2\theta\sin2\phi + A_6\sin2\theta\sin\phi + A_7\sin\theta\sin\phi$$8 vs. $$\frac{d\sigma}{d\Omega} \propto (1+\cos^2\theta) + \frac{A_0}{2}(1-3\cos^2\theta) + A_1\sin2\theta\cos\phi + \frac{A_2}{2}\sin^2\theta\cos2\phi + A_3\sin\theta\cos\phi + A_4\cos\theta + A_5\sin^2\theta\sin2\phi + A_6\sin2\theta\sin\phi + A_7\sin\theta\sin\phi$$9 initial-state contributions and probe possible small-$A_i$0 effects or higher-twist corrections in the deviation from leading-twist behavior.

7. Extended Applications and Outlook

The comprehensive approach employed in the CDF angular analysis is extensible to higher energies and alternative kinematic regimes. Measurements of all leading Drell–Yan angular coefficients at high mass and $A_i$1 provide critical inputs for the study of QCD spin correlations, the development of resummation techniques, and the extraction of electroweak parameters. The established agreement between data and theory signifies that perturbative QCD, supplemented by sophisticated simulation and fit techniques, remains a powerful framework for describing and predicting spin-dependent structures in neutral-current Drell–Yan processes (Collaboration et al., 2011).