EUROPEAN LABORATORY FOR PARTICLE PHYSICS
arXiv:1005.2855v1 [hep-ph] 17 May 2010
CERN–PH–EP / 2010–xxx
May 17, 2010
A new approach to
quarkonium polarization studies
Pietro Faccioli1) , Carlos Lourenço2) , and João Seixas1,3)
Abstract
Significant progress in understanding quarkonium production requires improved
polarization measurements, fully considering the intrinsic multidimensionality
of the problem. We propose a frame-invariant formalism which minimizes the
dependence of the measured result on the experimental acceptance, facilitates the
comparison with theoretical calculations, and provides a much needed control
over systematic effects due to detector limitations and analysis biases. This
formalism is a direct and generic consequence of the rotational invariance of the
dilepton decay distribution and is independent of any assumptions specific to
particular models of quarkonium production.
1)
2)
3)
Laboratório de Instrumentação e Fı́sica Experimental de Partı́culas (LIP),
Lisbon, Portugal
CERN, Geneva, Switzerland
Physics Department, Instituto Superior Técnico (IST), Lisbon, Portugal
Quarkonium polarization measurements should provide key information for the
understanding of quantum chromodynamics (QCD) [1], with the competing mechanisms dominating in the different theoretical approaches to quarkonium production
leading to very different predictions. In particular, non-relativistic QCD calculations [2], dominated by colour-octet components, predict that, at Tevatron energies
and asymptotically high pT , the directly produced J/ψ and ψ ′ mesons should be
almost fully transversely polarized (angular momentum component Jz = ±1) with
respect to their own momentum direction (helicity frame), while NLO calculations
of colour-singlet quarkonium production [3] indicate a strong longitudinal (Jz = 0)
polarization component.
The present experimental knowledge on quarkonium polarization is also contradictory and puzzling. The slightly longitudinal prompt-J/ψ polarization measured
by CDF [4] along the helicity axis (J/ψ momentum direction in the center-of-mass
of the collision system) is in clear disagreement with the expectations of all existing
models. The fact that these recent measurements disagree with the results previously
published by the same experiment [5] adds further confusion to the picture. Equally
disturbing is the lack of continuity between fixed-target and collider results [6]. Bottomonium polarization should be easier to interpret theoretically, given that the heavier bottom quark mass should satisfy the non-relativistic approximation much better
than the cc̄ case. However, the Υ(1S) Tevatron measurements [7], available in the
helicity frame and extending to pT values around 15 GeV/c, are contradictory: CDF
indicates unpolarized production; D0 indicates longitudinal polarization. The discrepancy between the two results cannot be justified by their different rapidity windows. At lower energy and pT , E866 [8] has shown yet a different polarization pattern:
the Υ(2S) and Υ(3S) states have maximal transverse polarization with respect to
the direction of motion of the colliding hadrons (Collins-Soper frame [9]), with no
significant dependence on transverse or longitudinal momentum, while the Υ(1S) is
only weakly polarized, indicating a dominant role of feed-down contributions.
To clarify this intricate situation, improved measurements are needed. So far,
most experiments have presented results based on a fraction of the physical information derivable from the data: only one polarization frame is used and only the polar
projection of the decay angular distribution is studied. These incomplete results
prevent model-independent physical conclusions. Moreover, such partial descriptions
of the observed physical processes reduce the chances of detecting possible biases
induced by insufficiently controlled systematic effects. The forthcoming LHC measurements, in particular, would benefit from an improved formalism. In Ref. [6] we
emphasized the importance of approaching the polarization measurement as a multidimensional problem, determining the full angular distribution in more than one
frame. In this letter we show the existence of a frame-independent relation among
the observable parameters of the dilepton decay angular distribution. The determination of invariant quantities facilitates the comparison between measurements,
and with theory, reducing the kinematic dependence of the results to their intrinsic
(and physically relevant) component. Invariant relations can also be used to perform
self-consistency checks which can expose unaccounted systematic effects or eventual
biases in the experimental analyses, a crucial advantage given the challenging nature
2
of quarkonium polarization measurements due, in particular, to the difficult subtraction of the spurious kinematic correlations induced by the detector acceptance.
We begin with preliminary considerations on the kinematics of the dilepton decay
of inclusively produced vector mesons. The most general observable distribution is
1
(1 + λϑ cos2 ϑ
(3 + λϑ )
+ λϕ sin2 ϑ cos 2ϕ + λϑϕ sin 2ϑ cos ϕ) ,
W (cos ϑ, ϕ) ∝
(1)
where ϑ and ϕ are the (polar and azimuthal) angles formed by the positive lepton
with, respectively, the polarization axis z and the production plane xz (containing the
directions of the colliding particles and of the decaying meson). All experimentally
definable polarization axes belong to the production plane. A transformation from
one observation frame (A) to another (B) is a rotation around the y axis by a δ angle,
the parameters changing as:
(A)
(A)
λϕ + Λ
λϑ − 3Λ
, λ(B)
,
ϕ =
1+Λ
1+Λ
(A)
(A)
(A)
λϑϕ cos 2δ − 21 (λϑ − λϕ ) sin 2δ
(B)
,
λϑϕ =
1+Λ
1 (A)
1 (A)
2
with Λ = (λϑ − λ(A)
λ sin 2δ .
ϕ ) sin δ −
2
2 ϑϕ
(B)
λϑ =
(2)
Since the magnitude of the “polar anisotropy”, λϑ , never exceeds 1 in any frame, we
deduce the frame-independent inequalities
|λϕ | ≤
1
(1 + λϑ ) ,
2
|λϑϕ | ≤
1
(1 − λϕ ) ,
2
(3)
which imply the bounds |λϕ | ≤ 1 and |λϑϕ | ≤ 1. More interestingly, we can see that
|λϕ | ≤ 0.5 when λϑ = 0 and must vanish when λϑ → −1.
In general, the transformation of the polarization parameters explicitly depends
on the production kinematics. Considering, for example, the Collins-Soper (CS) and
helicity (HX) frames, the angular terms in Eq. 2 are
sin2 δHX→CS = sin2 δCS→HX =
sin 2δHX→CS = − sin 2δCS→HX
p2T E 2
,
p2 (m2 + p2T )
2 m pT pL E
,
= 2 2
p (m + p2T )
(4)
where m, E, p, pT and pL are, respectively, the mass, the energy and the total, transverse and longitudinal momenta of the meson in the center-of-mass of the collision. As
a result, the observed quarkonium polarization has, in general, an “extrinsic”, framerelated, kinematic dependence, superimposed on the “intrinsic” physical dependence
due to the properties of the production processes and their varying mixture.
Such an extrinsic dependence can introduce artificial differences between the polarization results obtained by experiments probing different acceptance windows. Figure 1 shows how natural Υ polarizations λϑ = +1 and −1 in the CS frame (with
3
λHX
ϑ
λHX
ϑ
1
(a
1
0.5
0.5
0
-0.5
0
(d
-1
0
2
4
6
8
10
12
14
0
2
4
6
8
10
p [GeV/c]
12
14
p [GeV/c]
λHX
ϕ
T
λHX
ϕ
T
0.3
0
(e
-0.2
0.2
-0.4
-0.6
0.1
-0.8
(b
0
-1
0
2
4
6
8
10
12
14
0
2
4
6
8
10
p [GeV/c]
14
T
λHX
ϑϕ
λHX
ϑϕ
0.5
12
p [GeV/c]
T
0
(f
0.4
-0.2
0.3
-0.4
0.2
0.1
-0.6
(c
0
0
2
4
6
8
10
12
14
0
p [GeV/c]
2
4
6
8
10
12
14
p [GeV/c]
T
T
Figure 1: Kinematic dependence of the Υ(1S) decay angular distribution seen in
the HX frame, for natural polarizations λϑ = +1 (a-b-c) and λϑ = −1 (d-e-f) in the
CS frame. The curves correspond to different rapidity intervals; from the solid line:
|y| < 0.6 (CDF), |y| < 0.9 (ALICE), |y| < 1.8 (D0), |y| < 2.5 (ATLAS and CMS),
2 < |y| < 5 (LHCb). For simplicity the event populations were generated flat in
rapidity. The sign of λϑϕ depends on the definition of the y axis of the polarization
frame, here taken as sign(pL )(P~ ′ 1 × P~ ′ 2 )/|P~ ′1 × P~ ′ 2 |, where P~ ′ 1,2 are the momenta of
the colliding particles in the meson’s rest frame.
λϕ = λϑϕ = 0 and no intrinsic kinematic dependence) translate into different pT dependent polarizations measured in the HX frame in different rapidity acceptances.
In this simple case, a common choice of the “natural” frame (CS) by all experiments
would avoid such a misleading differentiation of results. In general, however, it may
be impossible to find one frame providing a simple representation of the quarkonium
polarization scenario. This is shown in Fig. 2, where we consider, for illustration,
that 60% of the Υ events have natural polarization λϑ = +1 in the CS frame and
the remaining fraction has λϑ = +1 in the HX frame. While the polarizations of the
two event subsamples are intrinsically independent of the production kinematics, in
neither frame will measurements performed in different transverse and longitudinal
momenta windows find identical results.
4
λCS
ϑ
λHX
ϑ
1
(a
1
0.5
0.5
0
0
(d
0
2
4
6
8
10
12
14
0
2
4
6
8
10
p [GeV/c]
12
14
p [GeV/c]
T
λϕ
λHX
ϕ
CS
T
0.3
(e
0.3
0.2
0.2
0.1
0.1
(b
0
0
0
2
4
6
8
10
12
14
0
2
4
6
8
10
p [GeV/c]
(c
14
T
λCS
ϑϕ
λHX
ϑϕ
0.5
12
p [GeV/c]
T
0
0.4
-0.1
0.3
-0.2
0.2
-0.3
0.1
-0.4
0
-0.5
(f
0
2
4
6
8
10
12
14
0
p [GeV/c]
2
4
6
8
10
12
14
p [GeV/c]
T
T
Figure 2: Kinematic dependence of the Υ(1S) decay angular distribution seen in the
HX (a-b-c) and CS (d-e-f) frames, when 60% (40%) of the events have full transverse
polarization in the CS (HX) frame. The curves represent measurements in different
acceptance ranges, as detailed in Fig. 1.
Also the comparison between experimental data and theory must consider the
experimental acceptance and efficiency. Experiments measure the net polarization
of the specific cocktail of quarkonium events accepted by the detector, trigger and
analysis cuts. If the polarization depends on the kinematics, the measured angular parameters depend on the effective population of collected events in the probed
phase space window. Two experiments covering the same kinematic interval may
find different average polarizations if they have different acceptance shapes in that
range. The problem can be solved by presenting the results in narrow intervals of
the probed phase space. Similarly, theoretical calculations of the average polarization
in a certain experiment should consider how the momentum distribution is distorted
by its acceptance. Alternatively, the prediction should avoid as much as possible
kinematic integrations or provide event-level information to be embedded into Monte
Carlo simulations of the experiments.
The general frame-transformation relations in Eq. 2 imply the existence of an
invariant quantity, definable in terms of λϑ , λϕ and λϑϕ , in one of the following
5
equivalent forms:
F{ci } =
(3 + λϑ ) + c1 (1 − λϕ )
.
c2 (3 + λϑ ) + c3 (1 − λϕ )
(5)
An account of the fundamental meaning of the frame-invariance of these quantities
can be found in Ref. [10]. We will consider here, specifically, the form
λ̃ ≡ F{−3,0,1} =
λϑ + 3λϕ
.
1 − λϕ
(6)
In the special case when the observed distribution is the superposition of n “elemen(i)
tary” distributions of the kind 1 + λϑ cos2 ϑ, with event weights f (i) , with respect to
n different polarization axes, λ̃ represents a weighted average of the n polarizations,
made irrespectively of the orientations of the corresponding axes:
λ̃ =
n
X
i=1
f (i)
3+
(i)
λ
(i) ϑ
λϑ
X
n
i=1
f (i)
(i)
3 + λϑ
.
(7)
The determination of an invariant quantity is immune to “extrinsic” kinematic dependencies induced by the observation perspective and is, therefore, less acceptancedependent than the anisotropy parameters λϑ , λϕ and λϑϕ . For instance, in the case
illustrated in Fig. 2, as well as in the simpler case of Fig. 1(a-b-c), any arbitrary
choice of the experimental observation frame will always yield the value λ̃ = 1, independently of kinematics. This particular case, where all contributing processes
are transversely polarized, is formally equivalent to the Lam-Tung relation [11], as
discussed in Ref. [10]. Analogously, the example represented in Fig. 1(d-e-f), or any
other case where all polarizations are longitudinal, yields λ̃ = −1.
The existence of frame-invariant parameters also provides a useful instrument for
experimental analyses. Checking, for example, that the same value of an invariant
quantity (Eq. 5) is obtained (within systematic uncertainties) in two distinct polarization frames is a non-trivial verification of the absence of unaccounted systematic
effects. In fact, detector geometry and/or data selection constraints strongly polarize
the reconstructed dilepton events. Background processes also affect the measured
polarization, if not well subtracted. The spurious anisotropies induced by detector
effects and background do not obey the frame transformation rules characteristic of
a physical J = 1 state. If not well corrected and subtracted, these effects will distort the shape of the measured decay distribution differently in different polarization
frames. In particular, they will violate the frame-independent relations between the
angular parameters. Any two physical polarization axes (defined in the rest frame of
the meson and belonging to the production plane) may be chosen to perform these
tests. The HX and CS frames are ideal choices at high pT , where they tend to be
orthogonal to each other (in Eq. 4, sin2 δ → 1 for pT ≫ m). At low pT , where the difference between the two frames vanishes, any of the two and its exact orthogonal may
be used to maximize the significance of the test. Given that λ̃ is “homogeneous” to
the anisotropy parameters, the difference λ̃(B) − λ̃(A) between the results obtained in
two frames provides a direct evaluation of the level of systematic errors not accounted
in the analysis.
6
We conclude with a summary of our messages. Choosing a given polarization
axis in an experimental analysis or theoretical calculation has a radical effect on
the magnitudes and signs of the polar and azimuthal anisotropies observed in the
dilepton decay distribution: all terms of the distribution must be determined to
provide unambiguously interpretable physical information. Rotational invariance imposes frame-invariant constraints on the polar and azimuthal anisotropy parameters
and, for any mixture of production mechanisms in a given kinematic condition, there
exists an invariant relation depending on one frame-independent parameter. Measurements of the anisotropy parameters are necessarily affected by the fact that the
transformation from one frame to another is an explicit function of the production
kinematics. These effects may result in a misleading interpretation of the measurements. Reporting polarization results in terms of invariant quantities facilitates the
comparison between different measurements, and with theory, reducing the kinematic
dependencies to their intrinsic (and physically most relevant) component. The invariant relation can also be used in the data analyses to perform self-consistency checks
which can expose unaccounted detector effects or eventual biases.
P.F. and J.S. acknowledge support from Fundação para a Ciência e a Tecnologia,
Portugal, under contracts SFRH/BPD/42343/2007 and CERN/FP/109343/2009.
References
[1] N. Brambilla et al. (QWG Coll.),
hep-ph/0412158, and references therein.
CERN Yellow Report 2005-005,
[2] E. Braaten, B.A. Kniehl and J. Lee, Phys. Rev. D 62, 094005 (2000), and
references therein.
[3] J.P. Lansberg, Eur. Phys. J. C 61 (2009) 693.
[4] A. Abulencia et al. (CDF Coll.), Phys. Rev. Lett. 99, 132001 (2007).
[5] T. Affolder et al. (CDF Coll.), Phys. Rev. Lett. 85, 2886 (2000).
[6] P. Faccioli, C. Lourenço, J. Seixas and H.K. Wöhri, Phys. Rev. Lett. 102, 151802
(2009).
[7] D. Acosta et al. (CDF Coll.), Phys. Rev. Lett. 88, 161802 (2002);
V.M. Abazov et al. (D0 Coll.), Phys. Rev. Lett. 101, 182004 (2008).
[8] C.N. Brown et al. (E866 Coll.), Phys. Rev. Lett. 86, 2529 (2001).
[9] J.C. Collins and D.E. Soper, Phys. Rev. D 16, 2219 (1977).
[10] P. Faccioli, C. Lourenço and J. Seixas, “Rotation-invariant relations in vector
meson decays into fermion pairs”, arXiv:1005.2601 [hep-ph].
[11] C.S. Lam and W.K. Tung, Phys. Rev. D 18, 2447 (1978);
C.S. Lam and W.K. Tung, Phys. Rev. D 21, 2712 (1980).
7