Search for
CP
violation in the decays
D
!
K
0
S
K
,
D
s
!
K
0
S
K
, and
D
s
!
K
0
S
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
G. Lynch,
5
H. Koch,
6
T. Schroeder,
6
D. J. Asgeirsson,
7
C. Hearty,
7
T. S. Mattison,
7
J. A. McKenna,
7
R. Y. So,
7
A. Khan,
8
V. E. Blinov,
9
A. R. Buzykaev,
9
V. P. Druzhinin,
9
V. B. Golubev,
9
E. A. Kravchenko,
9
A. P. Onuchin,
9
S. I. Serednyakov,
9
Yu. I. Skovpen,
9
E. P. Solodov,
9
K. Yu. Todyshev,
9
A. N. Yushkov,
9
D. Kirkby,
10
A. J. Lankford,
10
M. Mandelkern,
10
B. Dey,
11
J. W. Gary,
11
O. Long,
11
G. M. Vitug,
11
C. Campagnari,
12
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
A. J. Martinez,
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
P. Ongmongkolkul,
14
F. C. Porter,
14
A. Y. Rakitin,
14
R. Andreassen,
15
Z. Huard,
15
B. T. Meadows,
15
M. D. Sokoloff,
15
L. Sun,
15
P. C. Bloom,
16
W. T. Ford,
16
A. Gaz,
16
U. Nauenberg,
16
J. G. Smith,
16
S. R. Wagner,
16
R. Ayad,
17,
*
W. H. Toki,
17
B. Spaan,
18
K. R. Schubert,
19
R. Schwierz,
19
D. Bernard,
20
M. Verderi,
20
P. J. Clark,
21
S. Playfer,
21
D. Bettoni,
22a
C. Bozzi,
22a
R. Calabrese,
22a,22b
G. Cibinetto,
22a,22b
E. Fioravanti,
22a,22b
I. Garzia,
22a,22b
E. Luppi,
22a,22b
L. Piemontese,
22a
V. Santoro,
22a
R. Baldini-Ferroli,
23
A. Calcaterra,
23
R. de Sangro,
23
G. Finocchiaro,
23
P. Patteri,
23
I. M. Peruzzi,
23,
†
M. Piccolo,
23
M. Rama,
23
A. Zallo,
23
R. Contri,
24a,24b
E. Guido,
24a,24b
M. Lo Vetere,
24a,24b
M. R. Monge,
24a,24b
S. Passaggio,
24a
C. Patrignani,
24a,24b
E. Robutti,
24a
B. Bhuyan,
25
V. Prasad,
25
M. Morii,
26
A. Adametz,
27
U. Uwer,
27
H. M. Lacker,
28
T. Lueck,
28
P. D. Dauncey,
29
U. Mallik,
30
C. Chen,
31
J. Cochran,
31
W. T. Meyer,
31
S. Prell,
31
A. E. Rubin,
31
A. V. Gritsan,
32
N. Arnaud,
33
M. Davier,
33
D. Derkach,
33
G. Grosdidier,
33
F. Le Diberder,
33
A. M. Lutz,
33
B. Malaescu,
33
P. Roudeau,
33
M. H. Schune,
33
A. Stocchi,
33
G. Wormser,
33
D. J. Lange,
34
D. M. Wright,
34
J. P. Coleman,
35
J. R. Fry,
35
E. Gabathuler,
35
D. E. Hutchcroft,
35
D. J. Payne,
35
C. Touramanis,
35
A. J. Bevan,
36
F. Di Lodovico,
36
R. Sacco,
36
M. Sigamani,
36
G. Cowan,
37
D. N. Brown,
38
C. L. Davis,
38
A. G. Denig,
39
M. Fritsch,
39
W. Gradl,
39
K. Griessinger,
39
A. Hafner,
39
E. Prencipe,
39
R. J. Barlow,
40,
‡
G. D. Lafferty,
40
E. Behn,
41
R. Cenci,
41
B. Hamilton,
41
A. Jawahery,
41
D. A. Roberts,
41
C. Dallapiccola,
42
R. Cowan,
43
D. Dujmic,
43
G. Sciolla,
43
R. Cheaib,
44
P. M. Patel,
44,
§
S. H. Robertson,
44
P. Biassoni,
45a,45b
N. Neri,
45a
F. Palombo,
45a,45b
L. Cremaldi,
46
R. Godang,
46,
k
R. Kroeger,
46
P. Sonnek,
46
D. J. Summers,
46
X. Nguyen,
47
M. Simard,
47
P. Taras,
47
G. De Nardo,
48a,48b
D. Monorchio,
48a,48b
G. Onorato,
48a,48b
C. Sciacca,
48a,48b
M. Martinelli,
49
G. Raven,
49
C. P. Jessop,
50
J. M. LoSecco,
50
K. Honscheid,
51
R. Kass,
51
J. Brau,
52
R. Frey,
52
N. B. Sinev,
52
D. Strom,
52
E. Torrence,
52
E. Feltresi,
53a,53b
N. Gagliardi,
53a,53b
M. Margoni,
53a,53b
M. Morandin,
53a
A. Pompili,
53a
M. Posocco,
53a
M. Rotondo,
53a
G. Simi,
53a
F. Simonetto,
53a,53b
R. Stroili,
53a,53b
S. Akar,
54
E. Ben-Haim,
54
M. Bomben,
54
G. R. Bonneaud,
54
H. Briand,
54
G. Calderini,
54
J. Chauveau,
54
O. Hamon,
54
Ph. Leruste,
54
G. Marchiori,
54
J. Ocariz,
54
S. Sitt,
54
M. Biasini,
55a,55b
E. Manoni,
55a,55b
S. Pacetti,
55a,55b
A. Rossi,
55a,55b
C. Angelini,
56a,56b
G. Batignani,
56a,56b
S. Bettarini,
56a,56b
M. Carpinelli,
56a,56b,
{
G. Casarosa,
56a,56b
A. Cervelli,
56a,56b
F. Forti,
56a,56b
M. A. Giorgi,
56a,56b
A. Lusiani,
56a,56c
B. Oberhof,
56a,56b
E. Paoloni,
56a,56b
A. Perez,
56a
G. Rizzo,
56a,56b
J. J. Walsh,
56a
D. Lopes Pegna,
57
J. Olsen,
57
A. J. S. Smith,
57
F. Anulli,
58a
R. Faccini,
58a,58b
F. Ferrarotto,
58a
F. Ferroni,
58a,58b
M. Gaspero,
58a,58b
L. Li Gioi,
58a
M. A. Mazzoni,
58a
G. Piredda,
58a
C. Bu
̈
nger,
59
O. Gru
̈
nberg,
59
T. Hartmann,
59
T. Leddig,
59
C. Voß,
59
R. Waldi,
59
T. Adye,
60
E. O. Olaiya,
60
F. F. Wilson,
60
S. Emery,
61
G. Hamel de Monchenault,
61
G. Vasseur,
61
Ch. Ye
`
che,
61
D. Aston,
62
D. J. Bard,
62
J. F. Benitez,
62
C. Cartaro,
62
M. R. Convery,
62
J. Dorfan,
62
G. P. Dubois-Felsmann,
62
W. Dunwoodie,
62
M. Ebert,
62
R. C. Field,
62
B. G. Fulsom,
62
A. M. Gabareen,
62
M. T. Graham,
62
C. Hast,
62
W. R. Innes,
62
M. H. Kelsey,
62
P. Kim,
62
M. L. Kocian,
62
D. W. G. S. Leith,
62
P. Lewis,
62
D. Lindemann,
62
B. Lindquist,
62
S. Luitz,
62
V. Luth,
62
H. L. Lynch,
62
D. B. MacFarlane,
62
D. R. Muller,
62
H. Neal,
62
S. Nelson,
62
M. Perl,
62
T. Pulliam,
62
B. N. Ratcliff,
62
A. Roodman,
62
A. A. Salnikov,
62
R. H. Schindler,
62
A. Snyder,
62
D. Su,
62
M. K. Sullivan,
62
J. Va’vra,
62
A. P. Wagner,
62
W. F. Wang,
62
W. J. Wisniewski,
62
M. Wittgen,
62
D. H. Wright,
62
H. W. Wulsin,
62
V. Ziegler,
62
W. Park,
63
M. V. Purohit,
63
R. M. White,
63
J. R. Wilson,
63
A. Randle-Conde,
64
S. J. Sekula,
64
M. Bellis,
65
P. R. Burchat,
65
T. S. Miyashita,
65
E. M. T. Puccio,
65
M. S. Alam,
66
J. A. Ernst,
66
R. Gorodeisky,
67
N. Guttman,
67
D. R. Peimer,
67
A. Soffer,
67
S. M. Spanier,
68
J. L. Ritchie,
69
A. M. Ruland,
69
R. F. Schwitters,
69
B. C. Wray,
69
J. M. Izen,
70
X. C. Lou,
70
F. Bianchi,
71a,71b
D. Gamba,
71a,71b
S. Zambito,
71a,71b
L. Lanceri,
72a,72b
L. Vitale,
72a,72b
F. Martinez-Vidal,
73
A. Oyanguren,
73
P. Villanueva-Perez,
73
H. Ahmed,
74
J. Albert,
74
Sw. Banerjee,
74
F. U. Bernlochner,
74
H. H. F. Choi,
74
G. J. King,
74
R. Kowalewski,
74
M. J. Lewczuk,
74
I. M. Nugent,
74
J. M. Roney,
74
R. J. Sobie,
74
N. Tasneem,
74
T. J. Gershon,
75
P. F. Harrison,
75
T. E. Latham,
75
H. R. Band,
76
S. Dasu,
76
Y. Pan,
76
R. Prepost,
76
and S. L. Wu
76
PHYSICAL REVIEW D
87,
052012 (2013)
1550-7998
=
2013
=
87(5)
=
052012(10)
052012-1
Ó
2013 American Physical Society
(
B
A
B
AR
Collaboration)
1
Laboratoire d’Annecy-le-Vieux de Physique des Particules (LAPP), Universite
́
de Savoie,
CNRS/IN2P3, F-74941 Annecy-Le-Vieux, France
2
Departament ECM, Facultat de Fisica, Universitat de Barcelona, E-08028 Barcelona, Spain
3a
INFN Sezione di Bari, I-70126 Bari, Italy
3b
Dipartimento di Fisica, Universita
`
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 Universita
̈
t Bochum, Institut fu
̈
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
9
Budker Institute of Nuclear Physics SB RAS, Novosibirsk 630090, Russia
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
Fakulta
̈
t Physik, Technische Universita
̈
t Dortmund, D-44221 Dortmund, Germany
19
Technische Universita
̈
t Dresden, Institut fu
̈
r Kern und Teilchenphysik, D-01062 Dresden, Germany
20
Laboratoire Leprince-Ringuet, Ecole Polytechnique, CNRS/IN2P3, F-91128 Palaiseau, France
21
University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom
22a
INFN Sezione di Ferrara, I-44100 Ferrara, Italy
22b
Dipartimento di Fisica, Universita
`
di Ferrara, I-44100 Ferrara, Italy
23
INFN Laboratori Nazionali di Frascati, I-00044 Frascati, Italy
24a
INFN Sezione di Genova, I-16146 Genova, Italy
24b
Dipartimento di Fisica, Universita
`
di Genova, I-16146 Genova, Italy
25
Indian Institute of Technology Guwahati, Guwahati, Assam 781 039, India
26
Harvard University, Cambridge, Massachusetts 02138, USA
27
Universita
̈
t Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany
28
Humboldt-Universita
̈
t zu Berlin, Institut fu
̈
r Physik, Newtonstrasse 15, D-12489 Berlin, Germany
29
Imperial College London, London SW7 2AZ, United Kingdom
30
University of Iowa, Iowa City, Iowa 52242, USA
31
Iowa State University, Ames, Iowa 50011-3160, USA
32
Johns Hopkins University, Baltimore, Maryland 21218, USA
33
Laboratoire de l’Acce
́
le
́
rateur Line
́
aire, IN2P3/CNRS et Universite
́
Paris-Sud 11,
Centre Scientifique d’Orsay, B. P. 34, F-91898 Orsay Cedex, France
34
Lawrence Livermore National Laboratory, Livermore, California 94550, USA
35
University of Liverpool, Liverpool L69 7ZE, United Kingdom
36
Queen Mary, University of London, London E1 4NS, United Kingdom
37
University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 0EX, United Kingdom
38
University of Louisville, Louisville, Kentucky 40292, USA
39
Johannes Gutenberg-Universita
̈
t Mainz, Institut fu
̈
r Kernphysik, D-55099 Mainz, Germany
40
University of Manchester, Manchester M13 9PL, United Kingdom
41
University of Maryland, College Park, Maryland 20742, USA
42
University of Massachusetts, Amherst, Massachusetts 01003, USA
43
Laboratory for Nuclear Science, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
44
McGill University, Montre
́
al, Que
́
bec H3A 2T8, Canada
45a
INFN Sezione di Milano, I-20133 Milano, Italy
45b
Dipartimento di Fisica, Universita
`
di Milano, I-20133 Milano, Italy
46
University of Mississippi, University, Mississippi 38677, USA
47
Universite
́
de Montre
́
al, Physique des Particules, Montre
́
al, Que
́
bec H3C 3J7, Canada
48a
INFN Sezione di Napoli, I-80126 Napoli, Italy
48b
Dipartimento di Scienze Fisiche, Universita
`
di Napoli Federico II, I-80126 Napoli, Italy
49
NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, The Netherlands
50
University of Notre Dame, Notre Dame, Indiana 46556, USA
J. P. LEES
et al.
PHYSICAL REVIEW D
87,
052012 (2013)
052012-2
51
Ohio State University, Columbus, Ohio 43210, USA
52
University of Oregon, Eugene, Oregon 97403, USA
53a
INFN Sezione di Padova, I-35131 Padova, Italy
53b
Dipartimento di Fisica, Universita
`
di Padova, I-35131 Padova, Italy
54
Laboratoire de Physique Nucle
́
aire et de Hautes Energies, IN2P3/CNRS, Universite
́
Pierre et Marie Curie-Paris6,
Universite
́
Denis Diderot-Paris7, F-75252 Paris, France
55a
INFN Sezione di Perugia, I-06100 Perugia, Italy
55b
Dipartimento di Fisica, Universita
`
di Perugia, I-06100 Perugia, Italy
56a
INFN Sezione di Pisa, I-56127 Pisa, Italy
56b
Dipartimento di Fisica, Universita
`
di Pisa, I-56127 Pisa, Italy
56c
Scuola Normale Superiore di Pisa, I-56127 Pisa, Italy
57
Princeton University, Princeton, New Jersey 08544, USA
58a
INFN Sezione di Roma, I-00185 Roma, Italy
58b
Dipartimento di Fisica, Universita
`
di Roma La Sapienza, I-00185 Roma, Italy
59
Universita
̈
t Rostock, D-18051 Rostock, Germany
60
Rutherford Appleton Laboratory, Chilton, Didcot, Oxon OX11 0QX, United Kingdom
61
Centre de Saclay, CEA, Irfu, SPP, F-91191 Gif-sur-Yvette, France
62
SLAC National Accelerator Laboratory, Stanford, California 94309, USA
63
University of South Carolina, Columbia, South Carolina 29208, USA
64
Southern Methodist University, Dallas, Texas 75275, USA
65
Stanford University, Stanford, California 94305-4060, USA
66
State University of New York, Albany, New York 12222, USA
67
School of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel
68
University of Tennessee, Knoxville, Tennessee 37996, USA
69
University of Texas at Austin, Austin, Texas 78712, USA
70
University of Texas at Dallas, Richardson, Texas 75083, USA
71a
INFN Sezione di Torino, I-10125 Torino, Italy
71b
Dipartimento di Fisica Sperimentale, Universita
`
di Torino, I-10125 Torino, Italy
72a
INFN Sezione di Trieste, I-34127 Trieste, Italy
72b
Dipartimento di Fisica, Universita
`
di Trieste, I-34127 Trieste, Italy
73
IFIC, Universitat de Valencia-CSIC, E-46071 Valencia, Spain
74
University of Victoria, Victoria, British Columbia V8W 3P6, Canada
75
Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
76
University of Wisconsin, Madison, Wisconsin 53706, USA
(Received 14 December 2012; published 18 March 2013)
We report a search for
CP
violation in the decay modes
D
!
K
0
S
K
,
D
s
!
K
0
S
K
, and
D
s
!
K
0
S
using a data set corresponding to an integrated luminosity of
469 fb
1
collected with the
BABAR
detector at the PEP-II asymmetric energy
e
þ
e
storage rings. The decay rate
CP
asymmetries,
A
CP
,
are determined to be
ðþ
0
:
13
0
:
36
ð
stat
Þ
0
:
25
ð
syst
ÞÞ
%
,
ð
0
:
05
0
:
23
ð
stat
Þ
0
:
24
ð
syst
ÞÞ
%
, and
ðþ
0
:
6
2
:
0
ð
stat
Þ
0
:
3
ð
syst
ÞÞ
%
, respectively. These measurements are consistent with zero, and also
with the Standard Model prediction [
ð
0
:
332
0
:
006
Þ
%
for the
D
!
K
0
S
K
and
D
s
!
K
0
S
K
modes,
and
ðþ
0
:
332
0
:
006
Þ
%
for the
D
s
!
K
0
S
mode]. They are the most precise determinations to date.
DOI:
10.1103/PhysRevD.87.052012
PACS numbers: 11.30.Er, 13.25.Ft, 14.40.Lb
I. INTRODUCTION
The search for
CP
violation (
CPV
) in charm decays
provides a sensitive probe of physics beyond the Standard
Model (SM). Owing to its suppression within the SM, a
significant observation of direct
CPV
in charm decays
would indicate the possible presence of new physics effects
in the decay processes. In a previous article [
1
], we re-
ported a precise measurement of the
CP
asymmetry in the
D
!
K
0
S
mode, where the measured asymmetry was
found to be consistent with the value expected from indi-
rect
CPV
in the
K
0
system.
The LHCb and CDF Collaborations have recently re-
ported evidence for
CPV
in charm decays by measuring
the difference of
CP
asymmetries in the
D
0
!
K
þ
K
and
D
0
!
þ
channels [
2
,
3
], which is mainly sensitive to
direct
CPV
. The size of the world average direct
CP
*
Present address: The University of Tabuk, Tabuk 71491,
Saudi Arabia.
†
Also at Dipartimento di Fisica, Universita
`
di Perugia,
Perugia, Italy.
‡
Present address: The University of Huddersfield,
Huddersfield HD1 3DH, United Kingdom.
§
Deceased.
k
Present address: University of South Alabama, Mobile, AL
36688, USA.
{
Also at Universita
`
di Sassari, Sassari, Italy.
SEARCH FOR
CP
VIOLATION IN THE DECAYS
...
PHYSICAL REVIEW D
87,
052012 (2013)
052012-3
asymmetry difference,
ð
6
:
56
1
:
54
Þ
10
3
[
4
], sug-
gests either a significant enhancement of SM penguin
amplitudes or of new physics amplitudes (or both) in
charm decays [
5
]. Improved measurements of the
CP
asymmetries in the individual two-body modes, along
with measurements in other channels, are needed to deter-
mine the nature of the contributing amplitudes.
We present herein measurements of the decay rate
CP
asymmetry,
A
CP
, defined as
A
CP
¼
ð
D
þ
ð
s
Þ
!
f
Þ
ð
D
ð
s
Þ
!
f
Þ
ð
D
þ
ð
s
Þ
!
f
Þþ
ð
D
ð
s
Þ
!
f
Þ
;
(1)
in the decay modes
D
!
K
0
S
K
,
D
s
!
K
0
S
K
, and
D
s
!
K
0
S
. Previous measurements of
A
CP
in these
channels have been reported by the CLEO-c [
6
] and
Belle Collaborations [
7
]. As for the
A
CP
measurement
in
D
!
K
0
S
, we expect an
A
CP
asymmetry of
ð
0
:
332
0
:
006
Þ
%
[
8
] resulting from
CPV
in
K
0
K
0
mixing [
9
]. The sign of the
K
0
-induced asymmetry is
positive (negative) if a
K
0
(
K
0
) is present in the corre-
sponding tree-level Feynman diagram. Because it is iden-
tified by its
þ
decay, the intermediate state is a
coherent mix of
K
0
S
and
K
0
L
amplitudes. It has been shown
in Ref. [
10
] that the
K
0
S
K
0
L
interference term gives rise to
a measured
CP
asymmetry that depends on the range in
proper time over which the decay rates are integrated and
on the efficiency for the reconstruction of the intermediate
state as a function of its proper flight time.
For this analysis we employ a technique similar to that
used for our measurement of
CPV
in the
D
!
K
0
S
mode [
1
]. As a result, reference to our previous publication
is given for the description of some of the analysis details.
II. THE
BABAR
DETECTOR
AND EVENT SELECTION
The data used for these measurements were recorded at
or near the
ð
4
S
Þ
resonance by the
BABAR
detector at the
PEP-IIstorage rings and correspond to an integrated lumi-
nosity of
469 fb
1
. Charged particles are detected, and
their momenta measured, by a combination of a silicon
vertex tracker, consisting of five layers of double-sided
detectors, and a 40-layer central drift chamber, both oper-
ating in a 1.5 T axial magnetic field. Charged-particle
identification is provided by specific ionization energy
loss measurements in the tracking system and by the
measured Cherenkov angle from an internally reflecting
ring-imaging Cherenkov detector covering the central re-
gion of the detector. Electrons are detected by a CsI(Tl)
electromagnetic calorimeter. The
BABAR
detector, and the
coordinate system used throughout, are described in detail
in Refs. [
11
,
12
]. We validate the analysis procedure using
Monte Carlo (MC) simulation based on
GEANT4
[
13
]. The
MC samples include
e
þ
e
!
q
q
ð
q
¼
u;d;s;c
Þ
events,
simulated with JETSET [
14
] and
B
B
decays simulated
with the
EVTGEN
generator [
15
]. To avoid potential bias
in the measurements we finalize the event selection for
each channel, as well as the procedures for efficiency
correction, fitting, and the determination of the systematic
uncertainties and possible biases in the measurements,
prior to extracting the value of
A
CP
from the data.
Signal candidates are reconstructed by combining a
K
0
S
candidate, reconstructed in the decay mode
K
0
S
!
þ
,
with a charged pion or kaon candidate. A
K
0
S
candidate is
reconstructed from two oppositely charged tracks with an
invariant mass within a
10 MeV
=c
2
interval centered on
the nominal
K
0
S
mass [
8
], which is approximately
2
:
5
in
the measured
K
0
S
mass resolution. The
2
probability of the
þ
vertex fit must be greater than 0.1%. Motivated by
MC studies, we require the measured flight length of the
K
0
S
candidate to be at least 3 times greater than its uncer-
tainty, to reduce combinatorial background. A recon-
structed charged-particle track that has
p
T
400 MeV
=c
is selected as a pion or kaon candidate, where
p
T
is the
magnitude of the momentum in the plane perpendicular to
the
z
axis (transverse plane). In our measurement, we
require that a pion candidate not be identified as a kaon,
a proton, or an electron, and that a kaon candidate be
identified as a kaon, and not as a pion, a proton, or an
electron. Identification efficiencies and misidentification
rates for electron, pions, kaons, and protons with
2 GeV
=c
momentum in the laboratory frame are reported
in Table
I
. The criteria used to select pion or kaon candi-
dates are very effective in reducing the charge asymmetry
from track reconstruction and identification, as inferred
from studying the data control samples described below.
A vertex fit to the whole decay chain, constraining the
D
ð
s
Þ
production vertex to be within the
e
þ
e
interaction region,
is then performed [
16
]. We retain only
D
ð
s
Þ
candidates
having a
2
probability for this fit greater than 0.1%, and
an invariant mass
m
ð
K
0
S
h
Þ
,
h
¼
,
K
, within a
65 MeV
=c
2
interval centered on the nominal
D
ð
s
Þ
mass
[
8
], which is approximately equivalent to
8
in the
measured
D
ð
s
Þ
mass resolution.
TABLE I. Identification efficiencies and misidentification
rates for electron, pions, kaons, and protons with
2 GeV
=c
momentum in the laboratory frame. The values for kaons on
the third row refers to the identification criterion used to reject
kaons from the pion sample, while the values on the fourth row
to the criterion used in the kaon selection.
Misidentification rate [%]
Particle
Efficiency [%]
K
e
91
0.04
<
0
:
2
88
not applicable
1
K
(applied to
)
91
1
not applicable
K
(applied to
K
)
99
8
not applicable
p
80
0.2
0.2
J. P. LEES
et al.
PHYSICAL REVIEW D
87,
052012 (2013)
052012-4
We require further that the magnitude of the
D
s
candi-
date momentum in the
e
þ
e
center-of-mass system,
p
,be
between 2.6 and
5
:
0 GeV
=c
, in order to suppress combi-
natorial background from
B
B
events. For the
D
!
K
0
S
K
mode, the MC simulated sample shows that retaining
candidates with
p
between 2.0 and
5
:
0 GeV
=c
allows
signal candidates from
B
-meson decays, without introduc-
ing an excessive amount of combinatorial background.
Assuming that
CPT
is conserved, there is no contribution
to
A
CP
from
CP
violation in
B
meson decays from
Standard Model processes. Additional background rejec-
tion is obtained by requiring that the impact parameter of
the
D
ð
s
Þ
candidate with respect to the beam spot [
11
],
projected onto the transverse plane, be less than 0.3 cm,
and that the
D
ð
s
Þ
proper decay time,
t
xy
, be between
15
and 35 ps. The decay time is measured using
L
xy
, defined
as the distance of the
D
ð
s
Þ
decay vertex from the beam spot
projected onto the transverse plane.
In order to further optimize the sensitivity of the
A
CP
measurements, we construct a multivariate algorithm,
based on seven discriminating variables for each
D
ð
s
Þ
candidate:
t
xy
,
L
xy
,
p
, the momentum magnitude and
component in the transverse plane for the
K
0
S
candidate,
and also for the pion or kaon candidate. For the
D
!
K
0
S
K
and
D
s
!
K
0
S
K
modes the multivariate algorithm
with the best performance is a boosted decision tree [
17
],
while for the
D
s
!
K
0
S
mode the best algorithm is a
projective likelihood method [
17
]. The final selection cri-
teria, based on the outputs of the multivariate selectors, are
optimized using truth-matched signal and background can-
didates from the MC sample. For the optimization, we
maximize the
S=
ffiffiffiffiffiffiffiffiffiffiffiffiffi
S
þ
B
p
ratio, where
S
and
B
are the
numbers of signal and background candidates, respec-
tively, with invariant mass within
30 MeV
=c
2
of the
nominal
D
ð
s
Þ
mass, which is approximately
3
in the
measured mass resolution.
III. SIGNAL YIELD AND ASYMMETRY
EXTRACTION
For each mode the signal yield is extracted using a
binned maximum likelihood (ML) fit to the distribution
of the invariant mass
m
ð
K
0
S
h
Þ
for the selected
D
ð
s
Þ
candi-
dates. The total probability density function (PDF) is the
sum of signal and background components. The signal
PDF is modeled as a sum of two Gaussian functions for
the
D
ð
s
Þ
!
K
0
S
K
modes, and as a single Gaussian function
for the
D
s
!
K
0
S
mode. The background PDF is taken
as the sum of two components: a distribution describing the
invariant mass of misreconstructed charm meson decays,
and a combinatorial background modeling the mass
distribution from other sources. For the
D
!
K
0
S
K
(
D
s
!
K
0
S
) mode the charm background is mainly
from the tail of the invariant mass distribution for
D
s
!
K
0
S
K
(
D
!
K
0
S
) candidates. For the
D
s
!
K
0
S
K
mode, the misreconstructed charm background
originates mainly from
D
!
K
0
S
decays for which
the
is misidentified as a
K
. Assigning the wrong
mass to the pion shifts the reconstructed invariant mass,
and the resulting distribution is a broad peak with
mean value close to the
D
s
mass. For each mode, the
invariant mass distribution due to charm background is
modeled using a histogram PDF obtained from a MC
sample of simulated charm background decays. The com-
binatorial background is described by a first(second)-order
polynomial for the
D
s
!
K
0
S
mode (
D
!
K
0
S
K
and
D
s
!
K
0
S
K
modes). The fits to the
m
ð
K
0
S
h
Þ
distributions yield
ð
159
:
4
0
:
8
Þ
10
3
D
!
K
0
S
K
de-
cays,
ð
288
:
2
1
:
1
Þ
10
3
D
s
!
K
0
S
K
decays, and
ð
14
:
33
0
:
31
Þ
10
3
D
s
!
K
0
S
decays. The data and
the fit results are shown in Fig.
1
. All of the PDF parame-
ters are extracted from fits to the data.
For each channel, we determine
A
CP
by measuring the
signal yield asymmetry
A
defined as
]
2
c
)[GeV/
±
K
S
0
K
(
m
1.82
1.84
1.86
1.88
1.9
1.92
)
2
c
Event / (1 MeV/
0
2000
4000
6000
8000
10000
12000
14000
(a)
]
2
c
)[GeV/
±
K
S
0
K
(
m
1.92
1.94
1.96
1.98
2
2.02
)
2
c
Event / (1 MeV/
0
2000
4000
6000
8000
10000
12000
14000
16000
18000
20000
22000
(b)
]
2
c
)[GeV/
±
π
S
0
K
(
m
1.92
1.94
1.96
1.98
2
2.02
)
2
c
Event / (1 MeV/
0
500
1000
1500
2000
2500
(c)
FIG. 1 (color online). Invariant mass distribution for (a)
D
!
K
0
S
K
, (b)
D
s
!
K
0
S
K
, and (c)
D
s
!
K
0
S
candidates (points
with error bars). The solid curve shows the result of the fit to the data. The dashed curve represents the sum of all background
contributions, while the dotted curve indicates combinatorial background only.
SEARCH FOR
CP
VIOLATION IN THE DECAYS
...
PHYSICAL REVIEW D
87,
052012 (2013)
052012-5
A
¼
N
D
þ
ð
s
Þ
N
D
ð
s
Þ
N
D
þ
ð
s
Þ
þ
N
D
ð
s
Þ
;
(2)
where
N
D
þ
ð
s
Þ
(
N
D
ð
s
Þ
) is the number of
D
þ
ð
s
Þ
(
D
ð
s
Þ
) decays
determined from the fit to the invariant mass distribution.
The asymmetry
A
contains two contributions in addition to
A
CP
, namely the forward-backward (FB) asymmetry (
A
FB
),
and a detector-induced component. We measure
A
FB
together with
A
CP
using the selected data set, while we
correct the data for the detector-induced component using
coefficients derived from a control sample.
IV. CORRECTION OF DETECTOR-RELATED
ASYMMETRIES
We use a data-driven method, described in detail in
Ref. [
1
], to determine the charge asymmetry in track
reconstruction as a function of the magnitude of the track
momentum and its polar angle in the laboratory frame. The
method exploits the fact that
ð
4
S
Þ!
B
B
events provide a
sample evenly populated with positive and negative tracks,
free of any physics-induced asymmetries. The off-
resonance momentum distribution is subtracted from the
on-resonance one, to remove any contribution from
continuum, for which there is a FB asymmetry in the
center-of-mass frame. This sample is used to compute
the detector-related asymmetries in the reconstruction of
charged-particle tracks. Starting from a sample of
50
:
6fb
1
of data collected at the
ð
4
S
Þ
resonance and
an off-resonance data sample of
44
:
8fb
1
, we obtain a
large sample of charged-particle tracks and apply the same
charged pion or kaon track selection criteria used in the
reconstruction of the
D
ð
s
Þ
!
K
0
S
K
and
D
s
!
K
0
S
modes. Then, after subtracting the off-resonance contribu-
tion from the on-resonance sample, we obtain a sample of
more than
120
10
6
pion candidates, and
40
10
6
kaon
candidates, originating from
ð
4
S
Þ
decays. We use the full
off-resonance sample and an equivalent luminosity for the
on-resonance sample, because, due to the subtraction pro-
cedure, including additional data in the on-resonance sam-
ple does not improve the statistical error on the correction
ratios mentioned below. These candidates are then used to
compute the efficiency ratios for positive and negative
pions and kaons. The ratio values and their statistical errors
for pions and kaons are shown in Figs.
2
and
3
, respec-
tively. For the
D
ð
s
Þ
!
K
0
S
K
(
D
s
!
K
0
S
) modes, the
D
ð
s
Þ
(
D
s
) yields, in intervals of kaon (pion) momentum
and cosine of its polar angle,
cos
, are weighted with the
kaon (pion) efficiency ratios to correct for the detection
efficiency differences between
K
þ
and
K
(
þ
and
).
Momentum and cosine of its polar angle intervals are not
uniform in order to have similar statistics, and therefore
similar correction uncertainty, in each interval. Interval
sizes vary from (
0
:
05 GeV
=c
, 0.06) to (
4
:
4 GeV
=c
, 0.96),
where the first number is the momentum interval, and the
second its cosine of polar angle interval. The largest
correction is approximately 1% for pions and 2% for
kaons. After correcting the data for the detector-induced
component only
A
FB
and
A
CP
contribute to the measured
asymmetry
A
.
0.97
0.98
0.99
1
1.01
1.02
1.03
)
c
Momentum (GeV/
0123456
θ
cos
-1
-0.8
-0.6
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
0
0.002
0.004
0.006
0.008
0.01
0.012
0.014
0.016
)
c
Momentum (GeV/
0123456
θ
cos
-1
-0.8
-0.6
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
FIG. 3 (color online). (Top) The ratio between the detection
efficiency for
K
þ
and
K
, and (bottom) the corresponding
statistical errors. The values are computed using the numbers
of
K
þ
and
K
tracks in the selected control sample.
0.985
0.99
0.995
1
1.005
1.01
)
c
Momentum (GeV/
0123456
θ
cos
-1
-0.8
-0.6
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
0
0.001
0.002
0.003
0.004
0.005
0.006
0.007
)
c
Momentum (GeV/
0123456
θ
cos
-1
-0.8
-0.6
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
FIG. 2 (color online). (Top) The ratio between the detection
efficiency for
þ
and
, and (bottom) the corresponding
statistical errors. The values are computed using the numbers
of
þ
and
tracks in the selected control sample.
J. P. LEES
et al.
PHYSICAL REVIEW D
87,
052012 (2013)
052012-6
V. EXTRACTION OF
A
CP
AND
A
FB
Neglecting higher-order terms that contain
A
CP
and
A
FB
,
the resulting asymmetry can be expressed simply as the
sum of the two. Given that
A
FB
is an odd function of
cos
D
, where
D
is the polar angle of the
D
ð
s
Þ
candidate
momentum in the center-of-mass frame,
A
CP
and
A
FB
can
be written as a function of
j
cos
D
j
as follows:
A
CP
ðj
cos
D
jÞ¼
A
ðþj
cos
D
jÞþ
A
ðj
cos
D
jÞ
2
(3)
and
A
FB
ðj
cos
D
jÞ¼
A
ðþj
cos
D
jÞ
A
ðj
cos
D
jÞ
2
;
(4)
where
A
ðþj
cos
D
jÞ½
A
ðj
cos
D
jÞ
is the measured
asymmetry for the
D
ð
s
Þ
candidates in a positive (negative)
cos
D
interval.
A simultaneous ML fit to the
D
þ
ð
s
Þ
and
D
ð
s
Þ
invariant
mass distributions is carried out to extract the signal yield
asymmetry in each of ten equally spaced
cos
D
intervals,
starting with interval 1 at
½
1
:
0
;
0
:
8
. The PDF model
that describes the distribution in each subsample is the
same as that used in the fit to the full sample, but the
following parameters are allowed to float separately in
each subsample (referred to as split parameters): the yields
for signal, charm background and combinatorial candi-
dates; the asymmetries for signal and combinatorial can-
didates; the width, and the fraction of the Gaussian
function with the larger contribution to the signal PDF;
and the first-order coefficient of the polynomial that mod-
els the combinatorial background. For the
D
!
K
0
S
K
mode the yields for the charm background candidates in
intervals 1, 2, and 3 were fixed to 0 to obtain a fully
convergent fit. Since interval 10 contains the smallest
number of candidates, we use a single Gaussian function
to model the signal PDF for the
D
ð
s
Þ
!
K
0
S
K
modes. For
the
CP
asymmetry of charm background candidates we use
the same floating parameters as for the signal candidates,
because the largest source of
CP
asymmetry for both
samples is due to
CPV
in
K
0
K
0
mixing. For the
D
s
!
K
0
S
mode, where the primary charm background
channel,
D
!
K
0
S
, has the same magnitude but
opposite-sign asymmetry due to
K
0
K
0
mixing, we use
a separate parameter for the asymmetry of the charm
background candidates. To achieve a more stable fit, if
the fit results for a split parameter are statistically compat-
ible between two or more subsamples, the parameter is
forced to have the same floating value among those sub-
samples only. For the
D
s
!
K
0
S
mode the width of the
|
D
*
θ
|cos
0
0.2
0.4
0.6
0.8
1
(%)
CP
A
-2
-1
0
1
2
3
(a)
|
D
*
θ
|cos
0
0.2
0.4
0.6
0.8
1
(%)
CP
A
-1.5
-1
-0.5
0
0.5
1
(b)
|
D
*
θ
|cos
0
0.2
0.4
0.6
0.8
1
(%)
CP
A
-20
-15
-10
-5
0
5
10
(c)
|
D
*
θ
|cos
0
0.2
0.4
0.6
0.8
1
(%)
FB
A
-6
-5
-4
-3
-2
-1
0
(d)
|
D
*
θ
|cos
0
0.2
0.4
0.6
0.8
1
(%)
FB
A
-6
-5
-4
-3
-2
-1
0
(e)
|
D
*
θ
|cos
0
0.2
0.4
0.6
0.8
1
(%)
FB
A
-25
-20
-15
-10
-5
0
5
(f)
FIG. 4.
CP
asymmetry,
A
CP
, for (a)
D
!
K
0
S
K
, (b)
D
s
!
K
0
S
K
, and (c)
D
s
!
K
0
S
as a function of
j
cos
D
j
in the data
sample. The solid line represents the central value of
A
CP
and the gray band is the
1
interval, both obtained from a
2
minimization
assuming no dependence on
j
cos
D
j
. The corresponding forward-backward asymmetries,
A
FB
, are shown in (d)–(f).
SEARCH FOR
CP
VIOLATION IN THE DECAYS
...
PHYSICAL REVIEW D
87,
052012 (2013)
052012-7
first Gaussian function for the signal PDF is set to the same
floating value in intervals 1, 2, 3, and 4. The first-order
coefficient of the polynomial describing the combinatorial
background is set to the same floating value in intervals
4–8 (
D
!
K
0
S
K
), in intervals 4–8 (
D
s
!
K
0
S
K
),
and in intervals 2–7 (
D
s
!
K
0
S
). The final fit contains
70, 80, and 64 free parameters for the
D
!
K
0
S
K
,
D
s
!
K
0
S
K
, and
D
s
!
K
0
S
modes, respectively.
The
A
CP
and
A
FB
values for the five
j
cos
D
j
bins are
shown in Fig.
4
for the three decay modes. The weighted
average of the five
A
CP
values is
ð
0
:
16
0
:
36
Þ
%
for the
D
!
K
0
S
K
mode,
ð
0
:
00
0
:
23
Þ
%
for
D
s
!
K
0
S
K
,
and
ð
0
:
6
2
:
0
Þ
%
for
D
s
!
K
0
S
, where the errors are
statistical only.
We perform two tests to validate the analysis procedure
for each channel. The first involves generating 5000 toy
MC experiments with a statistics equal to data using the
PDF and the parameters obtained from the fit to data.
After extracting
A
CP
from each experiment, for the
D
!
K
0
S
K
and
D
s
!
K
0
S
K
modes, we deduce from the
mean of the
A
CP
pull distributions the presence of a small
bias in the fitted value of each fit parameter (the means
are
0
:
036
0
:
014
and
þ
0
:
041
0
:
014
, respectively).
To account for this effect we apply a correction to the
final values equal to
þ
0
:
013%
for the
D
!
K
0
S
K
mode, and
0
:
01%
for the
D
s
!
K
0
S
K
mode. The
A
CP
pull distributions show that the fit provides an accu-
rate estimate of the statistical error for all the modes. The
second test involves fitting a large number of MC events
from the full
BABAR
detector simulation. We measure
A
CP
from this MC sample to be consistent with the
generated value of zero.
VI. SYSTEMATICS
The main sources of systematic uncertainty are listed in
Table
II
for each decay mode, together with the overall
uncertainties. The primary sources of systematic uncer-
tainty are the detection efficiency ratios used to weight
the
D
ð
s
Þ
yields, and the contributions from misidentified
particles in the data control sample used to determine the
charge asymmetry in track reconstruction efficiency.
The technique used to remove the charge asymmetry due
to detector-induced effects produces a small systematic
uncertainty in the measurement of
A
CP
due to the statistical
error in the relative efficiency estimation. This systematic
uncertainty depends only on the type of charged particle
(pion or kaon) in the final state, and not on the initial state.
To estimate the systematic uncertainty on
A
CP
resulting
from this source, the relative charged-particle efficiency in
each interval of momentum and
cos
is randomly drawn
from a Gaussian distribution whose mean is the nominal
relative efficiency in that interval, and where the root-
mean-squared (rms) deviation is the corresponding
statistical error. For each mode, we generate 500 such
charged-particle relative-efficiency distributions and use
them to obtain 500
A
CP
values, following the procedure
described earlier to determine the nominal value of
A
CP
.
The rms deviation of these 500 values from the nominal
A
CP
is taken to be the systematic uncertainty. For the
D
ð
s
Þ
!
K
0
S
K
modes, the estimated systematic uncertainty
is 0.23%. For the
D
s
!
K
0
S
mode, we assign the same
systematic uncertainty, 0.06%, as that estimated for the
D
!
K
0
S
mode in Ref. [
1
].
The small fraction of misidentified particles in the ge-
neric track sample can introduce small biases in the esti-
mation of the efficiencies, and subsequently in the
A
CP
measurements. Because of the good agreement between
data and MC samples, we can use the simulated MC
candidates to measure the shift in the
A
CP
value from the
fit when the corrections are applied, and when they are not.
Again, this contribution depends only on the type of the
charged-particle track. Hence, for the
D
s
!
K
0
S
mode,
we assume the same shift obtained in Ref. [
1
], namely
þ
0
:
05%
. By fitting the
D
s
!
K
0
S
K
MC sample when
the corrections are applied, and again when not, we
obtain a shift of
þ
0
:
05%
and we assume this for both the
D
s
!
K
0
S
K
and
D
!
K
0
S
K
modes. For all the modes,
we shift the measured
A
CP
by this correction value and
then, conservatively, include the magnitude of this shift as
a contribution to the systematic uncertainty.
Using MC simulation, we evaluate an additional system-
atic uncertainty of
0
:
01%
due to a possible charge asym-
metry present in the control sample before applying the
TABLE II. Summary of the systematic uncertainty contributions for the
A
CP
measurement in each mode. The values are absolute
uncertainties, even though given as percentages. The total value corresponds to the sum in quadrature of the individual contributions.
Systematic uncertainty
D
!
K
0
S
K
[%]
D
s
!
K
0
S
K
[%]
D
s
!
K
0
S
[%]
Efficiency of particle-identification selectors
0.05
0.05
0.05
Statistics of the control sample
0.23
0.23
0.06
Misidentified tracks in the control sample
0.01
0.01
0.01
cos
D
interval size
0.04
0.02
0.27
K
0
K
0
regeneration
0.05
0.05
0.06
K
0
S
K
0
L
interference
0.015
0.014
0.008
Total
0.25
0.24
0.29
J. P. LEES
et al.
PHYSICAL REVIEW D
87,
052012 (2013)
052012-8
selection criteria. Another source of systematic uncertainty
is due to the choice of the
cos
D
interval size in the
simultaneous ML fit. The systematic uncertainty is
taken to be the largest absolute difference between the
nominal
A
CP
extracted using ten
cos
D
intervals and that
obtained when the fit is performed using either 8 or 12
intervals in
cos
D
. This is the dominant source of system-
atic uncertainty for the
D
s
!
K
0
S
mode, as shown in
Table
II
.
We also consider a possible systematic uncertainty due
to the regeneration of neutral kaons in the material of the
detector. The
K
0
and
K
0
mesons produced in the decay
processes can interact with the material in the tracking
volume before they decay. Following a method similar to
that described in Ref. [
18
], we compute the probability for
a
K
0
or a
K
0
meson to interact inside the
BABAR
tracking
system and estimate systematic uncertainties of 0.05%
(
D
ð
s
Þ
!
K
0
S
K
) and 0.06% (
D
s
!
K
0
S
).
Although the intermediate state is labeled as a
K
0
S
,we
apply a correction term to the measured
A
CP
to include the
effect of
K
0
S
K
0
L
interference in the intermediate state
[
19
]. This correction term depends on the proper time
range over which decay distributions are integrated, and
on the efficiency of the reconstruction of the
þ
final
state as a function of proper time. We compute the recon-
struction efficiency distribution as a function of proper
time using MC truth-matched
K
0
S
decays after the full
selection. Following the method in Ref. [
19
] we estimate
the asymmetry-correction term
A
CP
defined as
A
CP
¼
A
corr
CP
A
fit
CP
;
(5)
where
A
fit
CP
is the value obtained from the fit and
A
corr
CP
is the
corrected value. The correction terms are reported in
Table
III
and, to be conservative, we include their absolute
values as contributions to the systematic uncertainty
estimates. We also estimate the correction factor for the
D
!
K
0
S
mode using the
K
0
S
reconstruction efficiency
distribution after the selection detailed in Ref. [
1
] and
obtain the value
þ
0
:
002%
. All these corrections are rather
small, even compared to those estimated in a similar
analysis [
20
]. The smaller values of the corrections in the
present analysis are due to the improved efficiency for
K
0
S
mesons with short decay times that we obtain by applying
the requirement on the decay length divided by its uncer-
tainty, rather than on the decay length alone.
VII. CONCLUSION
In conclusion, we measure the direct
CP
asymmetry
A
CP
in the
D
!
K
0
S
K
,
D
s
!
K
0
S
K
, and
D
s
!
K
0
S
modes using approximately 159000, 288000, and 14000
signal candidates, respectively. The measured
A
CP
value
for each mode is reported in Table
III
, where the first
errors are statistical and the second are systematic. In the
last row of the table, we also report the
A
CP
values after
subtracting the expected
A
CP
contribution for each mode
due to
K
0
K
0
mixing. The results are consistent with
zero, and with the SM prediction, within 1 standard
deviation.
ACKNOWLEDGMENTS
We are grateful for the extraordinary contributions of
our PEP-II colleagues in achieving the excellent luminos-
ity and machine conditions that have made this work
possible. The success of this project also relies critically
on the expertise and dedication of the computing
organizations that support
BABAR
. The collaborating
institutions wish to thank SLAC for its support and
the kind hospitality extended to them. This work is
supported by the U.S. Department of Energy and
National Science Foundation, the Natural Sciences
and Engineering Research Council (Canada), the
Commissariat a
`
l’Energie Atomique and Institut National
de Physique Nucle
́
aire et de Physique des Particules
(France), the Bundesministerium fu
̈
r Bildung und
Forschung and Deutsche Forschungsgemeinschaft
(Germany), the Istituto Nazionale di Fisica Nucleare
(Italy), the Foundation for Fundamental Research on
Matter (Netherlands), the Research Council of Norway,
the Ministry of Education and Science of the Russian
Federation, Ministerio de Ciencia e Innovacio
́
n (Spain),
and the Science and Technology Facilities Council (United
Kingdom). Individuals have received support from the
Marie-Curie IEF program (European Union) and the A. P.
Sloan Foundation (USA).
TABLE III. Summary of the
A
CP
measurements. Where reported, the first uncertainty is statistical, and the second is systematic.
D
!
K
0
S
K
D
s
!
K
0
S
K
D
s
!
K
0
S
A
CP
value from the fit
ðþ
0
:
155
0
:
360
Þ
%
ð
0
:
00
0
:
23
Þ
%
ðþ
0
:
6
2
:
0
Þ
%
Correction for the bias from toy MC experiments
þ
0
:
013%
0
:
01%
not applied
Correction for the bias in the particle-identification
selectors
0
:
05%
0
:
05%
0
:
05%
Correction for the
K
0
S
K
0
L
interference (
A
CP
)
þ
0
:
015%
þ
0
:
014%
0
:
008%
A
CP
final value
ðþ
0
:
13
0
:
36
0
:
25
Þ
%
ð
0
:
05
0
:
23
0
:
24
Þ
%
ðþ
0
:
6
2
:
0
0
:
3
Þ
%
A
CP
contribution from
K
0
K
0
mixing
ð
0
:
332
0
:
006
Þ
%
ð
0
:
332
0
:
006
Þ
%
ðþ
0
:
332
0
:
006
Þ
%
A
CP
final value (charm only)
ðþ
0
:
46
0
:
36
0
:
25
Þ
%
ðþ
0
:
28
0
:
23
0
:
24
Þ
%
ðþ
0
:
3
2
:
0
0
:
3
Þ
%
SEARCH FOR
CP
VIOLATION IN THE DECAYS
...
PHYSICAL REVIEW D
87,
052012 (2013)
052012-9
[1] P. del Amo Sanchez
et al.
(
BABAR
Collaboration),
Phys.
Rev. D
83
, 071103(R) (2011)
.
[2] R. Aaij
et al.
(LHCb Collaboration),
Phys. Rev. Lett.
108
,
111602 (2012)
.
[3] T. Aaltonen
et al.
(CDF Collaboration),
Phys. Rev. Lett.
109
, 111801 (2012)
.
[4] Y. Amhis
et al.
(Heavy Flavor Averaging Group),
arXiv:1207.1158
; for a more recent update see
http://
www.slac.stanford.edu/xorg/hfag/
.
[5] G. Isidori, J. F. Kamenik, Z. Ligeti, and G. Perez,
Phys.
Lett. B
711
, 46 (2012)
; H.-Y. Cheng and C.-W. Chiang,
Phys. Rev. D
85
, 034036 (2012)
;
85
, 079903(E) (2012)
;B.
Bhattacharya, M. Gronau, and J. L. Rosner,
Phys. Rev. D
85
, 054014 (2012)
; G. F. Giudice, G. Isidori, and P.
Paradisi,
J. High Energy Phys. 04 (2012) 060
;E.
Franco, S. Mishima, and L. Silvestrini,
J. High Energy
Phys. 05 (2012) 140
; J. Brod, Y. Grossman, A. L. Kagan,
and J. Zupan,
J. High Energy Phys. 10 (2012) 161
.
[6] S. Dobbs
et al.
(CLEO Collaboration),
Phys. Rev. D
76
,
112001 (2007)
.
[7] B. R. Ko
et al.
(Belle Collaboration),
Phys. Rev. Lett.
104
,
181602 (2010)
.
[8] K. Nakamura (Particle Data Group),
J. Phys. G
37
, 075021
(2010)
.
[9] H. J. Lipkin and Z.-z. Xing,
Phys. Lett. B
450
, 405
(1999)
.
[10] Y. I. Azimov and A. A. Iogansen, Sov. J. Nucl. Phys.
33
,
205 (1981).
[11] B. Aubert
et al.
(
BABAR
Collaboration),
Nucl. Instrum.
Methods Phys. Res., Sect. A
479
, 1 (2002)
.
[12] W. Menges,
IEEE Nucl. Sci. Symp. Conf. Rec.
3
, 1470
(2006)
.
[13] S. Agostinelli
et al.
(Geant4 Collaboration),
Nucl.
Instrum. Methods Phys. Res., Sect. A
506
, 250 (2003)
.
[14] T. Sjostrand, S. Mrenna, and P. Z. Skands,
J. High Energy
Phys. 05 (2006) 026
.
[15] D. J. Lange,
Nucl. Instrum. Methods Phys. Res., Sect. A
462
, 152 (2001)
.
[16] W. D. Hulsbergen,
Nucl. Instrum. Methods Phys. Res.,
Sect. A
552
, 566 (2005)
.
[17] P. Speckmayer, A. Ho
̈
cker, J. Stelzer, and H. Voss,
J. Phys.
Conf. Ser.
219
, 032057 (2010)
.
[18] B. R. Ko, E. Won, B. Golob, and P. Pakhlov,
Phys. Rev. D
84
, 111501(R) (2011)
.
[19] Y. Grossman and Y. Nir,
J. High Energy Phys. 04 (2012)
002
.
[20] J. P. Lees
et al.
(
BABAR
Collaboration),
Phys. Rev. D
85
,
031102(R) (2012)
;
85
, 099904(E) (2012)
.
J. P. LEES
et al.
PHYSICAL REVIEW D
87,
052012 (2013)
052012-10