The $B^0_s$ and $B^0$ mixing frequencies, $\Delta m_s$ and $\Delta m_d$, are measured using a data sample corresponding to an integrated luminosity of 1.0 fb^{1} collected by the LHCb experiment in $pp$ collisions at a centre of mass energy of 7 TeV during 2011. Around 1.8x10^6 candidate events are selected of the type $B^0_{(s)} \to D^_{(s)} \mu^+$ (+ anything), where about half are from peaking and combinatorial backgrounds. To determine the B decay times, a correction is required for the momentum carried by missing particles, which is performed using a simulationbased statistical method. Associated production of muons or mesons allows us to tag the initialstate flavour and so to resolve oscillations due to mixing. We obtain \Delta m_s = (17.93 \pm 0.22 (stat) \pm 0.15 (syst)) ps^{1}, \Delta m_d = (0.503 \pm 0.011 (stat) \pm 0.013 (syst)) ps^{1}. The hypothesis of no oscillations is rejected by the equivalent of 5.8 standard deviations for $B^0_s$ and 13.0 standard deviations for $B^0$. This is the first observation of $B^0_s$ mixing to be made using only semileptonic decays.
Mass distributions for all selected signal candidates. Left, the $K^+K^\pi^+$ invariant mass, where the known mass of the $D^+_s$ has been subtracted. Right, the $D\mu$ normalized mass as defined in Eq. 2. Neutral candidates are those of the form $D^{\mp}\mu^\pm$, while doublecharged candidates are those of the form $D^{\pm}\mu^\pm$. The doublecharged candidates arise from several background sources, most of which are also present in the neutral sample. In the left plot, the neutral sample exhibits much larger $D$ mass peaks, indicative of the large $B$ signal component. 
Fig1.pdf [65 KiB] HiDef png [172 KiB] Thumbnail [145 KiB] *.C file 

Input to obtain the $k$factor correction from the fullysimulated $B^0_s$ sample. For each event the ratio of reconstructed to generated momentum, $p_{\textrm{rec}}/p_{\textrm{sim}}$ is plotted against the normalized $D\mu$ mass ($n$ in Eq. 2). The curve shows a fourthorder polynomial resulting from a fit to the mean of the distribution (in bins of $n$). 
Fig2.pdf [66 KiB] HiDef png [1 MiB] Thumbnail [635 KiB] *.C file 

Illustration of the decay time resolution obtained from a fully simulated $B^0$ signal sample. The left plots demonstrate the Gaussian fits (solid lines) using the full LHCb simulated data (filled), to determine the decay time resolution. Each measured (reconstructed and corrected) time, $t'$, is compared to the corresponding simulated decay time, $t$. The results are shown for several bins of $t'$. The dependence on decay time of the mean (bias, $\mu$) and width (standard deviation, $\sigma$) can be fitted with a quadratic or cubic function of either $t$ or $t'$. The right hand plot shows a quadratic fit to the widths. 
Fig3.pdf [74 KiB] HiDef png [340 KiB] Thumbnail [227 KiB] *.C file 

Distribution of measured $K^+K^\pi^+$ mass, where the known mass of the $D^+_s$ has been subtracted. Black points show the data, and the various lines overlay the result of the fit. The small step at $50$ MeV$c^{2}$ is the result of differences in tagging efficiency for the $B^0_s$ and $B^0$ hypotheses. 
Fig4.pdf [121 KiB] HiDef png [362 KiB] Thumbnail [292 KiB] *.C file 

Measured $B$ decaytime distribution, overlaid with projections of the fit, for backgroundonly regions. Top left: a region between the two signal peaks, $80$ to $20$ MeV$c^{2}$ (with respect to the known mass of the $D_s^+$), showing only low decay times. Top right: a region to the right of the signal peaks $20$ to $100$ MeV$c^{2}$, showing only low decay times. Bottom row: the same on an extended decaytime scale and logarithmic. The legend is the same as in Fig. 4. 
Fig5a[..].pdf [100 KiB] HiDef png [394 KiB] Thumbnail [348 KiB] *.C file 

Fig5b[..].pdf [99 KiB] HiDef png [397 KiB] Thumbnail [350 KiB] *.C file 

Fig5c[..].pdf [132 KiB] HiDef png [432 KiB] Thumbnail [342 KiB] *.C file 

Fig5d[..].pdf [128 KiB] HiDef png [394 KiB] Thumbnail [325 KiB] *.C file 

Measured $B$ decaytime distribution, overlaid with projections of the fit, for signal regions. Top left: for oddtags, high$n$ and a region of $\pm 20$ MeV$c^{2}$ around the $D^+_s$ mass peak, showing only low decay times, where $B^0_s$ oscillations can be clearly seen. Top right: for oddtags and all $n$ for a region of $\pm 20$ MeV$c^{2}$ around the $D^+$ mass peak, showing only low decay times. Bottom row: for both tags and all $n$ for regions of $\pm 20$ MeV$c^{2}$ around the $D^+_s$ (left) and $D^+$ (right) mass peaks. The legend is the same as in Fig. 4. 
Fig6a[..].pdf [98 KiB] HiDef png [445 KiB] Thumbnail [373 KiB] *.C file 

Fig6b[..].pdf [90 KiB] HiDef png [460 KiB] Thumbnail [383 KiB] *.C file 

Fig6c[..].pdf [122 KiB] HiDef png [360 KiB] Thumbnail [294 KiB] *.C file 

Fig6d[..].pdf [126 KiB] HiDef png [431 KiB] Thumbnail [340 KiB] *.C file 

Tagged (mixing) asymmetry, $(N_+N_)/(N_++N_)$, as a function of $B$ decay time. The left plot shows the asymmetry for events for a region of $\pm 20$ MeV$c^{2}$ around the $D^+_s$ mass peak, and the right plot shows the corresponding asymmetry around the $D^+$ mass peak. The black points show the data and the curves are projections of the fitted PDF. On the left plot the fast oscillations of $B^0_s$ are gradually washed out by the increasingly poor decaytime resolution. 
Fig7aleft.pdf [77 KiB] HiDef png [310 KiB] Thumbnail [293 KiB] *.C file 

Fig7b[..].pdf [65 KiB] HiDef png [226 KiB] Thumbnail [205 KiB] *.C file 

Result of using Fourier transforms to search for the $\Delta m_s$peak. The image on the left is constructed from bins of the $K^+K^\pi^+$ mass which are 25 MeV$c^{2}$ in width, analysed in steps of 5 MeV$c^{2}$ such that a smooth image is produced. The colour scale (bluegreenyellowred) is an arbitrary linear representation of the signal intensity; dark blue is used for zero and below. The vertical dashed line is drawn at $18.0$ ps$^{1}$. The apparent doublepeak structure is an artifact of this image. On the right a slice around the $D^+_s$ mass region shows only the peak as used to measure the central value and rms width. 
Fig8aleft.pdf [214 KiB] HiDef png [2 MiB] Thumbnail [651 KiB] *.C file 

Fig8b[..].pdf [26 KiB] HiDef png [73 KiB] Thumbnail [34 KiB] *.C file 

Animated gif made out of all figures. 
PAPER2013036.gif Thumbnail 
A selection of fitted parameter values, for which statistical uncertainties only are given. The $B^0_s$ signal fraction includes contributions from any detached $D^+_s$ production. When the omitted fractions (of combinatorial background components) are included, the total fraction sums to unity within each $n$ region separately. 
Table_1.pdf [46 KiB] HiDef png [144 KiB] Thumbnail [68 KiB] tex code 

A selection of fitted parameter values, for which statistical uncertainties only are given. The $B^0_s$ signal fraction includes contributions from any detached $D^+_s$ production. When the omitted fractions (of combinatorial background components) are included, the total fraction sums to unity within each $n$ region separately. 
Table_1.pdf [46 KiB] HiDef png [144 KiB] Thumbnail [68 KiB] tex code 

Sources of systematic uncertainty on $\Delta m_s$ and $\Delta m_d$. "Simulation" implies a combination of full LHCb simulation and pseudoexperiment studies. 
Table_2.pdf [54 KiB] HiDef png [94 KiB] Thumbnail [44 KiB] tex code 

A selection of fitted parameter values, for which statistical uncertainties only are given. The $B^0_s$ signal fraction includes contributions from any detached $D^+_s$ production. When the omitted fractions (of combinatorial background components) are included, the total fraction sums to unity within each $n$ region separately. 
Table_1.pdf [46 KiB] HiDef png [144 KiB] Thumbnail [68 KiB] tex code 

Sources of systematic uncertainty on $\Delta m_s$ and $\Delta m_d$. "Simulation" implies a combination of full LHCb simulation and pseudoexperiment studies. 
Table_2.pdf [54 KiB] HiDef png [94 KiB] Thumbnail [44 KiB] tex code 
Created on 17 August 2019.Citation count from INSPIRE on 23 August 2019.