The decaytimedependent $CP$ asymmetry in $B^{0}_{s}\to J/\psi K^{+} K^{}$ decays is measured using protonproton collision data, corresponding to an integrated luminosity of $1.9 \mathrm{fb^{1}}$, collected with the LHCb detector at a centreofmass energy of $13 \mathrm{TeV}$ in 2015 and 2016. Using a sample of approximately 117 000 signal decays with an invariant $K^{+} K^{}$ mass in the vicinity of the $\phi(1020)$ resonance, the $CP$violating phase $\phi_s$ is measured, along with the difference in decay widths of the light and heavy mass eigenstates of the $B^{0}_{s}$$\bar{B}^{0}_{s}$ system, $\Delta\Gamma_s$. The difference of the average $B^{0}_{s}$ and $B^{0}$ meson decay widths, $\Gamma_s\Gamma_d$, is determined using in addition a sample of $B^{0} \to J/\psi K^{+} \pi^{}$ decays. The values obtained are $\phi_s = 0.083\pm0.041\pm0.006 \mathrm{rad}$, $\Delta\Gamma_s = 0.077 \pm 0.008 \pm 0.003 \mathrm{ps^{1}}$ and $\Gamma_s\Gamma_d = 0.0041 \pm 0.0024 \pm 0.0015 \mathrm{ps^{1}}$, where the first uncertainty is statistical and the second systematic. These are the most precise single measurements of these quantities to date and are consistent with expectations based on the Standard Model and with a previous LHCb analysis of this decay using data recorded at centreofmass energies 7 and 8 TeV. Finally, the results are combined with recent results from $B^{0}_{s}\to J/\psi \pi^{+} \pi^{}$ decays obtained using the same dataset as this analysis, and with previous independent LHCb results.
Distribution of the invariant mass of $ B ^0_ s $ candidates, selected from simulated $ B ^0_ s \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K ^+ K ^ $ (green filled area), $\Lambda ^0_ b \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} p K ^ $ (solid red line) and $ B ^0 \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K ^+ \pi ^ $ (dotted blue line) decays. The distributions are weighted to correct differences in the kinematics and the resonance content between simulation and data. 
Fig1.pdf [15 KiB] HiDef png [247 KiB] Thumbnail [205 KiB] *.C file 

(a) Distribution of the invariant mass of selected $ B ^0_ s \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K ^+ K ^ $ decays. The signal component is shown by the longdashed red line, the background component by the dashed green line and the total fit function by the solid blue line. The background contribution due to $\Lambda ^0_ b \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} p K ^ $ decays is statistically subtracted. The contribution from $ B ^0 \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K ^+ K ^ $decays is not shown separately due to its small size. (b) Distribution of $ K ^+ K ^ $ invariant mass from selected $ B ^0_ s \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K ^+ K ^ $ decays. The background is subtracted using the $sPlot$ method. The dashed blue lines define the boundaries of the six $m(K^{+}K^{})$ bins that are used in the analysis. 
Fig2a.pdf [45 KiB] HiDef png [225 KiB] Thumbnail [192 KiB] *.C file 

Fig2b.pdf [24 KiB] HiDef png [167 KiB] Thumbnail [162 KiB] *.C file 

(a) Decaytime distribution of the prompt $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K ^+ K ^ $ calibration sample with the result of an unbinned maximumlikelihood fit overlaid in blue. The overall tripleGaussian resolution is represented by the dashed red line, while the two longlived and the wrongPV components are shown by the longdasheddotted and dashedmultipledotted brown and pink lines and the longdashed purple line, respectively. (b) Variation of the effective singleGaussian decaytime resolution, $\sigma_{\rm eff}$, as a function of the estimated percandidate decaytime uncertainty, $\delta_t$, obtained from the prompt $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K ^+ K ^ $ sample. The red line shows the result of a linear fit. The data points are positioned at the barycentre of each $\delta_t$ bin. The shaded histogram (see right $y$ axis) shows the distribution of $\delta_t$ in the backgroundsubtracted $ B ^0_ s \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K ^+ K ^ $ sample. 
Fig3a.pdf [50 KiB] HiDef png [245 KiB] Thumbnail [183 KiB] *.C file 

Fig3b.pdf [14 KiB] HiDef png [348 KiB] Thumbnail [212 KiB] *.C file 

Distribution of the invariant mass of selected (a) $ B ^0 \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K^+ \pi^$ and (b) $B^+ \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K^+$ decays used for the calibration and validation of the decaytime efficiency. The signal component is shown by the longdashed red line, the background component by the dashed green line and the total fit function by the solid blue line. 
Fig4a.pdf [43 KiB] HiDef png [198 KiB] Thumbnail [154 KiB] *.C file 

Fig4b.pdf [43 KiB] HiDef png [202 KiB] Thumbnail [161 KiB] *.C file 

Decaytime efficiency for the (a) 2015 unbiased, (b) 2015 biased, (c) 2016 unbiased and (d) 2016 biased $ B ^0_ s \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} \phi$ sample. The cubicspline function described in the text is shown by the blue line. For comparison, the black points show the efficiency when computed using histograms for each of the input component efficiencies. 
Fig5a.pdf [23 KiB] HiDef png [123 KiB] Thumbnail [116 KiB] *.C file 

Fig5b.pdf [23 KiB] HiDef png [122 KiB] Thumbnail [113 KiB] *.C file 

Fig5c.pdf [23 KiB] HiDef png [106 KiB] Thumbnail [95 KiB] *.C file 

Fig5d.pdf [22 KiB] HiDef png [106 KiB] Thumbnail [94 KiB] *.C file 

Normalised angular efficiency as a function of (a) $\cos\theta_K$, (b) $\cos\theta_\mu$ and (c) $\phi_h$, where in all cases the efficiency is integrated over the other two angles. The efficiency is evaluated using simulated $ B ^0_ s \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} \phi$ decays that have been weighted to match the kinematics and physics of $ B ^0_ s \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K ^+ K ^ $ decays in data, as described in the text. The points are obtained by dividing the angular distribution in the simulated sample by the distribution expected without any efficiency effect and the curves represent an even fourthorder polynomial parameterisation of each onedimensional efficiency. The figure is for illustration only as the angular efficiency is accounted for by normalisation weights in the signal PDF. 
Fig6.pdf [23 KiB] HiDef png [124 KiB] Thumbnail [72 KiB] *.C file 

Calibration of the OS tagger using $ B ^+ \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K ^+ $ decays. The black points show the average measured mistag probability, $\omega$, in bins of predicted mistag, $\eta$, the red line shows the calibration as described in the text and the yellow area the calibration uncertainty within one standard deviation. The shaded histogram shows the distribution, with arbitrary normalisation, of $\eta$ in the background subtracted $ B ^0_ s \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} \phi$ sample, summing over candidates tagged as $ B ^0_ s $ or $\overline{ B }{} {}^0_ s $. 
Fig7.pdf [26 KiB] HiDef png [218 KiB] Thumbnail [181 KiB] *.C file 

Distribution of the invariant mass of selected $ B ^0_ s \rightarrow D ^_ s \pi^+$ candidates (black points). The total fit function is shown as the solid blue line. The signal component is shown by the red longdashed line, the combinatorial background by the lightblue shortdashed line and other small background components are also shown as specified in the legend. Only the dominant backgrounds are shown. 
Fig8.pdf [41 KiB] HiDef png [240 KiB] Thumbnail [185 KiB] *.C file 

Variation of the effective singleGaussian decaytime resolution, $\sigma_{\rm eff}$, as a function of the estimated perevent decaytime uncertainty, $\delta_t$, obtained from the prompt $ D ^_ s \pi^+$ sample. The red line shows the result of a linear fit to the data and the yellow band its uncertainty within one standard deviation. 
Fig9.pdf [18 KiB] HiDef png [175 KiB] Thumbnail [117 KiB] *.C file 

(a) Distribution of the decay time for $ B ^0_ s \rightarrow D ^_ s \pi^+$ candidates tagged as mixed and unmixed with the projection of the fit result, which is described in the text. (b) Calibration of the SSK tagger using $ B ^0_ s \rightarrow D ^_ s \pi^{+}$ decays. The black points show the average measured mistag probability, $\omega$, in bins of predicted mistag, $\eta$, the red line shows the calibration obtained from the fit described in the text, and the yellow area the calibration uncertainty within one standard deviation. The shaded histogram shows the distribution of $\eta$ in the background subtracted $ B ^0_ s \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} \phi$ sample. 
Fig10a.pdf [24 KiB] HiDef png [257 KiB] Thumbnail [170 KiB] *.C file 

Fig10b.pdf [17 KiB] HiDef png [208 KiB] Thumbnail [160 KiB] *.C file 

Decaytime and helicityangle distributions for background subtracted $ B ^0_ s \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K ^+ K ^ $ decays (data points) with the onedimensional projections of the PDF at the maximumlikelihood point. The solid blue line shows the total signal contribution, which contains (longdashed red) $ C P$ even, (shortdashed green) $ C P$ odd and (dotteddashed purple) Swave contributions. Data and fit projections for the different samples considered (datataking year, trigger and tagging categories, $m(K^{+}K^{})$ bins) are combined. 
Fig11a.pdf [21 KiB] HiDef png [207 KiB] Thumbnail [170 KiB] *.C file 

Fig11b.pdf [19 KiB] HiDef png [224 KiB] Thumbnail [170 KiB] *.C file 

Fig11c.pdf [19 KiB] HiDef png [231 KiB] Thumbnail [181 KiB] *.C file 

Fig11d.pdf [21 KiB] HiDef png [221 KiB] Thumbnail [170 KiB] *.C file 

Regions of 68% confidence level in the $\phi_s$$\Delta\Gamma_s$ plane for the individual LHCb measurements and a combined contour (in blue). The $ B ^0_ s \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K ^+ K ^ $ (magenta) and $ B ^0_ s \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} \pi ^+ \pi ^ $ \cite{LHCbPAPER2019003} (red) contours show the Run 1 and Run 2 combined numbers. The $\phi_s$ \cite{CKMfitter2015} and $\Delta\Gamma_{s}$ \cite{Artuso:2015swg} predictions are indicated by the thin black rectangle. 
Fig12.pdf [103 KiB] HiDef png [398 KiB] Thumbnail [193 KiB] *.C file 

Animated gif made out of all figures. 
PAPER2019013.gif Thumbnail 
Calibration parameters for the OS and SSK taggers. Where given, the first uncertainty is statistical and the second is systematic. 
Table_1.pdf [36 KiB] HiDef png [62 KiB] Thumbnail [32 KiB] tex code 

Overall tagging performance for $ B ^0_ s \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K^{+}K^{}$. The uncertainty on $\epsilon_{\rm tag}D^{2}$ is obtained by varying the tagging calibration parameters within their statistical and systematic uncertainties summed in quadrature. 
Table_2.pdf [37 KiB] HiDef png [92 KiB] Thumbnail [43 KiB] tex code 

Correlation matrix including the statistical and systematic correlations between the parameters. 
Table_3.pdf [49 KiB] HiDef png [62 KiB] Thumbnail [30 KiB] tex code 

Correlation matrix for the results in Eq. \eqref{eq:comb1} taking into account correlated systematics between Run 1 and the 2015 and 2016 results. 
Table_4.pdf [49 KiB] HiDef png [65 KiB] Thumbnail [31 KiB] tex code 

Correlation matrix for the results in Eq. \eqref{eq:comb2} obtained taking into account correlated systematics between the considered analyses. 
Table_5.pdf [33 KiB] HiDef png [55 KiB] Thumbnail [26 KiB] tex code 

Values of the Swave parameters in each $m(K^{+}K^{})$ bin. The first uncertainty is statistical and the second systematic. 
Table_6.pdf [52 KiB] HiDef png [245 KiB] Thumbnail [111 KiB] tex code 
Created on 19 October 2019.