A study of $B^{0}_{s} \to \eta_{c} \phi$ and $B^{0}_{s} \to \eta_{c} \pi^{+}\pi^{}$ decays is performed using $pp$ collision data corresponding to an integrated luminosity of 3.0$ \rm fb^{1}$, collected with the LHCb detector in Run 1 of the LHC. The observation of the decay $B^{0}_{s} \to \eta_{c} \phi$ is reported, where the $\eta_{c}$ meson is reconstructed in the $p\bar p$, $K^+K^\pi^+\pi^$, $\pi^+\pi^\pi^+\pi^$ and $K^+K^K^+K^$ decay modes and the $\phi(1020)$ in the $K^+ K^$ decay mode. The decay $B^{0}_{s} \to J/\psi \phi$ is used as a normalisation channel. Evidence is also reported for the decay $B^{0}_{s} \to \eta_{c} \pi^{+}\pi^{}$, where the $\eta_{c}$ meson is reconstructed in the $p\bar p$ decay mode, using the decay $B^{0}_{s} \to J/\psi \pi^+ \pi^$ as a normalisation channel. The measured branching fractions are \begin{eqnarray*} {\mathcal B (B^{0}_{s} \to \eta_{c} \phi)} &=& \left(5.01 \pm 0.53 \pm 0.27 \pm 0.63 \right) \times 10^{4} , \nonumber \\ {\mathcal B (B^{0}_{s} \to \eta_{c} \pi^+ \pi^)} &=& \left(1.76 \pm 0.59 \pm 0.12 \pm 0.29 \right) \times 10^{4} , \end{eqnarray*} where in each case the first uncertainty is statistical, the second systematic and the third uncertainty is due to the limited knowledge of the external branching fractions.
Leading diagram corresponding to $ B ^0_ s \rightarrow \eta _ c \phi$ and $ B ^0_ s \rightarrow \eta _ c \pi ^+ \pi ^ $ decays, where the $\pi ^+ \pi ^ $ pair may arise from the decay of the $f_0(980)$ resonance. 
Fig1.pdf [8 KiB] HiDef png [95 KiB] Thumbnail [81 KiB] 

Distribution of $p\bar p$ invariantmass for $ B ^0_ s \rightarrow p\bar{p}\pi ^+ \pi ^ $ candidates obtained by the $_s{\cal P}lot$ technique. The solid black curve is the projection of the total fit result. The full blue, tightcrosshatched red and widecrosshatched black histograms show the $\eta _ c $ , $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu}$ and nonresonant $p\bar p$ contributions, respectively. The structure visible around $3.15 $GeV is found to be consistant with a statistical fluctuation. 
Fig2.pdf [27 KiB] HiDef png [835 KiB] Thumbnail [335 KiB] 

Onedimensional scan of $2 \Delta \ln L$ as a function of the $\eta _ c $ yield in the fit to the $p\bar{p}$ invariantmass distribution for $ B ^0_ s \rightarrow p\bar{p}\pi ^+ \pi ^ $ candidates. The dotted red and solid blue curves correspond to the result of the scan including statistical only and summed (statistical and systematic) uncertainties, respectively. 
Fig3.pdf [14 KiB] HiDef png [89 KiB] Thumbnail [67 KiB] 

Distribution of the $p\bar p$ and each of the $4h$ invariantmass spectra for $ B ^0_ s \rightarrow p\bar{p}(4h)\phi$ candidates obtained by the $_s{\cal P}lot$ technique. The solid black curve is the total result of the simultaneous fit. The full blue and tightcrosshatched red histograms show the $\eta _ c $ and $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu}$ contributions, respectively. In addition, for each of the $4h$ modes, the widecrosshatched and dotted black histograms show the $4h_{\rm NR}$ and $4h_{\rm bkg}$ contributions, respectively. The dashed blue curves show the interference between $ B ^0_ s \rightarrow \eta _ c (\rightarrow 4h)\phi$ and $ B ^0_ s \rightarrow (4h)_{\rm NR}\phi$ amplitudes. 
Fig4a.pdf [74 KiB] HiDef png [507 KiB] Thumbnail [250 KiB] 

Fig4b.pdf [78 KiB] HiDef png [505 KiB] Thumbnail [293 KiB] 

Fig4c.pdf [81 KiB] HiDef png [550 KiB] Thumbnail [301 KiB] 

Fig4d.pdf [80 KiB] HiDef png [455 KiB] Thumbnail [268 KiB] 

Distribution of the $p\bar{p}\pi ^+ \pi ^ $ invariant mass. Points with error bars show the data. The solid curve is the projection of the total fit result. The shortdashed blue, the dasheddoubledotted green, the dashedsingledotted yellow and mediumdashed red curves show the $ B ^0_ s \rightarrow p\bar{p}\pi ^+ \pi ^ $, $ B ^0 \rightarrow p\bar{p}\pi ^+ \pi ^ $, $ B ^0 \rightarrow p\bar{p} K ^+ \pi ^ $ and combinatorial background contributions, respectively. 
Fig5.pdf [21 KiB] HiDef png [306 KiB] Thumbnail [248 KiB] 

Distribution of the $p\bar{p} K ^+ K ^ $ (top) and $ K ^+ K ^ $ (bottom) invariant masses. Points with error bars show the data. The solid black curve is the projection of the total fit result. The shortdashed and dotted blue curves show the $ B ^0_ s \rightarrow p\bar{p}\phi$ and $ B ^0_ s \rightarrow p\bar{p} K ^+ K ^ $ contributions, respectively. The longdashed and mediumdashed red curves show the contributions of the combinatorial background with prompt $\phi$ and NR $ K ^+ K ^ $, respectively. The dasheddoubledotted green, dashedtripledotted green and dashedsingledottedyellow curves show the $ B ^0 \rightarrow p\bar{p}\phi$, $ B ^0 \rightarrow p\bar{p} K ^+ K ^ $ and $ B ^0 \rightarrow p\bar{p} K ^+ \pi ^ $ contributions, respectively. 
Fig6a.pdf [24 KiB] HiDef png [386 KiB] Thumbnail [273 KiB] 

Fig6b.pdf [20 KiB] HiDef png [294 KiB] Thumbnail [188 KiB] 

Distribution of the $2K2\pi K ^+ K ^ $ (top) and $ K ^+ K ^ $ (bottom) invariant masses. Points with error bars show the data. The solid black curve is the projection of the total fit result. The shortdashed and dotted blue curves show the $ B ^0_ s \rightarrow 2K2\pi\phi$ and $ B ^0_ s \rightarrow 2K2\pi K ^+ K ^ $ contributions, respectively. The longdashed and mediumdashed red curves show the contributions of the combinatorial background with prompt $\phi$ and NR $ K ^+ K ^ $, respectively. The dasheddoubledotted and dashedtripledotted green curves show the $ B ^0 \rightarrow 2K2\pi\phi$ and $ B ^0 \rightarrow 2K2\pi K ^+ K ^ $ contributions, respectively. 
Fig7a.pdf [30 KiB] HiDef png [344 KiB] Thumbnail [263 KiB] 

Fig7b.pdf [28 KiB] HiDef png [277 KiB] Thumbnail [207 KiB] 

Distribution of the $4\pi K ^+ K ^ $ (top) and $ K ^+ K ^ $ (bottom) invariant masses. Points with error bars show the data. The solid black curve is the projection of the total fit result. The shortdashed and dotted blue curves show the $ B ^0_ s \rightarrow 4\pi\phi$ and $ B ^0_ s \rightarrow 4\pi K ^+ K ^ $ contributions, respectively. The longdashed and mediumdashed red curves show the contributions of the combinatorial background with prompt $\phi$ and NR $ K ^+ K ^ $, respectively. The dasheddoubledotted and dashedtripledotted green curves show the $ B ^0 \rightarrow 4\pi\phi$ and $ B ^0 \rightarrow 4\pi K ^+ K ^ $ contributions, respectively. 
Fig8a.pdf [30 KiB] HiDef png [356 KiB] Thumbnail [274 KiB] 

Fig8b.pdf [29 KiB] HiDef png [285 KiB] Thumbnail [213 KiB] 

Distribution of the $4K K ^+ K ^ $ (top) and $ K ^+ K ^ $ (bottom) invariant masses. Points with error bars show the data. The solid black curve is the projection of the total fit result. The shortdashed and dotted blue curves show the $ B ^0_ s \rightarrow 4K\phi$ and $ B ^0_ s \rightarrow 4K K ^+ K ^ $ contributions, respectively. The longdashed and mediumdashed red curves show the contributions of the combinatorial background with prompt $\phi$ and NR $ K ^+ K ^ $, respectively. The dasheddoubledotted and dashedtripledotted green curves show the $ B ^0 \rightarrow 4K\phi$ and $ B ^0 \rightarrow 4K K ^+ K ^ $ contributions, respectively. 
Fig9a.pdf [29 KiB] HiDef png [315 KiB] Thumbnail [237 KiB] 

Fig9b.pdf [27 KiB] HiDef png [234 KiB] Thumbnail [171 KiB] 

Animated gif made out of all figures. 
PAPER2016056.gif Thumbnail 
Ratio of efficiencies between the normalisation and signal channels for each final state. 
Table_1.pdf [56 KiB] HiDef png [26 KiB] Thumbnail [12 KiB] tex code 

Yields of the different final states as obtained from the fit to the $p\bar{p}\pi ^+ \pi ^ $ invariantmass distribution and from the 2D fits in the $p\bar{p}(4h) K ^+ K ^ \times K ^+ K ^ $ invariantmass planes. Only statistical uncertainties are reported. The abbreviation "n/a" stands for "not applicable". 
Table_2.pdf [52 KiB] HiDef png [90 KiB] Thumbnail [48 KiB] tex code 

Results of the fit to the $p\bar p$ invariantmass spectra weighted for $ B ^0_ s \rightarrow p\bar{p}\pi ^+ \pi ^ $ candidates. Uncertainties are statistical only. The parameter $N_{\rm NR}$ corresponds to the yield of $B^0_s \rightarrow (p\bar p)_{\rm NR} \pi^+ \pi^$ candidates. The $\eta _ c $ yield does not appear since it is parametrised as a function of $\mathcal{B} ( B ^0_ s \rightarrow \eta _ c \pi ^+ \pi ^ )$, the measured value of which is reported in Eq. 22. 
[Error creating the table]  
Result of the simultaneous fit to the $p\bar{p}$ and $4h$ invariantmass spectra. Uncertainties are statistical only. The $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu}$ and $\eta _ c $ magnitudes do not appear since they are set to unity as reference and parametrised as a function of $\mathcal{B} ( B ^0_ s \rightarrow \eta _ c \phi)$, respectively. In the simultaneous fit, the $m_{\eta _ c }$ and $m_{ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} }$ parameters are shared across the four modes. The measured value of $\mathcal{B} ( B ^0_ s \rightarrow \eta _ c \phi)$ is reported in Eq. 23. The abbreviation "n/a" stands for "not applicable". 
Table_4.pdf [55 KiB] HiDef png [92 KiB] Thumbnail [49 KiB] tex code 

Fit fractions obtained from the parameters of the simultaneous fit to the $p\bar{p}$ and $4h$ invariantmass spectra. Uncertainties are statistical only. Due to interference between $ B ^0_ s \rightarrow \eta _ c (\rightarrow 4h) \phi$ and $ B ^0_ s \rightarrow (4h)_{\rm NR} \phi$ amplitudes, for the $4h$ final states the sum of fit fractions, $\sum_k {\rm FF}_k $, may be different from unity. The abbreviation "n/a" stands for "not applicable". 
Table_5.pdf [58 KiB] HiDef png [67 KiB] Thumbnail [37 KiB] tex code 

Summary of systematic uncertainties. The "Sum" of systematic uncertainties is obtained from the quadratic sum of the individual sources, except the external branching fractions, which are quoted separately. All values are in $\%$ of the measured branching fractions. The abbreviation "n/a" stands for "not applicable". 
Table_6.pdf [56 KiB] HiDef png [85 KiB] Thumbnail [41 KiB] tex code 

Statistical correlation matrix for the parameters from the simultaneous fit to the $p\bar{p}$ and $4h$ invariantmass spectra for $ B ^0_ s \rightarrow p\bar{p}\phi$ and $ B ^0_ s \rightarrow 4h\phi$ candidates. 
Table_7.pdf [53 KiB] HiDef png [217 KiB] Thumbnail [99 KiB] tex code 
Created on 18 January 2020.