The first observation of the decay $\kstarkstar$ is reported using 35\invpb of data collected by LHCb in protonproton collisions at a centreofmass energy of 7 TeV. A total of $49.8 \pm 7.5$ $B^0_s \rightarrow (K^+\pi^)(K^\pi^+)$ events are {observed within $\pm 50 \mevcc$ of the \Bs mass and $746 \mevcc < m_{K\pi}< 1046 \mevcc$, mostly coming from a resonant $\kstarkstar$ signal.} The branching fraction and the \CPaveraged \Kstarz longitudinal polarization fraction are measured to be {$\BR(\kstarkstar) = (2.81 \pm 0.46 ({\rm stat.}) \pm 0.45 ({\rm syst.}) \pm 0.34\, (f_s/f_d))\times10^{5}$} and $f_L = 0.31 \pm 0.12 ({\rm stat.}) \pm 0.04 ({\rm syst.})$.
Fit to the $ K ^+ \pi ^ K ^ \pi ^+ $ mass distribution of selected candidates. The fit model (dashed pink curve) includes a signal component that has two Gaussian components corresponding to the $ B ^0_s$ and $ B ^0$ decays. The background is described as an exponential component (dotted blue) plus the parametrization indicated in the text (dashdotted green). 
Fig1.pdf [19 KiB] HiDef png [276 KiB] Thumbnail [81 KiB] *.C file 

Background subtracted $K^+\pi^$ and $K^\pi^+$ combinations for selected candidates within a $\pm$50 $ {\mathrm{\,MeV\!/}c^2}$ window of the $ B ^0_s$ mass. The solid blue line shows the projection of the 2D fit model described in the text, indicating the $ K ^{*0}$ $\overline{ K }{} ^{*0}$ yield (dasheddotted red line) and a nonresonant component (blue dotted line), assumed to be a linear function times the twobody phase space. The dashed red line indicates the overall $ B ^0_s \rightarrow K ^{*0} X$ contribution. 
Fig2.pdf [17 KiB] HiDef png [244 KiB] Thumbnail [69 KiB] *.C file 

Fit to the mass distribution of selected $ B ^0 \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} \! K ^{*0} $ events. The dashed red curve is the Gaussian component for the $ B $ signal. The green dasheddotted line accounts for partially reconstructed $ B\rightarrow J/\psi X$ (see Eq. 2). The pink hatched region accounts for a possible $ B ^0_s \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} \phi$ contamination, parametrized as a sum of two CrystalBall functions [20]. The combinatorial background is parametrized as an exponential and indicated as a blue dotted line. 
Fig3.pdf [18 KiB] HiDef png [238 KiB] Thumbnail [74 KiB] *.C file 

$\cos\theta$ (above) and $\varphi$ (below) acceptance corrected distributions for events in the narrow window around the $ B ^0_s$ mass. The blue line is the projection of the fit model given by Eq. 3 for the measured values of the parameters $f_L$, $f_{\parallel}$ and $\delta_{\parallel}$. The dotted lines indicate $\pm 1\sigma$ variation of the $f_L$ central value. 
Fig4.pdf [19 KiB] HiDef png [179 KiB] Thumbnail [51 KiB] *.C file 

Animated gif made out of all figures. 
PAPER2011012.gif Thumbnail 
Fitted values of the model parameters for the mass spectrum, as described in the text. $N_s$ and $N_d$ are the number of events for the $ B ^0_s$ and $ B ^0$ signals, $\mu_s$ is the fitted mass value for the $ B ^0_s$ signal and $\sigma$ is the Gaussian width. The mass difference between $ B ^0_s$ and $ B ^0$ was fixed to its nominal value [4]. $N_{b}$ is the number of background events in the full mass range (49005800 $ {\mathrm{\,MeV\!/}c^2}$ ), and $c_{b}$ is the exponential parameter in the fit. $M_p$, $\sigma_p$ and $k_p$ are the parameters of Eq. (1). Finally, $f_{p}$ is the fraction of the background associated with Eq. (1). 
Table_1.pdf [82 KiB] HiDef png [158 KiB] Thumbnail [25 KiB] tex code 

Selection and trigger efficiencies obtained from simulation. The observed yield found for the signal and control channels in the full mass range are also indicated. The efficiency errors are statistical, derived from the size of the simulated samples. 
Table_2.pdf [52 KiB] HiDef png [42 KiB] Thumbnail [8 KiB] tex code 

Estimated systematic error sources in the $\cal B \left( B ^0_s \rightarrow K ^{*0} \!\overline{ K }{} ^{*0} \right)$ measurement. 
Table_3.pdf [48 KiB] HiDef png [130 KiB] Thumbnail [20 KiB] tex code 
Created on 16 February 2019.Citation count from INSPIRE on 21 February 2019.