Study of the CP asymmetry of B0 ->J/psi K0S decays in ALEPH

The decay B0 ->J/psi K0_S is reconstructed with J/psi ->e+ e- or mu+ mu- and K0_S ->pi+ pi-. From the full ALEPH dataset at LEP1 of about 4 million hadronic Z decays, 23 candidates are selected with an estimated purity of 71%. They are used to measure the CP asymmetry of this decay, given by sin 2beta in the Standard Model, with the result sin 2beta = 0.84 +0.82-1.04 +-0.16. This is combined with existing measurements from other experiments, and increases the confidence level that CP violation has been observed in this channel to 98%.


Introduction
In the Standard Model CP violation arises from a complex phase of the quark mixing matrix [1], and this can accommodate the observed CP violation in the K sector [2].Precise predictions can be made of relations between asymmetries expected in B decays, and their detailed study will provide an important test of the model.The first step, however, is to establish the existence of CP violation in B decays.This has been attempted using inclusive methods, where the expected asymmetry is small, O(10 −3 ), beyond the sensitivity of current experiments [3].The alternative is to use exclusive decays, where the asymmetries are predicted to be large, O (1), but the branching ratios are small.
The decay B 0 → J/ψ K 0 S is known as the gold-plated mode for such studies, due to its clean experimental signature and low theoretical uncertainty.The final state is a CP eigenstate, to which both B 0 and B 0 can decay.The interference between their direct and indirect decays via B 0 -B 0 mixing leads to a time-dependent CP asymmetry given by Here Γ(B 0 → J/ψ K 0 S ) represents the rate of particles that were produced as B 0 decaying to J/ψ K 0 S at proper time t, ∆m d is the oscillation frequency of the B 0 , and β is an angle of the "unitarity triangle" of the quark mixing matrix, given by the following combination of matrix elements: arg(−V cd V * cb /V td V * tb ) [4].Information can be obtained indirectly about the value of sin 2β, within the context of the Standard Model, from the combination of other measurements that constrain the matrix elements, such as ∆m d , charmless B decays and CP violation in the K sector.Many such fits have been made, typically preferring large positive values for sin 2β in the range 0.4-0.8[5]; a recent example gave sin 2β = 0.75 ± 0.09 [6].
The first published attempt at a direct measurement was made by OPAL [7].They selected 24 candidates with an estimated purity of 60% and reported a value outside the physical region, sin 2β = 3.2 +1.8 −2.0 ± 0.5, which was nevertheless interpreted as favouring large positive values.CDF published an analysis based on 395 candidates with a purity of about 40%, although half of the sample has poor proper-time determination [8].They measured sin 2β = 0.79 +0.41  −0.44 (statistical and systematic errors combined).The key to making such a measurement at LEP is to keep the efficiency as high as possible.Using the latest branching ratio B(B 0 → J/ψ K 0 ) = (8.9± 1.2) × 10 −4 [2], in the complete dataset of ALEPH, about 30 signal events are expected before the reconstruction efficiency is applied.The production state of the B 0 must also be determined (or "tagged").The precision on sin 2β scales as 1/(1 − 2w), where w is the mistag rate, the fraction of incorrectly tagged events.Using a neural-network technique to combine the information from many observables, the lowest possible mistag rate is aimed for, whilst providing a tag for every event.
In this paper, after a brief description of the ALEPH detector, the event selection is discussed.Details are given of the proper-time measurement, and the production-state tagging.The unbinned likelihood fit for the asymmetry is then presented, followed by a discussion of checks and systematic uncertainties.

Detector
A detailed description of the ALEPH detector can be found in [9] and its performance in [10].Charged particles are tracked in a two-layer silicon vertex detector with doublesided readout (r-φ and z), surrounded by a cylindrical drift chamber and a large time projection chamber (TPC), together measuring up to 33 space points along the trajectory.These detectors are immersed in a 1.5 T axial magnetic field, providing a transverse momentum resolution of ∆p/p = (6 × 10 −4 ) p at high momentum (for p in GeV/c) and a three-dimensional impact parameter resolution of 25 µm.The TPC also allows particle identification to be performed through the measurement of specific ionization (dE/dx).A finely segmented electromagnetic calorimeter of lead/wire-chamber sandwich construction surrounds the TPC.Estimators R T , R L and R I are formed for electron identification, for the transverse and longitudinal shower shape in the calorimeter and for the dE/dx in the TPC, respectively; they are calculated as the difference between the measured and expected value for electrons, divided by the expected uncertainty.The iron return yoke of the magnet is instrumented with streamer tubes to form a hadron calorimeter and is surrounded by two additional double layers of streamer tubes to aid muon identification.
The LEP1 data were recently reprocessed using improved reconstruction algorithms.A new pattern recognition algorithm for the vertex detector allows groups of nearby tracks to be analysed together, searching for hit assignments that minimize the overall χ 2 of the event.Information from the wires and pads of the TPC are also combined to improve the spatial and dE/dx resolution [11].
Monte Carlo simulated events are used to study both the signal and the background.The simulation is based on JETSET [12] and is described in detail in [13].To tune the selection cuts for background suppression, 6.5 million hadronic Monte Carlo events are used, corresponding to about 1.5 times the data statistics.In addition, a large sample of signal Monte Carlo events is used for the determination of the expected signal mass distribution, reconstruction efficiency, and for training of the neural network for tagging.

Event selection
Data taken by ALEPH in the years 1991-95 at the Z resonance are used, corresponding to 4.2 million hadronic Z decays.Hadronic events are selected in the data as described in [14].The production vertex position is reconstructed on an event-by-event basis using the constraint of the average beam-spot position [10].
First, the J/ψ → ℓ + ℓ − reconstruction is performed.The daughter tracks are required to have momentum greater than 2 GeV/c, distance of closest approach to the primary vertex of less than 2 cm transverse to and 10 cm along the beam axis, four or more hits in the TPC, and polar angle satisfying | cos θ | < 0.95.All oppositely-charged 2  S mass (a) for hadronic Monte Carlo, (b) for the data; the superimposed fits are described in the text, and the dashed lines indicate the signal region.In (a) the events that truly originated from a B 0 → J/ψ K 0 S decay are shaded, and the insert shows the mass distribution for a large sample of signal Monte Carlo events.pairs of such tracks, with opening angle satisfying cos θ ℓℓ > 0.85, are investigated for lepton identification.They must both be identified as muons or electrons using loose identification criteria.For muons, a pattern of hits in the HCAL consistent with a muon is required [15]; for electrons, cuts are made on the estimators: |R T | < 4, R L > −3 and R I > −4.The invariant mass of the lepton pair is required to be in the range 2.6-3.3GeV/c 2 , with the large window around the J/ψ mass (particularly on the low side) maintaining high efficiency in the presence of radiative decays or bremsstrahlung.
Next the K 0 S → π + π − reconstruction is performed, as described in [16].The distance of closest approach of the two daughters must be less than 5 mm, and the dE/dx for each daughter is required to be within three standard deviations of the expected value for a pion.The angle between the reconstructed directions of the K 0 S and J/ψ is required to satisfy cos θ ψK > 0.85, and the resultant of their momenta must satisfy p B > 22 GeV/c.The reconstructed K 0 S invariant mass is required to be within 15 MeV/c 2 of its nominal value.
A fit is made for the B 0 decay vertex, using the two lepton tracks and the K 0 S .The vertex is required to be successfully reconstructed, and the decay length is determined from the distance between the production and decay vertices, projected along the momentum vector of the B 0 candidate.The decay length is required to be greater than −1 mm, and its calculated error less than 1 mm.Finally, the J/ψ K 0 S invariant mass is calculated, assigning the nominal masses to the two particles.The resulting distribution is shown as an insert in Fig. 1 (a) for signal Monte Carlo.The peak is fitted with the sum of two Gaussian functions, the first accounting for 72% of the events with a width of 20 MeV/c 2 and the second with width 90 MeV/c 2 , offset to lower mass by 50 MeV/c 2 .The signal region is defined as 5.05-5.45GeV/c 2 , and the overall reconstruction efficiency is 28%.
When the event selection is applied to the hadronic Monte Carlo sample, the resulting . There are 20 events in the signal region, of which 11 are from background.The superimposed fit is the sum of an exponential shape to describe the background and the signal shape discussed above.The four fitted parameters are the background slope and normalization, and the signal mass and normalization.
The mass distribution that results when the event selection is applied to the data is shown in Fig. 1 (b).A clear signal is seen for the B 0 , with 23 events in the signal region.A fit similar to that of the hadronic Monte Carlo is made, giving a fitted B 0 mass of 5.26 ± 0.01 GeV/c 2 , slightly lower than, but consistent with, the world-average value of m B = 5.279 GeV/c 2 .The fit is used to assign an event-by-event background probability for each event in the signal region, which is used in the maximum likelihood fit for the CP asymmetry.Their sum corresponds to 6.6 background events, giving an average background fraction f bkg = 0.29 ± 0.06, where the uncertainty is estimated by varying the parametrizations within their statistical errors, and using alternative shapes for the background.The fitted shape and normalization of the background are consistent with those seen in the Monte Carlo.
After subtraction of the background, about 16 signal events remain, to be compared with the predicted number of 9 ± 2, calculated from the expected production rate and decay branching ratios.Figure 2 shows a particularly clean signal candidate, in which there are no other charged tracks in the signal hemisphere.

Proper-time determination
The decay length of the B 0 is determined from vertex reconstruction as described in the previous section.The decay-length resolution determined using the Monte Carlo simulation is reasonably described by a Gaussian distribution of width 200 µm, although there are small tails from events with poorly measured vertices.The uncertainty on the decay length is estimated by propagating the production and decay vertex errors, which are calculated in turn from the errors of the tracks used in their determination.The pull distribution, given by the difference of true and reconstructed decay lengths divided by the calculated decay-length uncertainty, is close to being normally distributed.It has a fitted Gaussian width of 1.1, indicating that the event-by-event estimate of decay-length uncertainty is reasonably accurate.
As the signal events are fully reconstructed, the momentum resolution is excellent.There is a Gaussian core of 1% relative error, with an overall RMS of 2.5%.This small uncertainty on the momentum measurement gives a significant contribution to the overall proper-time uncertainty only for events with long decay lengths, greater than about 1 cm.
The proper time is given by t = d B m B /p B .Its uncertainty is calculated as follows: The first term in parentheses is taken from the event-by-event measurement of the decay length d B and its error σ d , scaled by a factor 1.1 ± 0.2, where the uncertainty is taken into account for systematic studies; the second term is taken as (2 ± 2)%.Despite being conservative, these estimates of the uncertainty lead to a negligible effect on the measured asymmetry, as the characteristic scale of the proper-time development of the asymmetry is 1/∆m d ∼ 2 ps, much longer than the typical proper-time resolution of 0.1 ps.
Of the 11 background events in the signal region for the hadronic Monte Carlo, nine involve tracks from b-hadron decays.The apparent proper time and its error are determined for these background events, following the procedure discussed above, and a fit made for their effective lifetime τ bkg .The probability density function is taken as the convolution of an exponential lifetime distribution with a Gaussian resolution function, and the fit gives τ bkg = 1.4 +0.5 −0.3 ps.

Production-state tagging
The extraction of the CP asymmetry from the signal events requires the determination of whether they originated from a B 0 or B 0 at production.This is achieved by studying the properties of the rest of the event, excluding the lepton and pion pairs that come from the signal B 0 decay.For this purpose, two hemispheres are defined with respect to the thrust axis of the event, determined using both charged and neutral energy-flow objects [10].
These are used to separate information from the same and opposite sides of the event with respect to the B 0 candidate.Properties of the opposite side are used to determine the particle/antiparticle nature of the other b hadron that was produced in conjunction with the signal B 0 , and thus indirectly determine its production state.The same side carries information from the fragmentation process that produced the signal B 0 , which can also be used in the production-state determination.The calculation of the tag for the opposite side starts with the search for a secondary vertex due to b-hadron decay.This is achieved using a topological vertexing algorithm that combines information from all charged tracks in the hemisphere.Jets are reconstructed from the charged tracks and neutral energy-flow objects using the JADE algorithm [17], with a jet-resolution parameter of 0.02.The highest energy jet and the jet which forms the highest invariant mass with it are selected, and the secondary vertex is constrained to lie along the direction of the selected jet in its hemisphere, within errors.Each track is then assigned a relative probability P v that it originates from the secondary vertex.
Next, the b-hadron flight direction is estimated.Jets are reconstructed with a jetresolution parameter lowered to 0.0044, as this gives an improved estimate of the b-hadron direction [15].If more than one jet is found in the hemisphere considered, the b-jet candidate is chosen on the basis of the kinematic properties of its tracks and the presence of lepton candidates.The leading track of the b-jet candidate is then used as a seed for a second cone-based jet algorithm [18], which is taken as a first estimate of the b-hadron flight direction.In two-jet events, the thrust axis is chosen instead.The uncertainties on the reconstructed angles are parametrized from the simulation, as a function of the jet momentum.A second estimate of the direction is taken as the vector joining the primary and secondary vertices, and its error is parametrized as a function of the measured decay length.The two estimates are averaged using their parametrized errors, and the result taken as the b-hadron flight direction.
To construct the charge estimators, b-decay and fragmentation tracks must be distinguished.This is achieved by combining three variables for each track in the hemisphere: the rapidity and longitudinal momentum relative to the b-hadron flight direction, and the vertex assignment probability P v .They are combined in a neural network, along with a hemisphere-based indicator of the quality of the inclusive vertex finding, which acts as a "control variable".This does not directly discriminate between the b-decay and fragmentation tracks, but improves the overall performance of the neural network as the discriminating power of the other inputs vary as a function of the control variable.The single track separation output is shown in Fig. 3 (a), using neural networks trained separately for classes of tracks passing different quality criteria.The neural-network output value x sec is converted to a probability P sec that the track comes from the secondary vertex.The b-hadron momentum is then determined from the sum of charged energy assigned to the secondary vertex (based on P sec weights), the longitudinal fraction of neutral energy in the jet, and the missing energy in the hemisphere.
Nine charge-sensitive inputs are used for the opposite-side production-state tag: 1. Jet charge (κ = 0.5): where i runs over the charged tracks in the hemisphere, q i and p i are the charge and momentum of the track, d J is the b-hadron flight direction.
3. Charge of the track with highest longitudinal momentum in the hemisphere, calculated relative to the b-hadron flight direction.
4. Total charge: Q tot = i q i where i runs over the charged tracks in the b jet.

5.
Secondary vertex charge: Q vtx = i P i sec q i where i runs over all charged tracks in the hemisphere.

Weighted primary vertex charge:
where i runs over all charged tracks in the hemisphere and κ = 1.0.

7.
Weighted secondary vertex charge: as for Q pri but replacing (1 − P i sec ) with P i sec , and taking κ = 0.3.8. Decay kaon charge: kaons are identified using dE/dx, based upon the ratio of probabilities that the measured ionization is due to a kaon, relative to either a kaon or a pion.This variable is combined with the track momentum, longitudinal momentum and P sec value to select kaons from the b → c → s cascade using a further neural network.The reconstructed b-hadron momentum in the hemisphere is also used as a control variable.The charge of the track with the highest output from the neural network is used to sign the output value.9. Lepton charge: lepton identification is performed with tighter requirements than those described previously [11].If more than one lepton candidate is selected, that with the highest p T with respect to the jet axis is used.The lepton transverse momentum is signed by the track charge.Finally, six additional control variables are used: cos θ thrust , the charged track multiplicity, the spread in the measured values of P sec : σ sec = i P sec (1 − P sec ), the reconstructed b-hadron momentum, the reconstructed decay length, and the lepton momentum (if a lepton has been selected).These are combined along with the charge estimators using a neural network, which takes into account correlations between the variables.It provides a single output, x opp , which is shown in Fig. 3 (b).
Information on the production state from the same hemisphere as the B 0 candidate is limited to the charged tracks from fragmentation.Excluding the tracks from the B 0 , seven charge estimators are constructed using tracks coming from the primary vertex in the signal hemisphere.In a similar way to the opposite-side analysis, these include two jet charges with κ values of 0.5 and 1.0, the charge of the track with the highest longitudinal momentum with respect to the jet, one momentum-weighted and two longitudinalmomentum-weighted primary vertex charges with κ values of 1.0, 0.3 and 1.0, respectively, and finally the sum of track charges inside the jet.The output of the opposite-side neural network, x opp , is then combined with these same-side charge estimators and the following control variables: σ sec , cos θ thrust , and the same-side charged track multiplicity, using another neural network.The output of this event tag, x tag , is shown in Fig. 3 (c).
If events with x tag > 0.5 are taken to have a B 0 in the production state and those with x tag < 0.5 to have a B 0 , the fraction of incorrect tags (the average mistag rate) is 27.8%.The event-by-event value of the neural-network output x tag is used in the measurement of the asymmetry.Due to its peaked distribution, seen in Fig. 3 (c), the effective mistag rate is lower than the average value quoted above.Using an overlap integral technique [19] the effective mistag rate is found to be 24.1%,measured with an independent sample of simulated events, with an expected difference of (0.2 ± 0.2)% between the effective mistag values for B 0 and B 0 hemispheres.The opposite-side tag alone gives average and effective mistag rates of 31.4% and 27.6%, respectively.
Finally, the production-state tag q is calculated from the neural-network output x tag , correcting for the purity in each bin of Fig. 3 (c) using where F (x tag ) is the distribution of x tag for events produced as B 0 , shaded in the figure, and G(x tag ) is the distribution for B 0 ; q = +1 for events produced as B 0 and −1 for B 0 in the case of perfect tagging.The relationship between x tag and q, determined with signal Monte Carlo events, is consistent with linearity, as shown in Fig. 3 (d).It is parametrized as a straight line q = a tag + 2 b tag (x tag − 0.5), with fitted coefficients a tag = −0.01 ± 0.01 and b tag = 0.95 ± 0.01.

Asymmetry measurement
The production-state tag q is plotted against the proper time t in Fig. 4 (a), for the 23 events in the signal region of the data.Considering the events which have a clear tag result, |q| > 0.5, there are 4 "B 0 -like" events (positive q) and 9 "B 0 -like".Furthermore, the excess of events with negative q is more noticeable with increasing proper time, as would be expected for a negative CP asymmetry.
To extract a measurement of the asymmetry, an unbinned likelihood is calculated.The probability density function expected for the signal distribution is given by The B 0 lifetime τ d = 1.548±0.032ps and ∆m d = 0.472±0.017ps −1 are fixed to the central values of their world averages [2] (the uncertainties are taken into account in systematic error studies).
A convolution is made of this signal distribution and a Gaussian resolution function R(σ t ) with width given by the event-by-event proper-time resolution calculated in Proper time t (ps) Section 4. The result is illustrated in Fig. 4 (b), where the contributions from B 0 and B 0 are indicated as a function of proper time for fixed values of the asymmetry and resolution.
The probability density function for background events P bkg (t, q) is taken to have the same form as that of the signal but replacing τ d → τ bkg and sin 2β → a bkg , where the effective lifetime of the background τ bkg was determined in Section 4. The background asymmetry a bkg is taken to be zero, but is varied to study possible systematic effects.The likelihood of an event is then calculated as where f bkg is the event-by-event background fraction discussed in Section 3. A generalized likelihood function is used, adding a Poisson term to the combined likelihoods of the signal (b) Reconstructed asymmetry versus proper time, after correction for the average dilution in each time bin.The result of the fit for the asymmetry is shown superimposed; note that this is not fitted directly to the binned data shown. candidates: where n = 23 is the number of observed candidates, and ν is the number expected.A fit is made with two free parameters, sin 2β and ν, with the results sin 2β = 0.84 +0.82 −1.04 , ν = 23.3+5.1 −4.5 .A scan of the log-likelihood versus sin 2β is shown in Fig. 5 (a).The result of the measurement is displayed as a function of proper time in Fig. 5 (b).For the purpose of display the data are divided into three proper-time bins, and the average asymmetry is corrected for the average dilution (from mistagging and background) in each bin.Negative asymmetry, and thus positive sin 2β, is favoured.

Checks and systematic uncertainties
The statistical error from the fit has been checked using a fast Monte Carlo simulation.Event samples the same size as the data were generated according to Eq. 6 and 7, using the resolution function and tagging distributions discussed above.The value of sin 2β was set to that measured for the data, and all other parameters were fixed to their nominal values in the fit.The number of background events within the sample was varied according to Poisson statistics.The CP asymmetry of each sample was measured as if it were data, and the central value and errors recorded.This was repeated for many samples.The values reconstructed for the positive and negative errors of the data are consistent with the distributions of values seen in the simulated samples.Furthermore, the measured errors reasonably estimate the spread of the measured central value about its true input value: the pull is normally distributed.
The contributions to the systematic error are listed in Table 1.The lifetime and oscillation frequency of the signal were varied within their world-average uncertainties, and the effect on the measured asymmetry taken as a systematic error.For the lifetime, a check was made by fitting for the lifetime of selected signal Monte Carlo events, giving τ = 1.56 ± 0.02 ps, in agreement with the input value of 1.56 ps, indicating no evidence for a systematic bias.Similarly their asymmetry was measured to be −0.02± 0.03, in agreement with the input value of zero.The decay-length and momentum resolution parametrizations were varied as discussed in Section 4.
For the background, the uncertainty on the level was discussed in Section 3. The effective lifetime of the background was varied within the large uncertainty measured using the hadronic Monte Carlo, given in Section 4. To increase the statistics for the study of the CP asymmetry of the background, events in the side-band region 4.0 < m(J/ψ K 0 S ) < 5.0 GeV/c 2 were investigated.There are 140 such events in the data, with asymmetry 0.09 ± 0.23, indicating no significant effect.About one event is expected in the signal region from the decay B 0 → J/ψ K * 0 with K * 0 → K 0 S π 0 , predicted to be mainly CP even.This would correspond to a background asymmetry of less than 0.1.For the systematic error evaluation the asymmetry of the background, a bkg , was varied by ±0.2.
Uncertainties arising from the production-state tag can be considered as being due to knowledge of the overall mistag rate for signal events in data and the possible difference in the individual mistag rates for B 0 and B 0 events.Samples of events were isolated in data with a similar topology to that of the signal, where there is a b hadron of known charge with decay products that can be cleanly separated from fragmentation tracks.Three samples of B + candidates were selected for this purpose.The same neural network was used for tagging their production state as for B 0 → J/ψ K 0 S , except that the sign of the same-side charge estimators is reversed.This is required as a b quark combines with a d (u) quark to produce a B 0 (B + ), and the accompanying d (u) quarks have opposite sign.
1. B + → J/ψ K + decays: The signal selection is modified slightly, to select B + → J/ψ K + decays (and charge conjugate).The reconstructed J/ψ K + mass in the data is shown in Fig. 6 (a); a clear signal is seen, with 52 events in the signal region compared to an expected background of about 17.In this channel the efficiency is similar to that of J/ψ K 0 S , but the product branching ratio is about three times higher.Now the charge of the kaon indicates whether the signal was produced by B + or B − , and the resulting distributions of the neural-network output x tag are compared in Fig. 6 (b).From a study of side-band events, the background is found to have a similar tagging distribution to the signal.The average mistag rate is 26% in signal Monte Carlo and is measured to be (27±6)% in the data, giving confidence that the Monte Carlo simulation reproduces faithfully the features of the data used for tagging.The asymmetry and lifetime have also been measured for these events, with the results 0.09 ± 0.41 and 1.57 +0.23  −0.20 ps, respectively.The latter agrees with the world-average value of 1.653 ± 0.028 ps [2].calibration of the tag, the parametrization of q versus x tag .A degraded mistag rate leads to a reduced slope of the calibration, as illustrated in Fig. 3 (d).Similarly, a B 0 -B 0 tagging difference is seen as an offset to the calibration.The variations applied to the parametrization are listed in Table 1.
Adding the systematic error contributions in quadrature, the final measurement is sin 2β = 0.84 +0.82 −1.04 ± 0.16 , where the first error is statistical and the second systematic.

Conclusion
An analysis of B 0 → J/ψ K 0 S decays has been performed, with J/ψ → e + e − or µ + µ − and K 0 S → π + π − .A reconstruction efficiency of 28% is achieved.From the full dataset of ALEPH at LEP1 of 4.2 million hadronic Z decays, 23 candidates are selected with an estimated purity of 71%.They are used to measure the CP asymmetry of this decay, with the result sin 2β = 0.84 +0.84 −1.05 , where the uncertainty is dominated by the statistics.This result is compared with the other two published measurements in Fig. 8.The likelihood of each measurement is approximated as a Gaussian distribution on either side of the central value, with width equal to the sum in quadrature of statistical and systematic errors.The three log-likelihoods can be summed to give a combined result, i.e. neglecting any correlation between the systematic errors of the different experiments (expected to be small).The resulting likelihood is shown by the shaded distribution in the figure and corresponds to sin 2β = 0.88 +0.36  −0.39 .The integral of the likelihood for sin 2β > 0 is 77% for the ALEPH result alone, or 67% if the total integral is limited to the physical region |sin 2β| < 1.The corresponding fractions are 98.6% and 98.0% for the combination of the three analyses.Thus the confidence level that CP violation has been seen in this channel is increased to 98%, compared to 93% for CDF alone [8].
Preliminary results from the BABAR and BELLE experiments [21] are consistent with the combined result given above, with similar precision and slightly lower central values.

2 Figure 2 :
Figure2: Event display of a B 0 → J/ψ K 0 S candidate in the data, with J/ψ → e + e − : (a) fish-eye view of the projection transverse to the beam axis, (b) zoom into the vertex region for the same projection, with reconstructed vertices marked with ellipses, and the reconstructed neutrals marked with dashed lines.
Neural-network output x sec

Figure 3 :
Figure 3: Neural-network output distributions for signal Monte Carlo events: (a) separation of b-decay and fragmentation tracks, x sec , with the contribution from b-decay tracks shaded;(b) opposite-side production state, x opp , with the contribution from hemispheres containing b hadrons shaded; (c) overall event production state, x tag , with the contribution from B 0 decays shaded.(d) Calibration of the production-state tag q versus the neural-network output, with superimposed linear fit; the effect of degrading the mistag rate for systematic studies is indicated by the dashed line.

Figure 4 :
Figure 4: (a) Production-state tag q versus proper time t for the 23 events in the signal region of the data.The proper-time uncertainty is indicated for each point, and histograms of the two variables are shaded.The solid points indicate events with a background probability less than 30%.(b) Expected proper-time distribution of B 0 and B 0 events, for fixed values of the CP asymmetry and proper-time resolution: sin 2β = 0.7, σ t = 0.2 ps.

Figure 5 :
Figure 5: (a) Log-likelihood versus sin 2β from the fit, with respect to the minimum.

output x tag 6 Figure 6 :
Figure 6: (a) Reconstructed J/ψ K + mass, for the data.(b) Production-state tag neural-network output x tag , when applied to events from the signal region.
2. B + → D 0 ℓ + X decays: A selection is made of D 0 candidates, decaying to a kaon and a pion of opposite charge.Only those D 0 candidates are used which are found to lie in the same hemisphere as an identified lepton (e or µ) with transverse momentum greater than 0.5 GeV/c and charge opposite in sign to that of the pion from the D 0 decay.The estimated b purity of the sample is 87%, dominated by B ± decays.The lepton charge indicates whether the signal was produced by a B + or B − , and the resulting tag distributions are shown in Fig.7(a).The effective mistag

Figure 8 :
Figure 8: Likelihood distributions versus sin 2β for the three existing measurements (dashed lines) and their combination (shaded).

Table 1 :
Contributions to the systematic uncertainty on sin 2β.