An improved measurement of the decay $B^0_s\to\mu^+\mu^-$ and searches for the decays $B^0\to\mu^+\mu^-$ and $B^0_s\to\mu^+\mu^-\gamma$ are performed at the LHCb experiment using data collected in proton-proton collisions at $\sqrt{s} = 7, 8$ and $13$ TeV, corresponding to integrated luminosities of 1, 2 and 6 fb$^{-1}$, respectively. The $B^0_s\to\mu^+\mu^-$ branching fraction and effective lifetime are measured to be ${\mathcal{B}}(B^0_s\to\mu^+\mu^-)=\left(3.09^{+0.46+0.15}_{-0.43-0.11}\right)\times 10^{-9}$ and $\tau(B^0_s\to\mu^+\mu^-)=(2.07\pm 0.29\pm 0.03)$ ps, respectively, where the uncertainties include both statistical and systematic contributions. No significant signal for $B^0\to\mu^+\mu^-$ and $B^0_s\to\mu^+\mu^-\gamma$ decays is found and the upper limits $\mathcal{B}(B^0\to\mu^+\mu^-)<2.6\times 10^{-10}$ and $\mathcal{B}(B^0_s\to\mu^+\mu^-\gamma)<2.0\times 10^{-9}$ at 95 confidence level are determined, where the latter is limited to the range $m_{\mu\mu} > 4.9$ GeV$/c^2$. Additionally, the ratio between the $B^0\to\mu^+\mu^-$ and $B^0_s\to\mu^+\mu^-$ branching fractions is measured to be $\mathcal{R}_{\mu^+\mu^-}<0.095$ at 95 confidence level. The results are in agreement with the Standard Model predictions.
SM Feynman diagrams mediating (top) the $ B ^0_ s \rightarrow \mu ^+ \mu ^- $ and (bottom) the $ B ^0_ s \rightarrow \mu ^+ \mu ^- \gamma $ processes. Subpanels show (a) the so-called "penguin" diagram and (b) the "box" diagram for $ B ^0_ s \rightarrow \mu ^+ \mu ^- $ , and (c) an ISR contribution and (d) an FSR contribution to $ B ^0_ s \rightarrow \mu ^+ \mu ^- \gamma $ . |
[Failure to get the plot] | |
BDT distribution calibrated using corrected simulated $ B ^0_ s \rightarrow \mu ^+ \mu ^- $ decays (black circles) and combinatorial background from high dimuon-mass data sidebands (blue filled circles) in (left) Run 1 and (right) Run 2 data. Blue error bands represent the statistical uncertainty. |
[Failure to get the plot] | [Failure to get the plot] |
Mass distributions of selected (top) $ B ^0 \rightarrow K ^+ \pi ^- $ and (bottom) $ B ^0_ s \rightarrow K ^+ K ^- $ candidates in (left) Run 1 and (right) Run 2 data. The results of the fits used to determine the means of the $ B _{( s )}^0 \rightarrow \mu ^+ \mu ^- $ mass distributions are overlaid and the different components are detailed in the legends. |
Fig3_t[..].pdf [38 KiB] HiDef png [311 KiB] Thumbnail [226 KiB] *.C file |
|
Fig3_t[..].pdf [37 KiB] HiDef png [311 KiB] Thumbnail [226 KiB] *.C file |
||
Fig3_b[..].pdf [38 KiB] HiDef png [311 KiB] Thumbnail [233 KiB] *.C file |
||
Fig3_b[..].pdf [38 KiB] HiDef png [298 KiB] Thumbnail [219 KiB] *.C file |
||
Mass distributions of (top) $ { J \mskip -3mu/\mskip -2mu\psi } \rightarrow \mu^+\mu^-$, (centre) $\psi{(2S)} \rightarrow \mu^+\mu^-$, (bottom) $\varUpsilon(1S,2S,3S)\rightarrow \mu^+\mu^-$ candidates in (left) Run 1 and (right) Run 2 data. The result from the fit to determine the mass resolutions to each sample is overlaid, and the components are detailed in the legend. |
Fig4_t[..].pdf [23 KiB] HiDef png [228 KiB] Thumbnail [193 KiB] *.C file |
|
Fig4_t[..].pdf [23 KiB] HiDef png [229 KiB] Thumbnail [192 KiB] *.C file |
||
Fig4_c[..].pdf [22 KiB] HiDef png [233 KiB] Thumbnail [196 KiB] *.C file |
||
Fig4_c[..].pdf [22 KiB] HiDef png [213 KiB] Thumbnail [178 KiB] *.C file |
||
Fig4_b[..].pdf [25 KiB] HiDef png [327 KiB] Thumbnail [242 KiB] *.C file |
||
Fig4_b[..].pdf [25 KiB] HiDef png [333 KiB] Thumbnail [251 KiB] *.C file |
||
Fit to the measured mass resolutions of quarkonia resonances (blue dots) to obtain the mass resolution at the $ B ^0$ and $ B ^0_ s $ masses as indicated by the two dashed lines in (left) Run 1 and (right) Run 2 data. |
Fig5_left.pdf [73 KiB] HiDef png [146 KiB] Thumbnail [135 KiB] *.C file |
|
Fig5_right.pdf [73 KiB] HiDef png [144 KiB] Thumbnail [132 KiB] *.C file |
||
The calibrated BDT distribution for $ B ^0_ s \rightarrow \mu ^+ \mu ^- $ decays with $ A^{\mu\mu}_{\Delta\Gamma_s} =1$ (black) and $ B ^0 \rightarrow \mu ^+ \mu ^- $ decays (red) for (left) Run 1 and (right) Run 2, including the total uncertainty on the fraction per BDT region. The $ B ^0_ s \rightarrow \mu ^+ \mu ^- $ and $ B ^0 \rightarrow \mu ^+ \mu ^- $ distributions are determined on corrected simulated samples, as described in the text. |
Fig6_left.pdf [15 KiB] HiDef png [160 KiB] Thumbnail [125 KiB] *.C file |
|
Fig6_right.pdf [15 KiB] HiDef png [161 KiB] Thumbnail [126 KiB] *.C file |
||
The BDT distributions of (top) $ B ^0_ s \rightarrow \mu ^+ \mu ^- $ and $ B ^0_ s \rightarrow K ^+ K ^- $ decays and (bottom) $ B ^0 \rightarrow \mu ^+ \mu ^- $ and $ B ^0 \rightarrow K ^+ \pi ^- $ decays in (left) Run 1 and (right) Run 2 data, including the total uncertainty on the fraction per BDT region. For $ B ^0_ s \rightarrow \mu ^+ \mu ^- $ and $ B ^0 \rightarrow \mu ^+ \mu ^- $ decays, the distributions are determined on corrected simulated samples, as described in the text, and are shown in black. The $ B ^0_ s \rightarrow K ^+ K ^- $ and $ B ^0 \rightarrow K ^+ \pi ^- $ distributions are determined using fits to data as described in the text and are shown in red. |
Fig7_t[..].pdf [15 KiB] HiDef png [151 KiB] Thumbnail [116 KiB] *.C file |
|
Fig7_t[..].pdf [15 KiB] HiDef png [148 KiB] Thumbnail [115 KiB] *.C file |
||
Fig7_b[..].pdf [15 KiB] HiDef png [146 KiB] Thumbnail [113 KiB] *.C file |
||
Fig7_b[..].pdf [15 KiB] HiDef png [148 KiB] Thumbnail [116 KiB] *.C file |
||
Mass distribution of $ B ^+ \rightarrow { J \mskip -3mu/\mskip -2mu\psi } K ^+ $ candidates in data for different data-taking years. Superimposed is a fit to the distribution: the blue line shows the total fit, the red dashed line is the $ B ^+ \rightarrow { J \mskip -3mu/\mskip -2mu\psi } K ^+ $ component, the green dash-dotted line is the combinatorial background, the purple dash-three-dotted line is the $ B ^+ \rightarrow { J \mskip -3mu/\mskip -2mu\psi } \pi ^+ $ misidentified background. |
Fig8_t[..].pdf [36 KiB] HiDef png [233 KiB] Thumbnail [198 KiB] *.C file |
|
Fig8_t[..].pdf [37 KiB] HiDef png [238 KiB] Thumbnail [202 KiB] *.C file |
||
Fig8_c[..].pdf [36 KiB] HiDef png [238 KiB] Thumbnail [201 KiB] *.C file |
||
Fig8_c[..].pdf [37 KiB] HiDef png [239 KiB] Thumbnail [203 KiB] *.C file |
||
Fig8_b[..].pdf [37 KiB] HiDef png [239 KiB] Thumbnail [204 KiB] *.C file |
||
Fig8_b[..].pdf [37 KiB] HiDef png [240 KiB] Thumbnail [205 KiB] *.C file |
||
Mass distribution of $ B ^0 \rightarrow K ^+ \pi ^- $ candidates in data for different data-taking years, triggered independently of the signal. Superimposed is a fit to the distribution: the blue line shows the total fit, the red dashed line is the $ B ^0 \rightarrow K ^+ \pi ^- $ component, the magenta dashed line is the $\overline{ B } {}^0_ s \rightarrow K ^+ \pi ^- $ component, the green dashed line is the combinatorial background, and the brown dashed line is the partially reconstructed background component. |
Fig9_t[..].pdf [37 KiB] HiDef png [333 KiB] Thumbnail [251 KiB] *.C file |
|
Fig9_t[..].pdf [37 KiB] HiDef png [314 KiB] Thumbnail [239 KiB] *.C file |
||
Fig9_c[..].pdf [37 KiB] HiDef png [331 KiB] Thumbnail [251 KiB] *.C file |
||
Fig9_c[..].pdf [37 KiB] HiDef png [291 KiB] Thumbnail [223 KiB] *.C file |
||
Fig9_b[..].pdf [38 KiB] HiDef png [295 KiB] Thumbnail [224 KiB] *.C file |
||
Fig9_b[..].pdf [38 KiB] HiDef png [299 KiB] Thumbnail [229 KiB] *.C file |
||
Mass distribution of signal candidates (black dots) for (left) Run 1 and (right) Run 2 samples in regions of BDT. The result of the fit is overlaid (blue line) and the different components detailed in the legend. The solid bands represent the variation of the signal branching fractions within their total uncertainty. |
Fig10_[..].pdf [23 KiB] HiDef png [306 KiB] Thumbnail [211 KiB] *.C file |
|
Fig10_[..].pdf [22 KiB] HiDef png [296 KiB] Thumbnail [185 KiB] *.C file |
||
Fig10_[..].pdf [22 KiB] HiDef png [319 KiB] Thumbnail [192 KiB] *.C file |
||
Fig10_[..].pdf [23 KiB] HiDef png [321 KiB] Thumbnail [205 KiB] *.C file |
||
Fig10_[..].pdf [22 KiB] HiDef png [337 KiB] Thumbnail [202 KiB] *.C file |
||
Fig10_[..].pdf [23 KiB] HiDef png [318 KiB] Thumbnail [196 KiB] *.C file |
||
Fig10_[..].pdf [21 KiB] HiDef png [268 KiB] Thumbnail [161 KiB] *.C file |
||
Fig10_[..].pdf [22 KiB] HiDef png [309 KiB] Thumbnail [195 KiB] *.C file |
||
Fig10_[..].pdf [21 KiB] HiDef png [329 KiB] Thumbnail [187 KiB] *.C file |
||
Fig10_[..].pdf [24 KiB] HiDef png [453 KiB] Thumbnail [278 KiB] *.C file |
||
Two-dimensional representations of the branching fraction measurements for the decays (top) $ B ^0_ s \rightarrow \mu ^+ \mu ^- $ and $ B ^0 \rightarrow \mu ^+ \mu ^- $ , (bottom left) $ B ^0 \rightarrow \mu ^+ \mu ^- $ and $ B ^0_ s \rightarrow \mu ^+ \mu ^- \gamma $ and (bottom right) $ B ^0_ s \rightarrow \mu ^+ \mu ^- $ and $ B ^0_ s \rightarrow \mu ^+ \mu ^- \gamma $ . The $ B ^0_ s \rightarrow \mu ^+ \mu ^- \gamma $ branching fraction is limited to the range $m_{\mu\mu}>4.9\text{ Ge V /}c^2 $. The measured central values of the branching fractions are indicated with a blue dot. The profile likelihood contours for 68%, 95% and 99% CL regions of the result presented in this paper are shown as blue contours, while in the top plot the brown contours indicate the previous measurement [32] and the red cross shows the SM prediction. Figure reproduced from Ref. [37]. |
Fig11_top.pdf [16 KiB] HiDef png [364 KiB] Thumbnail [217 KiB] *.C file |
|
Fig11_[..].pdf [14 KiB] HiDef png [261 KiB] Thumbnail [163 KiB] *.C file |
||
Fig11_[..].pdf [15 KiB] HiDef png [248 KiB] Thumbnail [150 KiB] *.C file |
||
Results from the CL$_{\text{s}}$ scan used to obtain the limit on (left) $\cal B ( B ^0 \rightarrow \mu ^+ \mu ^- )$ and (right) $\cal B ( B ^0_ s \rightarrow \mu ^+ \mu ^- \gamma )$. The background-only expectation is shown by the red line and the 1$\sigma$ and 2$\sigma$ bands are shown as light blue and blue bands respectively. The observation is shown as the solid black line. The two dashed lines intersecting with the observation indicate the limits at 90% and 95% CL for the upper and lower line, respectively. |
Fig12_left.pdf [19 KiB] HiDef png [300 KiB] Thumbnail [199 KiB] *.C file |
|
Fig12_[..].pdf [20 KiB] HiDef png [302 KiB] Thumbnail [197 KiB] *.C file |
||
Dimuon mass distributions of $ B ^0_ s \rightarrow \mu ^+ \mu ^- $ candidates with the fit model used to perform the background subtraction for the measurement of the $ B ^0_ s \rightarrow \mu ^+ \mu ^- $ effective lifetime superimposed in the (left) low and (right) high BDT regions. |
Fig13_left.pdf [17 KiB] HiDef png [283 KiB] Thumbnail [241 KiB] *.C file |
|
Fig13_[..].pdf [16 KiB] HiDef png [271 KiB] Thumbnail [218 KiB] *.C file |
||
The functions used to model the decay-time efficiency in the (left) low and (right) high BDT regions in the fit for the $ B ^0_ s \rightarrow \mu ^+ \mu ^- $ effective lifetime. |
Fig14_left.pdf [14 KiB] HiDef png [162 KiB] Thumbnail [142 KiB] *.C file |
|
Fig14_[..].pdf [14 KiB] HiDef png [159 KiB] Thumbnail [141 KiB] *.C file |
||
(Top) Distribution of $ K ^+ K ^- $ mass with the fit models used to perform the background subtraction superimposed and (bottom) the background-subtracted decay-time distributions with the fit model used to determine the $ B ^0_ s \rightarrow K ^+ K ^- $ effective lifetime superimposed (bottom row). The distributions in the low and high BDT regions are shown in the left and right columns, respectively. |
Fig15_[..].pdf [18 KiB] HiDef png [282 KiB] Thumbnail [233 KiB] *.C file |
|
Fig15_[..].pdf [18 KiB] HiDef png [274 KiB] Thumbnail [224 KiB] *.C file |
||
Fig15_[..].pdf [15 KiB] HiDef png [191 KiB] Thumbnail [167 KiB] *.C file |
||
Fig15_[..].pdf [15 KiB] HiDef png [205 KiB] Thumbnail [180 KiB] *.C file |
||
(Top) Distribution of $ K ^+ \pi ^- $ mass with the fit models used to perform the background subtraction superimposed and (bottom) the background-subtracted decay-time distributions with the fit model used to determine the $ B ^0 \rightarrow K ^+ \pi ^- $ lifetime superimposed. The distributions in the low and high BDT regions are shown in the left and right columns, respectively. |
Fig16_[..].pdf [18 KiB] HiDef png [310 KiB] Thumbnail [247 KiB] *.C file |
|
Fig16_[..].pdf [18 KiB] HiDef png [290 KiB] Thumbnail [224 KiB] *.C file |
||
Fig16_[..].pdf [15 KiB] HiDef png [183 KiB] Thumbnail [152 KiB] *.C file |
||
Fig16_[..].pdf [15 KiB] HiDef png [191 KiB] Thumbnail [174 KiB] *.C file |
||
Background-subtracted decay-time distributions with the fit model used to extract the $ B ^0_ s \rightarrow \mu ^+ \mu ^- $ effective lifetime superimposed in the (left) low and (right) high BDT regions. |
Fig17_left.pdf [15 KiB] HiDef png [178 KiB] Thumbnail [162 KiB] *.C file |
|
Fig17_[..].pdf [15 KiB] HiDef png [174 KiB] Thumbnail [161 KiB] *.C file |
||
Animated gif made out of all figures. |
PAPER-2021-008.gif Thumbnail |
Luminosity-weighted signal mass shape parameter combinations per data set, including propagated uncertainties. Where appropriate, statistical and systematic uncertainties are added in quadrature. As the tail parameters determined in Run 1 and Run 2 are consistent, they are combined into common estimates. |
Table_1.pdf [61 KiB] HiDef png [81 KiB] Thumbnail [40 KiB] tex code |
|
Yields of the two normalisation channels with their combined statistical and systematic errors. |
Table_2.pdf [61 KiB] HiDef png [121 KiB] Thumbnail [61 KiB] tex code |
|
Efficiencies of reconstruction within the LHCb detector and selection for the signal and normalisation channels, averaged for the two running periods. The uncertainties include the statistical uncertainty from the simulated samples and the uncertainty of the tracking efficiency corrections. |
Table_3.pdf [67 KiB] HiDef png [89 KiB] Thumbnail [41 KiB] tex code |
|
Particle identification efficiencies for the signal and normalisation channels, averaged for the two running periods, where the first uncertainty is statistical and the second systematic. A data-simulation correction of the muon-system identification is included for channels with muons. For the $ B ^+ \rightarrow { J \mskip -3mu/\mskip -2mu\psi } K ^+ $ channel only the data-simulation correction part of the muon identification is reported, as no multivariate PID requirement is applied to this channel. |
Table_4.pdf [67 KiB] HiDef png [72 KiB] Thumbnail [33 KiB] tex code |
|
Trigger efficiencies for the signal and normalisation channels, averaged for the two data taking periods. The first uncertainty is statistical and the second systematic. |
Table_5.pdf [67 KiB] HiDef png [73 KiB] Thumbnail [35 KiB] tex code |
|
Efficiency on the signal channels of excluding the BDT region $\text{BDT}<0.25$, averaged for the two data taking periods. The uncertainties combine statistical and the second systematic uncertainties. |
Table_6.pdf [65 KiB] HiDef png [57 KiB] Thumbnail [26 KiB] tex code |
|
Single-event sensitivities, $\alpha( B ^+ )$, $\alpha( B ^0 )$ and $\alpha(\rm{Comb})$ for the three signal channels obtained for $ BDT>0.25$ with the two normalisation channels, $ B ^0 \rightarrow K ^+ \pi ^- $ and $ B ^+ \rightarrow { J \mskip -3mu/\mskip -2mu\psi } K ^+ $ , and combined, for Run 1 , Run 2 and the full data set. The first uncertainty is statistical and the second systematic. The expected yields assuming SM branching fractions, $N_{\rm{exp}}$, are also reported. The $ B ^0_ s \rightarrow \mu ^+ \mu ^- \gamma $ expected number does not include an uncertainty on the branching fraction. |
Table_7.pdf [69 KiB] HiDef png [124 KiB] Thumbnail [60 KiB] tex code |
|
Expected background yields per BDT region and for $ BDT >0.25$ with their total estimated uncertainties for Run 1 data. |
Table_8.pdf [70 KiB] HiDef png [53 KiB] Thumbnail [21 KiB] tex code |
|
Expected background yields per BDT region and for $ BDT >0.25$ with their total estimated uncertainties for Run 2 data. |
Table_9.pdf [70 KiB] HiDef png [48 KiB] Thumbnail [20 KiB] tex code |
|
Summary of the systematic uncertainties for the measurement of the $ B ^0_ s \rightarrow \mu ^+ \mu ^- $ effective lifetime. |
Table_10.pdf [61 KiB] HiDef png [50 KiB] Thumbnail [21 KiB] tex code |
Created on 18 May 2024.