of 8
Collins asymmetries in inclusive charged
KK
and
K
π
pairs produced
in
e
þ
e
annihilation
J. P. Lees,
1
V. Poireau,
1
V. Tisserand,
1
E. Grauges,
2
A. Palano,
3a,3b
G. Eigen,
4
B. Stugu,
4
D. N. Brown,
5
L. T. Kerth,
5
Yu. G. Kolomensky,
5
M. J. Lee,
5
G. Lynch,
5
H. Koch,
6
T. Schroeder,
6
C. Hearty,
7
T. S. Mattison,
7
J. A. McKenna,
7
R. Y. So,
7
A. Khan,
8
V. E. Blinov,
9a,9b,9c
A. R. Buzykaev,
9a
V. P. Druzhinin,
9a,9b
V. B. Golubev,
9a,9b
E. A. Kravchenko,
9a,9b
A. P. Onuchin,
9a,9b,9c
S. I. Serednyakov,
9a,9b
Yu. I. Skovpen,
9a,9b
E. P. Solodov,
9a,9b
K. Yu. Todyshev,
9a,9b
A. J. Lankford,
10
B. Dey,
11
J. W. Gary,
11
O. Long,
11
M. Franco Sevilla,
12
T. M. Hong,
12
D. Kovalskyi,
12
J. D. Richman,
12
C. A. West,
12
A. M. Eisner,
13
W. S. Lockman,
13
W. Panduro Vazquez,
13
B. A. Schumm,
13
A. Seiden,
13
D. S. Chao,
14
C. H. Cheng,
14
B. Echenard,
14
K. T. Flood,
14
D. G. Hitlin,
14
T. S. Miyashita,
14
P. Ongmongkolkul,
14
F. C. Porter,
14
M. Röhrken,
14
R. Andreassen,
15
Z. Huard,
15
B. T. Meadows,
15
B. G. Pushpawela,
15
M. D. Sokoloff,
15
L. Sun,
15
P. C. Bloom,
16
W. T. Ford,
16
A. Gaz,
16
J. G. Smith,
16
S. R. Wagner,
16
R. Ayad,
17
,
W. H. Toki,
17
B. Spaan,
18
D. Bernard,
19
M. Verderi,
19
S. Playfer,
20
D. Bettoni,
21a
C. Bozzi,
21a
R. Calabrese,
21a,21b
G. Cibinetto,
21a,21b
E. Fioravanti,
21a,21b
I. Garzia,
21a,21b
E. Luppi,
21a,21b
L. Piemontese,
21a
V. Santoro,
21a
A. Calcaterra,
22
R. de Sangro,
22
G. Finocchiaro,
22
S. Martellotti,
22
P. Patteri,
22
I. M. Peruzzi,
22
,
M. Piccolo,
22
A. Zallo,
22
R. Contri,
23a,23b
M. R. Monge,
23a,23b
S. Passaggio,
23a
C. Patrignani,
23a,23b
B. Bhuyan,
24
V. Prasad,
24
A. Adametz,
25
U. Uwer,
25
H. M. Lacker,
26
U. Mallik,
27
C. Chen,
28
J. Cochran,
28
S. Prell,
28
H. Ahmed,
29
A. V. Gritsan,
30
N. Arnaud,
31
M. Davier,
31
D. Derkach,
31
G. Grosdidier,
31
F. Le Diberder,
31
A. M. Lutz,
31
B. Malaescu,
31
P. Roudeau,
31
A. Stocchi,
31
G. Wormser,
31
D. J. Lange,
32
D. M. Wright,
32
J. P. Coleman,
33
J. R. Fry,
33
E. Gabathuler,
33
D. E. Hutchcroft,
33
D. J. Payne,
33
C. Touramanis,
33
A. J. Bevan,
34
F. Di Lodovico,
34
R. Sacco,
34
G. Cowan,
35
D. N. Brown,
36
C. L. Davis,
36
A. G. Denig,
37
M. Fritsch,
37
W. Gradl,
37
K. Griessinger,
37
A. Hafner,
37
K. R. Schubert,
37
R. J. Barlow,
38
G. D. Lafferty,
38
R. Cenci,
39
B. Hamilton,
39
A. Jawahery,
39
D. A. Roberts,
39
R. Cowan,
40
R. Cheaib,
41
P. M. Patel,
41
,*
S. H. Robertson,
41
N. Neri,
42a
F. Palombo,
42a,42b
L. Cremaldi,
43
R. Godang,
43
,**
D. J. Summers,
43
M. Simard,
44
P. Taras,
44
G. De Nardo,
45a,45b
G. Onorato,
45a,45b
C. Sciacca,
45a,45b
G. Raven,
46
C. P. Jessop,
47
J. M. LoSecco,
47
K. Honscheid,
48
R. Kass,
48
M. Margoni,
49a,49b
M. Morandin,
49a
M. Posocco,
49a
M. Rotondo,
49a
G. Simi,
49a,49b
F. Simonetto,
49a,49b
R. Stroili,
49a,49b
S. Akar,
50
E. Ben-Haim,
50
M. Bomben,
50
G. R. Bonneaud,
50
H. Briand,
50
G. Calderini,
50
J. Chauveau,
50
Ph. Leruste,
50
G. Marchiori,
50
J. Ocariz,
50
M. Biasini,
51a,51b
E. Manoni,
51a
A. Rossi,
51a
C. Angelini,
52a,52b
G. Batignani,
52a,52b
S. Bettarini,
52a,52b
M. Carpinelli,
52a,52b
,
††
G. Casarosa,
52a,52b
M. Chrzaszcz,
52a
F. Forti,
52a,52b
M. A. Giorgi,
52a,52b
A. Lusiani,
52a,52c
B. Oberhof,
52a,52b
E. Paoloni,
52a,52b
M. Rama,
52a
G. Rizzo,
52a,52b
J. J. Walsh,
52a
D. Lopes Pegna,
53
J. Olsen,
53
A. J. S. Smith,
53
F. Anulli,
54a
R. Faccini,
54a,54b
F. Ferrarotto,
54a
F. Ferroni,
54a,54b
M. Gaspero,
54a,54b
A. Pilloni,
54a,54b
G. Piredda,
54a
C. Bünger,
55
S. Dittrich,
55
O. Grünberg,
55
M. Hess,
55
T. Leddig,
55
C. Voß,
55
R. Waldi,
55
T. Adye,
56
E. O. Olaiya,
56
F. F. Wilson,
56
S. Emery,
57
G. Vasseur,
57
D. Aston,
58
D. J. Bard,
58
C. Cartaro,
58
M. R. Convery,
58
J. Dorfan,
58
G. P. Dubois-Felsmann,
58
W. Dunwoodie,
58
M. Ebert,
58
R. C. Field,
58
B. G. Fulsom,
58
M. T. Graham,
58
C. Hast,
58
W. R. Innes,
58
P. Kim,
58
D. W. G. S. Leith,
58
S. Luitz,
58
V. Luth,
58
D. B. MacFarlane,
58
D. R. Muller,
58
H. Neal,
58
T. Pulliam,
58
B. N. Ratcliff,
58
A. Roodman,
58
R. H. Schindler,
58
A. Snyder,
58
D. Su,
58
M. K. Sullivan,
58
J. Va
vra,
58
W. J. Wisniewski,
58
H. W. Wulsin,
58
M. V. Purohit,
59
J. R. Wilson,
59
A. Randle-Conde,
60
S. J. Sekula,
60
M. Bellis,
61
P. R. Burchat,
61
E. M. T. Puccio,
61
M. S. Alam,
62
J. A. Ernst,
62
R. Gorodeisky,
63
N. Guttman,
63
D. R. Peimer,
63
A. Soffer,
63
S. M. Spanier,
64
J. L. Ritchie,
65
R. F. Schwitters,
65
J. M. Izen,
66
X. C. Lou,
66
F. Bianchi,
67a,67b
F. De Mori,
67a,67b
A. Filippi,
67a
D. Gamba,
67a,67b
L. Lanceri,
68a,68b
L. Vitale,
68a,68b
F. Martinez-Vidal,
69
A. Oyanguren,
69
J. Albert,
70
Sw. Banerjee,
70
A. Beaulieu,
70
F. U. Bernlochner,
70
H. H. F. Choi,
70
G. J. King,
70
R. Kowalewski,
70
M. J. Lewczuk,
70
T. Lueck,
70
I. M. Nugent,
70
J. M. Roney,
70
R. J. Sobie,
70
N. Tasneem,
70
T. J. Gershon,
71
P. F. Harrison,
71
T. E. Latham,
71
H. R. Band,
72
S. Dasu,
72
Y. Pan,
72
R. Prepost,
72
and S. L. Wu
72
(
B
A
B
AR
Collaboration)
1
Laboratoire d
Annecy-le-Vieux de Physique des Particules (LAPP), Université de Savoie, CNRS/IN2P3,
F-74941 Annecy-Le-Vieux, France
2
Universitat de Barcelona, Facultat de Fisica, Departament ECM, E-08028 Barcelona, Spain
3a
INFN Sezione di Bari, I-70126 Bari, Italy
3b
Dipartimento di Fisica, Università di Bari, I-70126 Bari, Italy
4
University of Bergen, Institute of Physics, N-5007 Bergen, Norway
5
Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720, USA
6
Ruhr Universität Bochum, Institut für Experimentalphysik 1, D-44780 Bochum, Germany
7
University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
8
Brunel University, Uxbridge, Middlesex UB8 3PH, United Kingdom
9a
Budker Institute of Nuclear Physics, SB RAS, Novosibirsk 630090, Russia
9b
Novosibirsk State University, Novosibirsk 630090, Russia
9c
Novosibirsk State Technical University, Novosibirsk 630092, Russia
PHYSICAL REVIEW D
92,
111101(R) (2015)
1550-7998
=
2015
=
92(11)
=
111101(8)
111101-1
© 2015 American Physical Society
RAPID COMMUNICATIONS
10
University of California at Irvine, Irvine, California 92697, USA
11
University of California at Riverside, Riverside, California 92521, USA
12
University of California at Santa Barbara, Santa Barbara, California 93106, USA
13
University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, California 95064, USA
14
California Institute of Technology, Pasadena, California 91125, USA
15
University of Cincinnati, Cincinnati, Ohio 45221, USA
16
University of Colorado, Boulder, Colorado 80309, USA
17
Colorado State University, Fort Collins, Colorado 80523, USA
18
Technische Universität Dortmund, Fakultät Physik, D-44221 Dortmund, Germany
19
Laboratoire Leprince-Ringuet, Ecole Polytechnique, CNRS/IN2P3, F-91128 Palaiseau, France
20
University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom
21a
INFN Sezione di Ferrara, I-44122 Ferrara, Italy
21b
Dipartimento di Fisica e Scienze della Terra, Università di Ferrara, I-44122 Ferrara, Italy
22
INFN Laboratori Nazionali di Frascati, I-00044 Frascati, Italy
23a
INFN Sezione di Genova, I-16146 Genova, Italy
23b
Dipartimento di Fisica, Università di Genova, I-16146 Genova, Italy
24
Indian Institute of Technology Guwahati, Guwahati, Assam 781 039, India
25
Universität Heidelberg, Physikalisches Institut, D-69120 Heidelberg, Germany
26
Humboldt-Universität zu Berlin, Institut für Physik, D-12489 Berlin, Germany
27
University of Iowa, Iowa City, Iowa 52242, USA
28
Iowa State University, Ames, Iowa 50011-3160, USA
29
Physics Department, Jazan University, Jazan 22822, Saudi Arabia
30
Johns Hopkins University, Baltimore, Maryland 21218, USA
31
Laboratoire de l
Accélérateur Linéaire, IN2P3/CNRS et Université Paris-Sud 11, Centre Scientifique
d
Orsay, F-91898 Orsay Cedex, France
32
Lawrence Livermore National Laboratory, Livermore, California 94550, USA
33
University of Liverpool, Liverpool L69 7ZE, United Kingdom
34
Queen Mary, University of London, London, E1 4NS, United Kingdom
35
University of London, Royal Holloway and Bedford New College,
Egham, Surrey TW20 0EX, United Kingdom
36
University of Louisville, Louisville, Kentucky 40292, USA
37
Johannes Gutenberg-Universität Mainz, Institut für Kernphysik, D-55099 Mainz, Germany
38
University of Manchester, Manchester M13 9PL, United Kingdom
39
University of Maryland, College Park, Maryland 20742, USA
40
Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge,
Massachusetts 02139, USA
41
McGill University, Montréal, Québec, Canada H3A 2T8
42a
INFN Sezione di Milano, I-20133 Milano, Italy
42b
Dipartimento di Fisica, Università di Milano, I-20133 Milano, Italy
43
University of Mississippi, University, Mississippi 38677, USA
44
Université de Montréal, Physique des Particules, Montréal, Québec H3C 3J7, Canada
45a
INFN Sezione di Napoli, I-80126 Napoli, Italy
45b
Dipartimento di Scienze Fisiche, Università di Napoli Federico II, I-80126 Napoli, Italy
46
NIKHEF, National Institute for Nuclear Physics and High Energy Physics,
NL-1009 DB Amsterdam, The Netherlands
47
University of Notre Dame, Notre Dame, Indiana 46556, USA
48
Ohio State University, Columbus, Ohio 43210, USA
49a
INFN Sezione di Padova, I-35131 Padova, Italy
49b
Dipartimento di Fisica, Università di Padova, I-35131 Padova, Italy
50
Laboratoire de Physique Nucléaire et de Hautes Energies, IN2P3/CNRS, Université Pierre et Marie
Curie-Paris 6, Université Denis Diderot-Paris 7, F-75252 Paris, France
51a
INFN Sezione di Perugia, I-06123 Perugia, Italy
51b
Dipartimento di Fisica, Università di Perugia, I-06123 Perugia, Italy
52a
INFN Sezione di Pisa, I-56127 Pisa, Italy
52b
Dipartimento di Fisica, Università di Pisa, I-56127 Pisa, Italy
52c
Scuola Normale Superiore di Pisa, I-56127 Pisa, Italy
53
Princeton University, Princeton, New Jersey 08544, USA
54a
INFN Sezione di Roma, I-00185 Roma, Italy
54b
Dipartimento di Fisica, Università di Roma La Sapienza I-00185 Roma, Italy
55
Universität Rostock, D-18051 Rostock, Germany
J. P. LEES
et al.
PHYSICAL REVIEW D
92,
111101(R) (2015)
111101-2
RAPID COMMUNICATIONS
56
Rutherford Appleton Laboratory, Chilton, Didcot, Oxon OX11 0QX, United Kingdom
57
CEA, Irfu, SPP, Centre de Saclay, F-91191 Gif-sur-Yvette, France
58
SLAC National Accelerator Laboratory, Stanford, California 94309 USA
59
University of South Carolina, Columbia, South Carolina 29208, USA
60
Southern Methodist University, Dallas, Texas 75275, USA
61
Stanford University, Stanford, California 94305-4060, USA
62
State University of New York, Albany, New York 12222, USA
63
Tel Aviv University, School of Physics and Astronomy, Tel Aviv 69978, Israel
64
University of Tennessee, Knoxville, Tennessee 37996, USA
65
University of Texas at Austin, Austin, Texas 78712, USA
66
University of Texas at Dallas, Richardson, Texas 75083, USA
67a
INFN Sezione di Torino, I-10125 Torino, Italy
67b
Dipartimento di Fisica, Università di Torino, I-10125 Torino, Italy
68a
INFN Sezione di Trieste, I-34127 Trieste, Italy
68b
Dipartimento di Fisica, Università di Trieste, I-34127 Trieste, Italy
69
IFIC, Universitat de Valencia-CSIC, E-46071 Valencia, Spain
70
University of Victoria, Victoria, British Columbia V8W 3P6, Canada
71
Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
72
University of Wisconsin, Madison, Wisconsin 53706, USA
(Received 23 June 2015; published 22 December 2015)
We present measurements of Collins asymmetries in the inclusive process
e
þ
e
h
1
h
2
X
,
h
1
h
2
¼
KK
,
K
π
,
ππ
, at the center-of-mass energy of 10.6 GeV, using a data sample of
468
fb
1
collected by the
BABAR
experiment at the PEP-II
B
factory at SLAC National Accelerator Center. Considering hadrons in opposite
thrust hemispheres of hadronic events, we observe clear azimuthal asymmetries in the ratio of unlike sign to
like sign, and unlike sign to all charged
h
1
h
2
pairs, which increase with hadron energies. The
K
π
asymmetries are similar to those measured for the
ππ
pairs, whereas those measured for high-energy
KK
pairs are, in general, larger.
DOI:
10.1103/PhysRevD.92.111101
PACS numbers: 13.66.Bc, 13.87.Fh, 13.88.+e, 14.65.
q
The Collins effect
[1]
relates the transverse spin com-
ponent of a fragmenting quark to the azimuthal distribution
of final state hadrons about its flight direction. The chiral-
odd, transverse momentum-dependent Collins fragmenta-
tion function (FF) provides a unique probe of QCD, such as
factorization and evolution with the energy scale
Q
2
[2
5]
.
Additional interest has been sparked by the observation
of azimuthal asymmetries for pions and kaons in semi-
inclusive deep inelastic scattering experiments (SIDIS)
[6
10]
. These are sensitive to the product of a Collins
FF and a chiral-odd transversity parton distribution func-
tion (PDF), one of the three fundamental PDFs needed to
describe the spin content of the nucleon. Although these
observations require nonzero Collins FFs, independent
direct measurements of one of these chiral-odd functions
are needed to determine each of them.
In
e
þ
e
annihilation, one can measure the product of
two Collins FFs, and detailed measurements have been
made for pairs of charged pions
[11
13]
. No measurements
are available for
K
π
and
KK
pairs, which are sensitive to
different quark-flavor combinations, in particular the con-
tribution of the strange quark. Such measurements could be
combined with SIDIS data to simultaneously determine the
Collins FFs and transversity PDF for up, down, and strange
quarks
[14
19]
.
In this paper, we report the measurement of the Collins
effect (or Collins asymmetry) for inclusive production of
hadron pairs in the process
e
þ
e
q
̄
q
h
1
h
2
X
, where
h
1
;
2
¼
K

or
π

,
q
stands for light quarks
u
or
d
or
s
, and
X
for any combination of additional hadrons.
The probability that a transversely polarized quark (
q
)
with momentum direction
ˆ
k
and spin
S
q
fragments into a
hadron
h
carrying zero intrinsic spin with momentum
P
h
is
defined in terms of unpolarized
D
q
1
and Collins
H
q
1
fragmentation functions
[20]
:
D
q
h
ð
z;
P
hT
Þ¼
D
q
1
ð
z; P
2
hT
Þþ
H
q
1
ð
z; P
2
hT
Þ
ð
ˆ
k
×
P
hT
Þ
·
S
q
zM
h
;
ð
1
Þ
*
Deceased.
Present address: University of Tabuk, Tabuk 71491, Saudi
Arabia.
Also at Università di Perugia, Dipartimento di Fisica, I-06123
Perugia, Italy.
§
Present address: Laboratoire de Physique Nucléaire et de
Hautes Energies, IN2P3/CNRS, F-75252 Paris, France.
Present address: University of Huddersfield, Huddersfield
HD1 3DH, UK.
**
Present address: University of South Alabama, Mobile,
Alabama 36688, USA.
††
Also at Università di Sassari, I-07100 Sassari, Italy.
COLLINS ASYMMETRIES IN INCLUSIVE CHARGED
KK
...
PHYSICAL REVIEW D
92,
111101(R) (2015)
111101-3
RAPID COMMUNICATIONS
where
M
h
,
P
hT
, and
z
¼
2
E
h
=
ffiffiffi
s
p
are the hadron mass,
momentum transverse to
ˆ
k
, and fractional energy, respec-
tively, with
E
h
its total energy and
ffiffiffi
s
p
the
e
þ
e
center-of-
mass (c.m.) energy. The term including
H
1
introduces an
azimuthal modulation around the direction of the fragment-
ing quark, called Collins asymmetry.
In
e
þ
e
q
̄
q
events, the
q
and
̄
q
must be produced
back to back in the
e
þ
e
c.m. frame with their spin aligned.
For unpolarized
e
þ
and
e
beams at
BABAR
energies, the
q
and
̄
q
spins are polarized along either the
e
þ
or
e
beam
direction, so there is a large transverse component when the
angle between the
e
þ
e
and the
q
̄
q
axis is large. The
direction is unknown for a given event, but the correlation
can be exploited. Experimentally, the
q
and
̄
q
directions are
difficult to measure, but the event thrust axis
ˆ
n
[21,22]
approximates at leading order the
q
̄
q
axis, so an azimuthal
correlation between two hadrons in opposite thrust hemi-
spheres reflects the product of the two Collins functions.
Figure
1
shows the thrust reference frame (RF12)
[23]
.If
not otherwise specified, all kinematic variables are defined
in the
e
þ
e
c.m. frame. The Collins effect results in a
cosine modulation of the azimuthal angle
φ
12
¼
φ
1
þ
φ
2
of
the dihadron yields. Expressing the yield as a function of
φ
12
(after the integration over
P
hT
), and dividing by the
average bin content, we obtain the normalized rate
[11]
R
12
ð
φ
12
Þ¼
1
þ
sin
2
θ
th
1
þ
cos
2
θ
th
cos
φ
12
·
H
½
1

1
ð
z
1
Þ
̄
H
½
1

1
ð
z
2
Þ
D
½
0

1
ð
z
1
Þ
̄
D
½
0

1
ð
z
2
Þ
;
ð
2
Þ
where the sum over the involved quark flavors is implied,
θ
th
is defined in Fig.
1
,
z
1
ð
2
Þ
is the fractional energy of the
first (second) hadron, and the bar denotes the function for
the
̄
q
. Equation
(2)
involves only the moments of FF, which
are defined as
F
½
n

ð
z
i
Þ
Z
d
j
k
T
j
2

j
k
T
j
M
i

n
F
ð
z
i
;
j
k
T
j
2
Þ
;
ð
3
Þ
with
n
¼
0
, 1, and
j
k
T
j
the transverse momentum of the
quarks with respect to the hadrons they fragment into,
which, in this frame, is related to the measurement of the
transverse momenta of the two hadrons with respect to the
thrust axis.
Despite the simple form of the
R
12
normalized rate,
which involves only the product of moments of FFs, the
RF12 frame comes with several downsides, among others
of having to rely on Monte Carlo (MC) simulations when
using the thrust axis as a proxy for the leading-order
q
̄
q
axis. An alternative frame is the analogue of the Gottfried-
Jackson frame
[23,24]
which uses the momentum of one
hadron as a reference axis, and defines a single angle
φ
0
between the plane containing the two hadron momenta and
the plane defined by the beam and the reference axis. We
refer to this frame as RF0
[11,12]
. The corresponding
normalized yield in the
e
þ
e
c.m. system is
[23]
R
0
ð
2
φ
0
Þ¼
1
þ
sin
2
θ
2
1
þ
cos
2
θ
2
cos
2
φ
0
·
F
½ð
2
ˆ
h
·
k
T
ˆ
h
·
p
T
k
T
·
p
T
Þ
H
1
̄
H
1

ð
M
1
M
2
Þ
F
½
D
1
̄
D
1

;
ð
4
Þ
where
θ
2
is the angle between the hadron used as reference
and the beam axis
ˆ
h
is the unit vector in the direction of the
transverse momentum of the first hadron relative to the axis
defined by the second hadron, and
F
is used to denote the
convolution integral
F
½
X
̄
X

X
q
e
2
q
Z
d
2
k
T
d
2
p
T
δ
2
ð
p
T
þ
k
T
q
T
Þ
×
X
q
ð
z
1
;z
2
1
k
2
T
Þ
̄
X
q
ð
z
2
;z
2
2
p
2
T
Þ
;
ð
5
Þ
with
k
T
,
p
T
, and
q
T
the transverse momentum of the
fragmenting quark, antiquark, and virtual photon from
e
þ
e
annihilation, respectively, in the frame where the
two hadrons are collinear, and
X
ð
̄
X
Þ
D
1
ð
̄
D
1
Þ
or
H
1
ð
̄
H
1
Þ
. In this frame, specific assumptions on the
k
T
dependence of the involved functions are necessary to
explicitly evaluate the convolution integrals.
For this analysis we use a data sample of
468
fb
1
[25]
collected at the c.m. energy
ffiffiffi
s
p
10
.
6
GeV with the
BABAR
detector
[26,27]
at the SLAC National
Accelerator Laboratory. We use tracks reconstructed in
the silicon vertex detector and in the drift chamber (DCH)
and identified as pions or kaons in the DCH and in the
Cherenkov ring imaging detector (DIRC). Detailed MC
simulation is used to study detector effects and to estimate
contribution from various background sources. Hadronic
FIG. 1 (color online). Thrust reference frame (RF12). The
azimuthal angles
φ
1
and
φ
2
are the angles between the scattering
plane and the transverse hadron momenta
p
t
1
ð
t
2
Þ
around the thrust
axis
ˆ
n
. The polar angle
θ
th
is the angle between
ˆ
n
and the beam
axis. Note that the difference between
p
t
1
ð
t
2
Þ
and
P
hT
is that the
latter is calculated with respect to the
q
̄
q
axis.
J. P. LEES
et al.
PHYSICAL REVIEW D
92,
111101(R) (2015)
111101-4
RAPID COMMUNICATIONS
events are generated using the
J
etset
[28]
package and
undergo a full detector simulation based on G
EANT
4
[29]
.
We make a tight selection of hadronic events in order to
minimize biases due to detector acceptance and hard initial-
state photon radiation (ISR), as they can introduce fake
azimuthal modulations. Furthermore, final-state gluon (
q
̄
qg
)
radiation also leads to angular asymmetries to be taken into
account
[23]
. Requiring at least three charged tracks con-
sistentwiththe
e
þ
e
primaryvertexanda totalvisibleenergy
of the event in the laboratory frame
E
tot
>
11
GeV, we reject
e
þ
e
τ
þ
τ
and two-photon backgrounds, as well as ISR
(
q
̄
qg
) events with the photon (one jet) along the beam line.
About 10% of ISR photons are within our detector accep-
tance, and we reject events with a photon candidate with
energy above 2 GeV. We require an event thrust value
T>
0
.
8
to suppress
q
̄
qg
and
B
̄
B
events, and
j
cos
θ
th
j
<
0
.
6
so
that most tracks are within the detector acceptance.
We assign randomly the positive direction of the thrust
axis, and divide each event into two hemispheres by the
plane perpendicular to it. To ensure tracks are assigned to
the correct hemispheres, we require them to be within a 45°
angle of the thrust axis and to have
z>
0
.
15
.A
tight
identification algorithm is used to identify kaons (pions),
which is about 80% (90%) efficient and has misidentifi-
cation rates below 10% (5%). We select those pions and
kaons that lie within the DIRC acceptance region with a
polar angle in laboratory frame
0
.
45
rad
<
θ
lab
<
2
.
46
rad.
To minimize backgrounds, such as
e
þ
e
μ
þ
μ
γ
fol-
lowed by photon conversion, we require
z<
0
.
9
.
We construct all the possible pairs of selected tracks
reconstructed in opposite thrust hemispheres, and we
calculate the corresponding azimuthal angles
φ
1
,
φ
2
, and
φ
0
in the respective reference frames. In this way, we
identify three different samples of hadron pairs:
KK
,
K
π
,
and
ππ
. To reduce low-energy gluon radiation and the
contribution due to wrong hemispheres assignment, we
require
Q
t
<
3
.
5
GeV
=c
, where
Q
t
is the transverse
momentum of the virtual photon from
e
þ
e
annihilation
in the frame where the two hadrons are collinear
[23]
.
The analysis is performed in intervals of hadron frac-
tional energies with the following boundaries: 0.15, 0.2,
0.3, 0.5, 0.9, for a total of 16 two-dimensional
ð
z
1
;z
2
Þ
intervals.
For each of the three samples, we evaluate the normal-
ized yield distributions
R
12
and
R
0
for unlike (
U
), like (
L
),
and any charge combination (
C
) of hadron pairs as a
function of
φ
1
þ
φ
2
and
2
φ
0
, as shown in the left plot of
Fig.
2
for
KK
pairs, for example. These combinations of
charged hadrons contain different contributions of favored
and disfavored FFs, where a favored (disfavored) process
refers to the production of a hadron for which one (none) of
the valence quarks is of the same kind as the fragmenting
quark. In particular, by selecting
KK
pairs, we are able to
study the favored contribution
H
fav
s
of the strange quark,
not accessible when considering
ππ
pairs only.
The normalized distributions can be parametrized with a
cosine function:
R
i
α
¼
b
α
þ
a
i
α
cos
β
α
, where
α
¼
0
,
12
indicates the reference frames,
i
¼
U
,
L
,
C
the charge
combination of hadron pairs, and
β
12
ð
0
Þ
¼
φ
12
ð
2
φ
0
Þ
.
The
R
i
α
distributions are strongly affected by instrumen-
tal effects. In order to reduce the impact of the detector
acceptance, as well as any remaining effect from gluon
bremsstrahlung
[23]
, we construct two double ratios (DRs)
of normalized distributions,
R
U
α
=R
L
α
and
R
U
α
=R
C
α
. The two
ratios give access to the same physical quantities as the
independent
R
i
α
, that is the favored and disfavored FFs, but
in different combinations. We report the results for both
kinds of DRs, which are strongly correlated since they are
obtained by using the same data set. These are shown in the
right plot of Fig.
2
for
KK
pairs in RF12. At first order, the
double ratios are still parametrized by a function that is
linear in the cosine of the corresponding combination of
azimuthal angles:
R
ij
α
¼
R
i
α
R
j
α
B
ij
α
þ
A
ij
α
· cos
β
α
;
ð
6
Þ
with
B
and
A
free parameters, and
i
,
j
¼
U
,
L
,
C
. The
constant term
B
must be consistent with unity, while
A
contains the information about the favored and disfavored
Collins FFs.
We fit the binned
R
ij
α
distributions independently for
KK
,
K
π
, and
ππ
hadron pairs. Using the MC sample, we
evaluate the
K=
π
(mis)identification probabilities for the
16
ð
z
1
;z
2
Þ
intervals in each of the three samples. For
example, the probability
f
KK
KK
that a true
KK
pair is
reconstructed as
KK
pair is about 90% on average, slightly
decreasing at higher momenta, while the probability
f
KK
K
π
that a true
K
π
pair is identified as
KK
is about 10%, and
f
KK
ππ
is negligible.
The presence of background processes could introduce
azimuthal modulations not related to the Collins effect, and
modifies the measured asymmetry as follows:
2
φ
+
1
φ
-3
-2
-1
0
1
2
3
R
0.96
0.98
1
1.02
1.04
12
U
R
12
L
R
12
C
R
2
φ
+
1
φ
-3
-2
-1
0
1
2
3
12
UL
R
12
UC
R
FIG. 2 (color online). Distributions of normalized yields (left
plot) for unlike (
U
), like (
L
), and any charge combination (
C
)of
KK
pairs, and their double ratios (right plot) in RF12.
COLLINS ASYMMETRIES IN INCLUSIVE CHARGED
KK
...
PHYSICAL REVIEW D
92,
111101(R) (2015)
111101-5
RAPID COMMUNICATIONS
A
meas
KK
¼
F
KK
uds
·

X
nm
f
KK
nm
·
A
nm

þ
X
i
F
KK
i

X
nm
f
ð
KK
Þ
i
nm
·
A
i
nm

;
ð
7
Þ
with
nm
¼
KK
,
K
π
,
ππ
, and
i
¼
c
̄
c
,
B
̄
B
,
τ
þ
τ
. In Eq.
(7)
,
A
nm
are the true Collins asymmetries produced from the
fragmentation of light quarks in the three samples,
A
i
nm
is
the
i
th background asymmetry contribution, and
F
KK
uds
ð
i
Þ
are
the fractions of reconstructed kaon pairs coming from
uds
and background events, calculated from the respective MC
samples. By construction,
P
i
F
i
þ
F
uds
¼
1
. A similar
expression holds for
K
π
and
ππ
samples.
Previous studies
[11]
show that
e
þ
e
B
̄
B
and
τ
þ
τ
events have negligible
A
i
nm
,
F
B
̄
B
<
2%
, and
F
τ
þ
τ
signifi-
cantly different from zero only for the
ππ
sample at high
z
values. Since
F
c
̄
c
can be as large as 30%, and
A
c
̄
c
are
unknown, we determine
A
c
̄
c
nm
in Eq.
(7)
from samples
enhanced in
c
̄
c
by requiring the reconstruction of at least
one
D

meson from the decay
D

D
0
π

, with the
D
0
candidate reconstructed in the following four Cabibbo-
favored decay modes:
K
π
þ
,
K
π
þ
π
π
þ
,
K
0
s
π
þ
π
, and
K
π
þ
π
0
. These modes are assumed to provide a represen-
tative sample of
ππ
,
K
π
, and
KK
pairs to be used in the
correction, an assumption that is strengthened by the
observation that the background asymmetries for those
modes were found to be consistent. We solve the system of
equations for
A
meas
KK
,
A
meas
K
π
,
A
meas
ππ
, for the standard and
charm-enhanced samples, and we extract simultaneously
the Collins asymmetries
A
KK
,
A
K
π
, and
A
ππ
, corrected
for the contributions of the background and
K=
π
(mis)
identification. The dominant uncertainties related to this
procedure come from the limited statistics of the
D

-
enhanced sample and from the fractions
F
i
. The uncer-
tainties on the fractions are evaluated by data-MC com-
parison and amount to a few percent. All these uncertainties
are therefore included in the statistical error of the asym-
metries extracted from the system of Eq.
(7)
.
We test the DR method on the MC sample. Spin effects
are not simulated in MC, and so the DR distributions
should be uniform. However, when fitting the distributions
for reconstructed
KK
pairs with Eq.
(6)
, we measure a
cosine term in the full sample of
0
.
004

0
.
001
and
0
.
007

0
.
001
in the RF12 and RF0 frames, respectively,
indicating a bias. Smaller values are obtained for
K
π
and
ππ
pairs
[30]
. Studies performed on the MC samples, both
at generation level and after full simulation, demonstrate
that the main source of this bias is due to the emission of
ISR, which boosts the hadronic system and distorts the
angular distribution of the final state particles, resulting in
azimuthal modulations not related to the Collins effect.
This effect is more pronounced for
KK
pairs due to the
lower multiplicity with respect to the other two combina-
tions of hadrons. Assuming the bias, which is everywhere
smaller than the asymmetries measured in the data sample
in each bin, is additive, we subtract it from the background-
corrected asymmetry.
Using the
uds
MC sample, or light quark
e
þ
e
q
̄
q
MC events, we study the difference between measured and
true azimuthal asymmetries. The asymmetry is introduced
into the simulation by reweighting the events according to
the distribution
1

a
· cos
φ
gen
α
, where we use different
values of
a
ranging from 0 to 8% with positive (negative)
sign for
U
(
L
and
C
) hadron pairs, and
φ
gen
α
are the
azimuthal angles combinations calculated with respect to
the true
q
̄
q
axis in RF12, or the generated hadron
momentum in RF0. The reconstructed asymmetries in
RF12 are systematically underestimated for the three
samples of hadron pairs, as expected since we use the
thrust axis instead of the
q
̄
q
axis, while they are consistent
with the simulated ones in RF0, where only particle
identification and tracking reconstruction effects could
introduce possible dilution. Since we measure the same
dilution for
KK
,
K
π
, and
ππ
samples, the asymmetry is
corrected by rescaling
A
KK
,
A
K
π
, and
A
ππ
using the same
correction factor, which ranges from 1.3 to 2.3 increasing
with
z
, as shown in Fig.
3
. No corrections are needed for the
asymmetries measured in RF0. The uncertainties on the
correction factors are assigned as systematic contributions.
All systematic effects, if not otherwise specified, are
evaluated for each bin of
z
. The main contribution comes
from the MC bias. We compare the bias results from the
nominal selection, with those obtained by requiring differ-
ent cuts on
E
tot
, and/or by changing the detector acceptance
region for the hadrons. The largest variation of the bias is
combined in quadrature with the MC statistical error and
BIN
)
2
,z
1
(z
0
2
4
6
8
10
12
14
16
Correction factor
1
1.2
1.4
1.6
1.8
2
2.2
2.4
UL
12
A
UC
12
A
UL
0
A
UC
0
A
0.2
0.3
0.15
0.5
0.2
0.3
0.5
0.2
0.3
0.5
0.9
0.2
0.3
0.5
0.15
0.2
0.3
0.5
0.9
1
z
2
z
0.9/0.15
0.9/0.15
0.9/0.15
FIG. 3 (color online). Correction factors for the dilution of the
asymmetry due to the difference between the thrust and the
q
̄
q
axis. The open (full) markers, triangles and circles, refer to the
U=L
and
U=C
double ratios in the RF12 (RF0) frame, respec-
tively. The 16 (
z
1
,
z
2
) bins are shown on the
x
axis: in each
interval between the dashed lines,
z
1
is chosen in the following
ranges: [0.15, 0.2], [0.2, 0.3], [0.3, 0.5], and [0.5, 0.9], while
within each interval the points correspond to the four bins in
z
2
.
J. P. LEES
et al.
PHYSICAL REVIEW D
92,
111101(R) (2015)
111101-6
RAPID COMMUNICATIONS
taken as systematic uncertainty. The effects due to the
particle identification are evaluated using tighter and looser
selection criteria. The largest deviations with respect to the
nominal selection are taken as systematic uncertainties: the
average relative uncertainties are around 10%, 7%, and 5%
for the
KK
,
K
π
, and
ππ
pairs. Fitting the azimuthal
distributions using different bin sizes, we determine relative
systematic uncertainties, which are not larger than 5%,
1.9%, and 1% for the three samples. The systematic
uncertainty due to the
E
tot
cut is obtained by comparing
the measured asymmetries with those obtained with the
looser selection
E
tot
>
10
GeV. The average systematic
contribution is around 10% for the three samples in both
reference frames. We use different fitting functions with
additional higher harmonic terms. No significant changes in
the value of the cosine moments with respect to the standard
fits are found. As a cross-check of the double ratio method
we fit the difference of
R
i
distributions, and we compare the
two results. The difference between the two procedures is
negligible for
K
π
and
ππ
pairs, while it reaches 1% and 3%
for kaon pairs in RF12 and RF0, respectively. All the other
systematic contributions are negligible
[11]
.
The Collins asymmetries measured for the 16 two-
dimensional
ð
z
1
;z
2
Þ
bins, for reconstructed
KK
,
K
π
, and
ππ
hadron pairs, are shown in Fig.
4
for RF12 and RF0, and
are summarized in tables reported in the Supplemental
Material
[30]
. The asymmetries are corrected for the
background contributions and
K=
π
contamination follow-
ing Eq.
(7)
, the MC bias is subtracted, and the corrections
due to the dilution effects are applied. The total systematic
uncertainties are obtained by adding in quadrature the
individual contributions, and are represented by the bands
around the data points.
An increasing asymmetry with increasing hadron ener-
gies is visible for the
U=L
double ratio in both reference
frames. The largest effects, but with less precision, are
observed for
KK
pairs, for which
A
UL
12
is consistent with zero
at low
z
, and reaches 22% in the last
z
bin, while somewhat
smaller values are seen for
ππ
and
K
π
pairs. In particular, at
low
ð
z
1
;z
2
Þ
bins
A
UL
for
ππ
pairs is nonzero, in agreement
with the behavior observed in
[11]
. The small differences
between the two data sets are due to the different kinematic
region selected after the cut on cos
θ
th
. The
A
UC
asymmetry
is smaller than
A
UL
in all cases, and, for the
KK
pairs, the rise
of the asymmetry with the hadron energies is not evident. We
also note that the asymmetries for the
KK
pairs are larger
than the others when the
U=L
ratio is considered, while they
are at the same level, or lower, when they are extracted from
the
U=C
ratio.
In summary, we have studied for the first time in
e
þ
e
annihilation the Collins asymmetry for inclusive produc-
tion of
KK
and
K
π
pairs as a function of
ð
z
1
;z
2
Þ
in two
distinct reference frames. We measure the azimuthal
modulation of the double ratios
U=L
and
U=C
, which
are sensitive to the favored and disfavored Collins FFs for
light quarks. We simultaneously extract also the Collins
asymmetries for
ππ
pairs, which are found to be in
agreement with those obtained in previous studies
[11,13]
. The results reported in this paper and those
obtained from SIDIS experiments can be used in a global
analysis to extract the favored contribution of the strange
quark, and to improve the knowledge on the
u
and
d
fragmentation processes
[14
16]
.
We are grateful for the excellent luminosity and machine
conditions provided by our PEP-II2 colleagues, and for the
0 2 4 6 8 10121416
UL
12
A
0
0.05
0.1
0.15
0.2
KK
π
K
ππ
BIN
)
2
,z
1
(z
0246810121416
UC
12
A
0
0.02
0.04
0.06
0.08
0.2
0.3
0.15
0.5
0.2
0.3
0.5
0.2
0.3
0.5
0.9
0.2
0.3
0.5
0.15
0.2
0.3
0.5
0.9
1
z
2
z
0.9/0.15
0.9/0.15
0.9/0.15
0246810121416
UL
0
A
0
0.02
0.04
0.06
0.08
KK
π
K
ππ
BIN
)
2
,z
1
(z
0246810121416
UC
0
A
-0.01
0
0.01
0.02
0.03
0.04
0.2
0.3
0.15
0.5
0.2
0.3
0.5
0.2
0.3
0.5
0.9
0.2
0.3
0.5
0.15
0.2
0.3
0.5
0.9
1
z
2
z
0.9/0.15
0.9/0.15
0.9/0.15
FIG. 4 (color online). Comparison of
U=L
(top) and
U=C
(bottom) Collins asymmetries in RF12 (left) and RF0 (right) for
KK
,
K
π
,
and
ππ
pairs. The statistical and systematic uncertainties are represented by the bars and the bands around the points, respectively. The
16
ð
z
1
;z
2
Þ
bins are shown on the
x
axis: in each interval between the dashed lines,
z
1
is chosen in the following ranges: [0.15, 0.2], [0.2,
0.3], [0.3, 0.5], and [0.5, 0.9], while within each interval the points correspond to the four bins in
z
2
.
COLLINS ASYMMETRIES IN INCLUSIVE CHARGED
KK
...
PHYSICAL REVIEW D
92,
111101(R) (2015)
111101-7
RAPID COMMUNICATIONS
substantial dedicated effort from the computing organiza-
tions that support
BABAR
. The collaborating institutions
wish to thank SLAC for its support and kind hospitality.
This work is supported by DOE and NSF (U.S.), NSERC
(Canada), CEA and CNRS-IN2P3 (France), BMBF and
DFG (Germany), INFN (Italy), FOM (The Netherlands),
NFR (Norway), MES (Russia), MEC (Spain), and STFC
(U.K.). Individuals have received support from the Marie
Curie EIF (European Union) and the A. P. Sloan
Foundation.
[1] J. C. Collins,
Nucl. Phys.
B396
, 161 (1993)
.
[2] J. C. Collins and D. E. Soper,
Nucl. Phys.
B193
, 381 (1981)
.
[3] J. Collins, D. E. Soper, and G. Sterman,
Nucl. Phys.
B250
,
199 (1985)
.
[4] J. C. Collins and A. Metz,
Phys. Rev. Lett.
93
, 252001
(2004)
.
[5] P. Sun and F. Yuan,
Phys. Rev. D
88
, 114012 (2013)
.
[6] A. Airapetian
et al.
(HERMES Collaboration),
Phys. Rev.
Lett.
94
, 012002 (2005)
.
[7] A. Airapetian
et al.
(HERMES Collaboration),
Phys. Lett. B
693
, 11 (2010)
.
[8] C. Adolph
et al.
(COMPASS Collaboration),
Phys. Lett. B
717
, 376 (2012)
.
[9] M. Alekseev
et al.
(COMPASS Collaboration),
Phys. Lett.
B
673
, 127 (2009)
.
[10] A. Airapetian
et al.
(HERMES Collaboration),
Phys. Rev. D
87
, 012010 (2013)
.
[11] J. Lees
et al.
(
BABAR
Collaboration),
Phys. Rev. D
90
,
052003 (2014)
.
[12] R. Seidl, G. Perdekamp
et al.
(Belle Collaboration),
Phys.
Rev. D
78
, 032011 (2008)
.
[13] R. Seidl, G. Perdekamp
et al.
(Belle Collaboration),
Phys.
Rev. D
86
, 039905(E) (2012)
.
[14] A. Bacchetta, L. P. Gamberg, G. R. Goldstein, and A.
Mukherjee,
Phys. Lett. B
659
, 234 (2008)
.
[15] M. Anselmino, M. Boglione, U. D
Alesio, A. Kotzinian, F.
Murgia, A. Prokudin, and C. Türk,
Phys. Rev. D
75
, 054032
(2007)
.
[16] M. Anselmino, M. Boglione, U. D
Alesio, S. Melis,
F. Murgia, and A. Prokudin,
Phys. Rev. D
87
, 094019
(2013)
.
[17] R. L. Jaffe and X. Ji,
Phys. Rev. Lett.
67
, 552
(1991)
.
[18] J. Soffer,
Phys. Rev. Lett.
74
, 1292 (1995)
.
[19] A. Martin, F. Bradamante, and V. Barone,
Phys. Rev. D
91
,
014034 (2015)
.
[20] A. Bacchetta, U. D
Alesio, M. Diehl, and C. A. Miller,
Phys. Rev. D
70
, 117504 (2004)
.
[21] E. Farhi,
Phys. Rev. Lett.
39
, 1587 (1977)
.
[22] S. Brandt, C. Peyrou, R. Sosnowski, and A. Wroblewski,
Phys. Lett.
12
, 57 (1964)
.
[23] D. Boer,
Nucl. Phys.
B806
, 23 (2009)
.
[24] D. Boer, R. Jakob, and P. Mulders,
Nucl. Phys.
B504
, 345
(1997)
.
[25] J. P. Lees
et al.
(
BABAR
Collaboration),
Nucl. Instrum.
Methods Phys. Res., Sect. A
726
, 203 (2013)
.
[26] B. Aubert
et al.
(
BABAR
Collaboration),
Nucl. Instrum.
Methods Phys. Res., Sect. A
479
, 1 (2002)
.
[27] B. Aubert
et al.
(
BABAR
Collaboration),
Nucl. Instrum.
Methods Phys. Res., Sect. A
729
, 615 (2013)
.
[28] T. Sjöstrand,
arXiv:hep-ph/9508391.
[29] S. Agostinelli
et al.
(GEANT4 Collaboration),
Nucl. Ins-
trum. Methods Phys. Res., Sect. A
506
, 250 (2003)
.
[30] See Supplemental Material at
http://link.aps.org/
supplemental/10.1103/PhysRevD.92.111101
for tables
and figures.
J. P. LEES
et al.
PHYSICAL REVIEW D
92,
111101(R) (2015)
111101-8
RAPID COMMUNICATIONS