Search for the Rare Decays B 0 s → e + e − and B 0 → e + e −
Onderwater, C. J. G.; LHCb Collaboration
Published in:
Physical Review Letters DOI:
10.1103/PhysRevLett.124.211802
IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please check the document version below.
Document Version
Publisher's PDF, also known as Version of record
Publication date: 2020
Link to publication in University of Groningen/UMCG research database
Citation for published version (APA):
Onderwater, C. J. G., & LHCb Collaboration (2020). Search for the Rare Decays B 0 s → e + e − and B 0 → e + e −. Physical Review Letters, 124(21), [211802]. https://doi.org/10.1103/PhysRevLett.124.211802
Copyright
Other than for strictly personal use, it is not permitted to download or to forward/distribute the text or part of it without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license (like Creative Commons).
Take-down policy
If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim.
Downloaded from the University of Groningen/UMCG research database (Pure): http://www.rug.nl/research/portal. For technical reasons the number of authors shown on this cover page is limited to 10 maximum.
Search for the Rare Decays B
0s→ e
+e
−and B
0→ e
+e
− R. Aaijet al.*(LHCb Collaboration)
(Received 10 March 2020; accepted 22 April 2020; published 27 May 2020)
A search for the decays B0s → eþe− and B0→ eþe− is performed using data collected with the
LHCb experiment in proton-proton collisions at center-of-mass energies of 7, 8, and 13 TeV, corresponding to integrated luminosities of 1, 2, and2 fb−1, respectively. No signal is observed. Assuming no contribution from B0→ eþe− decays, an upper limit on the branching fractionBðB0s→ eþe−Þ < 9.4ð11.2Þ × 10−9 is obtained at 90(95)% confidence level. If no B0s → eþe− contribution is assumed, a limit of
BðB0→ eþe−Þ < 2.5ð3.0Þ × 10−9is determined at 90(95)% confidence level. These upper limits are more
than one order of magnitude lower than the previous values.
DOI:10.1103/PhysRevLett.124.211802
Searches for rare particle decays provide ideal probes for contributions from physics processes beyond the standard model (SM). Recent measurements of decays involving b → slþl− transitions (the inclusion of charge-conjugated processes is implied throughout this Letter) hint at deviations from SM predictions in lepton-flavor universality tests[1–6] and thus motivate measurements of decay rates into final states involving leptons. Following the observation of the decay B0s → μþμ− [7,8], the search for B0s → eþe− and B0→ eþe− decays provides an independent test of lepton-flavor universality. According to SM predictions (calculated from Ref.[9], neglecting QED corrections that are expected to be at the percent level), B0ðsÞ→eþe−decays have branch-ing fractions of BðB0s→eþe−Þ¼ð8.600.36Þ×10−14 and BðB0→eþe−Þ¼ð2.410.13Þ×10−15. With contributions beyond the SM, these branching fractions could be significantly larger, reaching values of Oð10−8Þ for BðB0
s→eþe−Þ and Oð10−10Þ for BðB0→eþe−Þ[10]. These values are close to the current experimental bounds of BðB0
s→eþe−Þ<2.8×10−7 and BðB0→eþe−Þ<8.3×10−8 at 90% confidence level (CL) [11], set by the CDF collaboration.
In this Letter, a search for B0s → eþe− and B0→ eþe− decays is presented using data collected with the LHCb experiment in proton-proton collisions at center-of-mass energies of 7 TeV in 2011, 8 TeV in 2012 and 13 TeV in 2015 and 2016, corresponding to integrated luminosities of 1, 2, and 2 fb−1, respectively. The signal yields are
determined from a fit to the data and normalized to those of the Bþ → J=ψKþ decay, where the J=ψ meson decays to eþe−, which has a precisely measured branching fraction [12]and a similar dielectron signature in the detector.
The LHCb detector [13,14] is a single-arm forward spectrometer covering the pseudorapidity range2 < η < 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. Different types of charged hadrons are distinguished using information from two ring-imaging Cherenkov detectors. Photons, electrons, and hadrons are identified by a calo-rimeter system consisting of scintillating-pad and pre-shower detectors, an electromagnetic 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[15], which consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, which applies a full event reconstruction. At the hardware trigger stage, events are required to have a high-energy deposit in the calorimeters associated with a signal electron candidate, or a muon candidate with high transverse momentum pT, or a photon, electron, or hadron candidate with high transverse energy from the decays of other particles from the pp collision. The software trigger requires a two-track secondary vertex with a significant displacement from any primary pp interaction vertex (PV). At least one charged particle must have high pT and be inconsistent with originating from a PV. A multivariate algorithm[16,17]is used in the trigger for the identification of secondary vertices consistent with the decay of a b *Full author list given at the end of the article.
Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI. Funded by SCOAP3.
hadron. Simulated samples are used to optimize the candidate selection, estimate selection efficiencies and describe the expected invariant-mass shapes of the signal candidates and background decays. In the simulation, pp collisions are generated usingPYTHIA[18]with a specific LHCb configuration[19]. Decays of unstable particles are described byEVTGEN[20], in which final-state radiation is generated using PHOTOS [21]. The interaction of the generated particles with the detector, and its response, are implemented using theGEANT4toolkit[22]as described in Ref.[23]. The simulation is corrected for data-simulation differences in B-meson production kinematics, detector occupancy, and isolation criteria[24] using Bþ → J=ψKþ and B0s→J=ψϕ decays, with J=ψ →eþe− andϕ → KþK−. Particle identification variables are calibrated using data from Bþ → J=ψKþ and D0→ K−πþ decays [25]. The calibration data are binned in momentum and pseudor-apidity of the particle as well as detector occupancy to account for possible differences in kinematics between the investigated decay and the calibration data.
The B0ðsÞ→ eþe− candidates are selected in events passing the trigger requirements by combining two tracks that are inconsistent with originating from any PV in the event and which form a good-quality secondary vertex. The tracks are also required to have a momentum larger than3 GeV=c and pT greater than500 MeV=c, and must be identified as electrons using information from the Cherenkov detectors and calorimeters. The dielectron candidate’s momentum must be aligned with the vector pointing from a PV (the associated PV) to the two-track vertex and have a considerable transverse component. The candidate must also have an invariant mass in the range½4166; 6566 MeV=c2.
The measured electron momenta are corrected for losses due to bremsstrahlung radiation by adding the momentum of photons consistent with being emitted upstream of the magnet [26]. Candidates in data and simulation are sepa-rated into three categories with either zero, one, or both electrons having a bremsstrahlung correction applied. To avoid experimenters’ bias, the narrowest dielectron invari-ant-mass region containing 90% of simulated B0s → eþe− decays, corresponding to a range of½4689; 5588 MeV=c2, was removed from the data set until the analysis procedure was finalized.
Candidates for the normalization mode, Bþ → J=ψKþ, are constructed similarly, but require an additional track consistent with being a kaon and originating from the same vertex as the dielectron candidate. The dielectron candidate must have an invariant mass in the range ½2450; 3176 MeV=c2, consistent with arising from a J=ψ meson decay. In addition, the reconstructed Bþ-candidate mass, when the dielectron candidate is constrained to the known J=ψ mass [12], must be above 5175 MeV=c2, suppressing partially reconstructed decays.
A boosted decision tree (BDT) algorithm[27–29]is used to separate B0ðsÞ → eþe−signal from random combinations of two electrons (combinatorial background). The BDT is trained separately for data taking periods 2011–2012 (Run 1) and 2015–2016 (Run 2) on simulated B0s→eþe− decays as signal proxy and dielectron candidates from data with a mass above5588 MeV=c2as background proxy. The split between the data taking periods is done to account for changes in the center-of-mass energies and trigger strategies, which significantly impact the data distributions and improve the BDT and the particle identification algorithms in Run 2. It is checked that the data behave consistently across the data-taking periods. The BDT input variables comprise of the following: kinematic information on the electron tracks and B candidate, information on the dis-placement of the electrons and B candidate from the associated PV, and isolation variables that quantify the compatibility of other tracks in the event with originating from the same decay as the B candidate[24,30]. Candidates with a BDT response compatible with that of the background are discarded, with the threshold chosen by maximizing the figure of merit ϵsignal=ð
ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
Nbackground
p
þ3=2Þ [31], where
ϵsignal is the signal efficiency and the expected background
yield in the signal region is Nbackground.
The final selected data set is separated by data-taking period and by category of bremsstrahlung correction. The branching fractionBðB0ðsÞ→ eþe−Þ is measured relative to that of the normalization channel via
BðB0 ðsÞ→ eþe−Þ ¼ NðB0ðsÞ→ eþe−Þ × α× BðBþ→ J=ψKþÞ × fdðsÞ fu −1 ; ð1Þ where α ≡εðBþ → J=ψKþÞ εðB0 ðsÞ→ eþe−Þ × 1 NðBþ→ J=ψKþÞ; ð2Þ εðB0
ðsÞ→eþe−Þ and εðBþ→J=ψKþÞ denote the efficiencies of the signal and normalization modes, and NðB0ðsÞ→eþe−Þ and NðBþ → J=ψKþÞ their yields. The normalization mode branching fraction (including that for the decay J=ψ → eþe−) is BðBþ→ J=ψKþÞ ¼ ð6.03 0.17Þ × 10−5, taken from Ref. [12]. The b-hadron fragmentation fraction ratio fd=fu is assumed to be unity, while fs=fu¼ 0.259 0.015 [32] is used for the Run 1 data and is scaled by 1.068 0.016 for the Run 2 data, according to Ref. [33], to account for center-of-mass energy differences. A measurement of fs=fu from Run 2 yields a consistent, but less precise, result[34].
The yield of the normalization mode is determined using an unbinned maximum-likelihood fit to the Kþ eþe− invariant mass separately for each year of data
taking and bremsstrahlung category. The fit model com-prises a Gaussian function with power-law tails[35]for the signal component, where the tail parameters are fixed from simulation, and an exponential function to describe combinatorial background. Summed over the bremsstrah-lung categories, the yield of the normalization mode is 20480 140 in the Run 1 data and 33080 180 in the Run 2 data.
The selection efficiencies εðB0ðsÞ → eþe−Þ and εðBþ → J=ψKþÞ are determined separately for each year of data taking and bremsstrahlung category using simulated decays that are weighted to better represent the data. Calibration data are used to evaluate particle-identification efficiencies [25]. Trigger efficiencies are also estimated from data, using the technique described in Ref.[36]. For simulated B0s → eþe− decays, the mean B0s lifetime[37]is assumed. The selection efficiency is assumed to be the same for both B0→ eþe−and B0s → eþe−decays, which is consistent with results from simulation. The normalization factors,α, are combined across the data-taking periods and, given in TableI, split by bremsstrahlung category (for the selection efficiency ratio between normalization and signal mode, see the Supplemental Material [38]).
In addition to the combinatorial background, back-grounds due to misidentification and partial reconstruction are present in the data. These backgrounds differ signifi-cantly between the categories of bremsstrahlung correction. Their invariant-mass shapes and relative contributions are evaluated using simulation. In the lower mass region, partially reconstructed backgrounds of the types B → Xeþe− and Bþ → D0ð→ Yþe−¯νeÞeþνe dominate, where X and Y represent hadronic systems. The main source of background in the B -mass region, however, stems from misidentified particles in the decays B0→ π−eþνeand B → hþh0−, where h and h0are hadrons. The latter has a peaking structure in the B -mass region. Backgrounds involving misidentified particles contribute mostly to categories in which at most one of the electrons has a bremsstrahlung correction applied. The contribution from combinatorial background is evaluated from same-sign lepton pairs in data and found to be small. The yields of the backgrounds are Gaussian constrained to their expected values, estimated from simulation using their known branching fractions [12].
The shape of the invariant mass of the B0s → eþe− and B0→ eþe− components is modeled using a Gaussian function with power-law tails, where the parameters are obtained from simulation and differ between each brems-strahlung category and year of data taking. The peak values and the widths of the functions are corrected for data-simulation differences by a factor determined from the normalization mode. The parameters of the B0s → eþe−and B0→ eþe− line shapes are fixed to the same values with the exception of the peak value, which is shifted according to the known B0s–B0 mass difference [12]. Due to the limited mass resolution, arising from imperfect bremsstrah-lung recovery, the line shapes from B0s→ eþe− and B0→ eþe− are highly overlapping. Therefore the branch-ing fraction of B0s → eþe− is obtained by performing a simultaneous fit to the dielectron invariant-mass distribu-tion of all six data sets while neglecting the contribudistribu-tion from B0→ eþe−, and vice versa. In these fits, the only shared parameters between categories are the branching fractions BðB0ðsÞ → eþe−Þ and BðBþ→ J=ψKþÞ, and the ratio of the fragmentation fractions fs=fu.
Systematic uncertainties are estimated separately for each data set. Dominant sources of systematic uncertainties in the normalization arise from the uncertainty on the fragmentation fraction ratio, the technique used to evaluate the trigger efficiencies, and the determination of particle-identification efficiencies; the systematic uncertainties from these sources extend to 5.8%, 5.3%, and 5.3% on the branching fractions, respectively. The uncertainty on BðBþ → J=ψKþÞ of 2.8% [12] is taken into account. A difference of up to 4.1% is found between the efficiency of the BDT selection on simulated Bþ→ J=ψKþ decays and Bþ → J=ψKþ decays in data, which is assigned as a systematic uncertainty. The fraction of candidates in each bremsstrahlung-correction category of the signal modes is taken from simulation. The difference between simulation and data is investigated using Bþ → J=ψKþ decays and its effect on the normalization, up to 4.0%, is taken as a systematic uncertainty. Systematic uncertainties on the invariant-mass resolution corrections are determined by repeating the correction procedure with pseudoexperiments obtained with the bootstrapping method[39], yielding up to 1.1%. A difference between the total selection efficiencies in the B0s → eþe− and B0→ eþe− channels of up to 2.5% is assigned as a systematic uncertainty on the B0→ eþe− normalization factor. Due to the presence of an additional kaon in the final state of the normalization mode, the track-reconstruction efficiency is different between the signal and normalization modes. An uncertainty of 1.1% is assigned to the branching fraction as a systematic uncertainty on the kaon reconstruction efficiency arising from the limited knowledge of the interactions in the detector material[40]. Finally, an uncertainty of 1.0% is assigned to account for small differences in detector occupancy between the signal
TABLE I. Normalization factors α for B0ðsÞ→ eþe−. The bremsstrahlung category denotes whether zero, one or both electrons are corrected for bremsstrahlung losses. The uncertain-ties include statistical uncertainuncertain-ties and uncertainuncertain-ties due to limited size of the simulated samples.
Bremsstrahlung category 2011–2012 ½10−5 2015–2016 ½10−5 No correction 2.85 0.24 1.84 0.08 One electron corrected 1.13 0.08 0.70 0.03 Both electrons corrected 1.73 0.20 1.04 0.06
and normalization mode arising from the trigger selection. The dominant sources of systematic uncertainties on the background composition are due to the imprecise knowl-edge of the branching fractions of the background compo-nents. The largest uncertainty of this type on the expected background yield in the B-mass region is 14%, determined from refitting the mass sidebands while varying the background components according to their uncertainties. Taking all correlations into account, overall single event sensitivities of½4.71 0.12ðstatÞ 0.33ðsystÞ × 10−10for B0s→eþe− and ½1.2710.034ðstatÞ0.063ðsystÞ×10−10 for B0→ eþe− are obtained.
The dielectron invariant-mass spectrum, summed over bremsstrahlung categories, is shown in Fig.1, with the result of the B0s → eþe−fit. The individual categories are shown in the Supplemental Material[38], as well as the distributions with the result of the B0→ eþe−fit. The measured branching fractions are BðB0s → eþe−Þ ¼ ð2.4 4.4Þ × 10−9 and BðB0→ eþe−Þ ¼ ð0.30 1.29Þ × 10−9, where the uncer-tainties include both statistical and systematic components. The results are in agreement with the background-only hypothesis.
Upper limits on the branching fractions are set using the CLs method [41], as implemented in the GAMMACOMBO framework[42,43]with a one-sided profile likelihood ratio [44]as test statistic. The likelihoods are computed from fits to the invariant-mass distributions. In the fits, the normali-zation factor, normalinormali-zation mode branching fraction, frag-mentation fraction ratio, and background yields are Gaussian constrained to their expected values within stat-istical and systematic uncertainties. Pseudoexperiments, in which the nuisance parameters are set to their fitted values from data, are used for the evaluation of the test statistic.
The expected and observed CLs
distribu-tions are shown in Fig. 2. The upper observed limits are BðB0s → eþe−Þ < 9.4ð11.2Þ × 10−9 and BðB0→ eþe−Þ < 2.5ð3.0Þ × 10−9 at 90(95)% confi-dence level. These are consistent with the expected upper limits of BðB0s → eþe−Þ < 7.0ð8.6Þ × 10−9 and BðB0→ eþe−Þ < 2.0ð2.5Þ × 10−9 at 90(95)% confidence level, obtained as the median of limits determined on background-only pseudoexperiments.
In conclusion, a search for the rare decays B0ðsÞ→ eþe− is performed using data from proton-proton collisions 4500 5000 5500 6000 6500 ] 2 c ) [MeV/ − e + e ( m 0 50 100 150 200 ) 2c Candidates / ( 120 MeV/ LHCb 2011−2012 data − e + e → s 0 B full model combinatorial decays e ν + e 0 D → + B decays − e + e X → B decays e ν − e + h → b X decays − h' + h → B 4500 5000 5500 6000 6500 ] 2 c ) [MeV/ − e + e ( m 0 200 400 ) 2 c Candidates / ( 120 MeV/ LHCb 2015−2016 data − e + e → s 0 B full model combinatorial decays e ν + e 0 D → + B decays − e + e X → B decays e ν − e + h → b X decays − h' + h → B
FIG. 1. Simultaneous fit to the dielectron invariant-mass distribution, withBðB0→ eþe−Þ fixed to zero. The sum of bremsstrahlung categories is shown for (left) Run 1 and (right) Run 2. The relative proportions of background contributions change between Run 1 and Run 2 due to different performances of the particle identification algorithms and BDT selections.
branching fraction − e + e → s 0 B 0 5 10 15 20 9 − 10 × S CL 0 0.2 0.4 0.6 0.8 1 Observed Expected σ 1 ± σ 2 ± 90.0% 95.0% LHCb branching fraction − e + e → 0 B 0 1 2 3 4 5 6 9 − 10 × S CL 0 0.2 0.4 0.6 0.8 1 Observed Expected σ 1 ± σ 2 ± 90.0% 95.0% LHCb
FIG. 2. CLsvalues as a function of the branching fractions of the decays (left) B0s→ eþe−and (right) B0→ eþe−. The red solid line
(black solid line with data points) corresponds to the distribution of the expected (observed) upper limits, and the light blue (dark blue) band contains the1σ ð2σÞ uncertainties on the expected upper limits. Thresholds corresponding to 90% and 95% confidence level are indicated with dashed lines. The observed values are plotted for branching fractions greater than the measured branching fraction in the data; the test statistic is defined to be nonzero only in that region.
recorded with the LHCb experiment, corresponding to a total integrated luminosity of5 fb−1. No excess of events is observed over the background. The resulting limits on the branching fractions areBðB0s → eþe−Þ < 9.4ð11.2Þ × 10−9 and BðB0→ eþe−Þ < 2.5ð3.0Þ × 10−9 at 90(95)% confi-dence level, when neglecting the contribution from the other decay. The mean B0s lifetime is assumed for B0s → eþe−decays. Assuming SM-like CP-odd (CP-even) B0s → eþe− decays, an increase (decrease) of 2.4% with respect to the quoted limit is found. The results improve the limits on these branching fractions [11] by more than one order of magnitude and constrain contributions beyond the SM, for example from scalar and pseudoscalar currents [10].
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 (Netherlands); MNiSW and NCN (Poland); MEN/IFA (Romania); MSHE (Russia); MinECo (Spain); SNSF and SER (Switzerland); NASU (Ukraine); STFC (United Kingdom); DOE NP and NSF (USA). We acknowledge the computing resources that are provided by CERN, IN2P3 (France), KIT and DESY (Germany), INFN (Italy), SURF (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 and OCEVU, and R´egion Auvergne-Rhône-Alpes (France); Key Research Program of Frontier Sciences of CAS, CAS PIFI, and the Thousand Talents Program (China); RFBR, RSF, and Yandex LLC (Russia); GVA, XuntaGal, and GENCAT (Spain); the Royal Society and the Leverhulme Trust (United Kingdom).
[1] R. Aaij et al. (LHCb Collaboration), Test of lepton univer-sality with B0→ K0lþl−decays,J. High Energy Phys. 08 (2017) 055.
[2] R. Aaij et al. (LHCb Collaboration), Search for Lepton-Universality Violation in Bþ→ Kþlþl− Decays, Phys. Rev. Lett. 122, 191801 (2019).
[3] R. Aaij et al. (LHCb Collaboration), Test of lepton univer-sality usingΛ0b→ pK−lþl−decays,J. High Energy Phys. 05 (2020) 040.
[4] J. P. Lees et al. (BABAR Collaboration), Measurement of branching fractions and rate asymmetries in the rare decays B → KðÞlþl−,Phys. Rev. D 86, 032012 (2012).
[5] A. Abdesselam et al. (Belle Collaboration), Test of lepton flavor universality in B → Klþl− decays, arXiv: 1908.01848.
[6] A. Abdesselam et al. (Belle Collaboration), Test of lepton flavor universality in B → Klþl− decays at Belle, arXiv:1904.02440.
[7] V. Khachatryan et al. (CMS and LHCb Collaborations), Observation of the rare B0s → μþμ− decay from the
com-bined analysis of CMS and LHCb data, Nature (London) 522, 68 (2015).
[8] R. Aaij et al. (LHCb Collaboration), Measurement of the B0s → μþμ−Branching Fraction and Effective Lifetime and
Search for B0→ μþμ− Decays, Phys. Rev. Lett. 118, 191801 (2017).
[9] M. Beneke, C. Bobeth, and R. Szafron, Power-enhanced leading-logarithmic QED corrections to Bq→ μþμ−,J. High
Energy Phys. 10 (2019) 232.
[10] R. Fleischer, R. Jaarsma, and G. Tetlalmatzi-Xolocotzi, In pursuit of New Physics with B0s;d→ lþl−,J. High Energy
Phys. 05 (2017) 156.
[11] T. Aaltonen et al. (CDF Collaboration), Search for the Decays B0s→ eþμ−and B0s → eþe−in CDF Run II,Phys.
Rev. Lett. 102, 201801 (2009).
[12] M. Tanabashi et al. (Particle Data Group), Review of particle physics, Phys. Rev. D 98, 030001 (2018), and 2019 update.
[13] A. A. Alves, Jr. et al. (LHCb Collaboration), The LHCb detector at the LHC, J. Instrum. 3, S08005 (2008). [14] R. Aaij et al. (LHCb Collaboration), LHCb detector
performance,Int. J. Mod. Phys. A 30, 1530022 (2015). [15] R. Aaij et al., The LHCb trigger and its performance in
2011,J. Instrum. 8, P04022 (2013).
[16] V. V. Gligorov and M. Williams, Efficient, reliable and fast high-level triggering using a bonsai boosted decision tree, J. Instrum. 8, P02013 (2013).
[17] T. Likhomanenko, P. Ilten, E. Khairullin, A. Rogozhnikov, A. Ustyuzhanin, and M. Williams, LHCb topological trigger reoptimization, J. Phys. Conf. Ser. 664, 082025 (2015).
[18] T. Sjöstrand, S. Mrenna, and P. Skands, PYTHIA 6.4 physics and manual, J. High Energy Phys. 05 (2006) 026; A brief introduction to PYTHIA 8.1,Comput. Phys. Commun. 178, 852 (2008).
[19] I. Belyaev et al., Handling of the generation of primary events in Gauss, the LHCb simulation framework,J. Phys. Conf. Ser. 331, 032047 (2011).
[20] D. J. Lange, The EvtGen particle decay simulation package, Nucl. Instrum. Methods Phys. Res., Sect. A 462, 152 (2001).
[21] P. Golonka and Z. Was, PHOTOS Monte Carlo: A precision tool for QED corrections in Z and W decays,Eur. Phys. J. C 45, 97 (2006).
[22] J. Allison et al. (GEANT4Collaboration), Geant4 develop-ments and applications, IEEE Trans. Nucl. Sci. 53, 270 (2006); S. Agostinelli et al. (Geant4 Collaboration), Geant4: A simulation toolkit,Nucl. Instrum. Methods Phys. Res., Sect. A 506, 250 (2003).
[23] M. Clemencic, G. Corti, S. Easo, C. R. Jones, S. Miglior-anzi, M. Pappagallo, and P. Robbe, The LHCb simulation application, Gauss: Design, evolution and experience, J. Phys. Conf. Ser. 331, 032023 (2011).
[24] R. Aaij et al. (LHCb Collaboration), Search for the rare decays B0s → μþμ−and B0→ μþμ−,Phys. Lett. B 708, 55
(2012).
[25] L. Anderlini et al., The PIDCalib package, Report No. LHCb-PUB-2016-021, 2016.
[26] R. Aaij et al. (LHCb Collaboration), Measurement of the B0→ K0eþe− branching fraction at low dilepton mass, J. High Energy Phys. 05 (2013) 159.
[27] Y. Freund and R. E. Schapire, A decision-theoretic gener-alization of on-line learning and an application to boosting, J. Comput. Syst. Sci. 55, 119 (1997).
[28] L. Breiman, J. H. Friedman, R. A. Olshen, and C. J. Stone, Classification and Regression Trees (Wadsworth International Group, Belmont, California, USA, 1984). [29] F. Pedregosa et al., Scikit-learn: Machine learning in
Python, J. Mach. Learn. Res. 12, 2825 (2011), and online athttp://scikit-learn.org/stable/.
[30] L. Gavardi, Search for lepton flavour violation in τ decays at the LHCb experiment, Master’s thesis, Universit degli studi di Milano-Bicocca, 2013, CERN-THESIS-2013-259.
[31] G. Punzi, Sensitivity of searches for new signals and its optimization, eConf C030908, MODT002 (2003). [32] R. Aaij et al. (LHCb Collaboration), Measurement
of the fragmentation fraction ratio fs=fdand its dependence
on B meson kinematics, J. High Energy Phys. 04 (2013) 001, fs=fd value updated in
LHCb-CONF-2013-011.
[33] R. Aaij et al. (LHCb Collaboration), Measurement of fs=fu Variation with Proton-Proton Collision Energy
and B-Meson Kinematics, Phys. Rev. Lett. 124, 122002 (2020).
[34] R. Aaij et al. (LHCb Collaboration), Measurement of b-hadron fractions in 13 TeV pp collisions,Phys. Rev. D 100, 031102(R) (2019).
[35] T. Skwarnicki, A study of the radiative cascade transitions between the Upsilon-prime and Upsilon resonances, Ph.D. thesis, Institute of Nuclear Physics, Krakow, 1986, DESY-F31-86-02.
[36] S. Tolk, J. Albrecht, F. Dettori, and A. Pellegrino, Data driven trigger efficiency determination at LHCb, Report No. LHCb-PUB-2014-039, 2014.
[37] K. A. Olive et al. (Particle Data Group), Review of particle physics,Chin. Phys. C 38, 090001 (2014).
[38] See the Supplemental Material at http://link.aps.org/ supplemental/10.1103/PhysRevLett.124.211802for the se-lection efficiency ratio of normalization and signal mode, the individual categories of the B0s → eþe− fit and the
distributions with the result of the B0→ eþe−fit. [39] B. Efron and R. J. Tibshirani, An Introduction to the
Bootstrap, Mono. Stat. Appl. Probab. (Chapman and Hall, London, 1993).
[40] R. Aaij et al. (LHCb Collaboration), Measurement of the track reconstruction efficiency at LHCb, J. Instrum. 10, P02007 (2015).
[41] A. L. Read, Presentation of search results: The CLs
tech-nique,J. Phys. G 28, 2693 (2002).
[42] R. Aaij et al. (LHCb Collaboration), Measurement of the CKM angleγ from a combination of LHCb results,J. High Energy Phys. 12 (2016) 087.
[43] M. Kenzie et al., GammaCombo: A statistical analysis framework for combining measurements, fitting datasets and producing confidence intervals,https://doi.org/10.5281/ zenodo.3371421.
[44] G. Cowan, K. Cranmer, E. Gross, and O. Vitells, Asymp-totic formulae for likelihood-based tests of new physics, Eur. Phys. J. C 71, 1554 (2011); Erratum,Eur. Phys. J. C 73, 2501 (2013).
R. Aaij,31 C. Abellán Beteta,49T. Ackernley,59 B. Adeva,45 M. Adinolfi,53 H. Afsharnia,9 C. A. Aidala,81S. Aiola,25 Z. Ajaltouni,9S. Akar,66P. Albicocco,22J. Albrecht,14F. Alessio,47M. Alexander,58A. Alfonso Albero,44G. Alkhazov,37 P. Alvarez Cartelle,60A. A. Alves Jr.,45S. Amato,2Y. Amhis,11L. An,21L. Anderlini,21G. Andreassi,48M. Andreotti,20 F. Archilli,16A. Artamonov,43M. Artuso,67K. Arzymatov,41E. Aslanides,10M. Atzeni,49B. Audurier,11S. Bachmann,16
J. J. Back,55S. Baker,60V. Balagura,11,a W. Baldini,20A. Baranov,41R. J. Barlow,61 S. Barsuk,11W. Barter,60 M. Bartolini,23,47,bF. Baryshnikov,78J. M. Basels,13G. Bassi,28 V. Batozskaya,35B. Batsukh,67 A. Battig,14A. Bay,48
M. Becker,14F. Bedeschi,28I. Bediaga,1A. Beiter,67 V. Belavin,41S. Belin,26 V. Bellee,48K. Belous,43I. Belyaev,38 G. Bencivenni,22E. Ben-Haim,12S. Benson,31A. Berezhnoy,39R. Bernet,49D. Berninghoff,16H. C. Bernstein,67 C. Bertella,47E. Bertholet,12A. Bertolin,27C. Betancourt,49F. Betti,19,c M. O. Bettler,54Ia. Bezshyiko,49 S. Bhasin,53
J. Bhom,33M. S. Bieker,14S. Bifani,52P. Billoir,12A. Bizzeti,21,dM. Bjørn,62M. P. Blago,47T. Blake,55F. Blanc,48 S. Blusk,67D. Bobulska,58V. Bocci,30O. Boente Garcia,45T. Boettcher,63A. Boldyrev,79A. Bondar,42,e N. Bondar,37,47
S. Borghi,61M. Borisyak,41 M. Borsato,16 J. T. Borsuk,33T. J. V. Bowcock,59C. Bozzi,20M. J. Bradley,60S. Braun,65 A. Brea Rodriguez,45M. Brodski,47J. Brodzicka,33A. Brossa Gonzalo,55D. Brundu,26E. Buchanan,53 A. Büchler-Germann,49A. Buonaura,49C. Burr,47A. Bursche,26A. Butkevich,40J. S. Butter,31 J. Buytaert,47 W. Byczynski,47S. Cadeddu,26H. Cai,72R. Calabrese,20,f L. Calero Diaz,22S. Cali,22R. Calladine,52M. Calvi,24,g
A. F. Campoverde Quezada,5 L. Capriotti,19,c A. Carbone,19,c G. Carboni,29R. Cardinale,23,b A. Cardini,26I. Carli,6 P. Carniti,24,gK. Carvalho Akiba,31A. Casais Vidal,45G. Casse,59M. Cattaneo,47G. Cavallero,47S. Celani,48R. Cenci,28,i
J. Cerasoli,10M. G. Chapman,53M. Charles,12Ph. Charpentier,47G. Chatzikonstantinidis,52M. Chefdeville,8 V. Chekalina,41C. Chen,3 S. Chen,26A. Chernov,33 S.-G. Chitic,47V. Chobanova,45S. Cholak,48M. Chrzaszcz,33 A. Chubykin,37V. Chulikov,37P. Ciambrone,22M. F. Cicala,55X. Cid Vidal,45G. Ciezarek,47F. Cindolo,19P. E. L. Clarke,57
M. Clemencic,47H. V. Cliff,54J. Closier,47J. L. Cobbledick,61V. Coco,47J. A. B. Coelho,11J. Cogan,10E. Cogneras,9 L. Cojocariu,36P. Collins,47T. Colombo,47 A. Contu,26N. Cooke,52 G. Coombs,58 S. Coquereau,44G. Corti,47 C. M. Costa Sobral,55B. Couturier,47D. C. Craik,63J. Crkovská,66A. Crocombe,55M. Cruz Torres,1,jR. Currie,57
C. L. Da Silva,66 E. Dall’Occo,14J. Dalseno,45,53C. D’Ambrosio,47A. Danilina,38P. d’Argent,47A. Davis,61 O. De Aguiar Francisco,47K. De Bruyn,47S. De Capua,61M. De Cian,48J. M. De Miranda,1L. De Paula,2M. De Serio,18,k P. De Simone,22J. A. de Vries,76C. T. Dean,66W. Dean,81D. Decamp,8L. Del Buono,12B. Delaney,54H.-P. Dembinski,14
A. Dendek,34V. Denysenko,49D. Derkach,79O. Deschamps,9 F. Desse,11F. Dettori,26,lB. Dey,7 A. Di Canto,47 P. Di Nezza,22S. Didenko,78H. Dijkstra,47 V. Dobishuk,51F. Dordei,26 M. Dorigo,28,mA. C. dos Reis,1 L. Douglas,58 A. Dovbnya,50K. Dreimanis,59M. W. Dudek,33L. Dufour,47P. Durante,47J. M. Durham,66D. Dutta,61M. Dziewiecki,16 A. Dziurda,33A. Dzyuba,37S. Easo,56U. Egede,69V. Egorychev,38S. Eidelman,42,eS. Eisenhardt,57S. Ek-In,48L. Eklund,58 S. Ely,67A. Ene,36E. Epple,66S. Escher,13J. Eschle,49S. Esen,31T. Evans,47A. Falabella,19J. Fan,3Y. Fan,5N. Farley,52
S. Farry,59D. Fazzini,11P. Fedin,38M. F´eo,47P. Fernandez Declara,47A. Fernandez Prieto,45F. Ferrari,19,c L. Ferreira Lopes,48F. Ferreira Rodrigues,2S. Ferreres Sole,31M. Ferrillo,49M. Ferro-Luzzi,47S. Filippov,40R. A. Fini,18 M. Fiorini,20,fM. Firlej,34K. M. Fischer,62C. Fitzpatrick,47T. Fiutowski,34F. Fleuret,11,aM. Fontana,47F. Fontanelli,23,b R. Forty,47V. Franco Lima,59M. Franco Sevilla,65M. Frank,47C. Frei,47D. A. Friday,58J. Fu,25,nQ. Fuehring,14W. Funk,47 E. Gabriel,57A. Gallas Torreira,45D. Galli,19,cS. Gallorini,27S. Gambetta,57Y. Gan,3M. Gandelman,2P. Gandini,25Y. Gao,4
L. M. Garcia Martin,46J. García Pardiñas,49B. Garcia Plana,45F. A. Garcia Rosales,11 L. Garrido,44D. Gascon,44 C. Gaspar,47D. Gerick,16E. Gersabeck,61M. Gersabeck,61T. Gershon,55D. Gerstel,10Ph. Ghez,8 V. Gibson,54 A. Gioventù,45P. Gironella Gironell,44L. Giubega,36C. Giugliano,20K. Gizdov,57V. V. Gligorov,12C. Göbel,70 E. Golobardes,44,h D. Golubkov,38A. Golutvin,60,78 A. Gomes,1,oP. Gorbounov,38I. V. Gorelov,39C. Gotti,24,g E. Govorkova,31 J. P. Grabowski,16R. Graciani Diaz,44 T. Grammatico,12L. A. Granado Cardoso,47E. Graug´es,44 E. Graverini,48G. Graziani,21A. Grecu,36R. Greim,31P. Griffith,20L. Grillo,61L. Gruber,47B. R. Gruberg Cazon,62C. Gu,3 E. Gushchin,40A. Guth,13Yu. Guz,43,47T. Gys,47P. A. Günther,16T. Hadavizadeh,62G. Haefeli,48C. Haen,47S. C. Haines,54 P. M. Hamilton,65Q. Han,7X. Han,16T. H. Hancock,62S. Hansmann-Menzemer,16N. Harnew,62T. Harrison,59R. Hart,31 C. Hasse,14M. Hatch,47J. He,5M. Hecker,60K. Heijhoff,31K. Heinicke,14A. M. Hennequin,47K. Hennessy,59L. Henry,46 J. Heuel,13A. Hicheur,68D. Hill,62M. Hilton,61P. H. Hopchev,48J. Hu,16J. Hu,71W. Hu,7W. Huang,5W. Hulsbergen,31 T. Humair,60R. J. Hunter,55M. Hushchyn,79D. Hutchcroft,59D. Hynds,31P. Ibis,14M. Idzik,34P. Ilten,52A. Inglessi,37 K. Ivshin,37R. Jacobsson,47S. Jakobsen,47 E. Jans,31B. K. Jashal,46A. Jawahery,65V. Jevtic,14F. Jiang,3 M. John,62 D. Johnson,47C. R. Jones,54B. Jost,47N. Jurik,62S. Kandybei,50M. Karacson,47J. M. Kariuki,53N. Kazeev,79M. Kecke,16 F. Keizer,54,47M. Kelsey,67M. Kenzie,55T. Ketel,32B. Khanji,47A. Kharisova,80K. E. Kim,67T. Kirn,13V. S. Kirsebom,48 S. Klaver,22K. Klimaszewski,35S. Koliiev,51A. Kondybayeva,78A. Konoplyannikov,38P. Kopciewicz,34R. Kopecna,16 P. Koppenburg,31M. Korolev,39I. Kostiuk,31,51O. Kot,51S. Kotriakhova,37L. Kravchuk,40R. D. Krawczyk,47M. Kreps,55
F. Kress,60S. Kretzschmar,13P. Krokovny,42,e W. Krupa,34W. Krzemien,35W. Kucewicz,33,p M. Kucharczyk,33 V. Kudryavtsev,42,eH. S. Kuindersma,31G. J. Kunde,66T. Kvaratskheliya,38D. Lacarrere,47G. Lafferty,61A. Lai,26 D. Lancierini,49J. J. Lane,61G. Lanfranchi,22C. Langenbruch,13O. Lantwin,49T. Latham,55F. Lazzari,28,qR. Le Gac,10 S. H. Lee,81R. Lef`evre,9A. Leflat,39,47O. Leroy,10T. Lesiak,33B. Leverington,16H. Li,71L. Li,62X. Li,66Y. Li,6Z. Li,67 X. Liang,67T. Lin,60R. Lindner,47V. Lisovskyi,14G. Liu,71X. Liu,3D. Loh,55A. Loi,26J. Lomba Castro,45I. Longstaff,58 J. H. Lopes,2G. Loustau,49G. H. Lovell,54Y. Lu,6D. Lucchesi,27,rM. Lucio Martinez,31Y. Luo,3A. Lupato,27E. Luppi,20,f
O. Lupton,55A. Lusiani,28,sX. Lyu,5S. Maccolini,19,c F. Machefert,11 F. Maciuc,36V. Macko,48 P. Mackowiak,14 S. Maddrell-Mander,53L. R. Madhan Mohan,53 O. Maev,37 A. Maevskiy,79D. Maisuzenko,37M. W. Majewski,34 S. Malde,62B. Malecki,47A. Malinin,77T. Maltsev,42,eH. Malygina,16G. Manca,26,lG. Mancinelli,10R. Manera Escalero,44
D. Manuzzi,19,c D. Marangotto,25,n J. Maratas,9,tJ. F. Marchand,8U. Marconi,19S. Mariani,21,21,47 C. Marin Benito,11 M. Marinangeli,48P. Marino,48J. Marks,16P. J. Marshall,59 G. Martellotti,30 L. Martinazzoli,47M. Martinelli,24,g
D. Martinez Santos,45F. Martinez Vidal,46A. Massafferri,1 M. Materok,13R. Matev,47A. Mathad,49Z. Mathe,47 V. Matiunin,38C. Matteuzzi,24K. R. Mattioli,81A. Mauri,49E. Maurice,11,aM. McCann,60L. Mcconnell,17A. McNab,61
R. McNulty,17J. V. Mead,59B. Meadows,64 C. Meaux,10G. Meier,14N. Meinert,74D. Melnychuk,35S. Meloni,24,g M. Merk,31A. Merli,25 M. Mikhasenko,47D. A. Milanes,73 E. Millard,55M.-N. Minard,8 O. Mineev,38L. Minzoni,20 S. E. Mitchell,57B. Mitreska,61D. S. Mitzel,47A. Mödden,14A. Mogini,12R. D. Moise,60T. Mombächer,14I. A. Monroy,73 S. Monteil,9 M. Morandin,27G. Morello,22M. J. Morello,28,sJ. Moron,34A. B. Morris,10A. G. Morris,55R. Mountain,67 H. Mu,3F. Muheim,57M. Mukherjee,7M. Mulder,47D. Müller,47K. Müller,49C. H. Murphy,62D. Murray,61P. Muzzetto,26 P. Naik,53T. Nakada,48R. Nandakumar,56T. Nanut,48I. Nasteva,2M. Needham,57N. Neri,25,nS. Neubert,16N. Neufeld,47 R. Newcombe,60T. D. Nguyen,48C. Nguyen-Mau,48,uE. M. Niel,11S. Nieswand,13N. Nikitin,39N. S. Nolte,47C. Nunez,81
A. Oblakowska-Mucha,34V. Obraztsov,43S. Ogilvy,58D. P. O’Hanlon,53 R. Oldeman,26,lC. J. G. Onderwater,75 J. D. Osborn,81A. Ossowska,33J. M. Otalora Goicochea,2 T. Ovsiannikova,38P. Owen,49 A. Oyanguren,46P. R. Pais,48
T. Pajero,28,28,47,sA. Palano,18 M. Palutan,22G. Panshin,80 A. Papanestis,56M. Pappagallo,57L. L. Pappalardo,20 C. Pappenheimer,64W. Parker,65C. Parkes,61G. Passaleva,21,47A. Pastore,18M. Patel,60 C. Patrignani,19,c A. Pearce,47 A. Pellegrino,31M. Pepe Altarelli,47S. Perazzini,19D. Pereima,38P. Perret,9L. Pescatore,48K. Petridis,53A. Petrolini,23,b
A. Petrov,77S. Petrucci,57M. Petruzzo,25,n B. Pietrzyk,8 G. Pietrzyk,48M. Pili,62D. Pinci,30J. Pinzino,47F. Pisani,19 A. Piucci,16V. Placinta,36 S. Playfer,57J. Plews,52M. Plo Casasus,45F. Polci,12M. Poli Lener,22M. Poliakova,67 A. Poluektov,10N. Polukhina,78,vI. Polyakov,67 E. Polycarpo,2G. J. Pomery,53S. Ponce,47A. Popov,43D. Popov,52 S. Poslavskii,43K. Prasanth,33L. Promberger,47 C. Prouve,45V. Pugatch,51A. Puig Navarro,49H. Pullen,62G. Punzi,28,i W. Qian,5J. Qin,5 R. Quagliani,12B. Quintana,8N. V. Raab,17R. I. Rabadan Trejo,10B. Rachwal,34J. H. Rademacker,53
M. Rama,28M. Ramos Pernas,45M. S. Rangel,2 F. Ratnikov,41,79G. Raven,32M. Reboud,8 F. Redi,48F. Reiss,12 C. Remon Alepuz,46Z. Ren,3 V. Renaudin,62S. Ricciardi,56D. S. Richards,56S. Richards,53K. Rinnert,59P. Robbe,11 A. Robert,12A. B. Rodrigues,48E. Rodrigues,59J. A. Rodriguez Lopez,73M. Roehrken,47A. Rollings,62V. Romanovskiy,43
M. Romero Lamas,45A. Romero Vidal,45J. D. Roth,81M. Rotondo,22M. S. Rudolph,67T. Ruf,47J. Ruiz Vidal,46 A. Ryzhikov,79J. Ryzka,34J. J. Saborido Silva,45N. Sagidova,37N. Sahoo,55B. Saitta,26,lC. Sanchez Gras,31 C. Sanchez Mayordomo,46R. Santacesaria,30C. Santamarina Rios,45 M. Santimaria,22E. Santovetti,29,w G. Sarpis,61 M. Sarpis,16A. Sarti,30C. Satriano,30,x A. Satta,29M. Saur,5D. Savrina,38,39 L. G. Scantlebury Smead,62S. Schael,13 M. Schellenberg,14M. Schiller,58H. Schindler,47M. Schmelling,15T. Schmelzer,14B. Schmidt,47O. Schneider,48
A. Schopper,47H. F. Schreiner,64M. Schubiger,31S. Schulte,48M. H. Schune,11R. Schwemmer,47B. Sciascia,22 A. Sciubba,22S. Sellam,68A. Semennikov,38A. Sergi,52,47N. Serra,49J. Serrano,10L. Sestini,27A. Seuthe,14P. Seyfert,47
D. M. Shangase,81 M. Shapkin,43L. Shchutska,48T. Shears,59L. Shekhtman,42,e V. Shevchenko,77E. Shmanin,78 J. D. Shupperd,67 B. G. Siddi,20R. Silva Coutinho,49L. Silva de Oliveira,2 G. Simi,27,r S. Simone,18,k I. Skiba,20 N. Skidmore,16T. Skwarnicki,67M. W. Slater,52 J. G. Smeaton,54A. Smetkina,38E. Smith,13I. T. Smith,57M. Smith,60
A. Snoch,31 M. Soares,19L. Soares Lavra,9 M. D. Sokoloff,64F. J. P. Soler,58B. Souza De Paula,2 B. Spaan,14 E. Spadaro Norella,25,nP. Spradlin,58F. Stagni,47M. Stahl,64S. Stahl,47P. Stefko,48O. Steinkamp,49,78 S. Stemmle,16 O. Stenyakin,43M. Stepanova,37H. Stevens,14S. Stone,67S. Stracka,28M. E. Stramaglia,48M. Straticiuc,36S. Strokov,80 J. Sun,26L. Sun,72Y. Sun,65P. Svihra,61K. Swientek,34A. Szabelski,35T. Szumlak,34M. Szymanski,47S. Taneja,61Z. Tang,3 T. Tekampe,14F. Teubert,47E. Thomas,47K. A. Thomson,59M. J. Tilley,60V. Tisserand,9 S. T’Jampens,8 M. Tobin,6 S. Tolk,47 L. Tomassetti,20,f D. Torres Machado,1 D. Y. Tou,12E. Tournefier,8 M. Traill,58 M. T. Tran,48E. Trifonova,78
C. Trippl,48A. Tsaregorodtsev,10 G. Tuci,28,iA. Tully,48N. Tuning,31A. Ukleja,35A. Usachov,31 A. Ustyuzhanin,41,79 U. Uwer,16A. Vagner,80V. Vagnoni,19A. Valassi,47G. Valenti,19M. van Beuzekom,31H. Van Hecke,66E. van Herwijnen,47
C. B. Van Hulse,17M. van Veghel,75 R. Vazquez Gomez,44P. Vazquez Regueiro,45C. Vázquez Sierra,31S. Vecchi,20 J. J. Velthuis,53M. Veltri,21,y A. Venkateswaran,67M. Veronesi,31M. Vesterinen,55J. V. Viana Barbosa,47D. Vieira,64 M. Vieites Diaz,48 H. Viemann,74X. Vilasis-Cardona,44,h G. Vitali,28A. Vitkovskiy,31A. Vollhardt,49 D. Vom Bruch,12 A. Vorobyev,37V. Vorobyev,42,eN. Voropaev,37R. Waldi,74J. Walsh,28J. Wang,3J. Wang,72J. Wang,6M. Wang,3Y. Wang,7 Z. Wang,49D. R. Ward,54H. M. Wark,59N. K. Watson,52D. Websdale,60A. Weiden,49C. Weisser,63B. D. C. Westhenry,53
D. J. White,61M. Whitehead,13 D. Wiedner,14G. Wilkinson,62M. Wilkinson,67I. Williams,54M. Williams,63 M. R. J. Williams,61T. Williams,52F. F. Wilson,56W. Wislicki,35M. Witek,33L. Witola,16G. Wormser,11S. A. Wotton,54 H. Wu,67K. Wyllie,47Z. Xiang,5D. Xiao,7Y. Xie,7H. Xing,71A. Xu,4J. Xu,5L. Xu,3M. Xu,7Q. Xu,5Z. Xu,4Z. Yang,3
Z. Yang,65Y. Yao,67L. E. Yeomans,59H. Yin,7J. Yu,7X. Yuan,67O. Yushchenko,43K. A. Zarebski,52M. Zavertyaev,15,v M. Zdybal,33M. Zeng,3D. Zhang,7L. Zhang,3S. Zhang,4W. C. Zhang,3,zY. Zhang,47A. Zhelezov,16Y. Zheng,5X. Zhou,5
Y. Zhou,5 X. Zhu,3 V. Zhukov,13,39J. B. Zonneveld,57and S. Zucchelli19,c (LHCb Collaboration)
1Centro Brasileiro de Pesquisas Físicas (CBPF), Rio de Janeiro, Brazil 2
Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil
3Center for High Energy Physics, Tsinghua University, Beijing, China 4
School of Physics State Key Laboratory of Nuclear Physics and Technology, Peking University, Beijing, China
5University of Chinese Academy of Sciences, Beijing, China 6
Institute Of High Energy Physics (IHEP), Beijing, China
7Institute of Particle Physics, Central China Normal University, Wuhan, Hubei, China 8
Universit´e Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, IN2P3-LAPP, Annecy, France
9Universit´e Clermont Auvergne, CNRS/IN2P3, LPC, Clermont-Ferrand, France 10
Aix Marseille Universit´e, CNRS/IN2P3, CPPM, Marseille, France
11Universit´e Paris-Saclay, CNRS/IN2P3, IJCLab, Orsay, France 12
LPNHE, Sorbonne Universit´e, Paris Diderot Sorbonne Paris Cit´e, CNRS/IN2P3, Paris, France
13I. Physikalisches Institut, RWTH Aachen University, Aachen, Germany 14
Fakultät Physik, Technische Universität Dortmund, Dortmund, Germany
15
Max-Planck-Institut für Kernphysik (MPIK), Heidelberg, Germany
16
Physikalisches Institut, Ruprecht-Karls-Universität Heidelberg, Heidelberg, Germany
17
School of Physics, University College Dublin, Dublin, Ireland
18
INFN Sezione di Bari, Bari, Italy
19
INFN Sezione di Bologna, Bologna, Italy
20
INFN Sezione di Ferrara, Ferrara, Italy
21
INFN Sezione di Firenze, Firenze, Italy
22
INFN Laboratori Nazionali di Frascati, Frascati, Italy
23
INFN Sezione di Genova, Genova, Italy
24
INFN Sezione di Milano-Bicocca, Milano, Italy
25
INFN Sezione di Milano, Milano, Italy
26
INFN Sezione di Cagliari, Monserrato, Italy
27
INFN Sezione di Padova, Padova, Italy
28
INFN Sezione di Pisa, Pisa, Italy
29
INFN Sezione di Roma Tor Vergata, Roma, Italy
30
INFN Sezione di Roma La Sapienza, Roma, Italy
31
Nikhef National Institute for Subatomic Physics, Amsterdam, Netherlands
32
Nikhef National Institute for Subatomic Physics and VU University Amsterdam, Amsterdam, Netherlands
33
Henryk Niewodniczanski Institute of Nuclear Physics Polish Academy of Sciences, Kraków, Poland
34
AGH—University of Science and Technology, Faculty of Physics and Applied Computer Science, Kraków, Poland
35National Center for Nuclear Research (NCBJ), Warsaw, Poland 36
Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest-Magurele, Romania
37Petersburg Nuclear Physics Institute NRC Kurchatov Institute (PNPI NRC KI), Gatchina, Russia 38
Institute of Theoretical and Experimental Physics NRC Kurchatov Institute (ITEP NRC KI), Moscow, Russia, Moscow, Russia
39Institute of Nuclear Physics, Moscow State University (SINP MSU), Moscow, Russia 40
Institute for Nuclear Research of the Russian Academy of Sciences (INR RAS), Moscow, Russia
41Yandex School of Data Analysis, Moscow, Russia 42
Budker Institute of Nuclear Physics (SB RAS), Novosibirsk, Russia
43Institute for High Energy Physics NRC Kurchatov Institute (IHEP NRC KI), Protvino, Russia, Protvino, Russia 44
ICCUB, Universitat de Barcelona, Barcelona, Spain
45Instituto Galego de Física de Altas Enerxías (IGFAE), Universidade de Santiago de Compostela, Santiago de Compostela, Spain 46
Instituto de Fisica Corpuscular, Centro Mixto Universidad de Valencia - CSIC, Valencia, Spain
47European Organization for Nuclear Research (CERN), Geneva, Switzerland 48
Institute of Physics, Ecole Polytechnique F´ed´erale de Lausanne (EPFL), Lausanne, Switzerland
49
Physik-Institut, Universität Zürich, Zürich, Switzerland
50
NSC Kharkiv Institute of Physics and Technology (NSC KIPT), Kharkiv, Ukraine
51
Institute for Nuclear Research of the National Academy of Sciences (KINR), Kyiv, Ukraine
52
53H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom 54
Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom
55Department of Physics, University of Warwick, Coventry, United Kingdom 56
STFC Rutherford Appleton Laboratory, Didcot, United Kingdom
57School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom 58
School of Physics and Astronomy, University of Glasgow, Glasgow, United Kingdom
59Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom 60
Imperial College London, London, United Kingdom
61Department of Physics and Astronomy, University of Manchester, Manchester, United Kingdom 62
Department of Physics, University of Oxford, Oxford, United Kingdom
63Massachusetts Institute of Technology, Cambridge, Massachusetts, USA 64
University of Cincinnati, Cincinnati, Ohio, USA
65University of Maryland, College Park, Maryland, USA 66
Los Alamos National Laboratory (LANL), Los Alamos, New Mexico, USA
67Syracuse University, Syracuse, New York, USA 68
Laboratory of Mathematical and Subatomic Physics, Constantine, Algeria [associated with Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil]
69
School of Physics and Astronomy, Monash University, Melbourne, Australia (associated with Department of Physics, University of Warwick, Coventry, United Kingdom)
70
Pontifícia Universidade Católica do Rio de Janeiro (PUC-Rio), Rio de Janeiro, Brazil [associated with Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil]
71
Guangdong Provencial Key Laboratory of Nuclear Science, Institute of Quantum Matter, South China Normal University, Guangzhou, China (associated with Center for High Energy Physics, Tsinghua University, Beijing, China)
72
School of Physics and Technology, Wuhan University, Wuhan, China (associated with Center for High Energy Physics, Tsinghua University, Beijing, China)
73
Departamento de Fisica, Universidad Nacional de Colombia, Bogota, Colombia (associated with LPNHE, Sorbonne Universit´e, Paris Diderot Sorbonne Paris Cit´e, CNRS/IN2P3, Paris, France)
74
Institut für Physik, Universität Rostock, Rostock, Germany (associated with Physikalisches Institut, Ruprecht-Karls-Universität Heidelberg, Heidelberg, Germany)
75
Van Swinderen Institute, University of Groningen, Groningen, Netherlands (associated with Nikhef National Institute for Subatomic Physics, Amsterdam, Netherlands)
76
Universiteit Maastricht, Maastricht, Netherlands (associated with Nikhef National Institute for Subatomic Physics, Amsterdam, Netherlands)
77
National Research Centre Kurchatov Institute, Moscow, Russia [associated with Institute of Theoretical and Experimental Physics NRC Kurchatov Institute (ITEP NRC KI), Moscow, Russia, Moscow, Russia]
78
National University of Science and Technology“MISIS”, Moscow, Russia [associated with Institute of Theoretical and Experimental Physics NRC Kurchatov Institute (ITEP NRC KI), Moscow, Russia, Moscow, Russia]
79
National Research University Higher School of Economics, Moscow, Russia (associated with Yandex School of Data Analysis, Moscow, Russia)
80
National Research Tomsk Polytechnic University, Tomsk, Russia [associated with Institute of Theoretical and Experimental Physics NRC Kurchatov Institute (ITEP NRC KI), Moscow, Russia, Moscow, Russia]
81
University of Michigan, Ann Arbor, USA (associated with Syracuse University, Syracuse, New York, USA)
aAlso at Laboratoire Leprince-Ringuet, Palaiseau, France. b
Also at Universit`a di Genova, Genova, Italy.
cAlso at Universit`a di Bologna, Bologna, Italy. d
Also at Universit`a di Modena e Reggio Emilia, Modena, Italy.
eAlso at Novosibirsk State University, Novosibirsk, Russia. f
Also at Universit`a di Ferrara, Ferrara, Italy.
gAlso at Universit`a di Milano Bicocca, Milano, Italy. h
Also at DS4DS, La Salle, Universitat Ramon Llull, Barcelona, Spain.
iAlso at Universit`a di Pisa, Pisa, Italy. j
Also at Universidad Nacional Autonoma de Honduras, Tegucigalpa, Honduras.
kAlso at Universit`a di Bari, Bari, Italy. l
Also at Universit`a di Cagliari, Cagliari, Italy.
mAlso at INFN Sezione di Trieste, Trieste, Italy. n
Also at Universit`a degli Studi di Milano, Milano, Italy.
oAlso at Universidade Federal do Triângulo Mineiro (UFTM), Uberaba-MG, Brazil. p
Also at AGH - University of Science and Technology, Faculty of Computer Science, Electronics and Telecommunications, Kraków, Poland.
qAlso at Universit`a di Siena, Siena, Italy. r
Also at Universit`a di Padova, Padova, Italy.
sAlso at Scuola Normale Superiore, Pisa, Italy. t
Also at MSU - Iligan Institute of Technology (MSU-IIT), Iligan, Philippines.
uAlso at Hanoi University of Science, Hanoi, Vietnam. v
Also at P.N. Lebedev Physical Institute, Russian Academy of Science (LPI RAS), Moscow, Russia.
wAlso at Universit`a di Roma Tor Vergata, Roma, Italy. x
Also at Universit`a della Basilicata, Potenza, Italy.
yAlso at Universit`a di Urbino, Urbino, Italy. z
