Precision measurement of the Bc+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {B}_c^{+} $$\end{document} meson mass

A precision measurement of the Bc+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {B}_c^{+} $$\end{document} meson mass is performed using proton- proton collision data collected with the LHCb experiment at centre-of-mass energies of 7, 8 and 13 TeV, corresponding to a total integrated luminosity of 9.0 fb−1. The Bc+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {B}_c^{+} $$\end{document} mesons are reconstructed via the decays Bc+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {B}_c^{+} $$\end{document}→ J/ψπ+, Bc+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {B}_c^{+} $$\end{document}→ J/ψπ+π−π+, Bc+→J/ψpp¯π+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {B}_c^{+}\to J/\psi p\overline{p}{\pi}^{+} $$\end{document}, Bc+→J/ψDs+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {B}_c^{+}\to J/\psi {D}_s^{+} $$\end{document}, Bc+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {B}_c^{+} $$\end{document}→ J/ψ D0K+ and Bc+→Bs0π+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {B}_c^{+}\to {B}_s^0{\pi}^{+} $$\end{document}. Combining the results of the individual decay channels, the Bc+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {B}_c^{+} $$\end{document} mass is measured to be 6274.47 ± 0.27 (stat) ± 0.17 (syst) MeV/c2. This is the most precise measurement of the Bc+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {B}_c^{+} $$\end{document} mass to date. The difference between the Bc+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {B}_c^{+} $$\end{document} and Bs0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {B}_s^0 $$\end{document} meson masses is measured to be 907.75 ± 0.37 (stat) ± 0.27 (syst) MeV/c.


Introduction
The B c meson family is unique in the Standard Model as its states contain two different heavy-flavour quarks, ab and a c quark. Quantum Chromodynamics (QCD) predicts that theb and c quarks are tightly bound in a compact system, with a rich spectroscopy of excited states. Studies of the B c mass spectrum can reveal information on heavy-quark dynamics and improve our understanding of the strong interaction. Due to the presence of two heavy-flavour quarks the mass spectrum of the B c states can be predicted with much better precision than many other hadronic systems. The mass spectrum of the B c family has been calculated with nonrelativistic quark potential models [1][2][3][4][5][6][7][8], nonperturbative phenomenological models [9,10], perturbative QCD [11,12], relativistic quark models [13][14][15][16][17], and lattice QCD [18][19][20][21][22][23]. The ground state of the B c meson family, denoted hereafter as B + c , decays only through the weak interaction, with a relatively long lifetime. The most accurate prediction of the B + c mass, M (B + c ) = 6278 ± 6 ± 4 MeV/c 2 [22], is obtained with unquenched lattice QCD.

JHEP07(2020)123
B + c → B 0 s π + . 1 The first two decays are chosen for their large signal yield, while the others have a low energy release. As the B 0 s mass is known with limited precision, the difference between the B + c and B 0 s masses, ∆M = M (B + c ) − M (B 0 s ), is also measured, such that improvements in the B + s mass measurement allow for a more precise B + c mass determination. The data sample corresponds to an integrated luminosity of 9.0 fb −1 , collected with the LHCb experiment in pp collisions at centre-of-mass energies of 7, 8 and 13 TeV. The integrated luminosity used in this analysis is at least three times the one used in previous LHCb measurements [26, 31-33] and the results of this paper supersede those earlier B + c mass measurements.

Detector and simulation
This LHCb detector [48,49] is a single-arm forward spectrometer covering the pseudorapidity range 2 < η < 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 [50], 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 [51,52] placed downstream of the magnet. The tracking system provides a measurement of the 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 momentum scale is calibrated using samples of B + → J/ψ K + and J/ψ → µ + µ − decays collected concurrently with the data sample used for this analysis [53,54]. The relative accuracy of this procedure is determined to be 3 × 10 −4 using samples of other fully reconstructed B, Υ and K 0 S -meson decays. 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 are distinguished using information from two ring-imaging Cherenkov detectors [55]. Photons, electrons and hadrons are identified by a calorimeter system consisting of a scintillating-pad and preshower detectors, an electromagnetic and a hadronic calorimeter. Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers [56]. The online event selection is performed by a trigger [57], which consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, which performs a full event reconstruction.
Simulated samples are used to model the effects of the detector acceptance, optimise signal selection and validate the analysis technique. In simulation, pp collisions are generated using Pythia 8 [58] with an LHCb specific configuration [59]. The production of B + c mesons is simulated using the dedicated generator BcVegPy [60]. Decays of hadrons are described by EvtGen [61], in which final-state radiation is generated using Photos 3 [62]. The interaction of the generated particles with the detector and its response are implemented using the Geant4 toolkit [63] as described in ref. [65]. 1 The inclusion of charge-conjugate modes is implied throughout this paper.

Event selection
The B + c candidates are reconstructed in the following decay modes: A pair of oppositely charged muons form J/ψ candidates. The D + s candidates are reconstructed via the D + s → K + K − π + and D + s → π + π − π + decays, while the D 0 is reconstructed using the D 0 → K − π + decay. The B 0 s candidates are reconstructed in the decay modes , and a multivariate classifier as used in ref.
[27] is employed to separate signal from combinatorial background. Then the B 0 s candidates are combined with an additional pion to reconstruct B + c candidates. All of the intermediate-state particles are required to have an invariant mass within three times the expected mass resolution around their known masses [46]. Muons, kaons, pions and protons are required to have good track-fit quality and high transverse momentum. The J/ψ and B + c candidates are required to have a good-quality vertex fit. A boosted decision tree [66][67][68] implemented within the TMVA [69] package optimises separation of the signal from combinatorial background for each decay mode. The classifiers are trained with simulated signal samples and a background proxy obtained from the upper mass sideband of the data, in the range [6.6, 7.0] GeV/c 2 . Kinematic variables that generically separate b-hadron decays from background are used in the training of the classifiers. The variables include the decay time, transverse momenta, vertex-fit quality of the B + c candidate, as well as variables related to the fact that the B + c meson is produced at the PV. The requirement on the classifiers is determined by maximising the signal significance S/ √ S + B, where S is the expected signal yield estimated using simulation, and B is the expected background yield evaluated in the upper sideband in data and extrapolated to the signal region.

Mass measurement
The B + c meson mass is determined in each decay mode by performing an unbinned maximum likelihood fit to the invariant mass distributions of the B + c candidates. The signal is described by a double-sided Crystal Ball (DSCB) function [70], while the background is described by an exponential function. The DSCB function comprises a Gaussian core with power-law tails to account for radiative effects. Parameters describing the radiative tails are determined from simulation.
The invariant mass of the B + c candidates is calculated from a kinematic fit [71], in which the B + c candidate is assumed to originate from its PV and the intermediatestate masses are constrained to their known values [46]. The PV of the B + c candidate is that with respect to which it has the smallest χ 2 IP . The χ 2 IP is defined as the difference in χ 2 of the PV fit with and without the particle in question. For B + c → B 0 s π + decays, the B 0 s mass is constrained to the value of 5366.89 ± 0.21 MeV/c 2 , which is an average of the measurements of the B 0 s mass performed by the LHCb collaboration [72][73][74][75].
and m(B 0 s ) are the reconstructed masses of B + c and B 0 s candidates. The mass difference ∆m is calculated with a kinematic fit [71], in which the B + c candidate is assumed to originate from the PV with the smallest χ 2 IP and the masses of the intermediate particles are constrained to their known values [46]. The fitting procedure for the mass difference is the same as for the mass fit. Figure 1 shows the invariant mass distributions and fit results for all B + c decay modes. Figure 2 shows the distributions of ∆m and fit results for the B + c → B 0 s (D − s π + )π + and B + c → B 0 s (J/ψ φ)π + decay modes. The lower limit of the mass window is chosen to exclude the partially reconstructed background while keeping sufficient left mass sideband. The signal yields, mass and resolution values as determined from fits to the individual mass distributions are given in table 1. For the B + c → B 0 s π + decays, the results of the fits to the ∆m distribution are reported in table 2.
The reconstructed invariant-mass distribution is distorted due to the missing energy from unreconstructed photons (bremsstrahlung) emitted by final-state particles. The resulting bias in the extracted B + c mass is studied with simulated samples for each decay channel, and is used to correct the mass obtained from the fit. Multiple scattering in detector material can decrease the observed opening angles among the B + c decay products, affecting the reconstructed B + c mass and decay length and thereby the selection efficiency. Such effect distorts the mass distribution after the event selection. The corresponding bias of the B + c mass measurement was studied with charmed hadrons (D + , D 0 , D + s , Λ + c ), and was found to be well reproduced by simulation [76]. A bias associated with the selection from simulated samples is assigned as a corresponding correction. The measured masses (M ) and mass difference (∆M ) are corrected for this bias (from -0.46 to 0.27 MeV/c 2 ) due to final-state radiation and the selection, and summarised in table 1 and 2.   Table 2. Signal yields, mass difference (∆M ) and resolution as obtained from fits shown in figure 2, together with the values corrected for the effects of final-state radiation and selection as described in the text. The uncertainties are statistical only.

Systematic uncertainties
To evaluate systematic uncertainties, the complete analysis is repeated varying assumed parameters, models and selection requirements. The observed differences in the B + c mass central values between the nominal result and the alternative estimates are considered as one standard-deviation uncertainties.
The systematic uncertainty of the B + c mass comprises uncertainties on the momentumscale calibration, energy loss corrections, signal and background models, the mass of the intermediate states and the uncertainty on the bias caused by the final-state radiation and selection.
The dominant source of systematic uncertainty arises due to the limited precision of the momentum-scale calibration. For each decay, this uncertainty is propagated to the B + c mass according to the energy release, which is the difference between the value of the B + c mass and the sum of the masses of its intermediate states. The amount of material traversed in the tracking system by a particle is known to 10% accuracy, which leads to an uncertainty on the estimated energy loss. This translates into a measured mass uncertainty of 0.03 MeV/c 2 for D 0 → K + K − π + π − decays [54]. The uncertainties on the B + c mass are scaled from that of the D 0 decay by the number of final-state particles. The uncertainties due to the limited size of simulated samples are taken as systematic uncertainties from the -6 -JHEP07(2020)123  Table 4. Summary of systematic uncertainties on the mass difference ∆M (in MeV/c 2 ) for the B 0 s (D − s π + )π + and B 0 s (J/ψ φ)π + decays. selection-induced bias on the B + c masses. The uncertainty on the masses of the intermediate states D + s , D 0 , B 0 s are propagated to the B + c mass measurement. The uncertainty related to the signal shape is estimated by using alternative signal models, including the sum of two Gaussian functions, a Hypatia function [77], the sum of a DSCB and a Gaussian function, and the sum of two DSCB functions. The differences of the fitted mass with final-state radiation corrections between the nominal and the alternative models are found to be smaller than 0.1 MeV/c 2 , which is taken as the corresponding systematic uncertainty. The uncertainty related to the background description is evaluated by using a first-order Chebyshev function instead of an exponential function.
The non-resonant contribution, for example the contribution of B + c → J/ψ π + π − π + decays to the B + c → J/ψ D + s (π + π − π + ) candidates, is found to be highly suppressed and have negligible effects on the mass measurement. The systematic uncertainties considered for the B + c mass and mass difference measurements are summarised in table 3 and 4, respectively.

Combination of the measurements
The combination of the B + c mass measurements is performed using the Best Linear Unbiased Estimate (BLUE) method [78][79][80]. In the combination, uncertainties arising from the momentum-scale calibration, energy loss corrections, and signal model are assumed to be 100% correlated, while all other sources of systematic uncertainty are assumed to be uncor-  Figure 3. Individual B + c mass measurements and their combination. The red (inner) cross-bars show the statistical uncertainties, and the blue (outer) cross-bars show the total uncertainties. related. The uncertainty on the momentum-scale calibration of the B 0 s mass (0.14 MeV/c 2 ) is assumed to be 100% correlated with that of the B + c mass. The individual mass measurements and the resulting combination are shown in figure 3. The individual measurements are consistent with each other. The breakdown of the combined systematic uncertainty is given in table 5. The weights of individual measurements returned by the BLUE method are listed in table 6. The weights are computed including all uncertainties. The measurement contributing most to the combination is obtained from the B + c → J/ψ D + s (K + K − π + ) decay. The negative weight for the B + c → J/ψ π + channel arises from the 100% correlation between the systematic uncertainties due to the momentum-scale calibration. This results in a larger statistical and smaller systematic uncertainty relative to an uncorrelated average.
The combination for the mass difference ∆M is shown in figure 4. The breakdown of the combined systematic uncertainty is given in table 5 Table 5. Breakdown of systematic uncertainties (in MeV/c 2 ) in the combination of the B + c mass and the mass difference ∆M . The total uncertainty is the sum in quadrature of the uncertainty of different sources.