A search for the leptonflavour violating decays $B^0_s \rightarrow e^\pm \mu^\mp$ and $B^0 \rightarrow e^\pm \mu^\mp$ is performed based on a sample of protonproton collision data corresponding to an integrated luminosity of 3 fb$^{1}$, collected with the LHCb experiment at centreofmass energies of 7 and 8 TeV. The observed yields are consistent with the backgroundonly hypothesis. Upper limits on the branching fractions are determined to be $\mathcal{B}(B^0_s \rightarrow e^\pm \mu^\mp) < 5.4 (6.3) \times 10^{9}$ and $\mathcal{B}(B^0 \rightarrow e^\pm \mu^\mp) < 1.0 (1.3) \times 10^{9}$ at $90\% (95\%)$ confidence level, which are the strongest limits on these decays to date.
Expected distribution of the BDT response for $ B ^0_{( s )} \rightarrow e^\pm\mu^\mp$ decays with recovered bremsstrahlung photons obtained from the $ B ^0 \rightarrow K ^+ \pi ^ $ control channel. The total uncertainty is shown as a light grey band. Each bin is normalised to its width. 
Fig1.pdf [13 KiB] HiDef png [68 KiB] Thumbnail [68 KiB] *.C file 

Invariant mass distributions of the two normalisation channels with fit functions superimposed: (left) $ B ^0 \rightarrow K ^+ \pi ^ $ and (right) $ B ^+ \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K ^+ $ . Pull distributions are shown below each plot. 
Fig2left.pdf [31 KiB] HiDef png [251 KiB] Thumbnail [212 KiB] *.C file 

Fig2right.pdf [75 KiB] HiDef png [281 KiB] Thumbnail [244 KiB] *.C file 

Distribution of the $ m_{e^{\pm}\mu^{\mp}}$ invariant mass of simulated $ B ^0_ s $ candidates with no (left) and one or more (right) recovered bremsstrahlung photons. The overlaid fit function is a modified Crystal Ball function with two tails on opposite sides. 
Fig3left.pdf [29 KiB] HiDef png [188 KiB] Thumbnail [172 KiB] *.C file 

Fig3right.pdf [32 KiB] HiDef png [182 KiB] Thumbnail [166 KiB] *.C file 

Distributions of the invariant mass of the $ B ^0_{( s )} \rightarrow e^\pm\mu^\mp$ candidates, $ m_{e^{\pm}\mu^{\mp}}$ , divided into bins of BDT response and two bremsstrahlung categories (left) without and (right) with bremsstrahlung photons recovered. The result of the fit is overlaid and the different components are detailed. The edges of the range that was examined only after finalising the selection and fit procedure are delimited by gray dashed vertical lines. This region includes 90% of the potential signal candidates. Given the result obtained from the fit, the $ B ^0 \rightarrow e^\pm\mu^\mp$ component is not visible in the plots. 
Fig4.pdf [46 KiB] HiDef png [523 KiB] Thumbnail [458 KiB] *.C file 

Results of the $\mathrm{CL_s}$ scan used to obtain the limit on (left) $\mathcal{B} ( B ^0 \rightarrow e^\pm\mu^\mp )$ and (right) $\mathcal{B} ( B ^0_ s \rightarrow e^\pm\mu^\mp )$. The backgroundonly expectation is shown by the dashed line and the $1\sigma$ and $2\sigma$ bands are shown as dark (green) and light (yellow) bands respectively. The observed limit is shown as the solid black line. 
Fig5left.pdf [14 KiB] HiDef png [114 KiB] Thumbnail [93 KiB] *.C file 

Fig5right.pdf [13 KiB] HiDef png [120 KiB] Thumbnail [98 KiB] *.C file 

Animated gif made out of all figures. 
PAPER2017031.gif Thumbnail 
Yields of normalisation channels obtained from fits to data. 
Table_1.pdf [42 KiB] HiDef png [16 KiB] Thumbnail [8 KiB] tex code 

Expected (assuming no signal) and observed upper limits for $\mathcal{B} ( B ^0_ s \rightarrow e^\pm\mu^\mp )$ and $\mathcal{B} ( B ^0 \rightarrow e^\pm\mu^\mp )$ at 95% (90%) CL. The upper limit on the $\mathcal{B} ( B ^0_ s \rightarrow e^\pm\mu^\mp )$ is evaluated under the assumption of pure heavy eigenstate contribution on the decay amplitude. 
Table_2.pdf [51 KiB] HiDef png [31 KiB] Thumbnail [15 KiB] tex code 
Created on 17 August 2019.Citation count from INSPIRE on 23 August 2019.