Search for the lepton-flavour violating decays

$$B_{\_}(s)^{0} \rightarrow e^{\pm} \mu^{\mp}$$

Roel Aaij, Bernardo Adeva, Marco Adinolfi, Ziad Ajaltouni, Simon Akar, Johannes Albrecht, Federico Alessio, Michael Alexander, Alejandro Alfonso Albero, Suvayu Ali, et al.

► To cite this version:

Roel Aaij, Bernardo Adeva, Marco Adinolfi, Ziad Ajaltouni, Simon Akar, et al.. Search for the lepton-flavour violating decays $B_{\_}(s)^{0} \rightarrow e^{\pm} \mu^{\mp}$. JHEP, 2018, 03, pp.078. 10.1007/JHEP03(2018)078. hal-01758630

HAL Id: hal-01758630

https://hal.science/hal-01758630

Submitted on 18 Jul 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Search for the lepton-flavour violating decays $B_s^0 \to e^\pm\mu^\mp$

LHCb collaboration†

## Abstract

A search for the lepton-flavour violating decays $B_s^0 \to e^\pm\mu^\mp$ and $B^0 \to e^\pm\mu^\mp$ is performed based on a sample of proton-proton collision data corresponding to an integrated luminosity of 3 fb⁻¹, collected with the LHCb experiment at centre-of-mass energies of 7 and 8 TeV. The observed yields are consistent with the background-only hypothesis. Upper limits on the branching fraction of the $B_s^0 \to e^\pm\mu^\mp$ decays are evaluated both in the hypotheses of an amplitude completely dominated by the heavy eigenstate and by the light eigenstate. The results are $\mathcal{B}(B_s^0 \to e^\pm\mu^\mp) < 6.3 (5.4) \times 10^{-9}$ and $\mathcal{B}(B_s^0 \to e^\pm\mu^\mp) < 7.2 (6.0) \times 10^{-9}$ at 95% (90%) confidence level, respectively. The upper limit on the branching fraction of the $B^0 \to e^\pm\mu^\mp$ decay is also evaluated, obtaining $\mathcal{B}(B^0 \to e^\pm\mu^\mp) < 1.3 (1.0) \times 10^{-9}$ at 95% (90%) confidence level. These are the strongest limits on these decays to date.

Published in JHEP 03 (2018) 078

© CERN on behalf of the LHCb collaboration, licence CC-BY-4.0.

<footnote>†Authors are listed at the end of this paper.</footnote>
# 1 Introduction

Processes that are suppressed or forbidden in the Standard Model (SM) are sensitive to potential contributions from new mediators, even if their masses are inaccessible to direct searches. Despite the fact that lepton-flavour violating (LFV) decays are forbidden within the SM, neutrino oscillation phenomena are proof that lepton flavour is not conserved in the neutral sector. However, LFV decays have not yet been observed, and their observation would be clear evidence of physics beyond the SM.

The study of LFV decays is particularly interesting in light of hints of lepton non-universality (LNU) effects in semileptonic decays [1] and $b \to s\ell\ell$ transitions [2, 3], which could be associated with LFV processes [4]. Possible explanations of these hints can be found in various scenarios beyond the SM, e.g. models with a new gauge $Z'$ boson [5] or leptoquarks [6, 7]. In these models, the branching fractions of the $B_s^0 \to e^\pm\mu^\mp$ and $B^0 \to e^\pm\mu^\mp$ decays<fnref n="1" /> can be enhanced up to $10^{-11}$. Other models also predict possible enhancement for $B_s^0 \to e^\pm\mu^\mp$ and $B^0 \to e^\pm\mu^\mp$ decays, e.g. heavy singlet Dirac neutrinos [8], supersymmetric models [9] and the Pati-Salam model [10]. The most stringent published limits on the branching fractions of these decays are currently $\mathcal{B}(B_s^0 \to e^\pm\mu^\mp) < 1.4 \times 10^{-8}$ and $\mathcal{B}(B^0 \to e^\pm\mu^\mp) < 3.7 \times 10^{-9}$ at 95% confidence level (CL) from the LHCb collaboration using data corresponding to $1 \text{ fb}^{-1}$ of integrated luminosity [11].

This article presents an analysis performed on a larger data sample, corresponding to an integrated luminosity of $3 \text{ fb}^{-1}$ of pp collisions collected at centre-of-mass energies of 7 and 8 TeV by the LHCb experiment in 2011 and 2012. In addition to a larger data sample, this analysis benefits from an improved selection and in particular a better performing multivariate classifier for signal and background separation. It supersedes the previous LHCb search for $B_s^0 \to e^\pm\mu^\mp$ and $B^0 \to e^\pm\mu^\mp$ decays [11].

Two normalisation channels are used: the $B^0 \to K^+\pi^-$ decay which has a similar topology to that of the signal, and the $B^+ \to J/\psi K^+$ decay, with $J/\psi \to \mu^+\mu^-$, which has an abundant yield and a similar purity and trigger selection. To avoid potential biases, $B_{(s)}^0 \to e^\pm\mu^\mp$ candidates in the signal region, $m_{e^\pm\mu^\mp} \in [5100, 5500]$ MeV/c², where $m_{e^\pm\mu^\mp}$ is the invariant mass of the $e^\pm\mu^\mp$ pair, were not examined until the selection and fitting procedure were finalised.

# 2 Detector and simulation

The LHCb detector [12, 13] is a single-arm forward spectrometer covering the pseudorapidity range $2 < \eta < 5$, designed for the study of particles containing b or c quarks. The detector includes a high-precision tracking system consisting of a silicon-strip vertex detector surrounding the pp interaction region, a large-area silicon-strip detector located upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip detectors and straw drift tubes placed downstream of the magnet. The tracking system provides a measurement of momentum, $p$, of charged particles with a relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200 GeV/c. The minimum distance of a track to a primary vertex (PV), the impact parameter (IP), is measured with a resolution of (15 + 29/$p_T$) µm, where $p_T$ is the component of the momentum transverse to the beam, in GeV/c. Different types of charged hadrons

<footnote n="1">¹Inclusion of charge conjugate processes is implied throughout the text.</footnote>
are distinguished using information from two ring-imaging Cherenkov detectors. Photons, electrons and hadrons are identified by a calorimeter system consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers.

The online event selection is performed by a trigger, which consists of a hardware stage, based on information from the muon and calorimeter systems, followed by a software stage that applies a full reconstruction of the event. The $B_{(s)}^0 \rightarrow e^{\pm}\mu^{\mp}$ candidates must fulfill the requirements of the electron or muon triggers. At the hardware stage, the electron trigger requires the presence of a cluster in the electromagnetic calorimeter with a transverse energy deposit, $E_T$, of at least 2.5 (3.0) GeV for 2011 (2012) data. The muon trigger selects muon candidates with $p_T$ higher than 1.5 (1.8) GeV/c for 2011 (2012) data. The software stage requires a two-track secondary vertex identified by a multivariate algorithm [14] to be consistent with the decay of a b hadron with at least two charged tracks, and at least one track with high $p_T$ and large IP with respect to any PV.

Simulated samples are used to evaluate geometrical, reconstruction and selection efficiencies for both signal and backgrounds, to train multivariate classifiers and to determine the shapes of invariant mass distributions of both signal and backgrounds. In the simulation, pp collisions are generated using PYTHIA [15] with a specific LHCb configuration [16]. Decays of hadronic particles are described by EVTGEN [17], in which final-state radiation is generated using PHOTOS [18]. The interaction of the generated particles with the detector, and its response, are simulated using the GEANT4 toolkit [19] as described in Ref. [20].

## 3 Selection

The $B_{(s)}^0 \rightarrow e^{\pm}\mu^{\mp}$ candidates in the events passing the trigger selection are constructed by combining pairs of tracks producing good quality secondary vertices that are separated from any PV in the downstream direction by a flight distance greater than 15 times its uncertainty. Only $B_{(s)}^0$ candidates with $p_T > 0.5$ GeV/c and a small impact parameter $\chi^2$, $\chi_{IP}^2$, are considered, where the $\chi_{IP}^2$ of a $B_{(s)}^0$ candidate is defined as the difference between the $\chi^2$ of the PV reconstructed with and without the considered candidate. The PV with the smallest $\chi_{IP}^2$ is associated to the $B_{(s)}^0$ candidate. The measured momentum of electron candidates is corrected for the loss of momentum due to bremsstrahlung. This correction is made by adding to the electron the momentum of photons consistent with being emitted from the electron before the magnet [21]. Since bremsstrahlung can affect the kinematic distribution of $B_{(s)}^0 \rightarrow e^{\pm}\mu^{\mp}$ candidates, the sample is split into two categories: candidates in which no photon is associated with the electron and candidates for which one or more photons are recovered. The fraction of electrons with recovered bremsstrahlung photons is about 60% for $B_{(s)}^0 \rightarrow e^{\pm}\mu^{\mp}$ decays. Only $B_{(s)}^0 \rightarrow e^{\pm}\mu^{\mp}$ candidates with $m_{e^{\pm}\mu^{\mp}} \in [4900, 5850]$ MeV/c$^2$ are retained to be further analysed.

Particles forming the $B_{(s)}^0 \rightarrow e^{\pm}\mu^{\mp}$ candidates are required to be well identified as an electron and a muon [22], using information from the Cherenkov detectors, the calorimeters and the muon stations. These identification criteria are optimised to keep high signal efficiency while maximising the rejection power for the two-body hadronic B decays, $B \rightarrow h^+h^-$, which are the major peaking backgrounds.
In order to reduce combinatorial background — combinations of two random tracks that can be associated to a common vertex — a loose requirement on the response of a multivariate classifier trained on simulated events is applied to the signal candidates. This classifier takes the following geometrical variables as input: the direction of the $B_{(s)}^0$ meson candidate; its impact parameter with respect to the assigned PV, defined as the PV with which it forms the smallest $\chi_{\text{IP}}^2$; the separation between the two outgoing leptonic tracks at their point of closest approach; and the minimum IP of each lepton particle with respect to any PV. In total 22 020 $B_{(s)}^0 \to e^{\pm}\mu^{\mp}$ candidates are selected, which are mainly comprised of combinatorial background that is made up of true electrons and muons.

The normalisation channels are selected with requirements as similar as possible to those used for the signal. The selection for $B^0 \to K^+\pi^-$ candidates is the same as for the $B_{(s)}^0 \to e^{\pm}\mu^{\mp}$ channel, except for the particle identification criteria which are changed into hadronic particle identification requirements. Similarly, the $B^+ \to J/\psi K^+$ candidate selection is also kept as similar as possible, applying the same selection used for the signal to the dimuon pair from the $J/\psi$, except for the particle identification requirements. Additionally, loose quality requirements are applied on the $B^+$ vertex and particle identification is required on both muons. Finally, a 60 MeV/$c^2$ mass window around the nominal $J/\psi$ mass and the requirement $1.4 < 1 + p_{J/\psi}/p_K < 20.0$ is used. The latter removes backgrounds that have at least one track that is misidentified and another that is not reconstructed, mainly $B \to J/\psi\pi^+X$, where X can be one or more particles.

# 4 BDT training and calibration

A Boosted Decision Tree (BDT) classifier is used to separate the $B_{(s)}^0 \to e^{\pm}\mu^{\mp}$ signal from the combinatorial background. The BDT is trained using a simulated sample of $B_s^0 \to e^{\pm}\mu^{\mp}$ events to describe the signal and a data sample of same-sign $e^{\pm}\mu^{\pm}$ candidates to describe the combinatorial background. The following input variables are used: the proper decay time of the $B_{(s)}^0$ candidate; the minimum $\chi_{\text{IP}}^2$ of the two leptons with respect to the assigned PV; the IP of the $B_{(s)}^0$ candidate with respect to its PV; the distance of closest approach between the two lepton tracks; the degree of isolation of the two tracks with respect to the other tracks in the same event [23]; the transverse momentum of the $B_{(s)}^0$ candidate; the cosine of the angle between the muon momentum in the $B_{(s)}^0$ candidate rest frame and the vector perpendicular to the $B_{(s)}^0$ candidate momentum and the beam axis; the flight distance of the $B_{(s)}^0$ candidate with respect to its PV; the $\chi^2$ of the $B_{(s)}^0$ candidate decay vertex; the maximum transverse momentum of the two decay products and their difference in pseudorapidity.

The BDT response is transformed such that it is uniformly distributed in the range [0,1] for the signal, while peaking at zero for the background. The linear correlation between the BDT response and the dilepton invariant mass is found to be around 4%.

Since the BDT is trained using only kinematic information of a two-body $B_{(s)}^0$ decay, its response is calibrated using $B^0 \to K^+\pi^-$ decays as a proxy. To avoid biases, $B^0 \to K^+\pi^-$ candidates are selected from candidates where the trigger decision did not depend on the presence of the $B^0$ decay products. Furthermore, the candidates are weighted to emulate the effect of the lepton triggers and the particle identification requirements. The number of $B^0 \to K^+\pi^-$ candidates in bins of BDT response is determined by fitting the $K^+\pi^-$ invariant mass distribution. As expected, the BDT response is found to be
<figure><img src="image_2.png" /><figcaption>Figure 1: Expected distribution of the BDT response for $B_{(s)}^0 \to e^\pm\mu^\mp$ decays with recovered bremsstrahlung photons obtained from the $B^0 \to K^+\pi^-$ control channel. The total uncertainty is shown as a light grey band. Each bin is normalised to its width.</figcaption></figure>

consistent with a uniform distribution across the range [0,1]. The distribution of the BDT response is also checked on a $B^0 \to K^+\pi^-$ simulated sample and a uniform distribution is obtained. Candidates with a value smaller than 0.25 are then excluded, as this region is highly contaminated by background, leaving a total of 476 signal candidates. The signal candidates are classified in a binned two-dimensional space formed by the BDT response and the two bremsstrahlung categories. The expected probability density function (PDF) of the BDT response for $B_{(s)}^0 \to e^\pm\mu^\mp$ decays with recovered bremsstrahlung photons is shown in Fig. 1.

Unrecovered bremsstrahlung photons emitted by signal electrons can affect the BDT response and are not accounted for in the calibration procedure since hadrons do not emit significant bremsstrahlung. The impact of bremsstrahlung on the BDT response distribution is evaluated using simulation and a correction is applied where no bremsstrahlung is recovered.

# 5 Normalisation

The $B_{(s)}^0 \to e^\pm\mu^\mp$ yields are obtained from a fit to the lepton-pair invariant mass distribution and translated into branching fractions according to

$$ \begin{aligned} \mathcal{B}(B_{(s)}^0 \to e^\pm\mu^\mp) &= \sum_i w^i \frac{B_{\text{norm}}^i}{N_{\text{norm}}^i} \frac{\varepsilon_{\text{norm}}^i}{\varepsilon_{\text{sig}}} \frac{f_q}{f_{d(s)}} \frac{\mathcal{L}_{\text{norm}}^i}{\mathcal{L}_{\text{sig}}} \times N_{B_{(s)}^0 \to e^\pm\mu^\mp} \\ &= \alpha_{B_{(s)}^0} \times N_{B_{(s)}^0 \to e^\pm\mu^\mp}, \end{aligned} \quad (1) $$
efficiencies are then combined with a weighted average over the properties of the electron and muon tracks of a $B_s^0 \to e^\pm\mu^\mp$ simulated sample.

Particle identification efficiencies are evaluated using calibration samples where the identity of one of the particles can be inferred by means uncorrelated to particle identification requirements. A tag-and-probe method is applied on $J/\psi \to \mu^+\mu^-$ and $J/\psi \to e^+e^-$ decay samples, where only one lepton, the tag, is required to be well identified and the identity of the other lepton is deduced. The single-track efficiencies, calculated as a function of kinematic variables, are then combined and averaged using the momentum distributions of the leptons in a $B_s^0 \to e^\pm\mu^\mp$ simulated sample.

The two normalisation factors $\alpha_{B_s^0}$ and $\alpha_{B^0}$ are determined to be $(2.48 \pm 0.17) \times 10^{-10}$ and $(6.16 \pm 0.23) \times 10^{-11}$. The total efficiencies for the $B^0 \to e^\pm\mu^\mp$, $B_s^0 \to e^\pm\mu^\mp$, $B^+ \to J/\psi K^+$ and $B^0 \to K^+\pi^-$ decays are respectively $(2.22 \pm 0.05)\%$, $(2.29 \pm 0.05)\%$, $(2.215 \pm 0.035)\%$ and $(0.360 \pm 0.021)\%$, where the efficiencies for $B_{(s)}^0 \to e^\pm\mu^\mp$ are for the full BDT and bremsstrahlung category range.

To validate the normalisation procedure, the ratio between the measured branching fractions of $B^0 \to K^+\pi^-$ and $B^+ \to J/\psi K^+$ is determined as

$$R_{\text{norm}} = \frac{N_{B^0 \to K^+ \pi^-} \times \varepsilon_{B^+ \to J/\psi K^+}}{N_{B^+ \to J/\psi K^+} \times \varepsilon_{B^0 \to K^+ \pi^-}} = 0.332 \pm 0.002 (\text{stat}) \pm 0.020 (\text{syst}), \quad (2)$$

where $\varepsilon_{B^+ \to J/\psi K^+}$ and $\varepsilon_{B^0 \to K^+ \pi^-}$ are the selection efficiencies for the $B^0 \to K^+\pi^-$ and $B^+ \to J/\psi K^+$ decays respectively. A correction of about 1% is applied in order to take into account the difference in luminosity between the two channels. The value obtained for $R_{\text{norm}}$ is in excellent agreement with the measured value of $0.321 \pm 0.013$ [27].

# 6 Backgrounds

In addition to the combinatorial background, the signal region is also potentially polluted by backgrounds from exclusive decays where one or more of the final-state particles are misidentified or not reconstructed. The potentially most dangerous of these backgrounds are hadronic $B \to h^+h'^-$ decays where both hadrons are misidentified as an electron-muon pair, resulting in peaking structures near the $B_s^0 \to e^\pm\mu^\mp$ signal mass. Other decays which could contribute, especially at low invariant masses, are $B_c^+ \to J/\psi \ell'^+\nu_\ell$ with $J/\psi \to \ell^+\ell^-$, $B^0 \to \pi^-\ell^+\nu_\ell$, $\Lambda_b^0 \to p\ell^-\bar{\nu}_\ell$ and $B^+ \to \pi^+J/\psi$ with $J/\psi \to \ell^+\ell^-$, where $\ell/\ell'^\pm = e^\pm$ or $\mu^\pm$. These decays do not peak under the signal but are potentially abundant. The expected number of candidates from each possible background decay that pass the signal selection is evaluated using simulation. The candidates are normalised to the number of $B^+ \to J/\psi K^+$ decays found in data as

$$N_X = N_{B^+ \to J/\psi K^+} \frac{f_q}{f_u} \frac{\mathcal{B}(X)}{\mathcal{B}(B^+ \to J/\psi K^+) \cdot \mathcal{B}(J/\psi \to \mu^+\mu^-)} \frac{\varepsilon(X)}{\varepsilon(B^+ \to J/\psi K^+)}, \quad (3)$$

where $N_X$ is the expected number of candidates from the X decay that fall into the $B_s^0 \to e^\pm\mu^\mp$ signal mass window; $f_q$ is the fragmentation fraction; $\mathcal{B}(X)$, $\mathcal{B}(B^+ \to J/\psi K^+)$ and $\mathcal{B}(J/\psi \to \mu^+\mu^-)$ are respectively the branching fractions of the decay under study, $B^+ \to J/\psi K^+$ and $J/\psi \to \mu^+\mu^-$ [27]; $\varepsilon(X)$ is the efficiency for each considered decay to pass the $B_s^0 \to e^\pm\mu^\mp$ selection; and $\varepsilon(B^+ \to J/\psi K^+)$ is the efficiency for $B^+ \to J/\psi K^+$ candidates to pass the respective selection.
The mass and BDT distributions of these background modes are evaluated using simulated samples, while the probabilities of misidentifying kaons, pions and protons as muons or electrons are determined from $D^{*+} \to D^0\pi^+$ with $D^0 \to K^-\pi^+$ and $\Lambda \to p\pi^-$ decays selected from data. The expected total number of $B \to h^+h'^-$ candidates is $0.11 \pm 0.02$ in the full BDT range, which is negligible. This yield estimation is cross-checked using data. A sample of $B \to h^+h'^-$ decays is selected by applying only a partial $B_{(s)}^0 \to e^{\pm}\mu^{\mp}$ selection: only the signal electron PID requirements are applied while the second particle is required to be identified as a pion. The application of these criteria still leaves a sizeable peak to be fit in data. The yield of decays identified as $B_{(s)}^0 \to e^{\pm}\pi^{\mp}$ is then modified to take into account the probability of a pion to be misidentified as a muon. After this correction the expected yield is compatible with the yield obtained using the simulation.

The expected yields of most of the other backgrounds are also found to be negligible. The only backgrounds which are relevant are $B^0 \to \pi^-\mu^+\nu_\ell$ and $\Lambda_b^0 \to p\ell^-\bar{\nu}_\ell$ for which $55 \pm 3$ and $82 \pm 39$ candidates, respectively, are expected in the full BDT range. The contributions from these two decays are included in the fit model.

# 7 Mass calibration

The invariant-mass distribution of $B_{(s)}^0 \to e^{\pm}\mu^{\mp}$ candidates is modelled by a modified Crystal Ball function [28] with two tails on opposite sides defined by two parameters each. The signal shape parameters are obtained from simulation, with data-driven scale factors applied to the core resolution to correct for possible data-simulation discrepancies. For this purpose, since there is no appropriate control channel with an electron and a muon in the final state, $J/\psi \to e^+e^-$ and $J/\psi \to \mu^+\mu^-$ decays are analysed comparing the mass resolution in data and simulation. The results are then combined to reproduce the effect on an $e^{\pm}\mu^{\mp}$ final state. Corrections to the widths of the mass are of the order of 10%. Since bremsstrahlung can significantly alter the mass shape by enhancing the tails, the fit model for $B_{(s)}^0 \to e^{\pm}\mu^{\mp}$ candidates is obtained separately for the two bremsstrahlung categories (see Fig. 3). The mass shape parameters are found to be independent of the particular BDT bin chosen and a single model for each bremsstrahlung category is therefore used.

# 8 Results

The data sample is split into two bremsstrahlung categories, which are further divided into seven subsets each depending on the BDT response covering the range from 0.25 to 1.0, with boundaries 0.25, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0. The region with BDT response lower than 0.25, which is mostly populated by combinatorial background, is excluded from the fit. The $B^0 \to e^{\pm}\mu^{\mp}$ and $B_s^0 \to e^{\pm}\mu^{\mp}$ yields are obtained from a single unbinned extended maximum likelihood fit performed simultaneously to the $m_{e^{\pm}\mu^{\mp}}$ distributions in each subset. The $B_{(s)}^0 \to e^{\pm}\mu^{\mp}$ fractional yields and the mass shape parameters in each category are Gaussian-constrained according to their expected values and uncertainties. The combinatorial background is modelled with an exponential function with independent yield and shape parameters in each subset. The exclusive backgrounds are included as separate components in the fit. Their mass shapes are modelled using nonparametric functions determined from simulation for each bremsstrahlung category. The overall yields
<figure><img src="image_3.png" /><figcaption>Figure 3: Distribution of the $m_{e^{\pm}\mu^{\mp}}$ invariant mass of simulated $B_s^0$ 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.</figcaption></figure>

<table><caption>Table 2: Expected (assuming no signal) and observed upper limits for $\mathcal{B}(B_s^0 \to e^{\pm}\mu^{\mp})$ and $\mathcal{B}(B^0 \to e^{\pm}\mu^{\mp})$ at 95% (90%) CL. The upper limit on the $\mathcal{B}(B_s^0 \to e^{\pm}\mu^{\mp})$ is evaluated under the assumption of pure heavy eigenstate contribution on the decay amplitude.</caption><thead><tr><td>channel</td><td>expected</td><td>observed</td></tr></thead><tbody><tr><td>$\mathcal{B}(B_s^0 \to e^{\pm}\mu^{\mp})$</td><td>$5.0 (3.9) \times 10^{-9}$</td><td>$6.3 (5.4) \times 10^{-9}$</td></tr><tr><td>$\mathcal{B}(B^0 \to e^{\pm}\mu^{\mp})$</td><td>$1.2 (0.9) \times 10^{-9}$</td><td>$1.3 (1.0) \times 10^{-9}$</td></tr></tbody></table>

and fractions of these backgrounds are Gaussian-constrained to their expected values. The result of this fit is shown in Fig. 4.

No significant excess of $B^0 \to e^{\pm}\mu^{\mp}$ or $B_s^0 \to e^{\pm}\mu^{\mp}$ decays is observed and upper limits on the branching fractions are set using the CL$_s$ method [29]. The ratio between the likelihoods in two hypotheses, signal plus background and background only, is used as the test statistic. The likelihoods are computed with nuisance parameters fixed to their nominal values. Pseudoexperiments, in which the nuisance parameters are varied according to their statistical and systematic uncertainties, are used for the evaluation of the test statistic. The resulting CL$_s$ scans are shown in Fig. 5 and upper limits at 95% and 90% confidence level are reported in Table 2.

Several systematic uncertainties can affect the evaluation of the limit on the $B_s^0 \to e^{\pm}\mu^{\mp}$ and $B^0 \to e^{\pm}\mu^{\mp}$ branching fractions through the normalisation formula in Eq. 1 and the fit model used to evaluate the signal yields. The systematic uncertainties are taken into account for the limit computation by constraining the respective nuisance parameters in the likelihood fit with a Gaussian distribution having the central value of the parameter as the mean and its uncertainty as the width. The nuisance parameters for the $B_{(s)}^0 \to e^{\pm}\mu^{\mp}$ yields are related to the calibration of the BDT response, the parameters of the signal shape, the estimated yields of the $B^0 \to \pi^{-}\mu^{+}\nu_{\ell}$ and $\Lambda_b^0 \to p\ell^{-}\bar{\nu}_{\ell}$ backgrounds and the fractional yield per bremsstrahlung category. For the limit on the $B_{(s)}^0 \to e^{\pm}\mu^{\mp}$ branching fractions, the nuisance parameters are in addition related to the signal efficiency, whose uncertainty is dominated by the systematic uncertainty on the trigger efficiencies, and the uncertainties on the efficiencies, branching fractions and yields of the normalisation channels. For the $B_s^0 \to e^{\pm}\mu^{\mp}$ branching fraction estimation, Eq. 1 also includes the hadronisation fraction
<figure><img src="image_4.png" /><figcaption>Figure 4: Distributions of the invariant mass of the $B_{(s)}^0 \to 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 \to e^{\pm}\mu^{\mp}$ component is not visible in the plots.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>Figure 5: Results of the CLs scan used to obtain the limit on (left) $\mathcal{B}(B^0 \to e^\pm\mu^\mp)$ and (right) $\mathcal{B}(B_s^0 \to e^\pm\mu^\mp)$. The background-only 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.</figcaption></figure>

$f_s/f_d$, which dominates the systematic uncertainty for the normalisation. The overall impact on the limits is evaluated to be below 5%.

The two $B_s^0$ mass eigenstates are characterised by a large lifetime difference. Depending on their contribution to the decay amplitude, the selection efficiency and the BDT shape can be affected. Given the negligible difference in lifetime for the $B^0$ system, this effect is not taken into account for the $B^0 \to e^\pm\mu^\mp$ limit evaluation. Two extreme cases can be distinguished: when only the heavy or the light eigenstate contributes to the total decay amplitude. For example, if the only contribution to the LFV $B_s^0 \to e^\pm\mu^\mp$ decay is due to neutrino oscillations, it is expected that the amplitude is dominated by the heavy eigenstate as for the $B_s^0 \to \mu^+\mu^-$ decay [23]. As the contribution to the total amplitude from the heavy and light eigenstate can have an effect on the acceptance, the limit on $\mathcal{B}(B_s^0 \to e^\pm\mu^\mp)$ is evaluated in the two extreme cases. The one reported in Table 2 and obtained from the CLs scan in Fig. 5, is evaluated assuming only a contribution from the heavy eigenstate. For the light eigenstate case the limit is found to be $\mathcal{B}(B_s^0 \to e^\pm\mu^\mp) < 7.2(6.0) \times 10^{-9}$ at 95% (90%) CL.

# 9 Summary

In summary, a search for the LFV decays $B_s^0 \to e^\pm\mu^\mp$ and $B^0 \to e^\pm\mu^\mp$ is performed using pp collision data collected at centre-of-mass energies of 7 and 8 TeV, corresponding to a total integrated luminosity of 3 fb$^{-1}$. No excesses are observed for these two modes and upper limits on the branching fractions are set to $\mathcal{B}(B_s^0 \to e^\pm\mu^\mp) < 6.3(5.4) \times 10^{-9}$ and $\mathcal{B}(B^0 \to e^\pm\mu^\mp) < 1.3(1.0) \times 10^{-9}$ at 95% (90%) CL, where only a contribution from the heavy eigenstate is assumed for the $B_s^0$ meson. If the $B_s^0$ amplitude is completely dominated by the light eigenstate, the upper limit on the branching fraction becomes $\mathcal{B}(B_s^0 \to e^\pm\mu^\mp) < 7.2(6.0) \times 10^{-9}$ at 95% (90%) CL. These results represent the best upper limits to date and are a factor 2 to 3 better than the previous results from LHCb [11].
# Acknowledgements

We express our gratitude to our colleagues in the CERN accelerator departments for the excellent performance of the LHC. We thank the technical and administrative staff at the LHCb institutes. We acknowledge support from CERN and from the national agencies: CAPES, CNPq, FAPERJ and FINEP (Brazil); MOST and NSFC (China); CNRS/IN2P3 (France); BMBF, DFG and MPG (Germany); INFN (Italy); NWO (The Netherlands); MNiSW and NCN (Poland); MEN/IFA (Romania); MinES and FASO (Russia); MinECo (Spain); SNSF and SER (Switzerland); NASU (Ukraine); STFC (United Kingdom); NSF (USA). We acknowledge the computing resources that are provided by CERN, IN2P3 (France), KIT and DESY (Germany), INFN (Italy), SURF (The Netherlands), PIC (Spain), GridPP (United Kingdom), RRCKI and Yandex LLC (Russia), CSCS (Switzerland), IFIN-HH (Romania), CBPF (Brazil), PL-GRID (Poland) and OSC (USA). We are indebted to the communities behind the multiple open-source software packages on which we depend. Individual groups or members have received support from AvH Foundation (Germany), EPLANET, Marie Skłodowska-Curie Actions and ERC (European Union), ANR, Labex P2IO, ENIGMASS and OCEVU, and Région Auvergne-Rhône-Alpes (France), RFBR and Yandex LLC (Russia), GVA, XuntaGal and GENCAT (Spain), Herchel Smith Fund, the Royal Society, the English-Speaking Union and the Leverhulme Trust (United Kingdom).

# References

[1] LHCb collaboration, R. Aaij et al., Measurement of the ratio of branching fractions $B(\bar{B}^0 \to D^{*+}\tau^{-}\bar{\nu}_{\tau})/B(\bar{B}^0 \to D^{*+}\mu^{-}\bar{\nu}_{\mu})$, Phys. Rev. Lett. **115** (2015) 111803, arXiv:1506.08614.

[2] LHCb collaboration, R. Aaij et al., Test of lepton universality with $B^0 \to K^{*0}\ell^+\ell^-$ decays, JHEP **08** (2017) 055, arXiv:1705.05802.

[3] LHCb collaboration, R. Aaij et al., Test of lepton universality using $B^+ \to K^+\ell^+\ell^-$ decays, Phys. Rev. Lett. **113** (2014) 151601, arXiv:1406.6482.

[4] D. Guadagnoli, Flavor anomalies on the eve of the Run-2 verdict, Mod. Phys. Lett. **A32** (2017) 1730006, arXiv:1703.02804.

[5] A. Crivellin et al., Lepton-flavour violating $B$ decays in generic $Z'$ models, Phys. Rev. **D92** (2015) 054013, arXiv:1504.07928.

[6] D. Bečirević, S. Fajfer, N. Košnik, and O. Sumensari, Leptoquark model to explain the $B$-physics anomalies, $R_K$ and $R_D$, Phys. Rev. **D94** (2016) 115021, arXiv:1608.08501.

[7] I. de Medeiros Varzielas and G. Hiller, Clues for flavor from rare lepton and quark decays, JHEP **06** (2015) 072, arXiv:1503.01084.

[8] A. Ilakovac, Lepton flavor violation in the standard model extended by heavy singlet Dirac neutrinos, Phys. Rev. **D62** (2000) 036010, arXiv:hep-ph/9910213.
[9] R. A. Diaz, R. Martinez, and C. E. Sandoval, *Improving bounds on flavor changing vertices in the two Higgs doublet model from $B^0-\bar{B}^0$ mixing*, Eur. Phys. J. C46 (2006) 403, arXiv:hep-ph/0509194.

[10] J. C. Pati and A. Salam, *Lepton Number as the Fourth Color*, Phys. Rev. D10 (1974) 275, Erratum-ibid. D11 (1975) 703.

[11] LHCb collaboration, R. Aaij et al., *Search for the lepton-flavour violating decays $B_s^0 \to e^\pm\mu^\mp$ and $B^0 \to e^\pm\mu^\mp$*, Phys. Rev. Lett. 111 (2013) 141801, arXiv:1307.4889.

[12] LHCb collaboration, A. A. Alves Jr. et al., *The LHCb detector at the LHC*, JINST 3 (2008) S08005.

[13] LHCb collaboration, R. Aaij et al., *LHCb detector performance*, Int. J. Mod. Phys. A30 (2015) 1530022, arXiv:1412.6352.

[14] V. V. Gligorov and M. Williams, *Efficient, reliable and fast high-level triggering using a bonsai boosted decision tree*, JINST 8 (2013) P02013, arXiv:1210.6861.

[15] T. Sjöstrand, S. Mrenna, and P. Skands, *A brief introduction to PYTHIA 8.1*, Comput. Phys. Commun. 178 (2008) 852, arXiv:0710.3820.

[16] I. Belyaev et al., *Handling of the generation of primary events in Gauss, the LHCb simulation framework*, J. Phys. Conf. Ser. 331 (2011) 032047.

[17] D. J. Lange, *The EvtGen particle decay simulation package*, Nucl. Instrum. Meth. A462 (2001) 152.

[18] P. Golonka and Z. Was, *PHOTOS Monte Carlo: A precision tool for QED corrections in Z and W decays*, Eur. Phys. J. C45 (2006) 97, arXiv:hep-ph/0506026.

[19] Geant4 collaboration, J. Allison et al., *Geant4 developments and applications*, IEEE Trans. Nucl. Sci. 53 (2006) 270; Geant4 collaboration, S. Agostinelli et al., *Geant4: A simulation toolkit*, Nucl. Instrum. Meth. A506 (2003) 250.

[20] M. Clemencic et al., *The LHCb simulation application, Gauss: Design, evolution and experience*, J. Phys. Conf. Ser. 331 (2011) 032023.

[21] LHCb collaboration, R. Aaij et al., *Measurement of the $B^0 \to K^{*0}e^+e^-$ branching fraction at low dilepton mass*, JHEP 05 (2013) 159, arXiv:1304.3035.

[22] F. Archilli et al., *Performance of the muon identification at LHCb*, JINST 8 (2013) P10020, arXiv:1306.0249.

[23] LHCb collaboration, R. Aaij et al., *Measurement of the $B_s^0 \to \mu^+\mu^-$ branching fraction and effective lifetime and search for $B^0 \to \mu^+\mu^-$ decays*, Phys. Rev. Lett. 118 (2017) 191801, arXiv:1703.05747.

[24] LHCb collaboration, R. Aaij et al., *Measurement of the fragmentation fraction ratio $f_s/f_d$ and its dependence on B meson kinematics*, JHEP 04 (2013) 001, arXiv:1301.5286, $f_s/f_d$ value updated in LHCb-CONF-2013-011.
[25] S. Tolk, J. Albrecht, F. Dettori, and A. Pellegrino, *Data driven trigger efficiency determination at LHCb*, LHCb-PUB-2014-039.

[26] L. Anderlini et al., *The PIDCalib package*, LHCb-PUB-2016-021.

[27] Particle Data Group, C. Patrignani et al., *Review of particle physics*, Chin. Phys. C40 (2016) 100001, and 2017 update.

[28] T. Skwarnicki, *A study of the radiative cascade transitions between the Upsilon-prime and Upsilon resonances*, PhD thesis, Institute of Nuclear Physics, Krakow, 1986, DESY-F31-86-02.

[29] A. L. Read, *Presentation of search results: The CLs technique*, J. Phys. G28 (2002) 2693.
LHCb collaboration

R. Aaij$^{40}$, B. Adeva$^{39}$, M. Adinolfi$^{48}$, Z. Ajaltouni$^{5}$, S. Akar$^{59}$, J. Albrecht$^{10}$, F. Alessio$^{40}$,
M. Alexander$^{53}$, A. Alfonso Albero$^{38}$, S. Ali$^{43}$, G. Alkhazov$^{31}$, P. Alvarez Cartelle$^{55}$,
A.A. Alves Jr$^{59}$, S. Amato$^{2}$, S. Amerio$^{23}$, Y. Amhis$^{7}$, L. An$^{3}$, L. Anderlini$^{18}$, G. Andreassi$^{41}$,
M. Andreotti$^{17,g}$, J.E. Andrews$^{60}$, R.B. Appleby$^{56}$, F. Archilli$^{43}$, P. d'Argent$^{12}$,
J. Arnau Romeu$^{6}$, A. Artamonov$^{37}$, M. Artuso$^{61}$, E. Aslanides$^{6}$, M. Atzeni$^{42}$, G. Auriemma$^{26}$,
M. Baalouch$^{5}$, I. Babuschkin$^{56}$, S. Bachmann$^{12}$, J.J. Back$^{50}$, A. Badalov$^{38,m}$, C. Baesso$^{62}$,
S. Baker$^{55}$, V. Balagura$^{7,b}$, W. Baldini$^{17}$, A. Baranov$^{35}$, R.J. Barlow$^{56}$, C. Barschel$^{40}$,
S. Barsuk$^{7}$, W. Barter$^{56}$, F. Baryshnikov$^{32}$, V. Batozskaya$^{29}$, V. Battista$^{41}$, A. Bay$^{41}$,
L. Beaucourt$^{4}$, J. Beddow$^{53}$, F. Bedeschi$^{24}$, I. Bediaga$^{1}$, A. Beiter$^{61}$, L.J. Bel$^{43}$, N. Beliy$^{63}$,
V. Bellee$^{41}$, N. Belloli$^{21,i}$, K. Belous$^{37}$, I. Belyaev$^{32,40}$, E. Ben-Haim$^{8}$, G. Bencivenni$^{19}$,
S. Benson$^{43}$, S. Beranek$^{9}$, A. Berezhnoy$^{33}$, R. Bernet$^{42}$, D. Beringhoff$^{12}$, E. Bertholet$^{8}$,
A. Bertolin$^{23}$, C. Betancourt$^{42}$, F. Betti$^{15}$, M.-O. Bettler$^{40}$, M. van Beuzekom$^{43}$, Ia. Bezshyiko$^{42}$,
S. Bifani$^{47}$, P. Billoir$^{8}$, A. Birnkraut$^{10}$, A. Bizzeti$^{18,u}$, M. Bjørn$^{57}$, T. Blake$^{50}$, F. Blanc$^{41}$,
S. Blusk$^{61}$, V. Bocci$^{26}$, T. Boettcher$^{58}$, A. Bondar$^{36,w}$, N. Bondar$^{31}$, I. Bordyuzhin$^{32}$,
A. Borgheresi$^{21,i}$, S. Borghi$^{56}$, M. Borisyak$^{35}$, M. Borsato$^{39}$, F. Bossu$^{7}$, M. Boubdir$^{9}$,
T.J.V. Bowcock$^{54}$, E. Bowen$^{42}$, C. Bozzi$^{17,40}$, S. Braun$^{12}$, T. Britton$^{61}$, J. Brodzicka$^{27}$,
D. Brundu$^{16}$, E. Buchanan$^{48}$, C. Burr$^{56}$, A. Bursche$^{16,f}$, J. Buytaert$^{40}$, W. Byczynski$^{40}$,
S. Cadeddu$^{16}$, H. Cai$^{64}$, R. Calabrese$^{17,g}$, R. Calladine$^{47}$, M. Calvi$^{21,i}$, M. Calvo Gomez$^{38,m}$,
A. Camboni$^{38,m}$, P. Campana$^{19}$, D.H. Campora Perez$^{40}$, L. Capriotti$^{56}$, A. Carbone$^{15,e}$,
G. Carboni$^{25,j}$, R. Cardinale$^{20,h}$, A. Cardini$^{16}$, P. Carniti$^{21,i}$, L. Carson$^{52}$, K. Carvalho Akiba$^{2}$,
G. Casse$^{54}$, L. Cassina$^{21}$, M. Cattaneo$^{40}$, G. Cavallero$^{20,40,h}$, R. Cenci$^{24,t}$, D. Chamont$^{7}$,
M.G. Chapman$^{48}$, M. Charles$^{8}$, Ph. Charpentier$^{40}$, G. Chatzikonstantinidis$^{47}$, M. Chefdeville$^{4}$,
S. Chen$^{16}$, S.F. Cheung$^{57}$, S.-G. Chitic$^{40}$, V. Chobanova$^{39,40}$, M. Chrzaszcz$^{42,27}$, A. Chubykin$^{31}$,
P. Ciambrone$^{19}$, X. Cid Vidal$^{39}$, G. Ciezarek$^{43}$, P.E.L. Clarke$^{52}$, M. Clemencic$^{40}$, H.V. Cliff$^{49}$,
J. Closier$^{40}$, J. Cogan$^{6}$, E. Cogneras$^{5}$, V. Cogoni$^{16,f}$, L. Cojocariu$^{30}$, P. Collins$^{40}$, T. Colombo$^{40}$,
A. Comerma-Montells$^{12}$, A. Contu$^{40}$, A. Cook$^{48}$, G. Coombs$^{40}$, S. Coquereau$^{38}$, G. Corti$^{40}$,
M. Corvo$^{17,g}$, C.M. Costa Sobral$^{50}$, B. Couturier$^{40}$, G.A. Cowan$^{52}$, D.C. Craik$^{58}$,
A. Crocombe$^{50}$, M. Cruz Torres$^{1}$, R. Currie$^{52}$, C. D'Ambrosio$^{40}$, F. Da Cunha Marinho$^{2}$,
E. Dall'Occo$^{43}$, J. Dalseno$^{48}$, A. Davis$^{3}$, O. De Aguiar Francisco$^{40}$, S. De Capua$^{56}$,
M. De Cian$^{12}$, J.M. De Miranda$^{1}$, L. De Paula$^{2}$, M. De Serio$^{14,d}$, P. De Simone$^{19}$, C.T. Dean$^{53}$,
D. Decamp$^{4}$, L. Del Buono$^{8}$, H.-P. Dembinski$^{11}$, M. Demmer$^{10}$, A. Dendek$^{28}$, D. Derkach$^{35}$,
O. Deschamps$^{5}$, F. Dettori$^{54}$, B. Dey$^{65}$, A. Di Canto$^{40}$, P. Di Nezza$^{19}$, H. Dijkstra$^{40}$,
F. Dordei$^{40}$, M. Dorigo$^{40}$, A. Dosil Suárez$^{39}$, L. Douglas$^{53}$, A. Dovbnya$^{45}$, K. Dreimanis$^{54}$,
L. Dufour$^{43}$, G. Dujany$^{8}$, P. Durante$^{40}$, R. Dzhelyadin$^{37}$, M. Dziewiecki$^{12}$, A. Dziurda$^{40}$,
A. Dzyuba$^{31}$, S. Easo$^{51}$, M. Ebert$^{52}$, U. Egede$^{55}$, V. Egorychev$^{32}$, S. Eidelman$^{36,w}$,
S. Eisenhardt$^{52}$, U. Eitschberger$^{10}$, R. Ekelhof$^{10}$, L. Eklund$^{53}$, S. Ely$^{61}$, S. Esen$^{12}$,
H.M. Evans$^{49}$, T. Evans$^{57}$, A. Falabella$^{15}$, N. Farley$^{47}$, S. Farry$^{54}$, D. Fazzini$^{21,i}$, L. Federici$^{25}$,
D. Ferguson$^{52}$, G. Fernandez$^{38}$, P. Fernandez Declara$^{40}$, A. Fernandez Prieto$^{39}$, F. Ferrari$^{15}$,
F. Ferreira Rodrigues$^{2}$, M. Ferro-Luzzi$^{40}$, S. Filippov$^{34}$, R.A. Fini$^{14}$, M. Fiorini$^{17,g}$, M. Firlej$^{28}$,
C. Fitzpatrick$^{41}$, T. Fiutowski$^{28}$, F. Fleuret$^{7,b}$, K. Fohl$^{40}$, M. Fontana$^{16,40}$, F. Fontanelli$^{20,h}$,
D.C. Forshaw$^{61}$, R. Forty$^{40}$, V. Franco Lima$^{54}$, M. Frank$^{40}$, C. Frei$^{40}$, J. Fu$^{22,q}$, W. Funk$^{40}$,
E. Furfaro$^{25,j}$, C. Färber$^{40}$, E. Gabriel$^{52}$, A. Gallas Torreira$^{39}$, D. Galli$^{15,e}$, S. Gallorini$^{23}$,
S. Gambetta$^{52}$, M. Gandelman$^{2}$, P. Gandini$^{22}$, Y. Gao$^{3}$, L.M. Garcia Martin$^{70}$,
J. García Pardiñas$^{39}$, J. Garra Tico$^{49}$, L. Garrido$^{38}$, P.J. Garsed$^{49}$, D. Gascon$^{38}$, C. Gaspar$^{40}$,
L. Gavardi$^{10}$, G. Gazzoni$^{5}$, D. Gerick$^{12}$, E. Gersabeck$^{12}$, M. Gersabeck$^{56}$, T. Gershon$^{50}$,
Ph. Ghez$^{4}$, S. Giani$^{41}$, V. Gibson$^{49}$, O.G. Girard$^{41}$, L. Giubega$^{30}$, K. Gizdov$^{52}$, V.V. Gligorov$^{8}$,
D. Golubkov$^{32}$, A. Golutvin$^{55}$, A. Gomes$^{1,a}$, I.V. Gorelov$^{33}$, C. Gotti$^{21,i}$, E. Govorkova$^{43}$,
J.P. Grabowski$^{12}$, R. Graciani Diaz$^{38}$, L.A. Granado Cardoso$^{40}$, E. Graugés$^{38}$, E. Graverini$^{42}$,
G. Graziani¹⁸, A. Grecu³⁰, R. Greim⁹, P. Griffith¹⁶, L. Grillo²¹, L. Gruber⁴⁰,
B.R. Gruberg Cazon⁵⁷, O. Grünberg⁶⁷, E. Gushchin³⁴, Yu. Guz³⁷, T. Gys⁴⁰, C. Göbel⁶²,
T. Hadavizadeh⁵⁷, C. Hadjivasiliou⁵, G. Haefeli⁴¹, C. Haen⁴⁰, S.C. Haines⁴⁹, B. Hamilton⁶⁰,
X. Han¹², T.H. Hancock⁵⁷, S. Hansmann-Menzemer¹², N. Harnew⁵⁷, S.T. Harnew⁴⁸, C. Hasse⁴⁰,
M. Hatch⁴⁰, J. He⁶³, M. Hecker⁵⁵, K. Heinicke¹⁰, A. Heister⁹, K. Hennessy⁵⁴, P. Henrard⁵,
L. Henry⁷⁰, E. van Herwijnen⁴⁰, M. Heß⁶⁷, A. Hicheur², D. Hill⁵⁷, C. Hombach⁵⁶,
P.H. Hopchev⁴¹, W. Hu⁶⁵, Z.C. Huard⁵⁹, W. Hulsbergen⁴³, T. Humair⁵⁵, M. Hushchyn³⁵,
D. Hutchcroft⁵⁴, P. Ibis¹⁰, M. Idzik²⁸, P. Ilten⁵⁸, R. Jacobsson⁴⁰, J. Jalocha⁵⁷, E. Jans⁴³,
A. Jawahery⁶⁰, F. Jiang³, M. John⁵⁷, D. Johnson⁴⁰, C.R. Jones⁴⁹, C. Joram⁴⁰, B. Jost⁴⁰,
N. Jurik⁵⁷, S. Kandybei⁴⁵, M. Karacson⁴⁰, J.M. Kariuki⁴⁸, S. Karodia⁵³, N. Kazeev³⁵,
M. Kecke¹², F. Keizer⁴⁹, M. Kelsey⁶¹, M. Kenzie⁴⁹, T. Ketel⁴⁴, E. Khairullin³⁵, B. Khanji¹²,
C. Khurewathanakul⁴¹, T. Kirn⁹, S. Klaver⁵⁶, K. Klimaszewski²⁹, T. Klimkovich¹¹, S. Koliiev⁴⁶,
M. Kolpin¹², I. Komarov⁴¹, R. Kopecna¹², P. Koppenburg⁴³, A. Kosmyntseva³²,
S. Kotriakhova³¹, M. Kozeiha⁵, L. Kravchuk³⁴, M. Kreps⁵⁰, F. Kress⁵⁵, P. Krokovny³⁶, w,
F. Kruse¹⁰, W. Krzemien²⁹, W. Kucewicz²⁷, l, M. Kucharczyk²⁷, V. Kudryavtsev³⁶, w,
A.K. Kuonen⁴¹, T. Kvaratskheliya³²,⁴⁰, D. Lacarrere⁴⁰, G. Lafferty⁵⁶, A. Lai¹⁶, G. Lanfranchi¹⁹,
C. Langenbruch⁹, T. Latham⁵⁰, C. Lazzeroni⁴⁷, R. Le Gac⁶, A. Leflat³³,⁴⁰, J. Lefrançois⁷,
R. Lefèvre⁵, F. Lemaitre⁴⁰, E. Lemos Cid³⁹, O. Leroy⁶, T. Lesiak²⁷, B. Leverington¹²,
P.-R. Li⁶³, T. Li³, Y. Li⁷, Z. Li⁶¹, T. Likhomanenko⁶⁸, R. Lindner⁴⁰, F. Lionetto⁴²,
V. Lisovskyi⁷, X. Liu³, D. Loh⁵⁰, A. Loi¹⁶, I. Longstaff⁵³, J.H. Lopes², D. Lucchesi²³, o,
M. Lucio Martinez³⁹, H. Luo⁵², A. Lupato²³, E. Luppi¹⁷, g, O. Lupton⁴⁰, A. Lusiani²⁴, X. Lyu⁶³,
F. Machefert⁷, F. Maciue³⁰, V. Macko⁴¹, P. Mackowiak¹⁰, S. Maddrell-Mander⁴⁸, O. Maev³¹,⁴⁰,
K. Maguire⁵⁶, D. Maisuzenko³¹, M.W. Majewski²⁸, S. Malde⁵⁷, B. Malecki²⁷, A. Malinin⁶⁸,
T. Maltsev³⁶, w, G. Manca¹⁶, f, G. Mancinelli⁶, D. Marangotto²², q, J. Maratas⁵, v,
J.F. Marchand⁴, U. Marconi¹⁵, C. Marin Benito³⁸, M. Marinangeli⁴¹, P. Marino⁴¹, J. Marks¹²,
G. Martellotti²⁶, M. Martin⁶, M. Martinelli⁴¹, D. Martinez Santos³⁹, F. Martinez Vidal⁷⁰,
D. Martins Tostes², L.M. Massacrier⁷, A. Massafferri¹, R. Matev⁴⁰, A. Mathad⁵⁰, Z. Mathe⁴⁰,
C. Matteuzzi²¹, A. Mauri⁴², E. Maurice⁷, b, B. Maurin⁴¹, A. Mazurov⁴⁷, M. McCann⁵⁵,⁴⁰,
A. McNab⁵⁶, R. McNulty¹³, J.V. Mead⁵⁴, B. Meadows⁵⁹, C. Meaux⁶, F. Meier¹⁰, N. Meinert⁶⁷,
D. Melnychuk²⁹, M. Merk⁴³, A. Merli²²,⁴⁰, q, E. Michielin²³, D.A. Milanes⁶⁶, E. Millard⁵⁰,
M.-N. Minard⁴, L. Minzoni¹⁷, D.S. Mitzel¹², A. Mogini⁸, J. Molina Rodriguez¹,
T. Mombächer¹⁰, I.A. Monroy⁶⁶, S. Monteil⁵, M. Morandin²³, M.J. Morello²⁴, t, O. Morgunova⁶⁸,
J. Moron²⁸, A.B. Morris⁵², R. Mountain⁶¹, F. Muheim⁵², M. Mulder⁴³, D. Müller⁵⁶, J. Müller¹⁰,
K. Müller⁴², V. Müller¹⁰, P. Naik⁴⁸, T. Nakada⁴¹, R. Nandakumar⁵¹, A. Nandi⁵⁷, I. Nasteva²,
M. Needham⁵², N. Neri²²,⁴⁰, S. Neubert¹², N. Neufeld⁴⁰, M. Neuner¹², T.D. Nguyen⁴¹,
C. Nguyen-Mau⁴¹, n, S. Nieswand⁹, R. Niet¹⁰, N. Nikitin³³, T. Nikodem¹², A. Nogay⁶⁸,
D.P. O'Hanlon⁵⁰, A. Oblakowska-Mucha²⁸, V. Obraztsov³⁷, S. Ogilvy¹⁹, R. Oldeman¹⁶, f,
C.J.G. Onderwater⁷¹, A. Ossowska²⁷, J.M. Otalora Goicochea², P. Owen⁴², A. Oyanguren⁷⁰,
P.R. Pais⁴¹, A. Palano¹⁴, d, M. Palutan¹⁹,⁴⁰, A. Papanestis⁵¹, M. Pappagallo¹⁴, d,
L.L. Pappalardo¹⁷, g, W. Parker⁶⁰, C. Parkes⁵⁶, G. Passaleva¹⁸,⁴⁰, A. Pastore¹⁴, d, M. Patel⁵⁵,
C. Patrignani¹⁵, e, A. Pearce⁴⁰, A. Pellegrino⁴³, G. Penso²⁶, M. Pepe Altarelli⁴⁰, S. Perazzini⁴⁰,
P. Perret⁵, L. Pescatore⁴¹, K. Petridis⁴⁸, A. Petrolini²⁰, h, A. Petrov⁶⁸, M. Petruzzo²², q,
E. Picatoste Olloqui³⁸, B. Pietrzyk⁴, M. Pikies²⁷, D. Pinci²⁶, F. Pisani⁴⁰, A. Pistone²⁰, h,
A. Piucci¹², V. Placinta³⁰, S. Playfer⁵², M. Plo Casasus³⁹, F. Polci⁸, M. Poli Lener¹⁹,
A. Poluektov⁵⁰, I. Polyakov⁶¹, E. Polycarpo², G.J. Pomery⁴⁸, S. Ponce⁴⁰, A. Popov³⁷,
D. Popov¹¹,⁴⁰, S. Poslavskii³⁷, C. Potterat², E. Price⁴⁸, J. Prisciandaro³⁹, C. Prouve⁴⁸,
V. Pugatch⁴⁶, A. Puig Navarro⁴², H. Pullen⁵⁷, G. Punzi²⁴, p, W. Qian⁵⁰, R. Quagliani⁷,⁴⁸,
B. Quintana⁵, B. Rachwal²⁸, J.H. Rademacker⁴⁸, M. Rama²⁴, M. Ramos Pernas³⁹,
M.S. Rangel², I. Raniuk⁴⁵, †, F. Ratnikov³⁵, G. Raven⁴⁴, M. Ravonel Salzgeber⁴⁰, M. Reboud⁴,
F. Redi⁵⁵, S. Reichert¹⁰, A.C. dos Reis¹, C. Remon Alepuz⁷⁰, V. Renaudin⁷, S. Ricciardi⁵¹,
S. Richards⁴⁸, M. Rih¹⁴⁰, K. Rinnert⁵⁴, V. Rives Molina³⁸, P. Robbe<fnref n="7" />, A. Robert<fnref n="8" />,
A.B. Rodrigues<fnref n="1" />, E. Rodrigues⁵⁹, J.A. Rodriguez Lopez⁶⁶, A. Rogozhnikov³⁵, S. Roiser⁴⁰,
A. Rollings⁵⁷, V. Romanovskiy³⁷, A. Romero Vidal³⁹, J.W. Ronayne<fnref n="13" />, M. Rotondo¹⁹,
M.S. Rudolph⁶¹, T. Ruf⁴⁰, P. Ruiz Valls⁷⁰, J. Ruiz Vidal⁷⁰, J.J. Saborido Silva³⁹,
E. Sadykhov³², N. Sagidova³¹, B. Saitta¹⁶,f, V. Salustino Guimaraes⁶²,
C. Sanchez Mayordomo⁷⁰, B. Sanmartin Sedes³⁹, R. Santacesaria²⁶, C. Santamarina Rios³⁹,
M. Santimaria¹⁹, E. Santovetti²⁵, j, G. Sarpis⁵⁶, A. Sarti¹⁹, k, C. Satriano²⁶, s, A. Satta²⁵,
D.M. Saunders⁴⁸, D. Savrina³²,³³, S. Schael<fnref n="9" />, M. Schellenberg<fnref n="10" />, M. Schiller⁵³, H. Schindler⁴⁰,
M. Schmelling<fnref n="11" />, T. Schmelzer<fnref n="10" />, B. Schmidt⁴⁰, O. Schneider⁴¹, A. Schopper⁴⁰, H.F. Schreiner⁵⁹,
M. Schubiger⁴¹, M.-H. Schune<fnref n="7" />, R. Schwemmer⁴⁰, B. Sciascia¹⁹, A. Sciubba²⁶, k,
A. Semennikov³², E.S. Sepulveda<fnref n="8" />, A. Sergi⁴⁷, N. Serra⁴², J. Serrano<fnref n="6" />, L. Sestini²³, P. Seyfert⁴⁰,
M. Shapkin³⁷, I. Shapoval⁴⁵, Y. Shcheglov³¹, T. Shears⁵⁴, L. Shekhtman³⁶, w, V. Shevchenko⁶⁸,
B.G. Siddi¹⁷, R. Silva Coutinho⁴², L. Silva de Oliveira<fnref n="2" />, G. Simi²³, o, S. Simone<fnref n="14" />, d, M. Sirendi⁴⁹,
N. Skidmore⁴⁸, T. Skwarnicki⁶¹, E. Smith⁵⁵, I.T. Smith⁵², J. Smith⁴⁹, M. Smith⁵⁵,
l. Soares Lavra<fnref n="1" />, M.D. Sokoloff⁵⁹, F.J.P. Soler⁵³, B. Souza De Paula<fnref n="2" />, B. Spaan<fnref n="10" />, P. Spradlin⁵³,
S. Sridharan⁴⁰, F. Stagni⁴⁰, M. Stahl<fnref n="12" />, S. Stahl⁴⁰, P. Stefko⁴¹, S. Stefkova⁵⁵, O. Steinkamp⁴²,
S. Stemmle<fnref n="12" />, O. Stenyakin³⁷, M. Stepanova³¹, H. Stevens<fnref n="10" />, S. Stone⁶¹, B. Storaci⁴²,
S. Stracka²⁴, p, M.E. Stramaglia⁴¹, M. Straticioc³⁰, U. Straumann⁴², J. Sun<fnref n="3" />, L. Sun⁶⁴,
W. Sutcliffe⁵⁵, K. Swientek²⁸, V. Syropoulos⁴⁴, T. Szumlak²⁸, M. Szymanski⁶³, S. T'Jampens<fnref n="4" />,
A. Tayduganov<fnref n="6" />, T. Tekampe<fnref n="10" />, G. Tellarini¹⁷, g, F. Teubert⁴⁰, E. Thomas⁴⁰, J. van Tilburg⁴³,
M.J. Tilley⁵⁵, V. Tisserand<fnref n="4" />, M. Tobin⁴¹, S. Tolk⁴⁹, L. Tomassetti¹⁷, g, D. Tonelli²⁴, F. Toriello⁶¹,
R. Tourinho Jadallah Aoude<fnref n="1" />, E. Tournefier<fnref n="4" />, M. Traill⁵³, M.T. Tran⁴¹, M. Tresch⁴²,
A. Trisovic⁴⁰, A. Tsaregorodtsev<fnref n="6" />, P. Tsopelas⁴³, A. Tully⁴⁹, N. Tuning⁴³,⁴⁰, A. Ukleja²⁹,
A. Usachov<fnref n="7" />, A. Ustyuzhanin⁵⁵, U. Uwer<fnref n="12" />, C. Vacca¹⁶, f, A. Vagner⁶⁹, V. Vagnoni¹⁵,⁴⁰,
A. Valassi⁴⁰, S. Valat⁴⁰, G. Valenti¹⁵, R. Vazquez Gomez⁴⁰, P. Vazquez Regueiro³⁹, S. Vecchi¹⁷,
M. van Veghel⁴³, J.J. Velthuis⁴⁸, M. Veltri¹⁸, r, G. Veneziano⁵⁷, A. Venkateswaran⁶¹,
T.A. Verlage<fnref n="9" />, M. Vernet<fnref n="5" />, M. Vesterinen⁵⁷, J.V. Viana Barbosa⁴⁰, B. Viaud<fnref n="7" />, D. Vieira⁶³,
M. Vieites Diaz³⁹, H. Viemann⁶⁷, X. Vilasis-Cardona³⁸, m, M. Vitti⁴⁹, V. Volkov³³,
A. Vollhardt⁴², B. Voneki⁴⁰, A. Vorobyev³¹, V. Vorobyev³⁶, w, C. Voß<fnref n="9" />, J.A. de Vries⁴³,
C. Vázquez Sierra³⁹, R. Waldi⁶⁷, C. Wallace⁵⁰, R. Wallace<fnref n="13" />, J. Walsh²⁴, J. Wang⁶¹,
D.R. Ward⁴⁹, H.M. Wark⁵⁴, N.K. Watson⁴⁷, D. Websdale⁵⁵, A. Weiden⁴², C. Weisser⁵⁸,
M. Whitehead⁴⁰, J. Wicht⁵⁰, G. Wilkinson⁵⁷, M. Wilkinson⁶¹, M. Williams⁵⁶, M.P. Williams⁴⁷,
M. Williams⁵⁸, T. Williams⁴⁷, F.F. Wilson⁵¹,⁴⁰, J. Wimberley⁶⁰, M. Winn<fnref n="7" />, J. Wishahi<fnref n="10" />,
W. Wislicki²⁹, M. Witek²⁷, G. Wormser<fnref n="7" />, S.A. Wotton⁴⁹, K. W straight⁵³, K. Wyllie⁴⁰, Y. Xie⁶⁵,
M. Xu⁶⁵, Z. Xu<fnref n="4" />, Z. Yang<fnref n="3" />, Z. Yang⁶⁰, Y. Yao⁶¹, H. Yin⁶⁵, J. Yu⁶⁵, X. Yuan⁶¹,
O. Yushchenko³⁷, K.A. Zarebski⁴⁷, M. Zavertyaev<fnref n="11" />, c, L. Zhang<fnref n="3" />, Y. Zhang<fnref n="7" />, A. Zhelezov<fnref n="12" />,
Y. Zheng⁶³, X. Zhu<fnref n="3" />, V. Zhukov³³, J.B. Zonneveld⁵², S. Zucchelli¹⁵.

<footnote n="1">¹Centro Brasileiro de Pesquisas Físicas (CBPF), Rio de Janeiro, Brazil</footnote>

<footnote n="2">²Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil</footnote>

<footnote n="3">³Center for High Energy Physics, Tsinghua University, Beijing, China</footnote>

<footnote n="4">⁴LAPP, Université Savoie Mont-Blanc, CNRS/IN2P3, Annecy-Le-Vieux, France</footnote>

<footnote n="5">⁵Clermont Université, Université Blaise Pascal, CNRS/IN2P3, LPC, Clermont-Ferrand, France</footnote>

<footnote n="6">⁶Aix Marseille Univ, CNRS/IN2P3, CPPM, Marseille, France</footnote>

<footnote n="7">⁷LAL, Université Paris-Sud, CNRS/IN2P3, Orsay, France</footnote>

<footnote n="8">⁸LPNHE, Université Pierre et Marie Curie, Université Paris Diderot, CNRS/IN2P3, Paris, France</footnote>

<footnote n="9">⁹I. Physikalisches Institut, RWTH Aachen University, Aachen, Germany</footnote>

<footnote n="10">¹⁰Fakultät Physik, Technische Universität Dortmund, Dortmund, Germany</footnote>

<footnote n="11">¹¹Max-Planck-Institut für Kernphysik (MPIK), Heidelberg, Germany</footnote>

<footnote n="12">¹²Physikalisches Institut, Ruprecht-Karls-Universität Heidelberg, Heidelberg, Germany</footnote>

<footnote n="13">¹³School of Physics, University College Dublin, Dublin, Ireland</footnote>

<footnote n="14">¹⁴Sezione INFN di Bari, Bari, Italy</footnote>
¹⁵Sezione INFN di Bologna, Bologna, Italy

¹⁶Sezione INFN di Cagliari, Cagliari, Italy

¹⁷Università e INFN, Ferrara, Ferrara, Italy

¹⁸Sezione INFN di Firenze, Firenze, Italy

¹⁹Laboratori Nazionali dell'INFN di Frascati, Frascati, Italy

²⁰Sezione INFN di Genova, Genova, Italy

²¹Università & INFN, Milano-Bicocca, Milano, Italy

²²Sezione di Milano, Milano, Italy

²³Sezione INFN di Padova, Padova, Italy

²⁴Sezione INFN di Pisa, Pisa, Italy

²⁵Sezione INFN di Roma Tor Vergata, Roma, Italy

²⁶Sezione INFN di Roma La Sapienza, Roma, Italy

²⁷Henryk Niewodniczanski Institute of Nuclear Physics Polish Academy of Sciences, Kraków, Poland

²⁸AGH - University of Science and Technology, Faculty of Physics and Applied Computer Science, Kraków, Poland

²⁹National Center for Nuclear Research (NCBJ), Warsaw, Poland

³⁰Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest-Magurele, Romania

³¹Petersburg Nuclear Physics Institute (PNPI), Gatchina, Russia

³²Institute of Theoretical and Experimental Physics (ITEP), Moscow, Russia

³³Institute of Nuclear Physics, Moscow State University (SINP MSU), Moscow, Russia

³⁴Institute for Nuclear Research of the Russian Academy of Sciences (INR RAN), Moscow, Russia

³⁵Yandex School of Data Analysis, Moscow, Russia

³⁶Budker Institute of Nuclear Physics (SB RAS), Novosibirsk, Russia

³⁷Institute for High Energy Physics (IHEP), Protvino, Russia

³⁸ICCUB, Universitat de Barcelona, Barcelona, Spain

³⁹Universidad de Santiago de Compostela, Santiago de Compostela, Spain

⁴⁰European Organization for Nuclear Research (CERN), Geneva, Switzerland

⁴¹Institute of Physics, Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland

⁴²Physik-Institut, Universität Zürich, Zürich, Switzerland

⁴³Nikhef National Institute for Subatomic Physics, Amsterdam, The Netherlands

⁴⁴Nikhef National Institute for Subatomic Physics and VU University Amsterdam, Amsterdam, The Netherlands

⁴⁵NSC Kharkiv Institute of Physics and Technology (NSC KIPT), Kharkiv, Ukraine

⁴⁶Institute for Nuclear Research of the National Academy of Sciences (KINR), Kyiv, Ukraine

⁴⁷University of Birmingham, Birmingham, United Kingdom

⁴⁸H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom

⁴⁹Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom

⁵⁰Department of Physics, University of Warwick, Coventry, United Kingdom

⁵¹STFC Rutherford Appleton Laboratory, Didcot, United Kingdom

⁵²School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom

⁵³School of Physics and Astronomy, University of Glasgow, Glasgow, United Kingdom

⁵⁴Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom

⁵⁵Imperial College London, London, United Kingdom

⁵⁶School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom

⁵⁷Department of Physics, University of Oxford, Oxford, United Kingdom

⁵⁸Massachusetts Institute of Technology, Cambridge, MA, United States

⁵⁹University of Cincinnati, Cincinnati, OH, United States

⁶⁰University of Maryland, College Park, MD, United States

⁶¹Syracuse University, Syracuse, NY, United States

⁶²Pontifícia Universidade Católica do Rio de Janeiro (PUC-Rio), Rio de Janeiro, Brazil, associated to ²

⁶³University of Chinese Academy of Sciences, Beijing, China, associated to ³

⁶⁴School of Physics and Technology, Wuhan University, Wuhan, China, associated to ³

⁶⁵Institute of Particle Physics, Central China Normal University, Wuhan, Hubei, China, associated to ³

⁶⁶Departamento de Fisica, Universidad Nacional de Colombia, Bogota, Colombia, associated to ⁸

⁶⁷Institut für Physik, Universität Rostock, Rostock, Germany, associated to ¹²

⁶⁸National Research Centre Kurchatov Institute, Moscow, Russia, associated to ³²
<sup>69</sup> *National Research Tomsk Polytechnic University, Tomsk, Russia, associated to* <sup>32</sup>

<sup>70</sup> *Instituto de Fisica Corpuscular, Centro Mixto Universidad de Valencia - CSIC, Valencia, Spain, associated to* <sup>38</sup>

<sup>71</sup> *Van Swinderen Institute, University of Groningen, Groningen, The Netherlands, associated to* <sup>43</sup>

<sup>a</sup> *Universidade Federal do Triângulo Mineiro (UFTM), Uberaba-MG, Brazil*

<sup>b</sup> *Laboratoire Leprince-Ringuet, Palaiseau, France*

<sup>c</sup> *P.N. Lebedev Physical Institute, Russian Academy of Science (LPI RAS), Moscow, Russia*

<sup>d</sup> *Università di Bari, Bari, Italy*

<sup>e</sup> *Università di Bologna, Bologna, Italy*

<sup>f</sup> *Università di Cagliari, Cagliari, Italy*

<sup>g</sup> *Università di Ferrara, Ferrara, Italy*

<sup>h</sup> *Università di Genova, Genova, Italy*

<sup>i</sup> *Università di Milano Bicocca, Milano, Italy*

<sup>j</sup> *Università di Roma Tor Vergata, Roma, Italy*

<sup>k</sup> *Università di Roma La Sapienza, Roma, Italy*

<sup>l</sup> *AGH - University of Science and Technology, Faculty of Computer Science, Electronics and Telecommunications, Kraków, Poland*

<sup>m</sup> *LIFAELS, La Salle, Universitat Ramon Llull, Barcelona, Spain*

<sup>n</sup> *Hanoi University of Science, Hanoi, Viet Nam*

<sup>o</sup> *Università di Padova, Padova, Italy*

<sup>p</sup> *Università di Pisa, Pisa, Italy*

<sup>q</sup> *Università degli Studi di Milano, Milano, Italy*

<sup>r</sup> *Università di Urbino, Urbino, Italy*

<sup>s</sup> *Università della Basilicata, Potenza, Italy*

<sup>t</sup> *Scuola Normale Superiore, Pisa, Italy*

<sup>u</sup> *Università di Modena e Reggio Emilia, Modena, Italy*

<sup>v</sup> *Iligan Institute of Technology (IIT), Iligan, Philippines*

<sup>W</sup> *Novosibirsk State University, Novosibirsk, Russia*

<sup>†</sup> *Deceased*