Study of B0_s anti-B0_s oscillations and B0_s lifetimes using hadronic decays of B0_s mesons

Oscillations of B0s mesons have been studied in samples selected from about 3.5 million hadronic Z decays detected by DELPHI between 1992 and 1995. One analysis uses events in the exclusive decay channels: B0s ->Ds- pi+ or Ds- a1+ and B0s ->anti-D0 K- pi+ or anti-D0 K- a1+, where the D decays are completely reconstructed. In addition, B0s anti-B0s oscillations have been studied in events with an exclusively reconstructed Ds accompanied in the same hemisphere by a high momentum hadron of opposite charge. Combining the two analyses, a limit on the mass difference between the physical B0s states has been obtained: Delta(m_B0s)>4.0 ps^{-1} at the 95% C.L. with a sensitivity of Delta(m_B0s) = 3.2 ps^{-1}. Using the latter sample of events, the B0s lifetime has been measured and an upper limit on the decay width difference between the two physical B0s states has been obtained: tau(B0s) = 1.53^{+0.16}_{-0.15}(stat.) +/- {0.07}(syst.) ps \Delta\Gamma(B0s)/\Gamma(B0s)<0.69 at the 95% C.L. The combination of these results with those obtained using Ds+- lepton-+ sample gives: Delta(m_B0s)>4.9 ps^{-1} at the 95% C.L. with a sensitivity of Delta(m_B0s) = 8.7 ps^{-1}. tau(B0s) = 1.46 +/- 0.11 ps and \Delta\Gamma(B0s)/\Gamma(B0s)<0.45 at the 95% C.L.

Oscillations of B 0 s mesons have been studied in samples selected from about 3.5 million hadronic Z decays detected by DELPHI between 1992 and 1995.One analysis uses events in the exclusive decay channels: B 0 s → D − s π + or D − s a + 1 and B 0 s → D 0 K − π + or D 0 K − a + 1 , where the D decays are completely reconstructed.In addition, B 0 s − B 0 s oscillations have been studied in events with an exclusively reconstructed D s accompanied in the same hemisphere by a high momentum hadron of opposite charge.Combining the two analyses, a limit on the mass difference between the physical B 0 s states has been obtained: ∆m B 0 s > 4.0 ps −1 at the 95% C.L. with a sensitivity of ∆m B 0 s = 3.2 ps −1 .Using the latter sample of events, the B 0 s lifetime has been measured and an upper limit on the decay width difference between the two physical B 0 s states has been obtained: 16  −0.15(stat.)± 0.07(syst.)ps ∆Γ B 0 s /Γ B 0 s < 0.69 at the 95% C.L. The combination of these results with those obtained using D ± s ℓ ∓ sample gives: ∆m B 0 s > 4.9 ps −1 at the 95% C.L. with a sensitivity of ∆m B 0 s = 8.7 ps −1 .τ B 0 s = 1.46 ± 0.11 ps and ∆Γ B 0 s /Γ B 0 s < 0.45 at the 95% C.L.
(Eur.Phys.J. C18(2000)229) 1 Introduction In this paper, the average lifetime of the B 0 s meson is measured and limits are derived on the oscillation frequency of the B 0 s -B 0 s system, ∆m B 0 s , and on the decay width difference, ∆Γ B 0 s , between mass eigenstates of this system.Starting with a B 0 s meson produced at time t=0, the probability, P(t), to observe a B 0 s or a B 0 s decaying at the proper time t can be written, neglecting effects from CP violation: where s ; L and H denote the light and heavy physical states, respectively; ∆Γ B 0 s and ∆m B 0 s are defined to be positive [1] and the plus (minus) signs in equation (1) refer to B 0 s (B 0 s ) decays.The oscillation period gives a direct measurement of the mass difference between the two physical states.The Standard Model predicts that ∆Γ B 0 s ≪ ∆m B 0 s and the previous expression simplifies to : for B 0 s → B 0 s and similarly: for B 0 s → B 0 s .The oscillation frequency, proportional to ∆m B 0 s , can be obtained from the fit of the time distributions given in relations (2) and (3), whereas expression (1), without distinguishing between the B 0 s and the B 0 s , can be used to determine the average lifetime and the difference between the lifetimes of the heavy and light mass eigenstates 1 .
B physics allows a precise determination of some of the parameters of the CKM matrix.All the nine elements can be expressed in term of 4 parameters that are, in Wolfenstein parameterisation [2], λ, A, ρ and η.The most uncertain parameters are ρ and η.
Several quantities which depend on ρ and η can be measured and, if the Standard Model is correct, they must define compatible values for the two parameters inside measurement errors and theoretical uncertainties.These quantities are ǫ K , the parameter introduced to measure CP violation in the K system, |V ub |/|V cb |, the ratio between the modulus of the CKM matrix elements corresponding to b→u and b→c transitions, and the mass difference ∆m B 0 d .In the Standard Model, B 0 q -B 0 q (q = d, s) mixing is a direct consequence of second order weak interactions.Having kept only the dominant top quark contribution, ∆m B 0 q can be expressed in terms of Standard Model parameters [3]: In this expression G F is the Fermi coupling constant; F (x t ), with , results from the evaluation of the second order weak "box" diagram responsible for the mixing and has a smooth dependence on x t ; η B is a QCD correction factor obtained at next-to-leading order in perturbative QCD [4].The dominant uncertainties in equation ( 4) come from the evaluation of the B meson decay constant f Bq and of the "bag" parameter B Bq .The two elements of the V CKM matrix are equal to: neglecting terms of order O(λ 4 ).
In the Wolfenstein parameterisation, |V ts | is independent of ρ and η.A measurement of ∆m B 0 s is thus a way to measure the value of the non-perturbative QCD parameters.Direct information on V td can be inferred by measuring ∆m B 0 d .Several experiments have accurately measured ∆m B 0 d , nevertheless this precision cannot be fully exploited to extract information on ρ and η because of the large uncertainty which originates in the evaluation of the non-perturbative QCD parameters.
An efficient constraint is the ratio between the Standard Model expectations for ∆m B 0 d and ∆m B 0 s , given by: ∆m A measurement of the ratio ∆m B 0 d /∆m B 0 s gives the same type of constraint, in the ρ − η plane, as a measurement of ∆m B 0 d , but is expected to be better under control from theory, since the ratios of the decay constants f Bq and of the bag parameters B Bq for B 0 d and B 0 s are better known than their individual values [5].
Using existing measurements which constrain ρ and η, except those on ∆m B 0 s , the distribution for the expected values of ∆m B 0 s can be obtained.It has been shown, for example, in the context of the Electroweak Standard Model and QCD assumptions, that ∆m B 0 s has to lie, at 68% C.L., between 8 ps −1 and 16 ps −1 and is expected to be smaller than 21 ps −1 at the 95% C.L. [6].
The search for B 0 s − B 0 s oscillations has been the subject of recent intense activity.No signal has been observed so far and the lower limit on the oscillation frequency comes from the combination of the results obtained at LEP, CDF and SLD experiments: ∆m B 0 s > 14.3 ps −1 at the 95% C.L. [7].The sensitivity of present measurements is estimated as 14.7 ps −1 .These results have been obtained by combining analyses which select various B 0 s decay channels 2 .Some analyses [8,9] use events containing simply a lepton emitted at large transverse momentum relative to the axis of the jet from which it emerges.In this case, the proper time is measured using an inclusive vertex algorithm to reconstruct the decay distance and the energy of the B hadron candidate.In other analyses, like D ± s ℓ ∓ [10,11], the identified lepton is accompanied, in the same hemisphere, by an exclusively reconstructed D s .In D ± s h ∓ analyses (see [12] and the present paper), such leptons are replaced by one or more charged hadrons.
Progress before the next generation accelerators is expected to come from improved analyses of these channels, but the sensitivity at high frequency is essentially limited by the damping of the reconstructed oscillation amplitude due to the limited resolution on the B 0 s proper time.In general, the proper time resolution σ t contains two terms and one of them is a time-dependent: where σ L is related to the decay distance resolution, σ p /p is the relative error on the momentum reconstruction and t is the proper decay time.The damping factor of the oscillation amplitude in this case is given by the following expression [13]: where x=1/ √ 2/∆m B 0 s 1/(σ p /p) and ERFC(x)=2/ √ π ∞ x e −t 2 dt.Exclusive (completely reconstructed) decays have a better time resolution for two reasons.As there is no missing particle in the decay, the B 0 s momentum is known with a good precision and therefore the contribution of the momentum uncertainty to the proper time resolution is negligible.Hence, σ t = σ L and the previous expression is simplified: In addition, the reconstructed channels are two-body or quasi two-body decays, with an opening angle of their decay products which is on average larger than in multi-body final states; this results in a better accuracy on their decay position determination.The B 0 s meson lifetime is expected to be equal to the B 0 d lifetime [14] within one percent.In the Standard Model, the ratio between the decay width and the mass differences in the B 0 -B 0 system is of the order of (m b /m t ) 2 , although large QCD corrections are expected.Explicit calculations to leading order in QCD correction, in the OPE formalism [1] predict: where the quoted error is dominated by the uncertainty related to hadronic matrix elements.Recent calculations [15] at next-to-leading order predict a lower value: 063 .An interesting approach consists in using the ratio between ∆Γ B 0 s and ∆m B 0 s [15]: to constrain the upper part of the ∆m B 0 s spectrum with an upper limit on ∆Γ B 0 s /Γ B 0 s .If, in future, the theoretical uncertainty can be reduced, this method can give an alternative way of determining ∆m B 0 s via ∆Γ B 0 s and, in conjunction with the determination of ∆m B 0 d , can provide an extra constraint on the ρ and η parameters.
After a description of the main features of the DELPHI detector, the event selection and the event simulation in Section 2, the two B 0 s analyses are described.Section 3 presents an analysis of B 0 s -B 0 s oscillations from a sample of 44 exclusively reconstructed B 0 s decays.Section 4 is dedicated to the D ± s h ∓ analysis, in which B 0 s → D s hX decays with fully reconstructed D s are selected.This analysis also includes a measurement of the B 0 s meson lifetime and sets a limit on the difference between the decay widths of the physical B 0 s states, ∆Γ B 0 s /Γ B 0 s .In Section 5 the combined limit on ∆m B 0 s is given.Finally in Section 6, these three results are combined with those obtained using a D ± s ℓ ∓ sample [10].

The DELPHI detector
The events used in this analysis have been recorded with the DELPHI detector at LEP operating at energies close to the Z peak.The DELPHI detector and its performance were described in detail elsewhere [16].In this section, components of the detector and their characteristics, which are the most relevant for these analyses are summarised.

Charged particle reconstruction
The detector elements used for tracking are the Vertex Detector (VD), the Inner Detector (ID), the Time Projection Chamber (TPC) and the Outer Detector (OD).The VD provides the high precision needed near the primary vertex.
For the data taken from 1991 to 1993, the VD consisted of three cylindrical layers of silicon detectors (at radii of 6.3, 9.0 and 10.9 cm) measuring points in the plane transverse to the beam direction (rφ coordinate3 ) in the polar angle range between 43 • and 137 • .In 1994, the first and the third layers were equipped with detector modules with double sided readout, providing a single hit precision of 7.6 µm in rφ, similar to that obtained previously, and 9 µm in z [17].For high momentum particles with associated hits in the VD, the impact parameter precision close to the interaction region is 20 µm in the rφ plane and 34 µm in the rz plane.Charged particle tracks are reconstructed with 95% efficiency and with a momentum resolution σ p /p < 1.5 × 10 −3 p (p in GeV/c) in the polar angle region between 25 • and 155 • .

Hadronic Z selection and event topology
Hadronic events from Z decays were selected by requiring a multiplicity of charged particles greater than four and a total energy of charged particles greater than 0.12 √ s, where √ s is the centre-of-mass energy and all particles were assumed to be pions; charged particles were required to have a momentum greater than 0.4 GeV/c and a polar angle between 20 • and 160 • .The overall trigger and selection efficiency is (95.0±0.1)%[18].A total of 3.5 million hadronic events was obtained from the 1992-1995 data.Each selected event was divided into two hemispheres separated by the plane transverse to the sphericity axis.A clustering analysis, based on the JETSET algorithm LUCLUS with default parameters [19], was used to define jets using both charged and neutral particles.These jets were used to compute the p out t of each particle of the event, as the transverse momentum of this particle with respect to the axis of the jet it belonged to, after having removed this particle from its jet.

Hadron identification
Hadron identification relied on the RICH detector and on the specific ionisation measurement performed by the TPC.
The RICH detector [16] used two radiators.A gas radiator separated kaons from pions between 3 and 9 GeV/c, where kaons gave no Cherenkov light whereas pions did, and between 9 and 16 GeV/c, using the measured Cherenkov angle.It also provided kaon/proton separation from 8 to 20 GeV/c.A liquid radiator, which has been fully operational since 1994, provided p/K/π separation in the momentum range between 0.7 and 8 GeV/c.
The specific energy loss per unit length (dE/dx) was measured in the TPC by using up to 192 sense wires.At least 30 contributing measurements were required to compute the truncated mean.In the momentum range between 3 and 25 GeV/c, this is fulfilled for 55% of the tracks and the dE/dx measurement has a precision of ±7%.
The combination of the two measurements, dE/dx and RICH, provides four levels of pion, kaon and proton tagging [20]: "very loose", "loose", "standard" and "tight" corresponding to different purities.The efficiency and purity of the tagging depend on the momentum of the particle.In the typical momentum range of the B decay products between 2 and 10 GeV/c, the average efficiency (purity) is about 75% (50%) for the "loose", 65% (60%), for the "standard" and 55% (70%) for the "tight" kaon tag.The "very loose" tag indicates that the particle is not identified as a pion.The purity is the proportion of genuine kaons in the selected sample.
Very recently a neural network algorithm has been developed in order to combine the different RICH identification packages optimally.The result of the neural network gives a tagging variable x net , which varies from −1 (pure background) to +1 (pure signal).In order to use the same definitions of the kaon tag through all the paper, the "very loose", "loose", "standard" and "tight" tags were defined for a tagging variable x net larger than −0.2, 0.0, 0.4 and 0.6, respectively.This second algorithm has been used in the D ± s h ∓ analysis (Section 4).

Primary and B vertices reconstruction
The average beam spot position evaluated using the events from different periods of data taking has been used as a constraint in the determination of the e + e − interaction point on an event-by-event basis [16].In 1994 and 1995 data, the position of the primary vertex transverse to the beam is determined with a precision of about 40 µm in the horizontal direction, and about 10 µm in the vertical direction.For 1992 and 1993 data, the uncertainties are larger by about 50%.
For both 1992-1993 and 1994-1995 data, the B decay length was estimated as L = L xy / sin θ B , where L xy is the measured distance between the primary vertex and the B decay vertex in the plane transverse to the beam direction and θ B is the polar angle of the B flight direction, estimated from the B decay products.The transverse decay length L xy was given the same sign as the scalar product of the B momentum with the vector joining the primary to the secondary vertex; this procedure gives signed decay length L.

Λ and K 0 s reconstruction
The Λ → pπ − and K 0 s → π + π − decays have been reconstructed using a kinematic fit.The distance in the rφ plane between the V 0 decay point and the primary vertex was required to be less than 90 cm.This condition meant that the decay products had track segments of at least 20 cm long in the TPC.The reconstruction of the V 0 vertex and selection cuts are described in detail in reference [16].Only K 0 s candidates passing the "tight" selection criteria defined in [16] were retained for this analysis.

b-tagging
The b-tagging package is described in reference [21].The impact parameters of the charged particle tracks, with respect to the primary vertex, were used to build the probability that all tracks come from this vertex.Due to the long B hadron lifetimes, the probability distribution is peaked near zero for events which contain beauty quark whereas it is flat for events containing only light quarks.

Event simulation
Simulated events were generated using the JETSET 7.3 program [19] with parameters tuned as in [22] and with an updated model of B decay branching fractions.B hadron semileptonic decays were simulated using the ISGW model [23].Generated events were passed through the full simulation of the DELPHI detector [16], and the resulting simulated raw data were processed through the same reconstruction and analysis programs as the real data.
, where the cascade photon or the neutral pion has not been reconstructed.All corresponding branching fractions are unmeasured.An estimation will be used in the following to evaluate the number of the expected B 0 s events.To evaluate the two-body B 0 s → D ( * )− s π + (a + 1 ) branching fractions, the equivalent decay channels for non-strange mesons were used.This evaluation is based on the numerical application of factorisation done in reference [24].
For each decay mode of the B 0 s → D ( * )− s M + , where M + is a π + or a ρ + , the branching fraction is estimated using the following formula: where the notation D ( * )− means D * (2010) − or D − meson.The superscript "exp" designates the experimentally measured values [25], and "th" the predictions from reference [24].The theoretical ratios in equation (10) are close to 1.
For the two-body decays containing D mesons with a 1 , the theoretical branching fractions are not available, and the same factors as for the ρ + channels have been assumed.
, which are described in the Section 3.4 of the current paper.The ratio of B 0 s to B 0 d branching fractions was taken equal to 1.1 (this factor corresponds to the average ratio of the four previous decay modes).Finally, the factor 3.3 was used as the ratio between the B 0 s decay modes containing a + 1 and π + mesons.This factor was taken as an average between the experimental values for the following B 0 d two-body decays: The theoretical calculation [26] gives for this ratio values between 3.4 and 3.5.Table 1 presents the evaluations of several branching fractions of interest for this paper.

Event Sample
Events were selected using the following decay channels: where the φ, K 0 s , f(980), K * 0 , K * − and a 1 are reconstructed in their charged decay channels: and D 0 mesons were reconstructed by considering charged particles in the same hemisphere with at least one VD hit.D meson candidates were accepted if their mass was within the intervals 1.93 − 2.01 GeV/c 2 for D s and 1.83 − 1.90 GeV/c 2 for D 0 .The D decay length was required to be positive and the χ 2 -probability of the fitted vertex to be larger than 10 −5 (10 −3 for the D 0 → Kπππ decay mode).Different selection criteria were used for different decay channels according to the optimisations derived from dedicated simulated samples.They are described in the following: The φ meson was reconstructed in the decay mode φ → K + K − by taking all possible pairs of oppositely charged particles if at least one of them was identified as a "very loose" kaon.The invariant mass of these pairs had to be within ±12 MeV/c 2 of the nominal φ mass value [25].The momenta of all three particles were required to be larger than 1 GeV/c.In the decay of the D s meson into a vector (φ) and a pseudoscalar meson (π), helicity conservation requires that the angle ψ, measured in the vector meson rest frame between the directions of its decay products and of the pseudoscalar meson, has a cos 2 ψ distribution.The value of | cos ψ| was required to be larger than 0.3.

•D −
s → φπ − π + π − The φ meson was reconstructed as in the previous channel.The momenta of all three pions had to be larger than 0.6 GeV/c.

•D −
s → f(980)π − The f(980) meson was reconstructed in the decay mode f(980) → π + π − by taking all possible pairs of oppositely charged particles classified as pions.This channel suffers from a large combinatorial background which was reduced by selecting candidates having an invariant π + π − mass within 15 MeV/c 2 of the nominal f(980) mass [25] and a total momentum larger than 8 GeV/c.The momenta of all three pions had to be larger than 1 GeV/c.
s meson reconstruction has been described in Section 2.5.In addition, the decay length of K 0 s candidates had to be positive and their momentum to be larger than 2.5 GeV/c.The momentum of the K − candidate, identified as a "very loose" kaon, had to be larger than 1 GeV/c.

•D −
s −→ K * 0 K − The K * 0 meson was reconstructed in the charged decay mode K * 0 → K + π − .The K + candidate was required to be identified as a "very loose" kaon.The momenta of both particles had to be larger than 1 GeV/c and the invariant mass of the pair had to be within ±40 MeV/c 2 of the nominal K * 0 mass [25].The value of | cos ψ| (see the D − s → φπ − selection) had to be larger than 0.3.For the K − candidate, the combinatorial background is higher than for the kaon coming from K * 0 resonance.The momentum of the K − candidate from D − s had to exceed 2.5 GeV/c and it had to be identified as a "loose" kaon.
The K * 0 meson was reconstructed as previously, but with an invariant mass within ±60 MeV/c 2 of the nominal K * 0 mass [25].This mass interval was chosen larger than in the previous K * 0 selection, because of the additional constraint from the K * − mass.The K * − meson was reconstructed in the decay mode K * − → K 0 s π − .The K 0 s meson reconstruction was discussed previously.The invariant mass of K * − candidates had to be within ±60 MeV/c 2 of the nominal K * − mass [25].
The D 0 meson decay to K + π − has been reconstructed by combining a kaon candidate, identified with the "loose" tag, with an oppositely charged pion with momentum larger than 1 GeV/c.The D 0 meson decay to K + π − π + π − has been reconstructed by combining a kaon candidate, identified with the "standard" tag, with three pions, each of them having a momentum larger than 0.5 GeV/c.In order to reduce the combinatorial background, for both decay modes, the kaon momentum was required to be larger than 2.5 GeV/c and the D 0 momentum candidate to be larger than 10 GeV/c.
Selected D 0 and D − s mesons were used to reconstruct B 0 s candidates by fitting a common vertex for the D The π momentum had to be larger than 4 GeV/c.The a 1 candidates were reconstructed using the combination of three pions with momenta larger than 0.8 GeV/c and with an invariant mass situated within the interval 0.95-1.50GeV/c 2 .At least one of the two π + π − combinations was required to have an invariant mass lying within ±150 MeV/c 2 of the nominal ρ mass [25].The a 1 momentum had to be larger than 5 GeV/c (6 GeV/c for the D 0 → K + π − π + π + channel).For all candidates, the D 0 and D − s meson decay distance, relative to the B vertex had to be positive.Events with an estimated error on the B 0 s decay distance larger than 250 µm and those having a vertex χ 2 -probability smaller than 10 −3 were removed.In order to reduce the combinatorial background from charm and light quarks, the b-tagging probability for the whole event and for the hemisphere opposite to the reconstructed B meson had to be smaller than 0.1.Additional selections which depend on the B 0 s decay channel were applied, mainly for D 0 decays which suffer from a larger combinatorial background than D − s candidates: The momentum of the B 0 s had to be larger than 22 GeV/c.For B 0 s → D − s a + 1 candidates, only the combination with the largest B 0 s momentum has been kept.
decays, the B 0 s momentum had to be larger than 27 GeV/c.Because of a high combinatorial background for B 0 s → D 0 K − a + 1 decays, the B 0 s momentum had to be larger than 29 GeV/c in the D 0 → K + π − channel and larger than 33 GeV/c in the In each event, only the candidate with the largest B 0 s momentum was kept.For the K − candidate, identified as "standard", the momentum had to be larger than 2 GeV/c.The source of the D 0 and K − meson pair can be an excited D meson state.The constraint on the mass value of such a state allows more soft K − mesons to be considered.The selection on the K − momentum was reduced to 1.5 GeV/c, if the invariant mass of the D 0 K − pair was compatible with the mass of the orbitally excited D sJ (2573) state (the most probable value of J is 2 [25]).One event was detected, which is compatible with the decay channel B 0 s → D sJ (2573  1 and the corresponding reconstruction efficiencies.These efficiencies vary from (0.2 ± 0.1)% for the The fitted or expected number of signal events and the fraction of combinatorial background are given inside a mass window corresponding to ±3σ(±2σ) for the main(satellite) peaks.In real data the number of events in these mass windows are 11 (33).The B 0 s mass in the simulation is 5.370 GeV/c 2 .In the last column the fraction of physical background, discussed in detail in Section 3.3, is given.
The mass spectrum, obtained by summing up the contributions from the different channels is shown in Figure 1.The data are indicated by points with error bars and the fit result is shown by the solid line.The mass distribution in the signal region was fitted using two Gaussian functions of different widths to account for the exclusive B 0 s signal (main narrow peak) and for the presence of partially reconstructed B 0 s decays (satellite wider peak).According to the expected branching fractions, it was assumed that the dominant contributions for the satellite peak come from the following decay channels: , where the cascade photon or neutral pion (or both) were not reconstructed.In the simulation, it was verified that the mass distribution of the satellite peak can be described by a Gaussian function as shown in Figure 2. The central values and the widths of the Gaussian functions (main and satellite peaks) in Figure 1 were left free to vary in the fit.The combinatorial background was fitted using an exponential function with a slope fixed according to the simulation.This slope was verified in several ways.Using 2.9 million bb and 5.3 million qq simulated events, the expected mass spectrum was obtained after having removed the B 0 s signal contribution in the mass region of the satellite and main peaks.The sum of the contributions from b events, charm events and light-quark events is shown in Figure 1 and is in agreement with that obtained in data.Two further checks were performed, using events selected in the side-bands of the D − s and D 0 candidates and from wrong sign combinations.The slopes of both distributions are in agreement with those obtained by fitting the data and the simulation.
The fit to the mass distribution yielded a signal of 8±4 B 0 s decays in the main peak and 15 ± 8 events in the satellite peak.The probability that the background has fluctuated to give the observed number of events in the main peak is 3 × 10 −4 .Table 2 gives the characteristics of the observed signals and the comparison with simulation.Specific decay channels from B 0 s , B 0 d , B + and Λ 0 b were simulated with statistics ranging from ten to several thousand times the expected rates in real data.These samples were used to determine efficiencies and reflection probabilities discussed below.

Reflections from B 0 d and Λ 0 b decays
For several B 0 s decay channels a possible physical background originates from nonstrange B decays, when one of the pions or proton is misidentified as a kaon (kinematic reflections).The main decay channels are B In the mass region of the main peak, 0.32 ± 0.13 events are expected to originate from kinematic reflections (with a 19% contribution from Λ 0 b decays).These estimates were obtained using dedicated Monte Carlo samples of B 0 d and Λ 0 b events passed through the full reconstruction algorithms and satisfying the selection criteria.The corresponding uncertainties come from the limited knowledge of the assumed branching fractions and from the simulated events statistics.
To study the contribution in the satellite peak from kinematic reflections, the same decay processes were considered with channels in which the D meson is accompanied by a ρ.This gives 1.3 ± 0.7 events.
In order to look for possible signals coming from kinematic reflections in real data, B 0 s candidates were considered in turn as B 0 d or Λ 0 b hadrons by changing the kaon into a pion or an antiproton, respectively.The mass distributions obtained are similar to those expected from genuine simulated B 0 s mesons and do not show any accumulation of events in the B 0 d or Λ 0 b mass regions.

Reconstruction of non-strange B meson decays
Non-strange B mesons, decaying into a D 0 and a small number of pions, were reconstructed in order to verify the B 0 s reconstruction algorithms.The B + → D 0 π + , B 0 d → D * (2010) − π + and B 0 d → D * (2010) − a + 1 decay channels were studied and their mass distributions are shown in Figures 3a and 3b.They were fitted using a Gaussian function for the signal and using an exponential function for the combinatorial background.An additional Gaussian function was used in Figure 3a to account for the signal coming from the following decay channels: 2007) 0 ρ + , where the π 0 and/or γ from the D * (2007) 0 and/or ρ decays were not reconstructed.The main selection criteria which were imposed are rather similar to those used in the B 0 s analysis.The information relevant to these reconstructed channels is given in Table 3.The numbers of observed events are in agreement with expectations.

Channel
B meson mass Characteristics of the non-strange B decay mode and comparison with the simulation.The error on the expected number of events comes from the errors in the branching fractions and reconstruction efficiencies.The branching fraction of the decay mode in the last row is estimated in this paper.The number of observed events in the main peak can be translated into a branching fraction: The quoted error is completely dominated by statistics.The CLEO Collaboration has given an upper limit [25] only on the decay channel with a pion pair: ).This limit is consistent with those which can be extracted from the measured channels in this paper: For this evaluation it was assumed (see Section 3.1) that the decay channel with a + 1 is produced with a rate which is (3.3 ± 0.9) times larger than the channel with π + .

3.5
Study of B 0 s -B 0 s oscillations The combined tagging variable is defined as: The variable x tag varies between -1 and 1.Large and positive values of x tag correspond to a high probability that the produced quark is a b rather than a b in a given hemisphere.
A set of nine discriminant variables has been selected for this analysis.One set (three variables) is determined in the hemisphere which contains the B 0 s meson, the other set (five variables) in the hemisphere opposite to the B 0 s meson, and one variable is common to both hemispheres.Details concerning the definition of these variables and the method are given in reference [10].
An event was classified as a mixed or as an unmixed candidate according to the sign, Q D , of the D s electric charge for decay channels containing a D s meson, or according to the sign of the kaon for D 0 K − π + decay channels, relative to the sign of the x tag variable.Mixed candidates are defined by requiring x tag •Q D < 0, and unmixed candidates correspond to x tag • Q D > 0. The probability, ǫ tag b , for tagging correctly the b or b quark, from the measurement of x tag , was evaluated by using a dedicated simulated event sample and was found to be (74.5±0.5)%.
The corresponding probability for events in the combinatorial background was obtained using real data candidates selected in the side-bands of the D signal: the probabilities to classify these events as mixed or unmixed candidates are called ǫ mix comb and ǫ unmix comb , respectively.For the reflection events the analogous probabilities are called ǫ mix ref , ǫ unmix ref and their values were taken from the simulation.

Proper time resolution
For each event, the B 0 s proper decay time was obtained from the measured decay length (L B 0 s ) and the estimate of the B 0 s momentum (p B 0 s ).The measured position of the D − s or D 0 decay vertex, the momentum, and the corresponding measurement errors, were used to reconstruct a D − s or a D 0 particle.A candidate B 0 s decay vertex was obtained by intercepting the trajectory of this particle with the other charged particle tracks which are supposed to come from the B 0 s decay vertex.For the main peak, the B 0 s momentum is precisely known since all decay products are reconstructed.For the satellite peak, this is no longer true, because there are one or two undetected neutral(s) (π 0 and/or γ) in the B 0 s decay.As discussed in Section 3.2, the satellite peak is assumed to be composed of D The first four decay channels give D − s π + (a + 1 ) and the last two channels give D 0 K − π + (a + 1 ) final states.Kinematic mass constrained fits (imposing the mass of B 0 s and of the intermediate states with corresponding errors) were performed assuming always the presence of a γ coming from the D * − s or D * (2007) 0 meson decays.This is a good approximation also for decay channels containing a ρ + meson.These mass constrained fits were performed on events lying in the mass region corresponding to ±2σ of the fitted mass of the satellite peak.
Except for the combinatorial background contribution, the predicted proper time distributions were obtained by convoluting the theoretical functions with resolution functions evaluated from simulated events.Due to the different decay length resolutions (for the different Vertex Detector configurations), proper time resolutions were used for data samples taken in 1992-1993 and 1994-1995 separately.The proper time resolution, R Bs (t−t i ), was evaluated from the difference between the generated time (t) and the reconstructed time (t i ), fitting this distribution by the sum of two Gaussian functions: where G 1 (t − t i , σ 1 ) and G 2 (t − t i , σ 2 ) are Gaussian functions with the corresponding resolutions σ 1 and σ 2 .In the general case, the σ i resolutions depend on the momentum uncertainty (see equation ( 2)).But in this analysis, as the B 0 s momentum is known with a good precision, the contribution of the momentum uncertainty to the proper time resolution is negligible.f 2 is the fraction of the second Gaussian function, which, by convention, is that with the worse resolution.The values of the corresponding parameters, obtained from simulated events, are given in Table 4.The time resolution for events in the satellite peak is only slightly worse than for those belonging to the main peak.

Main peak Decay channels
Years σ 1 (ps) σ 2 (ps) f 2 all B 0 s channels 1992-1993 0.068 0.18 0.27 all B 0 s channels 1994-1995 0.065 0.12 0.42 Satellite peak The time distribution P comb (t i ) for the combinatorial background under the main peak is obtained from real data by fitting the time distribution of wrong-sign and right-sign events lying in the side-bands of the B 0 s mass distribution.It was verified, in the simulation, that the time distribution for these classes of events is similar to that obtained for events lying under the B 0 s mass peak.The time distribution P comb (t i ) for the combinatorial background under the satellite peak was taken directly from the simulation, since it has a dependence on the measured B 0 s mass, due to the procedure used to reconstruct the B momentum for these events.The time distribution P ref (t i ) for the reflection was also taken from the simulation.

Fitting procedure
The oscillation analysis is performed in the framework of the amplitude method [13], which consists in measuring, for each value of the frequency ∆m B 0 s , an amplitude A and its error σ(A).The parameter A is introduced in the time evolution of pure B 0 s or B 0 s states, so that the value A = 1 corresponds to a genuine signal for oscillation.The time-dependent probability that B 0 s is detected as a B 0 s or B 0 s is then: where the plus (minus) signs refer to B 0 s (B 0 s ) decays.The 95% C.L. excluded region for ∆m B 0 s is obtained by evaluating the probability that, in at most 5% of the cases, a real signal having an amplitude equal to unity would give an observed amplitude smaller than the one measured.This corresponds to the condition: In the amplitude approach, the error σ(A(∆m B 0 s )) is related to the probability to exclude a given value of ∆m B 0 s .The sensitivity is defined as the value of ∆m B 0 s which would just be excluded, if the value of A, as measured in the experiment, was zero for all values of ∆m B 0 s , i.e. it is the expected 95% confidence limit that an experiment of the same statistics and resolution would be expected to achieve (a real experiment would set a higher or lower limit, because of fluctuations in the values of A as a function of ∆m B 0 s ).The probability distributions for mixed and unmixed4 events are: where f Bs , f ref and f comb are the relative fractions of the B 0 s , reflection and combinatorial background events, respectively, which satisfy the condition f Bs +f ref +f comb =1.The expressions for the different probability densities are: s mixing probability: • Reflection mixing probability: • Combinatorial background mixing probability: The parameters ǫ tag b , ǫ mix ref and ǫ mix comb were defined in Section 3.5.1.The amplitude analysis is then performed using all events selected in a mass region between 4.83 and 5.46 GeV/c 2 .The variation of the background level as a function of the reconstructed mass was included in this analysis on an event-by-event basis.Figure 5 shows the variation of the measured amplitude as a function of ∆m B 0 s .With the present level of the statistics, this analysis provides a negligible low limit.On the other hand, its most important feature, despite the low statistics, is the relatively small error on the amplitude at high values of ∆m B 0 s with respect to the more inclusive analyses.Due to the limited statistics, the error on A (σ(A)) can be asymmetric.The error has been symmetrized by taking the larger value.Figure 6a shows the variation of the measured error on the amplitude as a function of ∆m B 0 s .In Figure 6b the results from the D ± s ℓ ∓ [10] and from the D ± s h ∓ (see the second part of the current paper) are compared with the exclusive B 0 s analysis.It should be noted that, due to the better proper time resolution, the resolution σ(A) of the B 0 s exclusive analysis increases more slowly.The ratio of the corresponding errors of B 0 s exclusive and D ± s ℓ ∓ analyses is about 5 (2) at the low (high) values of ∆m B 0 s .The behaviour of σ(A) was investigated using a "toy" Monte Carlo generated with the same characteristics as those measured in real data.The individual toy experiments show a similar behaviour of σ(A), as in the real data, as a function of ∆m B 0 s .For each value of ∆m B 0 s , the distribution of σ(A) and its variance was obtained.The central value and the region corresponding to a ±2σ(A) variation are also shown in Figure 6a.

Study of systematic uncertainties
Systematic uncertainties were evaluated by varying the parameters, which were kept constant in the fit according to their measured or expected errors.
• Systematics from the tagging purity: a variation of ±3% on the expected tagging purity for the signal was considered, following the results given in [10].• Systematics from the background level and kinematic reflection: the levels of background and of kinematic reflection were varied separately for the main peak and for the satellite peak according to the statistical uncertainties given in Table 2. Since the level of the combinatorial background was used as a function of the reconstructed mass on an event-by-event basis, the measured central mass positions and the corresponding widths were also varied by ±1σ around their fitted values.• Systematics from the expected resolution on the B decay proper time: the widths σ 1 and σ 2 described in Section 3.5.2and given in Table 4 were varied by ±10% [10,27,28].• Finally, uncertainties in the values of the B 0 s meson lifetime (τ B 0 s = 1.46 ± 0.06 ps [7]) and in the parameterisation of the proper time of the combinatorial background were taken into account.
In Figure 5, which presents the variation of the measured amplitude as a function of ∆m B 0 s , the shaded area shows the contribution from systematics.

D ± s h ∓ analysis with fully reconstructed D s
Events with an exclusively reconstructed D s accompanied by one or more large momentum hadrons were used to perform a second oscillation analysis.This channel is similar to the D ± s ℓ ∓ final state [10] but, instead of a charged lepton, it uses charged hadrons.It provides a larger number of events but suffers from an ambiguity in the choice of the hadrons and from a reduced B 0 s purity of the selected sample.This approach has already been used in DELPHI to measure the B 0 s lifetime [27].

Event selection
The D s meson is selected in two decay channels: candidates were reconstructed by making combinations of three charged particles located in the same event hemisphere, with a momentum larger than 1 GeV/c and with reconstructed tracks associated to at least one VD hit.The value of | cos ψ| (see Section 3.2) was required to be larger than 0.4 and 0.6 in the D − s → φπ − and D − s → K * 0 K − decay channels, respectively.In order to reduce the combinatorial background from charm and light quarks, the b-tagging probability for the whole event was required to be smaller than 0.5.
In this analysis the neural network algorithm for hadron identification described in Section 2.3 was used.For the D − s → φπ − decay mode, the invariant mass of φ candidates was required to be within ±12 MeV/c 2 of the nominal φ mass and the φ momentum was required to be larger than 5 GeV/c.If the K + K − invariant mass was within ±4 MeV/c 2 of the nominal φ mass, both kaon candidates were required to be identified as "very loose" kaons, otherwise to be identified as "loose" kaons.
For the D − s → K * 0 K − decay mode, the invariant mass of the K * 0 candidates was required to be within ±60 MeV/c 2 of the nominal K * 0 mass value and the K * 0 momentum was required to exceed 5 GeV/c.The momentum of the K − candidate from D − s was required to exceed 3 GeV/c.To suppress the physical background from the D − → K + π − π − kinematic reflection, the K − candidate was required to be identified as a "standard" kaon.For K * 0 → K + π − decays, "loose" identified K + were also accepted.In order to suppress the combinatorial background, the K − candidate was required to be identified as "tight" kaon, if the invariant mass of the K * 0 candidates was out of ±20 MeV/c 2 of the nominal K * 0 mass value, and the value of | cos ψ| was required to be larger than 0.8.
The selection of a hadron accompanying the D s candidate is based on an impact parameter technique.A sample of tracks coming predominantly from B hadron decays was preselected by using their impact parameters and the corresponding errors, both with respect to the primary vertex (Im p , σ p ) and to the secondary D s decay vertex (Im s , σ s ).The hadron was then searched for amongst the preselected particles in the event, by requiring that its charge was opposite to the D s charge and that it had the largest momentum.The efficiency of the hadron selection was about 80% and, among the selected hadrons, (84 ± 4)% came from a B vertex.Details on the track preselection as well as on the hadron selection are given in reference [27].
The B decay vertex was reconstructed by fitting the selected hadron and the D s candidate to a common vertex.The χ 2 -probability of the fitted B 0 s vertex has been required to be larger than 10 −3 .In order to increase the resolution on the measured decay length, only reconstructed events with a decay length error smaller than 0.07 cm were kept.
If the previous procedure failed, a new attempt was made, using an inclusive algorithm which allowed several hadrons to be attached to the D s candidate.This algorithm is based on the difference in the rapidity distributions for particles coming from fragmentation and from B decays.The fragmentation particles, on average, have lower rapidity [29] than B decay products.The charged particles were ordered in increasing values of the rapidity and were attached in turn to the secondary D s vertex.Up to three particles were accepted with their total charge equal to ±1 or 0. The rapidity was calculated as 0.5 log ((E + P L )/(E − P L )), where E is the energy of the particle (assumed to be a pion) and P L its longitudinal momentum with respect to the thrust axis of the event.
Only particles with a momentum greater than 1 GeV/c were accepted.In addition, for the tracks satisfying the condition Im p /σ p < 3, it was required that Im s /σ s < Im p /σ p .Events with a decay length error smaller than 0.07 cm and a χ 2 -probability of the fitted B 0 s vertex larger than 10 −3 were kept.Table 5: Number of D s mesons reconstructed in the φπ − and K * 0 K − decay channels after selection of the accompanying hadron(s).The fraction of combinatorial background, given in parentheses, has been evaluated using a mass interval of ±2σ centred on the fitted D s mass.
The selected sample, specified as D ± s h ∓ in the following, contained about 30% of such "multi-hadron" events.Among them, about 20% have one, 50% two and 30% three selected hadrons.The multi-hadron and single-hadron events were treated in a similar way. Figure 7 shows the D s signals after selection of the accompanying hadron(s).The mass distributions were fitted with two Gaussian functions of equal widths to account for the D − s and D − signals and with an exponential function for the combinatorial background.All parameters were allowed to vary in the fit.Table 5 gives the number of observed events in the D s signal, after background subtraction, and the fraction of combinatorial background.The D s signal region corresponds to a mass interval of ±2σ centred on the fitted D s mass.
The reconstructed number of D − s → φπ − events using 1994-1995 data is about 1.5 larger than those obtained using 1992-1993 data.This factor reflects the difference in statistics between the two data sets.Such "statistical scaling" does not apply to the D − s → K * 0 K − decay mode, where the particle identification, which was better for the 1994-1995 data set (see Section 2), plays a more important role.
Four events are in common with the exclusively reconstructed B 0 s sample: one event in the main peak and three in the satellite peak.These events were removed from the D ± s h ∓ sample for the oscillation analysis.

Sample composition
The D ± s h ∓ sample contained a large physical background due to D s from non-strange B hadron decays and from cc fragmentation.Four sources of events, originating from B decays, were considered: two from B 0 s and two from non-B 0 s mesons, namely, B decaying to one or two charmed mesons comprising at least one D s .The relative fractions of these sources were calculated using five input parameters: [25]; • the probability in the non-strange B meson, that a D − s is produced at the lower vertex: Br(B u,d → D − s X) [32]; • the probability to have two charmed hadrons in a b-decay: Br(b →DDX) [33].
The last two probabilities were assumed to be the same for all B species.To estimate the first two branching fractions , the averaged production rate of D s from all B species, Br(b → D ± s X)×Br(D ± s → φπ ± ), measured by the ALEPH, DELPHI and OPAL collaborations [30], and the equivalent quantity Br(B u,d → D ± s X) × Br(D ± s → φπ ± ), measured at the Υ(4S) by the CLEO and ARGUS collaborations [31], were used.The following fractions for the different B decays contributing to the D ± s h ∓ signal were evaluated: The contribution F cc from direct charm was estimated from the measurement of D s production in charm events at LEP [30], taking into account the Z partial widths into b and c quark pairs: Finally, the relative proportion of the combinatorial background f bkg was taken directly from the fit of the real data mass distributions (see Table 5).

4.3
Discriminant variables to increase the B 0 s purity To increase the effective purity in B 0 s of the selected sample, five variables were used which allow the separation of the signal and background components.These variables are: the D s mass, the D s momentum, the value of | cos ψ|, the χ 2 -probability of the fitted D s decay vertex and the value of the b-tagging variable measured in the hemisphere opposite to that of the D s meson.
The relative components in the selected sample, defined in the previous section, were calculated on an event-by-event basis: F bkg (v i )/tot where v i indicates the i-th discriminant variable, F bc , F b , F c and F bkg are the probability density functions for the b and c together, b, c and the combinatorial background events, respectively.The relative charm contribution is f cc = F cc (1 − f bkg ).In these expressions, the total normalisation factor is: All discriminant variables, except for the b-tagging, were used to separate b and c events together from the combinatorial background (bkg).The b-tagging variable was used to distinguish separately the b from c and from combinatorial background events.The distributions of discriminant variables are shown in Figure 8 for events selected in the signal region, after having subtracted the corresponding distributions of background events obtained using events within the side-bands of the D s signal.The corresponding distributions for background events are also shown in Figure 8.For the b-tagging variable, comparison with simulation is shown for the sum of the b and c events.The agreement between real data and simulated distributions is satisfactory.
The use of this procedure is equivalent to increase the B 0 s purity by 20%.

4.4
Measurement of the B 0 s lifetime

Proper time resolution
For each event, the B 0 s proper decay time was obtained from the measured decay length (L B 0 s ) and the estimate of the B 0 s momentum (p B 0 s ).The technique used is described in [27].The predicted decay time distributions were obtained by convoluting the theoretical distributions with resolution functions evaluated from simulated events.Different parameterisations were used for the two Vertex Detector configurations installed in 1992-1993 and 1994-1995.The proper time resolution was obtained from the (t-t i ) distribution of the difference between the generated (t) and the reconstructed (t i ) proper times.The following distributions were considered: • R B,1D (t − t i ) is the resolution function for B decays with only one charmed meson (D s ) in the final state.It is parameterised as the sum of three Gaussian functions.The width of the second Gaussian is taken to be proportional to the width of the first one.The third Gaussian describes the component with a selected hadron coming from the primary vertex; the fraction (f 3 ) of these events decreases exponentially as a function of the proper time: where t 2 ; f 2 and f 3 are the fractions of the second and third Gaussian functions, respectively; f 3 is defined as The values for the decay length resolutions, σ Lj , the momentum resolutions, σ pj /p, the relative fractions, f j , and the coefficients s j , are given in Table 6.
• R B,2D (t − t i ) is the resolution function for B decays with two D mesons in the final state.In this case, the selected hadron often does not originate directly from the B vertex, but from the second D vertex, hence resulting in a worse resolution.This resolution function is parameterised in a similar way as R B,1D (t − t i ) and the values of the corresponding parameters are shown in Table 6.

Likelihood fit
The B 0 s mean lifetime and the time distribution of the combinatorial background were fitted simultaneously, using selected events lying within a mass interval of ±2σ centred on the measured D s mass (2953 events) and events lying in the D s mass side-band (3373 events) between 2.1 and 2.3 GeV/c 2 .The probability density function for the measured proper time, t i , can be written as: The different probability densities are expressed as convolutions of the physical probability densities with the appropriate resolution functions: • For the signal: where t is the true proper time.• For the physical background coming from all B meson decays: where q runs over the various b-hadron species contributing to this background and f Bq,1D , f Bq,2D are their corresponding fractions.• For the combinatorial background, the following function was used: Three distributions were considered: an exponential for poorly measured events having negative t (with lifetime τ − ), an exponential for the flying background (with lifetime τ + ) and a central Gaussian for the non-flying one.The seven parameters (f − , f + , τ + , τ − , σ − , σ + , σ 0 ) were allowed to vary in the fit.• For "charm" candidates, the function P cc (t i ) has the same form as P bkg (t i ) and has been parameterised using simulated events.In this case all the parameters were fixed in the fit.The B 0 s lifetime fit was performed in the proper time interval between -4 ps and 12 ps and the result of the fit, shown in Figure 9, is: Table 7: Sources of systematic uncertainties on the B 0 s lifetime.
The systematic error due to the uncertainties in the relative fractions of the different D s sources corresponds to a ±1σ variation of the fractions f used in the likelihood fit, excluding f bkg which is studied separately.
The estimate of systematics related to the evaluation of the B 0 s purity, on an eventby-event basis, was obtained in the following way.The distributions of the different quantities contributing to the discriminant variable (Figure 8) for signal and background events were re-weighted with a linear function in order to maximise the agreement between data and the simulation.The B s lifetime distribution was refitted with new probability distributions and the difference between the corresponding values of the fitted lifetime taken as a systematic error.
Uncertainties on the determination of the resolution of the proper time receive two contributions: one from errors on the decay distance evaluation and the other from errors on the measurement of the B 0 s momentum.The systematic error coming from uncertainties on the time resolution was evaluated by varying the widths σ L and σ p of the resolution function by ±10% [10,27,28].Finally, simulated B 0 s events, generated with a lifetime of 1.6 ps and satisfying the same selection criteria as the real data, have a fitted lifetime of 1.64 ± 0.04 ps.The B 0 s lifetime value obtained in the Section 4.4.2 was not corrected but the statistical error of this comparison (±0.04 ps) was included in the systematic error.The measured B 0 s lifetime was found to be: τ B 0 s = 1.53 +0.16 −0.15 (stat.)± 0.07(syst.)ps.A similar analysis by the ALEPH collaboration [12] gave a consistent result: τ B 0 s = 1.47 ± 0.14(stat.)± 0.08(syst.)ps.

Lifetime difference between B 0 s mass eigenstates
The B 0 s (or B 0 s ) mesons are superpositions of the two mass eigenstates: Neglecting CP violation, the time probability density is then given by: where The first term of equation (20) refers to final states of identified beauty flavor, as in the D s ℓ case [10], it does include D s +light mesons final states.The second term corresponds to D ( * )+ s D ( * )− s final states, which are dominantly (98%) CP even eigenstates [34].The value of x cp was taken as the ratio of F Bs,2D to (F Bs,1D +F Bs,2D ): x cp = 0.  [14]).The log-likelihood function described in Section 4.4.2 was modified by replacing the probability density for B 0 s by ( 20) convoluted with the appropriate resolution functions.It was minimized in the (τ B 0 s , ∆Γ B 0 s /Γ B 0 s ) plane and the likelihood difference with respect to its minimum ∆L (Figure 10-a) was computed: ).
The upper limit on ∆Γ B 0 s /Γ B 0 s , calculated from P(∆Γ B 0 s /Γ B 0 s ), is: ∆Γ B 0 s /Γ B 0 s < 0.67 at the 95% C.L.This limit takes into account both statistical uncertainties and the systematic coming from the uncertainty 5  where x i sys are the n parameters considered in the systematic uncertainty and f (x i sys ) are the corresponding probability densities.
Only three systematics were considered here: the relative fraction of combinatorial background f bkg , the B 0 s purity of the selected sample and the x cp fraction.Other systematic uncertainties are expected to be small as they are in the lifetime measurement.
The ∆Γ B 0 s /Γ B 0 s probability distribution, obtained with the inclusion of the systematics, is shown in Figure 10-c, the most probable value for ∆Γ B 0 s /Γ B 0 s is 0.35 and the upper limit at the 95% confidence level is: ∆Γ B 0 s /Γ B 0 s < 0.69 at the 95% C.L. It should be noted that the world average of the B 0 s lifetime cannot be used as a constraint in such an analysis, since it depends on ∆Γ B 0 s and on Γ B 0 s .Moreover, this dependence is also different for different decay channels.In the D ± s h ∓ case the expression of the average B 0 s lifetime is given by: 4.6 Study of B 0 s -B 0 s oscillations 4.6.1 Tagging procedure The tagging algorithm was explained in Section 3.5.1.An event is classified as a mixed or an unmixed candidate according to the relative signs of the D s electric charge, Q D , and of the tagging purity variable, x tag .
Mixed candidates have x tag • Q D < 0, and unmixed ones x tag • Q D > 0. The probability, ǫ tag b , of tagging the b or the b quark correctly from the measurement of x tag was evaluated using a dedicated simulated event sample.The average tagging purity of the x tag variable, given by the simulation for true B 0 s → D − s h + X decays, is (71.4±0.4)%.The purity is lower than that obtained in the D s ℓ sample [10] because not all B s charged decay products are reconstructed in the present analysis.It was verified that the tagging purity is the same for different B hadron species and varies by less than about ±2% whether the B 0 s has oscillated or not.This effect is taken into account in the systematics.The corresponding probability distribution for events in the combinatorial background was obtained using data candidates selected in the side-bands of the D s signal: the probabilities of classifying these events as mixed or as unmixed candidates are called ǫ mix bkg and ǫ unmix bkg , respectively.For the charm events, the analogous probabilities are called ǫ mix cc , ǫ unmix cc and their values were taken from the simulation.

Fitting procedure
From the expected proper time distributions and the tagging probabilities, the probability functions for mixed and unmixed event candidates are +f ef f cc P mix cc (t i ) + f ef f bkg P mix bkg (t i ), where t i is the reconstructed proper time.The analytical expressions for the different probability densities are given in the following, with t being the true proper time: • B 0 s signal mixing probability: • Physical background mixing probabilities: • Mixing probability for charm component: • Combinatorial background mixing probability: The oscillation analysis was performed in the framework of the amplitude method [13] as described in Section 3.5.3.Considering only statistical uncertainties, the limit is: with a corresponding sensitivity equal to 3.1 ps −1 .At ∆m B 0 s = 10 ps −1 , the error on the amplitude is 2.3 (see Figure 11).

Study of systematic uncertainties
Systematic uncertainties were evaluated by varying the parameters which were kept constant in the fit, according to their measured or expected errors.
• Systematics from the tagging probability: a conservative variation of ±3% on the expected tagging probability for the signal and for the other three processes contributing to the D ± s h ∓ candidates was used.The same variation is assumed for the tagging purity for the charm and combinatorial background events.The central values of these purities were fixed to the simulated ones.
• Systematics from the B 0 s purity: the same procedure already applied to the lifetime measurement has been used.
• Systematics from the resolution on the B decay proper time: the same procedure already applied to the lifetime measurement was used.In addition, the systematic error due to the variation of the proper time distribution of the combinatorial background was considered: the parameters used to define the background shape in the lifetime fit were varied according to their fitted errors.
The inclusion of systematic uncertainties lowers the sensitivity to 2.7 ps −1 and the 95% C.L limit becomes ∆m B 0 s > 4.1 ps −1 .The analogous analysis has been performed by the ALEPH collaboration [12], which set a limit at ∆m B 0 s > 3.9 ps −1 at the 95% C.L., with a better sensitivity of 4.1 ps −1 .

5
Combined limit on ∆m B 0 s using exclusive B 0 s and D ± s h ∓ events This paper presents two analyses on ∆m B 0 s using exclusively reconstructed B 0 s mesons and D ± s h ∓ events.Their results have been combined (Figure 12), taking into account correlations between systematic uncertainties in the two amplitude measurements [7].A limit at the 95% confidence level is obtained: with a corresponding sensitivity equal to 3.2 ps −1 (with statistical errors only, the limit would be ∆m B 0 s > 4.0 ps −1 at the 95% C.L. with a sensitivity of 4.4 ps −1 ). Figure 12 also shows the error on the amplitude for different values of ∆m B 0 s .

6
Combination of the DELPHI ∆m B 0 s , τ B 0 s and ∆Γ B 0 s /Γ B 0 s analyses DELPHI has performed three analyses on ∆m B 0 s using D ± s ℓ ∓ candidates [10], D ± s h ∓ events and exclusively reconstructed B 0 s mesons.They were combined, taking into account correlations between the event samples and between systematic uncertainties in the different amplitude measurements (Figure 13a).This gives the following limit for ∆m B 0 s : ∆m B 0 s > 4.9 ps −1 at the 95% C.L. with a sensitivity of ∆m B 0 s = 8.7 ps −1 The exclusion probability for this limit is 88%.The variation, with ∆m B 0 s of the error on the amplitude is given in Figure 13b.
The results of two analyses on τ B 0 s and ∆Γ B 0 s /Γ B 0 s using D ± s ℓ ∓ and D ± s h ∓ events provide the following DELPHI results:

Conclusion
Using about 3.5 million hadronic Z decays registered by DELPHI between 1992 and 1995, two samples of events have been selected.The first one consists of 44 reconstructed B 0 s events: 11 candidates (including 30% of background) are completely reconstructed and 33 candidates (including 55% of background) are partially reconstructed (π 0 and/or γ are not detected).This analysis used twelve different decay channels of the B 0 s meson and is a first attempt to use such events for the oscillation studies.Due to the excellent proper time resolution, this sample gives some contribution in the high ∆m B 0 s region.The second sample contains 2953 D ± s h ∓ candidates (including 60% of background) with completely reconstructed D − s mesons in the φπ − and K * 0 K − decay channels.Using the D ± s h ∓ sample, three studies have been performed.The B 0 s lifetime has been measured and a limit on the fractional width difference between the two physical B 0 s states has been obtained:         The distributions on the left are for events selected in the signal region after having subtracted the corresponding distributions of background events, which have been obtained using events situated in the side-bands of the D s signal.The corresponding distributions for background events are shown on the right.The points with error bars correspond to the data and the histograms are simulated events.For the | cos ψ| distribution, only the φπ + decay mode is shown because the cut on this variable was set at 0.6 for K * 0 K − and at 0.4 for φπ + channel.In addition, the D s mass was used as the fifth discriminating variable: the signal has a Gaussian and the background has an exponential distribution.

FinallyFigure 4
Figure 4 shows the mass distribution for these two decay channels.It was fitted with an exponential function for the background and two Gaussian functions for the signals corresponding to the main and satellite peaks, with respective widths equal to about 34 MeV/c 2 and 60 MeV/c 2 , as obtained from the simulation.The position of the wider Gaussian, corresponding to B 0 d → D * (2007) 0 π − π + and B 0 d → D * (2007) 0 π − a + 1 signals with unreconstructed π 0 or γ, and the slope of the background exponential function were also taken from the simulation.The fit to the mass distribution yielded a signal of 8 ± 4 B 0 d decays in the main peak.The probability that the background has fluctuated to give the observed number of events is 3 × 10 −2 .The number of observed events in the main peak can be translated into a branching fraction: Br(B 0 d → D 0 π − π + + D 0 π − a + 1 ) = (3.6 ± 1.9) × 10 −3 .The quoted error is completely dominated by statistics.The CLEO Collaboration has given an upper limit[25] only on the decay channel with a pion pair: Br(B 0 d → D 0 π − π + ) < 1.6 × 10 −3 (at the

3. 5 . 1
Algorithm for tagging b or b quark at production time The signature of the initial production of a b (b) quark in the jet containing the B 0 s or B 0 s candidate was determined using a combination of different variables sensitive to the initial quark state.For each individual variable X i , the probability density functions f b (X i ) (f b (X i )) for b (b) quarks were built and the ratio

2 )Figure 1 :Figure 2 :
Figure 1: B 0s mass spectrum for the candidates selected in the twelve decay channels described in Section 3.2.The data are indicated by the points with error bars and the result of the fit has been superimposed.Details on this fit are given in the text.The histograms represent the expected contribution from beauty events (after having removed the exclusively reconstructed B 0 s decay channels), from charm events and from light quark events.

2 )Figure 4 :
Figure 4: B 0 d mass spectrum for the sum of B 0 d → D 0 π − π + and B 0 d → D 0 π − a + 1 decays selected in the four decay channels described in Section 3.4.The data are indicated by the points with error bars and the result of the fit has been superimposed.The histograms represent the expected contribution from beauty events (after having removed the exclusively reconstructed B 0 d decay channels), from charm events and from light quark events.The widths of the signals have been fixed according to the values found in the simulation.Details on the fit are given in the text.

Figure 5 :
Figure 5: Exclusive B 0 s analysis: variation of the oscillation amplitude A as a function of ∆m B 0 s .The lower continuous line corresponds to A + 1.645 σ(A) where σ(A) includes statistical uncertainties only, while the shaded area shows the contribution from systematics.The dashed-dotted line corresponds to the sensitivity curve.The lines at A=0 and A=1 are also given.The points with error bars are real data.

Figure 6 :
Figure 6: σ(A) as a function of ∆m B 0 s .a) Exclusive B 0 s analysis: the full curve shows the result from the data.The dashed curve shows the average result from 100 toy experiments with the same statistics as in data and the shaded area gives the ±2σ(A) region around this average.The systematic effects are not included.b) Comparison of the σ(A) as function of ∆m B 0 s for three analyses: the full curve shows the result from the exclusive B 0 s analysis, the dashed and dotted curves show the result from D ± s h ∓ and D ± s ℓ ∓ analyses, respectively.

Figure 8 :
Figure8: Distributions of the variables used to increase the B s purity.The distributions on the left are for events selected in the signal region after having subtracted the corresponding distributions of background events, which have been obtained using events situated in the side-bands of the D s signal.The corresponding distributions for background events are shown on the right.The points with error bars correspond to the data and the histograms are simulated events.For the | cos ψ| distribution, only the φπ + decay mode is shown because the cut on this variable was set at 0.6 for K * 0 K − and at 0.4 for φπ + channel.In addition, the D s mass was used as the fifth discriminating variable: the signal has a Gaussian and the background has an exponential distribution.

Figure 10 :
Figure 10: D ± s h ∓ analysis: 68%, 95%, 99% C.L. contours of the negative log-likelihood in the plane ∆Γ B 0 s /Γ B 0 s − τ B 0 s a) without and b) with τ B 0 s = τ B 0 d constraint.The point indicates the minimum.c) Probability density distribution for ∆Γ Bs /Γ Bs ; the two lightly shaded regions at the 68% C.L. and dark one at the 95% C.L. are also shown.

Figure 11 :
Figure 11: D ± s h ∓ analysis.Left plot: variation of the oscillation amplitude A as a function of ∆m B 0 s .The lower continuous line corresponds to A + 1.645 σ(A) where σ(A) includes statistical uncertainties only, while the shaded area shows the contribution from systematics.The dashed-dotted line corresponds to the sensitivity curve.The lines at A=0 and A=1 are also given.The points with error bars are real data.Right plot: variation of the error on the amplitude as a function of ∆m B 0 s , including systematic uncertainties.

Figure 12 :
Figure 12: Combination of the D ± s h ∓ and exclusive B 0 s analyses.Left plot: variation of the oscillation amplitude A as a function of ∆m B 0 s .The lower continuous line corresponds to A + 1.645 σ(A) where σ(A) includes statistical uncertainties only, while the shaded area shows the contribution from systematics.The dashed-dotted line corresponds to the sensitivity curve.The lines at A=0 and A=1 are also given.The points with error bars are real data.Right plot: variation of the error on the amplitude as a function of ∆m B 0 s including systematic uncertainties.It should be noted that the D ± s h ∓ analysis dominates at low values of ∆m B 0 s and the exclusive B 0 s analysis dominates at large values of ∆m B 0 s .

Figure 13 :
Figure 13: Combined DELPHI analysis: a) variation of the oscillation amplitude A as a function of ∆m B 0 s .The lower continuous line corresponds to A + 1.645 σ(A) where σ(A) includes statistical uncertainties only, while the shaded area shows the contribution from systematics.The dashed-dotted line corresponds to the sensitivity curve.The lines at A=0 and A=1 are also given.The points with error bars are real data.b) variation of the error on the amplitude as a function of ∆m B 0 s .

Table 1 :
Values of the B 0 s branching fractions used in the current analysis.

Table 2 :
Characteristics of the B 0 s signals and comparison with the simulation.The expected numbers of events were calculated taking into account the different branching fractions as given in Table [25]expected mass difference of 708.9 ± 1.8 MeV/c 2 is in good agreement with the observed one, taking into account the full width of this state which is 15 +5 −4 MeV/c 2[25].

Table 4 :
Proper time resolution from simulation for 1992-1993 and 1994-1995 data sets and for the different decay channels.It has been parameterised by the sum of two Gaussian functions.

Table 6 :
Fitted values of the parameters used in the resolution functions R B,1D (t − t i ) and R B,2D (t − t i ).

Table 7 .
16−0.15 (stat.)ps.4.4.3Systematic uncertainties on the B 0 s lifetime The contributions to systematic uncertainties on the B 0 s lifetime measurement are summarised in ) .
The uncertainty due to the theoretical prediction of the equality of the τ B 0 s and τ B 0 d lifetimes is negligible with respect to the present error on τ B 0 d .