New algorithms for identifying the flavour of \({\mathrm {B}} ^0\) mesons using pions and protons

The European Physical Journal C, Apr 2017

Two new algorithms for use in the analysis of \(pp\) collision are developed to identify the flavour of \({\mathrm {B}} ^0\) mesons at production using pions and protons from the hadronization process. The algorithms are optimized and calibrated on data, using \({{\mathrm {B}} ^0} \!\rightarrow D^{-}\pi ^{+}\) decays from \(pp\) collision data collected by LHCb at centre-of-mass energies of 7 and 8 TeV . The tagging power of the new pion algorithm is 60% greater than the previously available one; the algorithm using protons to identify the flavour of a \({\mathrm {B}} ^0\) meson is the first of its kind.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-017-4731-y.pdf

New algorithms for identifying the flavour of \({\mathrm {B}} ^0\) mesons using pions and protons

Eur. Phys. J. C New algorithms for identifying the flavour of B0 mesons using pions and protons LHCb Collaboration 0 0 CERN , 1211 Geneva 23 , Switzerland Two new algorithms for use in the analysis of pp collision are developed to identify the flavour of B0 mesons at production using pions and protons from the hadronization process. The algorithms are optimized and calibrated on data, using B0 ? D?? + decays from pp collision data collected by LHCb at centre-of-mass energies of 7 and 8 TeV . The tagging power of the new pion algorithm is 60% greater than the previously available one; the algorithm using protons to identify the flavour of a B0 meson is the first of its kind. 1 Introduction Violation of CP symmetry in the B system was observed for the first time in the interference between mixing and decay processes [1]. Any measurement of a decay-time-dependent asymmetry requires the determination of the flavour of the B meson at production. For B mesons produced in pp collisions, this information is obtained by means of several flavour-tagging algorithms that exploit the correlations between B flavour and other particles in the event. Algorithms determining the flavour content of B meson by using particles associated to its production are called sameside (SS) taggers. As an example, in the production of B0 mesons from excited charged B mesons decaying via strong interaction to B0? +, the pion charge identifies the initial flavour of the B0 meson.1 A charge correlation can also arise from the hadronization process of the b quark. When a b and a d quark hadronize as a B0 meson, it is likely that the corresponding d quark ends up in a charged pion (ud), or in an antiproton (u u d). The B0 meson and the pion or antiproton are produced in nearby regions of phase space. Other algorithms used at LHCb, called opposite-side (OS) taggers [2,3], attempt to identify the flavour of the other b hadron produced in the same event. 1 The inclusion of charge-conjugate processes is implied throughout the paper, unless otherwise noted. e-mail: A simple cut-based SS algorithm selecting pions was successfully used by LHCb for tagging B0 ? J/? KS0 decays [4] in the measurement of sin 2?, and an SS kaon tagger [5] based on a neural network was used to determine the flavour of Bs0 mesons in measurements of the CP-violating phase ?s [6? 8]. This paper presents two new SS algorithms exploiting the charge correlation of pions and protons with B0 mesons, denoted SS? and SS p. This is the first time that protons are used for flavour tagging. The two algorithms are combined into a single tagger, SScomb. Both algorithms are based on multivariate selections and are optimized, calibrated and validated using B0 ? D?? + and B0 ? K +?? decays collected by LHCb in Run 1. The performance of a flavour-tagging algorithm is measured by its tagging efficiency ?tag, mistag fraction ?, dilution D, and tagging power ?eff , defined as ?tag = R +R +W W+ U , ? = R +W W , D = 1 ? 2?, ?eff = ?tag D2, where R, W , and U are the numbers of correctly-tagged, incorrectly-tagged, and untagged B0 signal candidates. The tagging power determines the sensitivity to the measurement of a decay-time-dependent CP asymmetry [9], as it quantifies the effective reduction in the sample size of flavour-tagged B0 candidates. It is the figure of merit used to optimize the algorithms. Each algorithm provides a decision on the flavour of the B0 candidate and an estimate of the probability ? that this decision is incorrect. The probability is used to determine a weight applied to the B0 candidate, in order to maximize the tagging power of a sample of B0 mesons in a time-dependent analysis. The probabilities provided by the two SS taggers are used to combine their decisions into the SScomb decision, which can be further combined with the decision of other taggers [2,3]. The expected relationship between the flavour of charged and neutral B mesons and the charge of the tagging particle is reported in Table 1. For a B+ meson the same correlation as for a B0 meson holds in the case of protons, but with opposite Table 1 Expected correlation between the flavour of a B meson and the hadronization products K ? charge in the case of pions. In addition, the tagging kaons carry the same charge as pions, while they are neutral for a B0. Since misidentified hadrons affect the tagging efficiency and the mistag fraction of charged and neutral mesons in different ways, B+ decays cannot be reliably used for the tuning and calibration of the SS taggers. As a consequence, B0 decays are used, and a time-dependent analysis is required to determine the mistag fraction. 2 Detector The LHCb detector [10,11] 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, 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. Regular reversal of the magnet polarity allows a quantitative assessment of detector-induced charge asymmetries. 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/ pT) ?m, where pT is the component of the momentum transverse to the beam, in GeV/c. Particularly relevant for this analysis is the identification of the different species of charged hadrons, which mainly relies on the information of two ring-imaging Cherenkov detectors. The first one covers the low and intermediate momentum region 2?40 GeV/c over the full spectrometer angular acceptance of 25?300 mrad. The second Cherenkov detector covers the high momentum region 15?100 GeV/c over the angular range 15?120 mrad [12]. 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 [13], 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 muon with high pT or a hadron, photon or electron with high transverse energy in the calorimeters. The software trigger requires a two-, three- or four-track secondary vertex detached from the PV. A multivariate algorithm [14] is used for the identification of secondary vertices consistent with the decay of a b hadron. Samples of simulated events are used to model the signal mass and decay-time distributions. In the simulation, pp collisions are generated using Pythia [15,16] with a specific LHCb configuration [17]. Decays of hadronic particles are described by EvtGen [18], in which final-state radiation is generated using Photos [19]. The interaction of the generated particles with the detector, and its response, are implemented using the Geant4 toolkit [20,21] as described in Ref. [22]. 3 Development of the same-side taggers The SS? and SS p algorithms are developed following similar strategies. A sample of B0 mesons decaying into the flavourspecific final state D?? +, with D? candidates reconstructed in the final state K +????, is selected using requirements similar to those presented in Ref. [23]. The sample is collected from pp collisions at ?s = 8 TeV , corresponding to an integrated luminosity of 2fb?1. Tagging pion or proton candidates, with their charge correlated with the B0 flavour, are selected by means of a set of loose selection requirements and a multivariate classifier, as described below. The B0 ? D?? + candidates are separated randomly into three disjoint subsamples of equal size. The first sample is used for training the multivariate classifiers, the second is used for determining the probability of an incorrect tagging decision, and the third is used to evaluate the calibration of the mistag probability. The correctness of a tagging decision is evaluated by comparing the charge of the tagging particle with the B0 decay flavour as determined by the reconstructed final state. Those B0 candidates that have oscillated before decaying enter the training process with an incorrectly assigned production flavour. In the training phase the dilution is reduced by requiring the decay time of the reconstructed B0 mesons to be smaller than 2.2 ps. This value was optimized with simulated events and reduces the fraction of oscillated candidates to about 11%, keeping 66% of the original sample. The signal and background components of the B0 sample are determined by an unbinned maximum likelihood fit to the D??+ mass distribution of the selected candidates in the region [5.2, 5.5] GeV/c2. The signal is described by a Johnson?s SU distribution [24], while the combinatorial Fig. 1 Mass distribution of B0 ? D?? + candidates with fit projections overlaid. Data points (black dots) correspond to the B0 candidates selected in the 2 fb?1 data sample collected at ?s = 8 TeV . The solid blue curve represents the total fit function which is the sum of signal (red dashed) and combinatorial background (green dash-dotted) background is modelled by the sum of an exponential function and a constant. All parameters are free to vary in the fit. A small component of B0 ? D? K + decays (?1.2% as estimated from simulation), with the kaon misidentified as a pion, is neglected in the fit. The number of signal candidates in the full 2fb?1 sample, estimated by the mass fit and shown in Fig. 1, is 300 370 ? 674. The fit to the mass distribution is used to assign event-by-event weights (sWeights), using the sPlot technique [25]. The weights are subsequently used to subtract the background contribution when training the SS? and SS p classifiers and in the fits to the B0 decay-time distribution. The loose selection requirements reduce the multiplicity of pion (proton) candidates to 2.3 (1.7) per B0 ? D?? + signal candidate, and are reported in Table 2. Only tracks with hits in all tracking detectors are considered as tagging candidates. The following observables are used: the ? 2/ndf of the track fit, where ndf is the number of degrees of freedom, the track transverse momentum ptrack, the ratio between T the track impact parameter with respect to the PV associated to the B0 meson and the error on this variable IP/?IP, the ratio between the track impact parameter with respect to any other PV in the event and its error IPPU/?IPPU , the difference between the logarithms of the likelihood of the proton and pion hypothesis log L p ? log L? , or kaon and pion hypothesis log L K ? log L? . The likelihoods for the various mass hypothesis are determined using the track and the Cherenkov angle information, as described in Ref. [26]. For particles passing the loose selection criteria the efficiency to identify a pion is 89% with a kaon misidentification probability of 2%, while the efficiency to identify a proton is 92% with a pion misidentification probability of 5%. Since mutually exclusive particle identification criteria are imposed, a given tagging track is identified either as a pion or as a proton. If more than one PV is reconstructed in the event, the PV Table 2 Loose selection requirements for the SS? and SS p algorithms. The variables used as input for the BDT classifiers are indicated by <3 >0.4 ? <4 ? >5 ? associated to the B0 meson is the one which has the smallest increase in the vertex-fit ? 2 when adding the B0 meson to the PV. Additional requirements are introduced on the system formed by the tagging particle and the reconstructed B0 meson. They are applied to the total transverse momentum of the system pTtot, the difference between the pseudorapidity of the B0 candidate and the tagging particle ?, the azimuthal angle ? between the B0 candidate and the tagging particle, and the difference between the invariant mass of the system and the mass of the B0 and of the tagging particle Q = m(B0 + h) ? m(B0) ? m(h), where h denotes the hadron, ? or p. The vertex formed by the B0 meson and the tagging particle is required to have the ? 2 of vertex fit 2 ?B0?track, less than 100. The multivariate classifiers used for the selection of the tagging particles are boosted decision trees (BDT) [27] using the AdaBoost [28] method to enhance and to increase the stability with respect to statistical fluctuations. This choice has been shown to be optimal with respect to the achievable tagging power. The classifiers take most of the above observables as input, as specified in Table 2. In addition the BDTs use the following variables: the momentum of the tagging particle ptrack, the transverse momentum of the B0 candidate pTB0 , the separation of tagging particle and the B0 candidate R = ?2 + ?2, and the number of tracks contributing to the PV fit P Vtracks. The sWeights are used to subtract the contribution of background B0 candidates in the training of the classifiers. The charge of the tagging particle determines 0?1 0?1 Fig. 3 Measured average mistag fraction in bins of (left) SS? and (right) SS p BDT output. The plots are obtained with the test sample of background-subtracted B0 ? D?? + candidates. The green-shaded area shows the confidence range within ?1? the flavour of the B0 candidate. In case of multiple tagging particle candidates per B0 candidate, the tagging particle with the highest BDT output value is chosen. The BDT outputs, ?BDT, are shown in Fig. 2. The global separation between signal and background is small, but enough to provide useful information to determine the flavour of the B0 candidate, as shown below. 4 Evaluation and calibration of mistag probability 4.1 The SS? and SSp taggers The BDT output is transformed into an estimate of the mistag probability through linear regression. The decay-time distribution of all tagged B0 candidates is considered and the dilution due to mixing is decoupled by means of a full timedependent analysis. Tagged B0 candidates are divided into eight bins of the BDT output and for each bin the probability of an incorrect tagging decision is determined from an unbinned maximum likelihood fit to the distribution of the measured decay time t of the candidates, using the sWeights. The probability density function (PDF) for the signal is described as S(t, q) = N a(t ) e?t /?d (1 + q(1 ? 2?) cos( md t )) ?R(t ? t ), where t represents the true decay time, N is a normalization factor, ? is the average mistag fraction in the bin, q is the mixing state (q = +1 when the flavour at production and the flavour at decay are the same, q = ?1 otherwise), R(t ? t ) is the decay-time resolution and a(t ) is the decay-time acceptance. The B0 lifetime ?d , and the mixing frequency md , are fixed in the fit to their known values [29]. Equation 2 is obtained under the assumption of zero width difference d and neglecting the production and detection asymmetries between B0 and B0. The decay-time resolution is modelled by a Gaussian function with a fixed width of 50 fs, as determined from simulation. The decay-time acceptance a(t ), is described by a parametric function based on cubic splines [30] whose nodes have fixed position and whose parameters are determined from data. Figure 3 shows the measured average mistag rate per subsample, interpolated with a third-order polynomial that represents ? as a function of ?BDT, for the SS? and SS p taggers. This polynomial parametrization is then used to determine the mistag probability ?(?BDT) of a B0 candidate. Tagging particles with ?(?BDT) > 0.5 are rejected. With the third subsample of B0 candidates, it is checked that the estimated mistag probability corresponds to the true value by measuring the mistag fraction ? with an unbinned likelihood fit to the decay-time distribution of the B0 candidates. Possible differences between the mistag probability of B0 and B0 mesons may arise from the different interaction crosssections of hadrons and antihadrons in the detector material and from differences in detection efficiencies of positive and negative hadrons. They are taken into account in the decaytime fit by defining the variables where ? B0 and ? B0 are the mistag fractions related to B0 and B0. Assuming a linear relation between the measured and estimated mistag fractions, the calibration functions are written as ? B0 (?) = p0B0 + p1B0 (? ? ? ), ? B0 (?) = p0B0 + p1B0 (? ? ? ), where piB0 and piB0 (with i = 0, 1) are the calibration parameters. The average calibration parameters and the differences between the B0 and B0 parameters are defined as pi = ( piB0 + piB0 )/2, pi = piB0 ? piB0 . The use of the arithmetic mean ? of the ? distribution aims at decorrelating p0 and p1. A perfect calibration corresponds to p0 = ? and p1 = 1. A difference in the number of reconstructed and tagged B0 and B0 mesons arises from several possible sources. Two of these sources are considered in the fit by introducing an asymmetry in the detection efficiency of the final state particles, defined as Adet = Table 3 Calibration parameters for the SS?, SS p and SScomb taggers where the first uncertainties are statistical and the second are systematic Atag = With these additional inputs, the PDF becomes S(t , q) = N a(t ) e?t /?d (Ccosh + Ccos cos( ?R(t ? t ). The coefficients Ccosh and Ccos are + (1 ? Aprod ? Atag) and an asymmetry of the tagging efficiencies, defined as ? (1 ? Aprod ? Atag) where r is the B meson flavour at decay (r = +1 for B0 ? D?? +, r = ?1 for B0 ? D+? ?) and d is the tagging decision (d = +1 for ? + ( p), d = ?1 for ? ? ( p)). These coefficients also take into account the production asymmetry, Aprod = NNBB00 ?+NNBB00 , and the asymmetry in mixing, or flavour-specific asymmetry, asdl. These two asymmetries cannot be distinguished from the tagging and detection asymmetries and are fixed in the fit. The production asymmetry is fixed to the value measured in Ref. [31], Aprod = (?0.58 ? 0.70)%, while asdl is fixed to the world average asdl = (?0.15 ? 0.17)% [32]. The effect of their uncertainties on the calibration parameters is included in the systematic uncertainty. The calibration parameters for the two taggers obtained in the fit to the calibration sample of B0 ? D?? + decays are reported in Table 3. The correlations between the calibration parameters are below 10%, except for the asymmetry ? p0 p1 p0 p1 Atag Table 4 Tagging efficiencies and tagging power of the SS?, SS p and SScomb algorithms. The SScomb efficiencies are shown splitting the sample in candidates tagged exclusively by SS? or SS p, or by both. As explained in the text, there is a large overlap between the SS? and SS p taggers of the tagging efficiencies, which has a correlation of about 16% with p0 and p1 and about 64% with Adet. For the SS? tagger, Atag, p0 and p1 are zero within one standard deviation, showing no significant difference in tagging behaviour between B0 and B0 decays. For the SS p tagger, it is found that p0 < 0, as a consequence of the higher interaction cross-section of anti-protons with matter compared to protons. A similar effect is reported for kaon taggers [5]. The fit result of the detection asymmetry is comparable for the two taggers ( AdSeSt? = (?0.87 ? 0.48)%, AdSeStp = (?0.66 ? 0.62)%) and in agreement with that found in Ref. [33]. The systematic uncertainties on the parameters will be described in Sect. 5. After calibration, the total tagging power of the sample is calculated as where si is the sWeight of the candidate i , N and Ntag are the numbers of total and tagged candidates, having mistag probability ?i , and the average mistag fraction ?(?i ) is calculated using Eqs. 3 and 4. Candidates with a mistag probability larger than 0.5 are considered untagged and are removed from the sum in the numerator, effectively setting ?(?i ) = 0.5. The tagging performances for the SS? and SS p taggers are reported in Table 4. The fit of the decay-time distribution is repeated after dividing events into bins of predicted mistag probability. The distribution of ? and the dependence of the measured mistag fraction on ? are shown in Fig. 4 with the linear fits superimposed, demonstrating the expected linearity. In Figs. 5 and N unmix?N mix 6 the time-dependent mixing asymmetries A = N unmix+N mix are shown for each of the five bins. 4.2 The SScomb tagger Even though a given tagging particle can be selected by only one of the SS? or the SS p taggers, both taggers may find a candidate track in the same event. About 50% of the candidates tagged by SS? are also tagged by SS p, and 80% of the candidates tagged by SS p are also tagged by SS?. When both taggers provide a decision, they are combined into a single decision. Since the correlation between the SS? and SS p decisions, and between their mistag probabilities, is found to be small, it is neglected when combining them using the following formulae p(b) = p(b) = where p(b) and p(b) are the probabilities that the signal B meson contains a b or a b quark respectively, and di is the tagging decision of the tagger i = SS?, SS p. The normalized probabilities are P(b) = 1 ? P(b). For P(b) > P(b) the combined tagging decision is d = +1 and the final mistag probability is ? = P(b). Otherwise, the combined tagging decision and the mistag probability are d = ?1 and ? = P(b). The combination procedure, which assumes no correlation, is validated by checking the combined mistag probability a posteriori. Assuming a linear relation between the predicted mistag probability and the true mistag fraction, the calibration parameters in the overlapping sample give ( p0 ? ? ) = 0.010 ? 0.005 and ( p1 ? 1) = 0.01 ? 0.08. The calibration is repeated on the sample of all B0 candidates tagged by the SScomb tagger, and the calibration parameters derived from the unbinned likelihood fit with the PDF of Eq. 8, reported in Table 3, demonstrate its validity. The performance of SScomb is reported in Table 4. The total tagging power obtained by the combined algorithm is (2.11?0.11)%, a relative increase of 25% compared to that provided by the SS? tagger alone. A higher tagging power can be obtained from the combination of the SScomb tagger with the OS tagger. The OS tagger is the combination of various OS tagging algorithms using electrons and muons from semileptonic decays of b hadrons, kaons from b ? c ? s decay chains and the inclusive reconstruction of a secondary vertex of the decay products of the opposite side b hadron. The SS and OS taggers are found to be uncorrelated, so their combination follows the same procedure as the combination of SS? and SS p into SScomb. The calibration of the combined mistag probability is verified a posteriori with a fit of the decay-time distribution of the B0 candidates. For B0 ? D?? + decays, the total tagging efficiency and the total tagging power are (84.48 ? 0.26)% and (5.14 ? 0.15)%, respectively. On the same sample, the s t i un0.12 y rra 0.1 t i rab0.08 s t i n u 0.1 y r a itr0.08 b r a0.06 Fig. 4 (Top left) distribution of the mistag probability ?SS? and (top right) measured mistag fraction ? as a function of ?SS?. (Bottom left) distribution of the mistag probability ?SSp and (bottom right) measured use of the OS tagger only provides a tagging efficiency and a tagging power of (37.95 ? 0.15)% and (3.52 ? 0.17)%, respectively. 5 Validation and systematic uncertainties A possible dependence of the calibration parameters of the SS taggers on properties of the event sample is checked by repeating the calibration after splitting the data according to the data-taking conditions (magnet polarity), global event properties (total number of reconstructed tracks, number of primary vertices) or according to the kinematic properties of the B0 meson (transverse momentum, pseudorapidity and azimuthal angle). The average mistag probability has a weak dependence on the number of tracks in the event. On the other hand, it decreases as a function of the transverse momentum since the number of random tracks decreases at high p B . T The tagging efficiency is nearly constant for pions, while the requirement on proton identification reduces the number of proton candidates at high pTB . A similar dependence is present versus the pseudorapidity of the B0 meson. Since the average mistag fraction and the p0 parameter decrease with increasmistag fraction ? as a function of ?SSp. The green-shaded area shows the 68% confidence range ing pTB0 , the calibration remains valid in all subsamples, with variations below two standard deviations. The portability of the mistag calibration, from the training data sample to other data samples and other B0 decay modes, is validated using an independent sample of B0 ? D?? + decays collected at ?s = 7 TeV (corresponding to an integrated luminosity of 1fb?1) and a sample of B0 ? K +?? decays collected at ?s = 8 TeV (corresponding to an integrated luminosity of 2fb?1). The same selection criteria and fitting procedure as described above are used for the B0 ? D?? + validation sample at ?s = 7 TeV . The calibration parameters for the SS?, SS p, and SScomb taggers determined from an unbinned maximum likelihood fit to the decay-time distribution are compatible with those derived in the 8 TeV sample. Consistent values of tagging power are found for all taggers. The selection criteria and the mass model for the B0 ? K +?? candidates are described in Ref. [34]. The decaytime acceptance is parametrized using cubic splines with six nodes, whose positions are fixed and whose coefficients are free in the fit. The decay-time resolution is described by a Gaussian function with parameters determined from simulation. The parameters shown in Table 5 demonstrate a good portability of the mistag calibration, with p0 ? ? ? 0 and try 0.6 e m0.4 m y sa 0.2 g inx 0 i M?0.2 try 0.6 e m0.4 m y sa 0.2 g inx 0 i M?0.2 p1 ? 1 ? 0 as expected. A lower tagging power is measured in this channel, giving (1.06 ? 0.09)%, (0.42 ? 0.06)%, and (1.37 ? 0.13)% for SS?, SS p and SScomb, respectively, as expected from the lower average pT of the selected B0 candidates. Several sources of systematic uncertainties on the calibration parameters are studied and the associated uncertainties are reported in Table 6. Uncertainties related to the mass model and background unfolding procedure are assessed by repeating the calibration replacing the sWeights derived in the fit to the mass distribution of all B0 candidates by the sWeights derived after restricting the sample to tagged B0 candidates. In a second test, the signal mass model is replaced by a Hypatia function [35] convolved with a Gaussian function. The sum in quadrature of the variations of the calibration parameters observed in the two tests is taken as uncertainty on the mass model. Uncertainties related to the decay-time acceptance model are assessed by changing the number of nodes in the cubic splines from six to nine and are found to be negligible. A negligible uncertainty is associated to the decay-time resolution model. The mistag model uncertainties are assessed by comparing the calibration parameters derived in the nominal fit and those derived in fits with the mistag probability binned in categories. Five, seven and nine bins are tested and the largest observed variation of the parameters is taken as try 0.6 e m0.4 m y sa 0.2 g inx 0 i M?0.2 Fig. 6 Mixing asymmetry in bins of mistag probability using the SS p tagger a systematic uncertainty. Differences between the results of the two implementations of the time-dependent fit are due to the dependence of the mistag probability on the decay time. Pseudoexperiments are generated where the mistag probability has the same dependence on time as in data and are fitted with the two approaches. The difference in parameters is similar to or smaller than that observed in data. Uncertainties related to neglecting d and possible CP violation in the B0 ? D?? + decays in the decay-time fit, are studied by performing pseudoexperiments in which changes associated with the parameter under study are incorporated in the generation and neglected in the subsequent fit. Terms proportional to the relevant CP parameters are added to the PDF in Eq. 8 and the values of the parameters are taken from Table 6 Systematic uncertainties on the calibration parameters of SS?, SS p and SScomb taggers. The total systematic uncertainty is the squared sum of all contributions. A dash indicates a value negligible with respect to the quoted precision Table 7 Systematic uncertainties related to the decay-time model. A dash indicates a value negligible with respect to the quoted precision Ref. [32]. The associated systematic uncertainties are taken to be the changes in the calibration parameters with respect to perfect calibration ( p0 = ? , p1=1), used in the generation. Uncertainties related to the variation of Aprod and asdl, which are fixed in the decay-time fit, are evaluated with pseudoexperiments where the parameters are varied within their uncertainties. The uncertainties are determined in the SS? configuration and attributed to both taggers. A breakdown of the systematic uncertainties related to the decay-time model is shown in Table 7. 6 Conclusion Two new same-side algorithms are developed to determine the production flavour of B0 mesons using pions and protons from the hadronization process. This is the first time that protons are used to identify the flavour of a B0 meson. The algorithms are optimized and calibrated on data using B0 ? D?? + decays. The calibration parameters of the taggers are reported in Table 3. The efficiency and mistag probability of the taggers depend on the kinematic properties of the B0 decay mode under study. Estimated mistag probabilities match the true mistag fraction throughout the phase space. The new SS? tagger provides a tagging power that is greater by 60% relative to the previous algorithm using pions, employed in Ref. [4]. Adding the combination of the two new algorithms to the existing OS taggers provides a relative increase of the total tagging power of about 40%. 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); NSFC (China); CNRS/IN2P3 (France); BMBF, DFG and MPG (Germany); INFN (Italy); FOM and 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), IFINHH (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), Conseil G?n?ral de Haute-Savoie, Labex ENIGMASS and OCEVU, R?gion Auvergne (France), RFBR and Yandex LLC (Russia), GVA, XuntaGal and GENCAT (Spain), Herchel Smith Fund, The Royal Society, Royal Commission for the Exhibition of 1851 and the Leverhulme Trust (United Kingdom). Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP3. LHCb Collaboration 1 Centro Brasileiro de Pesquisas F?sicas (CBPF), Rio de Janeiro, Brazil 2 Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil 3 Center for High Energy Physics, Tsinghua University, Beijing, China 4 LAPP, Universit? Savoie Mont-Blanc, CNRS/IN2P3, Annecy-Le-Vieux, France 5 Clermont Universit?, Universit? Blaise Pascal, CNRS/IN2P3, LPC, Clermont-Ferrand, France 6 CPPM, Aix-Marseille Universit?, CNRS/IN2P3, Marseille, France 7 LAL, Universit? Paris-Sud, CNRS/IN2P3, Orsay, France 8 LPNHE, Universit? Pierre et Marie Curie, Universit? Paris Diderot, CNRS/IN2P3, Paris, France 9 I. Physikalisches Institut, RWTH Aachen University, Aachen, Germany a Universidade Federal do Tri?ngulo Mineiro (UFTM), Uberaba-MG, Brazil b Laboratoire Leprince-Ringuet, Palaiseau, France c P.N. Lebedev Physical Institute, Russian Academy of Science (LPI RAS), Moscow, Russia d Universit? di Bari, Bari, Italy e Universit? di Bologna, Bologna, Italy f Universit? di Cagliari, Cagliari, Italy g Universit? di Ferrara, Ferrara, Italy h Universit? di Genova, Genova, Italy i Universit? di Milano Bicocca, Milano, Italy j Universit? di Roma Tor Vergata, Roma, Italy k Universit? di Roma La Sapienza, Roma, Italy l AGH-University of Science and Technology, Faculty of Computer Science, Electronics and Telecommunications, Krak?w, Poland m LIFAELS, La Salle, Universitat Ramon Llull, Barcelona, Spain n Hanoi University of Science, Hanoi, Viet Nam o Universit? di Padova, Padova, Italy p Universit? di Pisa, Pisa, Italy q Universit? degli Studi di Milano, Milano, Italy r Universit? di Urbino, Urbino, Italy s Universit? della Basilicata, Potenza, Italy t Scuola Normale Superiore, Pisa, Italy u Universit? di Modena e Reggio Emilia, Modena, Italy v Iligan Institute of Technology (IIT), Iligan, Philippines w Novosibirsk State University, Novosibirsk, Russia ?Deceased 1. Particle Data Group , K. A. Olive et al., Review of particle physics. CP violation in the quark sector (pages 223 - 234 ). Chin. Phys . C. 38 , 090001 ( 2014 ). doi:10.1088/ 1674 - 1137 /38/9/090001 2. LHCb collaboration , R. Aaij et al., Opposite-side flavour tagging of B mesons at the LHCb experiment . Eur. Phys. J. C . 72 , 2022 ( 2012 ). doi:10.1140/epjc/s10052- 012 - 2022 -1. arXiv: 1202 .4979 3. LHCb collaboration , R. Aaij et al., B flavour tagging using charm decays at the LHCb experiment . JINST 10 , P10005 ( 2015 ). doi:10. 1088/ 1748 - 0221 /10/10/P10005. arXiv: 1507 .07892 4. LHCb collaboration , R. Aaij et al., Measurement of CP violation in B0 ? J /? K S0 decays. Phys. Rev. Lett . 115 , 031601 ( 2015 ). doi:10.1103/PhysRevLett.115.031601. arXiv:1503.07089 5. LHCb collaboration , R. Aaij et al., Neural-network-based same side kaon tagging algorithm calibrated with Bs0 ? Ds?? + and Bs? 2 ( 5840 )0 ? B+ K ? decays . JINST 11 , P05010 ( 2016 ). doi:10. 1088/ 1748 - 0221 /11/05/P05010. arXiv: 1602 .07252 6. LHCb collaboration , R. Aaij et al., Precision measurement of CP violation in Bs0 ? J /? K + K ? decays. Phys. Rev. Lett . 114 , 041801 ( 2015 ). doi:10.1103/PhysRevLett.114.041801. arXiv:1411.3104 7. LHCb collaboration , R. Aaij et al., Measurement of the CPviolating phase ?s in B?s0 ? J /? ? +? ? decays. Phys. Lett. B . 736 , 186 ( 2014 ). doi:10.1016/j.physletb. 2014 .06.079. arXiv:1405.4140 8. LHCb collaboration , R. Aaij et al., Measurement of the CPviolating phase ?s in B?s0 ? Ds+ Ds? decays . Phys. Rev. Lett . 113 , 211801 ( 2014 ). doi:10.1103/PhysRevLett.113.211801. arXiv:1409.4619 9. A.J. Bevan , B. Golob , Th. Mannel , S. Prell , B.D. Yabsley (eds.), The Physics of the B Factories. Eur. Phys. J.C . 74 , 3026 ( 2014 ). doi:10.1140/epjc/s10052- 014 - 3026 -9. arXiv: 1406 .6311 10. LHCb collaboration, A.A. Alves Jr . et al., The LHCb detector at the LHC . JINST 3, S08005 ( 2008 ). doi:10.1088/ 1748 - 0221 /3/08/ S08005 11. LHCb collaboration, R. Aaij et al., LHCb detector performance . Int. J. Mod. Phys. A . 30 , 1530022 ( 2015 ). doi:10.1142/ S0217751X15300227. arXiv: 1412 .6352 12. M. Adinolfi et al., Performance of the LHCb RICH detector at the LHC . Eur. Phys. J. C . 73 , 2431 ( 2013 ). doi:10.1140/epjc/ s10052- 013 - 2431 -9. arXiv: 1211 .6759 13. R. Aaij et al., The LHCb trigger and its performance in 2011 . JINST 8, P04022 ( 2013 ). doi:10.1088/ 1748 - 0221 /8/04/P04022. arXiv: 1211 .3055 14. V.V. Gligorov , M. Williams , Efficient, reliable and fast high-level triggering using a bonsai boosted decision tree . JINST 8, P02013 ( 2013 ). doi:10.1088/ 1748 - 0221 /8/02/P02013. arXiv: 1210 .6861 15. T. Sj?strand , S. Mrenna , P. Skands , PYTHIA 6.4 physics and manual. JHEP 05 , 026 ( 2006 ). doi:10.1088/ 1126 - 6708 /2006/05/026. arXiv:hep-ph/ 0603175 16. T. Sj?strand , S. Mrenna , P. Skands , A brief introduction to PYTHIA 8.1. Comput . Phys. Commun . 178 , 852 ( 2008 ). doi:10.1016/j.cpc. 2008 .01.036. arXiv:0710.3820 17. I. Belyaev et al., Handling of the generation of primary events in Gauss, the LHCb simulation framework . J. Phys. Conf. Ser . 331 , 032047 ( 2011 ). doi:10.1088/ 1742 - 6596 /331/3/032047 18. D.J. Lange , The EvtGen particle decay simulation package . Nucl. Instrum. Methods A . 462 , 152 ( 2001 ). doi:10.1016/ S0168-9002(01)00089- 4 19. P. Golonka , Z. Was , PHOTOS Monte Carlo: a precision tool for QED corrections in Z and W decays . Eur. Phys. J. C . 45 , 97 ( 2006 ). doi:10.1140/epjc/s2005- 02396 -4. arXiv:hep-ph/ 0506026 20. Geant4 collaboration , J. Allison et al., Geant4 developments and applications . IEEE Trans. Nucl. Sci . 53 , 270 ( 2006 ). doi:10.1109/ TNS. 2006 .869826 21. Geant4 collaboration, S. Agostinelli et al., Geant4: a simulation toolkit . Nucl. Instrum. Methods A . 506 , 250 ( 2003 ). doi:10.1016/ S0168-9002(03)01368- 8 22. M. Clemencic et al., The LHCb simulation application , Gauss: design, evolution and experience. J. Phys. Conf. Ser . 331 , 032023 ( 2011 ). doi:10.1088/ 1742 - 6596 /331/3/032023 23. LHCb collaboration, R. Aaij et al., Measurement of the B?s0 meson lifetime in Ds+? ? decays . Phys. Rev. Lett . 113 , 172001 ( 2014 ). doi:10.1103/PhysRevLett.113.172001. arXiv:1407.5873 24. N.L. Johnson , Systems of frequency curves generated by methods of translation . Biometrika 36 , 149 ( 1949 ). doi:10.1093/biomet/36. 1- 2 . 149 25. M. Pivk , F.R. Le Diberder , sPlot: a statistical tool to unfold data distributions . Nucl. Instrum. Methods A . 555 , 356 ( 2005 ). doi:10. 1016/j.nima. 2005 .08.106. arXiv:physics/0402083 26. A.A. Alves Jr . et al., Performance of the LHCb muon system . JINST 8, P02022 ( 2013 ). doi:10.1088/ 1748 - 0221 /8/02/P02022. arXiv: 1211 .1346 27. L. Breiman , J.H. Friedman , R.A. Olshen , C.J. Stone , Classification and regression trees (Wadsworth international group , Belmont , 1984 ) 28. R.E. Schapire , Y. Freund , A decision-theoretic generalization of on-line learning and an application to boosting . J. Comput. Syst. Sci. 55 , 119 ( 1997 ). doi:10.1006/jcss.1997.1504 29. Particle Data Group , K.A. Olive et al., Review of particle physics. Chin. Phys. C . 38 , 090001 ( 2014 ). doi:10.1088/ 1674 - 1137 /38/9/ 090001 (and 2015 update) 30. C. de Boor, A Practical Guide to Splines, revised edn. (Springer, New York, 2001 ) 31. LHCb collaboration, R. Aaij et al., Measurement of the B? 0-B0 and B ?s0-Bs0 production asymmetries in pp collisions at ?s = 7 TeV . Phys. Lett . B. 739 , 218 ( 2014 ). doi:10.1016/j.physletb. 2014 . 10.005. arXiv:1408.0275 32. Heavy Flavor Averaging Group , Y. Amhis et al., Averages of bhadron , c-hadron, and ? -lepton properties as of summer 2014 . arXiv: 1412 .7515. Updated results and plots available at http:// www.slac.stanford.edu/xorg/hfag/ 33. LHCb collaboration, R. Aaij et al., Measurement of the semileptonic CP asymmetry in B0-B? 0 mixing . Phys. Rev. Lett . 114 , 041601 ( 2015 ). doi:10.1103/PhysRevLett.114.041601. arXiv:1409.8586 34. LHCb collaboration, R. Aaij et al., First measurement of timedependent CP violation in Bs0 ? K + K ? decays . JHEP 10 , 183 ( 2013 ). doi:10.1007/ JHEP10(2013)183 . arXiv:1308.1428 35. D. Martinez-Santos , F. Dupertuis , Mass distributions marginalized over per-event errors . Nucl. Instrum. Methods A . 764 , 150 ( 2014 ). doi:10.1016/j.nima. 2014 .06.081. arXiv:1312.5000


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-017-4731-y.pdf

New algorithms for identifying the flavour of \({\mathrm {B}} ^0\) mesons using pions and protons, The European Physical Journal C, 2017, DOI: 10.1140/epjc/s10052-017-4731-y