Physics Letters B 762 (2016) 484–492 Contents lists available at ScienceDirect Physics Letters B www.elsevier.com/locate/physletb Measurement of the B 0s → J /ψ η lifetime The LHCb Collaboration a r t i c l e i n f o Article history: Received 21 July 2016 Received in revised form 30 September 2016 Accepted October 2016 Available online 11 October 2016 Editor: M Doser This paper is dedicated to the memory of our friend and colleague Ailsa Sparkes a b s t r a c t Using a data set corresponding to an integrated luminosity of fb−1 , collected by the LHCb experiment in pp collisions at centre-of-mass energies of and TeV, the effective lifetime in the B 0s → J /ψ η decay mode, τeff , is measured to be τeff = 1.479 ± 0.034 (stat) ± 0.011 (syst) ps Assuming CP conservation, τeff corresponds to the lifetime of the light B 0s mass eigenstate This is the first measurement of the effective lifetime in this decay mode © 2016 The Author Published by Elsevier B.V This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/) Funded by SCOAP3 Introduction Studies of B 0s –B s mixing provide important tests of the Standard Model (SM) of particle physics In the SM, mixing occurs via box diagrams Extensions to the SM may introduce additional CP-violating phases that alter the value of the B 0s –B 0s mixing weak phase, φs , from that of the SM [1] The B 0s system exhibits a sizeable difference in the decay widths L and H , where L and H refer to the light and heavy B 0s mass eigenstates, respectively The effective lifetime, τeff , of a B 0s meson decay mode is measured by approximating the decay time distribution, determined in the B 0s rest system, by a single exponential function For final states that can be accessed by both B 0s and B 0s mesons the effective lifetime depends on their CP components and is also sensitive to φs [2,3] In this analysis τeff is determined for the CP-even B 0s → J /ψ η decay mode As φs is measured to be small [4,5] the mass eigenstates are also CP eigenstates to better than per mille level and τeff measured in B 0s → J /ψ η decays is equal, to good approximation, to the lifetime of the light B 0s mass eigenstate, τL ∝ L−1 In the SM τL is predicted to be 1.43 ± 0.03 ps [6] Measurements − of τL have previously been reported by LHCb in the B 0s → D + s Ds and B 0s → K + K − decay modes [7,8] The latter is dominated by penguin diagrams, which could arise within and beyond the SM and gives rise to direct CP violation This then leads to a differ− ent τeff , when compared to measurements in the B 0s → D + s Ds and B 0s → J /ψ η decays which are mediated by tree diagrams Improved precision on the effective lifetimes τL and τH will enable more stringent tests of the consistency between direct measurements of the decay width difference s = L − H measured in B 0s → J /ψφ decays and those inferred using effective lifetimes The measurement of the effective B 0s → J /ψ η lifetime presented in this Letter uses fb−1 of data collected in pp collisions at centre-of-mass energies of TeV and TeV during 2011 and 2012 using the LHCb detector The J /ψ meson is reconstructed via the dimuon decay mode and the η meson via the diphoton decay mode The presence of only two charged particles in the final state minimizes systematic uncertainties related to the tracking system Detector and simulation The LHCb detector [9,10] is a single-arm forward spectrometer covering the pseudorapidity range < η < 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 (TT) located upstream of a dipole magnet with a bending power of about Tm, and three stations of siliconstrip detectors and straw drift tubes placed downstream of the magnet The tracking system provides a measurement of momentum, p, of charged particles with a relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200 GeV/c Large samples of J /ψ → μ+ μ− and B + → J /ψ K + decays, collected concurrently with the data set used here, were used to calibrate the momentum scale of the spectrometer to a precision of 0.03% [11] 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 Photons, electrons and hadrons are identified by a calorimeter system consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter The calorimeter response is calibrated using samples of π → γ γ decays For this http://dx.doi.org/10.1016/j.physletb.2016.10.006 0370-2693/© 2016 The Author Published by Elsevier B.V This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/) Funded by SCOAP3 The LHCb Collaboration / Physics Letters B 762 (2016) 484–492 analysis a further calibration was made using the decay η → γ γ , which results in a precision of 0.07% on the neutral energy scale 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 [12], which consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, where a full event reconstruction is made Candidate events are required to pass the hardware trigger, which selects muon and dimuon candidates with high p T based upon muon system information The subsequent software trigger is composed of two stages The first performs a partial event reconstruction and requires events to have two well-identified oppositely charged muons with an invariant mass larger than 2.7 GeV/c The second stage performs a full event reconstruction Events are retained for further processing if they contain a displaced J /ψ → μ+ μ− candidate The decay vertex is required to be well separated from each reconstructed PV of the proton–proton interaction by requiring the distance between the PV and the J /ψ decay vertex divided by its uncertainty to be greater than three This introduces a nonuniform efficiency for b-hadron candidates that have a decay time less than 0.1 ps Simulated pp collisions are generated using Pythia [13] with a specific LHCb configuration [14] Decays of hadronic particles are described by EvtGen [15], in which final-state radiation is generated using Photos [16] The interaction of the generated particles with the detector, and its response, are implemented using the Geant4 toolkit [17] as described in Ref [18] Selection A two-step procedure, is used to optimize the selection of B 0s → J /ψ η decay candidates These studies use simulated data samples together with the high mass sideband of the data (5650 < m( J /ψ η) < 5850 MeV/c ), which is not used in the subsequent determination of τeff In a first step, loose selection criteria are applied that reduce background significantly whilst retaining high signal efficiency Subsequently, a multivariate selection (MVA) is used to reduce further the combinatorial background This is optimized using pseudoexperiments to obtain the best precision on the measured B 0s lifetime The selection starts from a pair of oppositely charged particles, identified as muons, that form a common decay vertex Combina2 of the muon torial background is suppressed by requiring that χIP candidates to all reconstructed PVs to be larger than four To ensure a high reconstruction efficiency the muon candidates are required to have a pseudorapidity between 2.0 and 4.5 The invariant mass of the dimuon candidate must be within 50 MeV/c of the known J /ψ mass [19] In addition, the trigger requirement that the J /ψ decay length divided by its uncertainty is greater than three is reapplied Photons are selected from neutral clusters reconstructed in the electromagnetic calorimeter [10] that have a transverse energy in excess of 300 MeV and a confidence level to be a photon, Pγ , greater than 0.009 The latter requirement has an efficiency of 98% for the simulated signal sample whilst removing 23% of the background in the high mass sideband To suppress combinatorial background, if either of the photons when combined with any other photon candidate in the event has an invariant mass within 25 MeV/c of the known π meson mass [19] the candidate is rejected The quantity χIP is defined as the difference between the structed with and without the considered particle χ of the PV recon- 485 Candidate η → γ γ decays are selected from diphoton combinations with an invariant mass within 70 MeV/c of the known η mass [19] and with a transverse momentum larger than GeV/c The decay angle between the photon momentum in the η rest frame and the direction of Lorentz boost from the laboratory frame to the η rest frame, θη∗ , is required to satisfy cos θη∗ < 0.8 The J /ψ and η candidates are combined to form candidate B 0(s) mesons The average number of PVs in each event is 1.8 (2.0) for the 2011 (2012) dataset When multiple PVs are reconstructed, the one with the minimum χIP to the B 0(s) candidate is chosen A kinematic fit is performed to improve the invariant mass resolution [20] In this fit the momentum vector of the B 0(s) candidate is constrained to point to the PV and the intermediate resonance masses are constrained to their known values The reduced χ of this fit, χ /ndf, is required to be less than five The measured B 0(s) decay time must be larger than 0.3 ps and less than 10 ps If more than one PV is reconstructed in an event the properties of the unassociated vertices are studied Any candidate for which there is a second PV with χIP < 50 is rejected This requirement has an efficiency of 98% that is almost flat as a function of the decay time and reduces background due to incorrect association of the B 0(s) candidate to a PV Finally, as in Ref [21], the position of the PV along the beam-line is required to be within 10 cm of the nominal interaction point, where the standard deviation of this variable is approximately cm This criterion leads to a 10% reduction in signal yield but defines a fiducial region where the reconstruction efficiency is uniform The second step of the selection process is based on a neural network [22], which is trained using the simulated signal sample and the high-mass sideband of the data for background Seven variables that show good agreement between data and simulation and that not significantly bias the B 0(s) decay time distribution are used to train the neural net: the χ /ndf of the kinematic fit; the p T of the B 0(s) and η mesons; the minimum p T of the two pho- tons; cos θη∗ ; the minimum Pγ of the two photons and the total hit multiplicity in the TT sub-detector The requirement on the MVA output was chosen to minimize the statistical uncertainty on the fitted τeff using a sample of 100 pseudoexperiments The chosen value removes 94% of background candidates whilst retaining 69% of the signal candidates After applying these requirements 2% of events contain multiple candidates from which only one, chosen at random, is kept Fit model The effective lifetime is determined by performing a twodimensional maximum likelihood fit to the unbinned distributions of the B 0(s) candidate invariant mass and decay time t= m·l p , where l is the candidate decay length, p the candidate momentum and m the reconstructed invariant mass of the candidate The fit is performed for candidates with 5050 < m( J /ψ η) < 5650 MeV/c and 0.3 < t < 10 ps The fit model has four components: the B 0s → J /ψ η signal, background from the B → J /ψ η decay, background from partially reconstructed B 0s → J /ψ η X decays, and combinatorial background In the fit, the decay-time distribution of each component is convolved with a Gaussian resolution function whose width is fixed to the standard deviation of the decay-time resolution in simulated data A decay-time acceptance function accounts for the dependence of the signal efficiency on several effects The procedure 486 The LHCb Collaboration / Physics Letters B 762 (2016) 484–492 Table Acceptance parameters due to the selection requirements ( A sel ) The correlation coefficients are ρc0 c1 = 0.51, ρc0 c2 = 0.62 and ρc1 c2 = 0.95 Parameter Value c0 (6.5 ± 0.4) × 10−3 ps−1 (6.6 ± 0.3) ps−1 1.50 ± 0.04 c1 c2 Table Values of the β and γ factor fitting the quadratic form The first uncertainty is statistical and the second from the propagation of the uncertainty on the efficiency versus the distance of closest approach obtained with the B + → J /ψ K + calibration sample The correlation coefficient between β and γ is 0.8 Sample β [%] γ [%] 2011 data 0.01 0.39 ± 0.06− +0.07 0.004 0.115 ± 0.021− +0.001 2012 0.001 0.93 ± 0.080+ −0.01 0.006 0.051 ± 0.023− +0.006 used to model the decay-time acceptance is described in detail in Ref [21] The overall acceptance, A tot , is factorised into the product of the selection ( A sel ), trigger ( A trig ) and vertex ( A β ) acceptance functions, determined as described below The effect of the selection requirements, dominated by the cut on the displacement of the muons from the PV, is studied using simulation and parameterised with the form A sel = − c0t + (c t )−c2 , where t is the decay time, and c , c and c are parameters determined from the simulation and summarized in Table In the second level of the software trigger a cut is applied on the decay length significance of the J /ψ candidate, which biases the decay time distribution The trigger efficiency, A trig , is measured separately for the 2011 and 2012 datasets using events that are selected by a dedicated prescaled trigger in which this requirement is removed It increases approximately linearly from 98% at t = 0.3 ps to 100% ps The resulting acceptance shape is parameterised in bins of decay time with linear interpolation between the bins Finally, the reconstruction efficiency of the vertex detector decreases as the distance of closest approach of the decay products to the pp beam-line increases This effect is studied using B + → J /ψ K + decays where the kaon is reconstructed without using vertex detector information [21] and parameterised with the form Aβ = − βt − γ t2, where the parameters β and γ are determined separately for the 2011 and 2012 data The obtained values are summarized in Table Fig shows the overall acceptance curves obtained for the 2011 and 2012 datasets The shape of A tot is mainly determined by A sel , whose uncertainty is dominated by the size of the simulation sample The overall acceptance correction is relatively small Fitting the simulated data with and without the correction τeff changes by 13 fs The invariant mass distribution for the B 0s → J /ψ η signal is parameterised by a Student’s t-distribution The Bukin [23] and JohnsonSU [24] functions are considered for systematic variations In the fit to the data, the shape parameters of this distribution are fixed to the simulation values The decay time distribution for this component is modelled with an exponential function convolved with the detector resolution and multiplied by the detector acceptance, as discussed above Fig Total acceptance function, A tot for 2011 data (black dashed line) and 2012 data (solid red) The second component in the fit accounts for the B → J /ψ η decay As the invariant mass resolution is approximately 48 MeV/c this overlaps with the B 0s signal mode Its mass distribution is modelled, analogously to the B 0s component, with a Student’s t-distribution, with resolution parameters fixed to values determined in the simulation The mass difference between the B 0s and B mesons, and the B lifetime, are fixed to their known central values: m( B 0s ) − m( B ) = 87.29 ± 0.26 MeV/c [25] and τ ( B ) = 1.519 ± 0.005 ps [19] and the uncertainty propagated to the systematic error Similarly, the relative yield of the B and B 0s components, f r , is fixed to (7.3 ± 0.8)% calculated from the average of the branching fractions measurements made by the Belle [26, 27] and LHCb collaborations [28], and the measured fragmentation fractions [29–31] Combinatorial background is modelled by a first order Chebyshev polynomial in mass and the sum of two exponentials in decay time In the fit to the data the lifetime of the shorter lived component is fixed to the value found in the fit to the sideband As a systematic variation of the mass model, an exponential function is considered Background from partially reconstructed decays of b hadrons is studied using a simulated bb sample Using this sample an additional background component, due to partially reconstructed B 0s → J /ψ η X decays, is identified Background from this source lies at invariant masses below 5100 MeV/c and has a lifetime of 1.33 ± 0.10 ps This component is modelled by a Novosibirsk function [32] in mass and an exponential in time All parameters for this component apart from the yield are fixed to the simulation values in the fit to the data The fit has eight free parameters: the yield of the B 0s → J /ψ η component (N B s ), the combinatorial background yield (N comb ), the partially reconstructed background yield (N partial ), the B 0s mass, the lifetime of the signal component (τeff ), the coefficient of the combinatorial background component in mass (acomb ), the longer lived background lifetime (τcomb ) and the fraction of the shortlived background ( f comb ) Independent fits are performed for the 2011 and 2012 data and a weighted average of the two lifetime values is made The correctness of the fit procedure is validated using the full simulation and pseudoexperiments No significant bias is found and the uncertainties estimated by the fit are found to be accurate Results Fig shows the fit projections in mass and decay time for the 2011 and 2012 data The corresponding fit results are summarized in Table The fitted signal yields of the two years scale according to the known integrated luminosity and b-hadron production cross-section There is some tension in the relative yield of the partially reconstructed background between the two years However, this parameter is almost uncorrelated with τeff and this tension The LHCb Collaboration / Physics Letters B 762 (2016) 484–492 487 Fig Mass and decay time distributions for the 2011 dataset (top row) and 2012 dataset (bottom row) The fit model described in the text is superimposed (red line) The partially reconstructed component is shown in solid yellow (dark grey), the combinatorial background in solid green (light grey) and the B component as open blue The pull, i.e the difference between the observed and fitted value divided by the uncertainty, is shown below each of the plots (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.) Table Parameters of the fit to B 0(s) → J /ψ η candidates for the 2011 and 2012 datasets Uncertainties are statistical only Fit parameter Fitted value 2011 N B 0s m B [MeV/c ] s τeff [ps] N comb N partial 960 ± 42 5365.6 ± 1.8 1.485 ± 0.060 1898 ± 64 81 ± 26 2012 2061 ± 60 5369.6 ± 1.3 1.476 ± 0.041 3643 ± 89 345 ± 39 −0.37 ± 0.05 −0.31 ± 0.03 f comb 0.52 ± 0.03 0.49 ± 0.02 τcomb [ps] 0.97 ± 0.06 0.82 ± 0.04 acomb has no impact on the result The average of the fitted values of is τeff τeff = 1.479 ± 0.034 ps, where the uncertainty is statistical The main source of systematic uncertainty is due to the modelling of the decay time acceptance function (Section 4) Varying the parameters of the acceptance function within their correlated uncertainties, a variation of the fitted lifetime of 10 fs is found, which is assigned as a systematic uncertainty Uncertainties on A sel due to the parameterisation of this effect are evaluated to be negligible by replacing the functional form with a histogram The statistical and systematic uncertainties on A β are evaluated by repeating the fit and varying the parameterisation within its uncertainties The statistical uncertainty on A trig is propagated by generating an ensemble of histograms with each bin varied within its statistical uncertainty Systematic uncertainties on A trig are estimated to be small by varying the binning of the histogram and considering an alternative analytic form In simulation studies the efficiency of the MVA is found to be independent of the decay time within uncertainties Conservatively, allowing for a linear dependence, an uncertainty of 1.7 fs is assigned The influence of the decay time resolution is estimated by increasing its value from 51 to 70 fs This variation covers the variation of the resolution with decay time and any possible discrepancy in the resolution between data and simulation The change in τeff from this variation is negligible The impact of the uncertainties in f r , the B 0s –B mass splitting, and the B lifetime are evaluated by repeating the fit procedure varying these parameters within their quoted uncertainties Further uncertainties arise from the modelling of the time distributions of the background components In the default fit the lifetime of the short-lived component is fixed to the value found in a fit to the mass sideband Removing this constraint changes the result by fs, which is assigned as a systematic uncertainty The uncertainty due to the fixed lifetime of the partially reconstructed component is found to be negligible Uncertainties arising from the modelling of the signal and background mass distributions are evaluated using the discrete profiling method described in Ref [33] and found to be negligible Further small uncertainties arise due to the limited knowledge of the length scale of the detector along the beam axis (z-scale), the charged particle momentum scale and the neutral particle energy scale The stability of the result has been tested against a number of possible variations, such as the fitted invariant mass range, the requirement on the IP of the muons, the MVA requirement and analysing the sample according to the number of reconstructed PVs No significant change in the final result is found and hence no further systematic uncertainty is assigned All the uncertainties are summarized in Table Adding them in quadrature leads to a total systematic uncertainty of 11.1 fs which 488 The LHCb Collaboration / Physics Letters B 762 (2016) 484–492 Table Systematic uncertainties on the lifetime measurement Uncertainties less than 0.1 fs are indicated by a dash Source Uncertainty [fs] A sel A β (stat) A β (syst) A trig (stat) A trig (syst) MVA Time resolution fr B 0s –B mass difference B lifetime Releasing τback Varying τpartial Mass model Momentum scale z-scale 10.0 2.0 0.1 0.6 0.6 1.7 – 1.2 – 0.2 4.0 – – – 0.3 Total 11.1 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 FANO (Russia); MinECo (Spain); SNSF and SER (Switzerland); NASU (Ukraine); STFC (United Kingdom); NSF (USA) We acknowledge the computing resources that are provided by CERN, IN2P3 (France), KIT and DESY (Germany), INFN (Italy), SURF (The Netherlands), PIC (Spain), GridPP (United Kingdom), RRCKI and Yandex LLC (Russia), CSCS (Switzerland), IFIN-HH (Romania), CBPF (Brazil), PL-GRID (Poland) and OSC (USA) We are indebted to the communities behind the multiple open source software packages on which we depend Individual groups or members have received support from AvH Foundation (Germany), EPLANET, Marie Skłodowska-Curie Actions and ERC (European Union), 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) Appendix A Supplementary material Supplementary material related to this article can be found online at http://dx.doi.org/10.1016/j.physletb.2016.10.006 References [1] LHCb collaboration, R Aaij, et al., A Bharucha, et al., Implications of LHCb measurements and future prospects, Eur Phys J C 73 (2013) 2373, arXiv: 1208.3355 [2] R Fleischer, R Knegjens, Effective lifetimes of B s decays and their constraints Fig Summary of measurements of τL The yellow band corresponds to the 2015 HFAG central value and uncertainty (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.) is dominated by the size of the simulation sample used to determine the acceptance and to validate the analysis procedure on the B 0s –B s mixing parameters, Eur Phys J C 71 (2011) 1789, arXiv: 1109.5115 [3] R Fleischer, R Knegjens, G Ricciardi, Exploring CP violation and η –η mixing with the B 0s,d → J /ψ η( ) systems, Eur Phys J C 71 (2011) 1798, arXiv: 1110.5490 [4] LHCb collaboration, R Aaij, et al., Precision measurement of CP violation in B 0s → J /ψ K + K − decays, Phys Rev Lett 114 (2015) 041801, arXiv:1411.3104 [5] LHCb collaboration, R Aaij, et al., Measurement of the CP-violating phase φs in Summary B s → J /ψ π + π − decays, Phys Lett B 736 (2014) 186, arXiv:1405.4140 [6] A Lenz, Theoretical update of B-mixing and lifetimes, in: 2012 Electroweak Interactions and Unified Theories, Moriond, 2012, arXiv:1205.1444 Using data collected by LHCb, the effective lifetime in the B 0s → J /ψ η decay mode is measured to be τeff = 1.479 ± 0.034 (stat) ± 0.011 (syst) ps In the limit of CP conservation, τeff is equal to the lifetime of the light B 0s mass eigenstate τ L The present measurement is consistent with, and has similar precision to, the effective lifetime de− + − termined using the B 0s → D + s D s decay mode [7], τeff ( D s D s ) = 1.379 ± 0.026 (stat) ± 0.017 (syst) ps and also with the value measured in the B 0s → K + K − mode [8], τeff ( K + K − ) = 1.407 ± 0.016 (stat) ± 0.007 (syst) ps, where penguin diagrams are expected to be more important Averaging the tree level measurements gives τeff = 1.42 ± 0.02 ps in good agreement with the expectations of the Standard Model [6], τL = 1.43 ± 0.03 ps and the value quoted by HFAG [34] from measurements made in the B 0s → J /ψφ mode, τL = 1.420 ± 0.006 ps The values from these different measurements are compared in Fig + [7] LHCb collaboration, R Aaij, et al., Measurement of the B s → D − s D s and B s → D− D+ effective lifetimes, Phys Rev Lett 112 (2014) 111802, arXiv:1312.1217 s [8] LHCb collaboration, R Aaij, et al., Effective lifetime measurements in the B 0s → + − K K , B → K + π − and B 0s → π + K − decays, Phys Lett B 736 (2014) 446, arXiv:1406.7204 [9] LHCb collaboration, A.A Alves Jr., et al., The LHCb detector at the LHC, J Instrum (2008) S08005 [10] LHCb collaboration, R Aaij, et al., LHCb detector performance, Int J Mod Phys A 30 (2015) 1530022, arXiv:1412.6352 [11] LHCb collaboration, R Aaij, et al., Measurements of the b0 , b− , and b− baryon masses, Phys Rev Lett 110 (2013) 182001, arXiv:1302.1072 [12] R Aaij, et al., The LHCb trigger and its performance in 2011, J Instrum (2013) P04022, arXiv:1211.3055 [13] T Sjöstrand, S Mrenna, P Skands, PYTHIA 6.4 physics and manual, J High Energy Phys 05 (2006) 026, arXiv:hep-ph/0603175; T Sjöstrand, S Mrenna, P Skands, A brief introduction to PYTHIA 8.1, Comput Phys Commun 178 (2008) 852, arXiv:0710.3820 [14] I Belyaev, et al., Handling of the generation of primary events in Gauss, the LHCb simulation framework, J Phys Conf Ser 331 (2011) 032047 [15] D.J Lange, The EvtGen particle decay simulation package, Nucl Instrum Methods, Sect A 462 (2001) 152 [16] P Golonka, Z Was, PHOTOS Monte Carlo: a precision tool for QED corrections in Z and W decays, Eur Phys J C 45 (2006) 97, arXiv:hep-ph/0506026 The LHCb Collaboration / Physics Letters B 762 (2016) 484–492 [17] Geant4 collaboration, J Allison, et al., Geant4 developments and applications, IEEE Trans Nucl Sci 53 (2006) 270; Geant4 collaboration, S Agostinelli, et al., Geant4: a simulation toolkit, Nucl Instrum Methods, Sect A 506 (2003) 250 [18] M Clemencic, et al., The LHCb simulation application, Gauss: design, evolution and experience, J Phys Conf Ser 331 (2011) 032023 [19] Particle Data Group, K.A Olive, et al., Review of particle physics, Chin Phys C 38 (2014) 090001, and 2015 update [20] W.D Hulsbergen, Decay chain fitting with a Kalman filter, Nucl Instrum Methods, Sect A 552 (2005) 566, arXiv:physics/0503191 [21] LHCb collaboration, R Aaij, et al., Measurements of the B + , B , B 0s meson and baryon lifetimes, J High Energy Phys 04 (2014) 114, arXiv:1402.2554 b [22] A Hoecker, et al., TMVA: toolkit for multivariate data analysis, PoS ACAT (2007) 040, arXiv:physics/0703039 [23] BABAR collaboration, J.P Lees, et al., Branching fraction measurements of the color-suppressed decays B → D (∗)0 π , D (∗)0 η , D (∗)0 ω , and D (∗)0 η and mea0 surement of the polarization in the decay B → D ∗0 ω , Phys Rev D 84 (2011) 112007, arXiv:1107.5751; Phys Rev D 87 (2013) 039901 (Erratum) [24] N.L Johnson, Systems of frequency curves generated by methods of translation, Biometrika 36 (1–2) (1949) 149 [25] LHCb collaboration, R Aaij, et al., Observation of the decay B s → ψ(2S ) K + π − , Phys Lett B 747 (2015) 484, arXiv:1503.07112 489 [26] Belle collaboration, M.C Chang, et al., Observation of the decay B → J /ψ η , Phys Rev Lett 98 (2007) 131803, arXiv:hep-ex/0609047 [27] Belle collaboration, M.C Chang, et al., Measurement of B → J /ψ η( ) and constraint on the η –η mixing angle, Phys Rev D 85 (2012) 091102, arXiv: 1203.3399 [28] LHCb collaboration, R Aaij, et al., Study of η –η mixing from measurement of B 0(s) → J /ψ η( ) decay rates, J High Energy Phys 01 (2015) 024, arXiv: 1411.0943 [29] LHCb collaboration, R Aaij, et al., Measurement of b hadron production fractions in TeV pp collisions, Phys Rev D 85 (2012) 032008, arXiv:1111.2357 [30] LHCb collaboration, R Aaij, et al., Measurement of the fragmentation fraction ratio f s / f d and its dependence on B meson kinematics, J High Energy Phys 04 (2013) 001, arXiv:1301.5286 [31] LHCb collaboration, Updated average f s / f d b-hadron production fraction ratio for TeV pp collisions, LHCb-CONF-2013-011 [32] Belle collaboration, H Ikeda, et al., A detailed test of the CsI(Tl) calorimeter for BELLE with photon beams of energy between 20-MeV and 5.4-GeV, Nucl Instrum Methods, Sect A 441 (2000) 401 [33] P.D Dauncey, et al., Handling uncertainties in background shapes: the discrete profiling method, J Instrum 10 (2015) P04015, arXiv:1408.6865 [34] Heavy Flavor Averaging Group, Y Amhis, et al., Averages of b-hadron, 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/ LHCb Collaboration R Aaij 39 , B Adeva 38 , M Adinolfi 47 , Z Ajaltouni , S Akar , J Albrecht 10 , F Alessio 39 , M Alexander 52 , S Ali 42 , G Alkhazov 31 , P Alvarez Cartelle 54 , A.A Alves Jr 58 , S Amato , S Amerio 23 , Y Amhis , L An 40 , L Anderlini 18 , G Andreassi 40 , M Andreotti 17,g , J.E Andrews 59 , R.B Appleby 55 , O Aquines Gutierrez 11 , F Archilli , P d’Argent 12 , J Arnau Romeu , A Artamonov 36 , M Artuso 60 , E Aslanides , G Auriemma 26 , M Baalouch , I Babuschkin 55 , S Bachmann 12 , J.J Back 49 , A Badalov 37 , C Baesso 61 , W Baldini 17 , R.J Barlow 55 , C Barschel 39 , S Barsuk , W Barter 39 , V Batozskaya 29 , B Batsukh 60 , V Battista 40 , A Bay 40 , L Beaucourt , J Beddow 52 , F Bedeschi 24 , I Bediaga , L.J Bel 42 , V Bellee 40 , N Belloli 21,i , K Belous 36 , I Belyaev 32 , E Ben-Haim , G Bencivenni 19 , S Benson 39 , J Benton 47 , A Berezhnoy 33 , R Bernet 41 , A Bertolin 23 , F Betti 15 , M.-O Bettler 39 , M van Beuzekom 42 , S Bifani 46 , P Billoir , T Bird 55 , A Birnkraut 10 , A Bitadze 55 , A Bizzeti 18,u , T Blake 49 , F Blanc 40 , J Blouw 11 , S Blusk 60 , V Bocci 26 , T Boettcher 57 , A Bondar 35 , N Bondar 31,39 , W Bonivento 16 , A Borgheresi 21,i , S Borghi 55 , M Borisyak 67 , M Borsato 38 , F Bossu , M Boubdir , T.J.V Bowcock 53 , E Bowen 41 , C Bozzi 17,39 , S Braun 12 , M Britsch 12 , T Britton 60 , J Brodzicka 55 , E Buchanan 47 , C Burr 55 , A Bursche , J Buytaert 39 , S Cadeddu 16 , R Calabrese 17,g , M Calvi 21,i , M Calvo Gomez 37,m , P Campana 19 , D Campora Perez 39 , L Capriotti 55 , A Carbone 15,e , G Carboni 25,j , R Cardinale 20,h , A Cardini 16 , P Carniti 21,i , L Carson 51 , K Carvalho Akiba , G Casse 53 , L Cassina 21,i , L Castillo Garcia 40 , M Cattaneo 39 , Ch Cauet 10 , G Cavallero 20 , R Cenci 24,t , M Charles , Ph Charpentier 39 , G Chatzikonstantinidis 46 , M Chefdeville , S Chen 55 , S.-F Cheung 56 , V Chobanova 38 , M Chrzaszcz 41,27 , X Cid Vidal 38 , G Ciezarek 42 , P.E.L Clarke 51 , M Clemencic 39 , H.V Cliff 48 , J Closier 39 , V Coco 58 , J Cogan , E Cogneras , V Cogoni 16,39,f , L Cojocariu 30 , G Collazuol 23,o , P Collins 39 , A Comerma-Montells 12 , A Contu 39 , A Cook 47 , S Coquereau , G Corti 39 , M Corvo 17,g , C.M Costa Sobral 49 , B Couturier 39 , G.A Cowan 51 , D.C Craik 51 , A Crocombe 49 , M Cruz Torres 61 , S Cunliffe 54 , R Currie 54 , C D’Ambrosio 39 , E Dall’Occo 42 , J Dalseno 47 , P.N.Y David 42 , A Davis 58 , O De Aguiar Francisco , K De Bruyn , S De Capua 55 , M De Cian 12 , J.M De Miranda , L De Paula , M De Serio 14,d , P De Simone 19 , C.-T Dean 52 , D Decamp , M Deckenhoff 10 , L Del Buono , M Demmer 10 , D Derkach 67 , O Deschamps , F Dettori 39 , B Dey 22 , A Di Canto 39 , H Dijkstra 39 , F Dordei 39 , M Dorigo 40 , A Dosil Suárez 38 , A Dovbnya 44 , K Dreimanis 53 , L Dufour 42 , G Dujany 55 , K Dungs 39 , P Durante 39 , R Dzhelyadin 36 , A Dziurda 39 , A Dzyuba 31 , N Déléage , S Easo 50 , U Egede 54 , V Egorychev 32 , S Eidelman 35 , S Eisenhardt 51 , U Eitschberger 10 , R Ekelhof 10 , L Eklund 52 , Ch Elsasser 41 , S Ely 60 , S Esen 12 , H.M Evans 48 , T Evans 56 , A Falabella 15 , N Farley 46 , S Farry 53 , R Fay 53 , D Fazzini 21,i , D Ferguson 51 , V Fernandez Albor 38 , F Ferrari 15,39 , F Ferreira Rodrigues , M Ferro-Luzzi 39 , S Filippov 34 , R.A Fini 14 , M Fiore 17,g , M Fiorini 17,g , M Firlej 28 , C Fitzpatrick 40 , T Fiutowski 28 , F Fleuret 7,b , K Fohl 39 , M Fontana 16 , F Fontanelli 20,h , D.C Forshaw 60 , R Forty 39 , V Franco Lima 53 , M Frank 39 , C Frei 39 , J Fu 22,q , E Furfaro 25,j , C Färber 39 , 490 The LHCb Collaboration / Physics Letters B 762 (2016) 484–492 A Gallas Torreira 38 , D Galli 15,e , S Gallorini 23 , S Gambetta 51 , M Gandelman , P Gandini 56 , Y Gao , J García Pardiđas 38 , J Garra Tico 48 , L Garrido 37 , P.J Garsed 48 , D Gascon 37 , C Gaspar 39 , L Gavardi 10 , G Gazzoni , D Gerick 12 , E Gersabeck 12 , M Gersabeck 55 , T Gershon 49 , Ph Ghez , S Gianì 40 , V Gibson 48 , E Gillies 51 , O.G Girard 40 , L Giubega 30 , K Gizdov 51 , V.V Gligorov , D Golubkov 32 , A Golutvin 54,39 , A Gomes 1,a , I.V Gorelov 33 , C Gotti 21,i , M Grabalosa Gándara , R Graciani Diaz 37 , L.A Granado Cardoso 39 , E Graugés 37 , E Graverini 41 , G Graziani 18 , A Grecu 30 , P Griffith 46 , L Grillo 21 , B.R Gruberg Cazon 56 , O Grünberg 65 , E Gushchin 34 , Yu Guz 36 , T Gys 39 , C Göbel 61 , T Hadavizadeh 56 , C Hadjivasiliou , G Haefeli 40 , C Haen 39 , S.C Haines 48 , S Hall 54 , B Hamilton 59 , X Han 12 , S Hansmann-Menzemer 12 , N Harnew 56 , S.T Harnew 47 , J Harrison 55 , M Hatch 39 , J He 62 , T Head 40 , A Heister , K Hennessy 53 , P Henrard , L Henry , J.A Hernando Morata 38 , E van Herwijnen 39 , M Heß 65 , A Hicheur , D Hill 56 , C Hombach 55 , W Hulsbergen 42 , T Humair 54 , M Hushchyn 67 , N Hussain 56 , D Hutchcroft 53 , M Idzik 28 , P Ilten 57 , R Jacobsson 39 , A Jaeger 12 , J Jalocha 56 , E Jans 42 , A Jawahery 59 , M John 56 , D Johnson 39 , C.R Jones 48 , C Joram 39 , B Jost 39 , N Jurik 60 , S Kandybei 44 , W Kanso , M Karacson 39 , J.M Kariuki 47 , S Karodia 52 , M Kecke 12 , M Kelsey 60 , I.R Kenyon 46 , M Kenzie 39 , T Ketel 43 , E Khairullin 67 , B Khanji 21,39,i , C Khurewathanakul 40 , T Kirn , S Klaver 55 , K Klimaszewski 29 , S Koliiev 45 , M Kolpin 12 , I Komarov 40 , R.F Koopman 43 , P Koppenburg 42 , A Kozachuk 33 , M Kozeiha , L Kravchuk 34 , K Kreplin 12 , M Kreps 49 , P Krokovny 35 , F Kruse 10 , W Krzemien 29 , W Kucewicz 27,l , M Kucharczyk 27 , V Kudryavtsev 35 , A.K Kuonen 40 , K Kurek 29 , T Kvaratskheliya 32,39 , D Lacarrere 39 , G Lafferty 55,39 , A Lai 16 , D Lambert 51 , G Lanfranchi 19 , C Langenbruch , B Langhans 39 , T Latham 49 , C Lazzeroni 46 , R Le Gac , J van Leerdam 42 , J.-P Lees , A Leat 33,39 , J Lefranỗois , R Lefèvre , F Lemaitre 39 , E Lemos Cid 38 , O Leroy , T Lesiak 27 , B Leverington 12 , Y Li , T Likhomanenko 67,66 , R Lindner 39 , C Linn 39 , F Lionetto 41 , B Liu 16 , X Liu , D Loh 49 , I Longstaff 52 , J.H Lopes , D Lucchesi 23,o , M Lucio Martinez 38 , H Luo 51 , A Lupato 23 , E Luppi 17,g , O Lupton 56 , A Lusiani 24 , X Lyu 62 , F Machefert , F Maciuc 30 , O Maev 31 , K Maguire 55 , S Malde 56 , A Malinin 66 , T Maltsev 35 , G Manca , G Mancinelli , P Manning 60 , J Maratas 5,v , J.F Marchand , U Marconi 15 , C Marin Benito 37 , P Marino 24,t , J Marks 12 , G Martellotti 26 , M Martin , M Martinelli 40 , D Martinez Santos 38 , F Martinez Vidal 68 , D Martins Tostes , L.M Massacrier , A Massafferri , R Matev 39 , A Mathad 49 , Z Mathe 39 , C Matteuzzi 21 , A Mauri 41 , B Maurin 40 , A Mazurov 46 , M McCann 54 , J McCarthy 46 , A McNab 55 , R McNulty 13 , B Meadows 58 , F Meier 10 , M Meissner 12 , D Melnychuk 29 , M Merk 42 , A Merli 22,q , E Michielin 23 , D.A Milanes 64 , M.-N Minard , D.S Mitzel 12 , J Molina Rodriguez 61 , I.A Monroy 64 , S Monteil , M Morandin 23 , P Morawski 28 , A Mordà , M.J Morello 24,t , J Moron 28 , A.B Morris 51 , R Mountain 60 , F Muheim 51 , M Mulder 42 , M Mussini 15 , D Müller 55 , J Müller 10 , K Müller 41 , V Müller 10 , P Naik 47 , T Nakada 40 , R Nandakumar 50 , A Nandi 56 , I Nasteva , M Needham 51,∗ , N Neri 22 , S Neubert 12 , N Neufeld 39 , M Neuner 12 , A.D Nguyen 40 , C Nguyen-Mau 40,n , S Nieswand , R Niet 10 , N Nikitin 33 , T Nikodem 12 , A Novoselov 36 , D.P O’Hanlon 49 , A Oblakowska-Mucha 28 , V Obraztsov 36 , S Ogilvy 19 , R Oldeman 48 , C.J.G Onderwater 69 , J.M Otalora Goicochea , A Otto 39 , P Owen 41 , A Oyanguren 68 , P.R Pais 40 , A Palano 14,d , F Palombo 22,q , M Palutan 19 , J Panman 39 , A Papanestis 50 , M Pappagallo 14,d , L.L Pappalardo 17,g , C Pappenheimer 58 , W Parker 59 , C Parkes 55 , G Passaleva 18 , A Pastore 14,d , G.D Patel 53 , M Patel 54 , C Patrignani 15,e , A Pearce 55,50 , A Pellegrino 42 , G Penso 26,k , M Pepe Altarelli 39 , S Perazzini 39 , P Perret , L Pescatore 46 , K Petridis 47 , A Petrolini 20,h , A Petrov 66 , M Petruzzo 22,q , E Picatoste Olloqui 37 , B Pietrzyk , M Pikies 27 , D Pinci 26 , A Pistone 20 , A Piucci 12 , S Playfer 51 , M Plo Casasus 38 , T Poikela 39 , F Polci , A Poluektov 49,35 , I Polyakov 60 , E Polycarpo , G.J Pomery 47 , A Popov 36 , D Popov 11,39 , B Popovici 30 , C Potterat , E Price 47 , J.D Price 53 , J Prisciandaro 38 , A Pritchard 53 , C Prouve 47 , V Pugatch 45 , A Puig Navarro 40 , G Punzi 24,p , W Qian 56 , R Quagliani 7,47 , B Rachwal 27 , J.H Rademacker 47 , M Rama 24 , M Ramos Pernas 38 , M.S Rangel , I Raniuk 44 , G Raven 43 , F Redi 54 , S Reichert 10 , A.C dos Reis , C Remon Alepuz 68 , V Renaudin , S Ricciardi 50 , S Richards 47 , M Rihl 39 , K Rinnert 53,39 , V Rives Molina 37 , P Robbe 7,39 , A.B Rodrigues , E Rodrigues 58 , J.A Rodriguez Lopez 64 , P Rodriguez Perez 55 , A Rogozhnikov 67 , S Roiser 39 , V Romanovskiy 36 , A Romero Vidal 38 , J.W Ronayne 13 , M Rotondo 23 , M.S Rudolph 60 , T Ruf 39 , P Ruiz Valls 68 , J.J Saborido Silva 38 , E Sadykhov 32 , N Sagidova 31 , B Saitta 16,f , V Salustino Guimaraes , C Sanchez Mayordomo 68 , B Sanmartin Sedes 38 , R Santacesaria 26 , C Santamarina Rios 38 , The LHCb Collaboration / Physics Letters B 762 (2016) 484–492 491 M Santimaria 19 , E Santovetti 25,j , A Sarti 19,k , C Satriano 26,s , A Satta 25 , D.M Saunders 47 , D Savrina 32,33 , S Schael , M Schellenberg 10 , M Schiller 39 , H Schindler 39 , M Schlupp 10 , M Schmelling 11 , T Schmelzer 10 , B Schmidt 39 , O Schneider 40 , A Schopper 39 , K Schubert 10 , M Schubiger 40 , M.-H Schune , R Schwemmer 39 , B Sciascia 19 , A Sciubba 26,k , A Semennikov 32 , A Sergi 46 , N Serra 41 , J Serrano , L Sestini 23 , P Seyfert 21 , M Shapkin 36 , I Shapoval 17,44,g , Y Shcheglov 31 , T Shears 53 , L Shekhtman 35 , V Shevchenko 66 , A Shires 10 , B.G Siddi 17 , R Silva Coutinho 41 , L Silva de Oliveira , G Simi 23,o , S Simone 14,d , M Sirendi 48 , N Skidmore 47 , T Skwarnicki 60 , E Smith 54 , I.T Smith 51 , J Smith 48 , M Smith 55 , H Snoek 42 , M.D Sokoloff 58 , F.J.P Soler 52 , D Souza 47 , B Souza De Paula , B Spaan 10 , P Spradlin 52 , S Sridharan 39 , F Stagni 39 , M Stahl 12 , S Stahl 39 , P Stefko 40 , S Stefkova 54 , O Steinkamp 41 , O Stenyakin 36 , S Stevenson 56 , S Stoica 30 , S Stone 60 , B Storaci 41 , S Stracka 24,t , M Straticiuc 30 , U Straumann 41 , L Sun 58 , W Sutcliffe 54 , K Swientek 28 , V Syropoulos 43 , M Szczekowski 29 , T Szumlak 28 , S T’Jampens , A Tayduganov , T Tekampe 10 , G Tellarini 17,g , F Teubert 39 , C Thomas 56 , E Thomas 39 , J van Tilburg 42 , V Tisserand , M Tobin 40 , S Tolk 48 , L Tomassetti 17,g , D Tonelli 39 , S Topp-Joergensen 56 , F Toriello 60 , E Tournefier , S Tourneur 40 , K Trabelsi 40 , M Traill 52 , M.T Tran 40 , M Tresch 41 , A Trisovic 39 , A Tsaregorodtsev , P Tsopelas 42 , A Tully 48 , N Tuning 42 , A Ukleja 29 , A Ustyuzhanin 67,66 , U Uwer 12 , C Vacca 16,39,f , V Vagnoni 15,39 , S Valat 39 , G Valenti 15 , A Vallier , R Vazquez Gomez 19 , P Vazquez Regueiro 38 , S Vecchi 17 , M van Veghel 42 , J.J Velthuis 47 , M Veltri 18,r , G Veneziano 40 , A Venkateswaran 60 , M Vernet , M Vesterinen 12 , B Viaud , D Vieira , M Vieites Diaz 38 , X Vilasis-Cardona 37,m , V Volkov 33 , A Vollhardt 41 , B Voneki 39 , D Voong 47 , A Vorobyev 31 , V Vorobyev 35 , C Voß 65 , J.A de Vries 42 , C Vázquez Sierra 38 , R Waldi 65 , C Wallace 49 , R Wallace 13 , J Walsh 24 , J Wang 60 , D.R Ward 48 , H.M Wark 53 , N.K Watson 46 , D Websdale 54 , A Weiden 41 , M Whitehead 39 , J Wicht 49 , G Wilkinson 56,39 , M Wilkinson 60 , M Williams 39 , M.P Williams 46 , M Williams 57 , T Williams 46 , F.F Wilson 50 , J Wimberley 59 , J Wishahi 10 , W Wislicki 29 , M Witek 27 , G Wormser , S.A Wotton 48 , K Wraight 52 , S Wright 48 , K Wyllie 39 , Y Xie 63 , Z Xing 60 , Z Xu 40 , Z Yang , H Yin 63 , J Yu 63 , X Yuan 35 , O Yushchenko 36 , M Zangoli 15 , K.A Zarebski 46 , M Zavertyaev 11,c , L Zhang , Y Zhang , Y Zhang 62 , A Zhelezov 12 , Y Zheng 62 , A Zhokhov 32 , V Zhukov , S Zucchelli 15 Centro Brasileiro de Pesquisas Físicas (CBPF), Rio de Janeiro, Brazil Universidade Federal Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil Center for High Energy Physics, Tsinghua University, Beijing, China LAPP, Université Savoie Mont-Blanc, CNRS/IN2P3, Annecy-Le-Vieux, France Clermont Université, Université Blaise Pascal, CNRS/IN2P3, LPC, Clermont-Ferrand, France CPPM, Aix-Marseille Université, CNRS/IN2P3, Marseille, France LAL, Université Paris-Sud, CNRS/IN2P3, Orsay, France LPNHE, Université Pierre et Marie Curie, Université Paris Diderot, CNRS/IN2P3, Paris, France I Physikalisches Institut, RWTH Aachen University, Aachen, Germany 10 Fakultät Physik, Technische Universität Dortmund, Dortmund, Germany 11 Max-Planck-Institut für Kernphysik (MPIK), Heidelberg, Germany 12 Physikalisches Institut, Ruprecht-Karls-Universität Heidelberg, Heidelberg, Germany 13 School of Physics, University College Dublin, Dublin, Ireland 14 Sezione INFN di Bari, Bari, Italy 15 Sezione INFN di Bologna, Bologna, Italy 16 Sezione INFN di Cagliari, Cagliari, Italy 17 Sezione INFN di Ferrara, Ferrara, Italy 18 Sezione INFN di Firenze, Firenze, Italy 19 Laboratori Nazionali dell’INFN di Frascati, Frascati, Italy 20 Sezione INFN di Genova, Genova, Italy 21 Sezione INFN di Milano Bicocca, Milano, Italy 22 Sezione INFN di Milano, Milano, Italy 23 Sezione INFN di Padova, Padova, Italy 24 Sezione INFN di Pisa, Pisa, Italy 25 Sezione INFN di Roma Tor Vergata, Roma, Italy 26 Sezione INFN di Roma La Sapienza, Roma, Italy 27 Henryk Niewodniczanski Institute of Nuclear Physics Polish Academy of Sciences, Kraków, Poland 28 AGH – University of Science and Technology, Faculty of Physics and Applied Computer Science, Kraków, Poland 29 National Center for Nuclear Research (NCBJ), Warsaw, Poland 30 Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest-Magurele, Romania 31 Petersburg Nuclear Physics Institute (PNPI), Gatchina, Russia 32 Institute of Theoretical and Experimental Physics (ITEP), Moscow, Russia 33 Institute of Nuclear Physics, Moscow State University (SINP MSU), Moscow, Russia 34 Institute for Nuclear Research of the Russian Academy of Sciences (INR RAN), Moscow, Russia 35 Budker Institute of Nuclear Physics (SB RAS) and Novosibirsk State University, Novosibirsk, Russia 36 Institute for High Energy Physics (IHEP), Protvino, Russia 492 37 The LHCb Collaboration / Physics Letters B 762 (2016) 484–492 ICCUB, Universitat de Barcelona, Barcelona, Spain Universidad de Santiago de Compostela, Santiago de Compostela, Spain European Organization for Nuclear Research (CERN), Geneva, Switzerland Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland Physik-Institut, Universität Zürich, Zürich, Switzerland Nikhef National Institute for Subatomic Physics, Amsterdam, The Netherlands Nikhef National Institute for Subatomic Physics and VU University Amsterdam, Amsterdam, The Netherlands NSC Kharkiv Institute of Physics and Technology (NSC KIPT), Kharkiv, Ukraine Institute for Nuclear Research of the National Academy of Sciences (KINR), Kyiv, Ukraine University of Birmingham, Birmingham, United Kingdom H.H Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom Department of Physics, University of Warwick, Coventry, United Kingdom STFC Rutherford Appleton Laboratory, Didcot, United Kingdom School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom School of Physics and Astronomy, University of Glasgow, Glasgow, United Kingdom Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom Imperial College London, London, United Kingdom School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom Department of Physics, University of Oxford, Oxford, United Kingdom Massachusetts Institute of Technology, Cambridge, MA, United States University of Cincinnati, Cincinnati, OH, United States University of Maryland, College Park, MD, United States Syracuse University, Syracuse, NY, United States Pontifícia Universidade Católica Rio de Janeiro (PUC-Rio), Rio de Janeiro, Brazil w University of Chinese Academy of Sciences, Beijing, China x Institute of Particle Physics, Central China Normal University, Wuhan, Hubei, China x Departamento de Fisica, Universidad Nacional de Colombia, Bogota, Colombia y Institut für Physik, Universität Rostock, Rostock, Germany z National Research Centre Kurchatov Institute, Moscow, Russia aa Yandex School of Data Analysis, Moscow, Russia aa Instituto de Fisica Corpuscular (IFIC), Universitat de Valencia-CSIC, Valencia, Spain ab Van Swinderen Institute, University of Groningen, Groningen, The Netherlands ac 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 * a b c d e f Corresponding author E-mail address: matthew.needham@cern.ch (M Needham) Universidade Federal Triângulo Mineiro (UFTM), Uberaba-MG, Brazil Laboratoire Leprince-Ringuet, Palaiseau, France P.N Lebedev Physical Institute, Russian Academy of Science (LPI RAS), Moscow, Russia Università di Bari, Bari, Italy Università di Bologna, Bologna, Italy g Università di Cagliari, Cagliari, Italy 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 LIFAELS, La Salle, Universitat Ramon Llull, Barcelona, Spain Hanoi University of Science, Hanoi, Viet Nam Università di Padova, Padova, Italy Università di Pisa, Pisa, Italy Università degli Studi di Milano, Milano, Italy Università di Urbino, Urbino, Italy Università della Basilicata, Potenza, Italy Scuola Normale Superiore, Pisa, Italy Università di Modena e Reggio Emilia, Modena, Italy Iligan Institute of Technology (IIT), Iligan, Philippines Associated to Universidade Federal Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil Associated to Center for High Energy Physics, Tsinghua University, Beijing, China Associated to LPNHE, Université Pierre et Marie Curie, Université Paris Diderot, CNRS/IN2P3, Paris, France Associated to Physikalisches Institut, Ruprecht-Karls-Universität Heidelberg, Heidelberg, Germany Associated to Institute of Theoretical and Experimental Physics (ITEP), Moscow, Russia m n o p q r s t u v w x y z aa ab ac Associated to ICCUB, Universitat de Barcelona, Barcelona, Spain Associated to Nikhef National Institute for Subatomic Physics, Amsterdam, The Netherlands ... the sum of two exponentials in decay time In the fit to the data the lifetime of the shorter lived component is fixed to the value found in the fit to the sideband As a systematic variation of the. .. p T of the B 0(s) and η mesons; the minimum p T of the two pho- tons; cos θη∗ ; the minimum Pγ of the two photons and the total hit multiplicity in the TT sub-detector The requirement on the. .. Table Values of the β and γ factor fitting the quadratic form The first uncertainty is statistical and the second from the propagation of the uncertainty on the efficiency versus the distance of closest