# Measurements of the $B^{+}$, $B^{0}$, $B_{s}^{0}$ meson and $\Lambda b^{0}$ baryon lifetimes

R. Aaij, B. Adeva, M. Adinolfi, A. Affolder, Ziad Zj Ajaltouni, J. Albrecht, F.

Alessio, Alexander M. Hoole, S. Ali, G. Alkhazov, et al.

► To cite this version:

R. Aaij, B. Adeva, M. Adinolfi, A. Affolder, Ziad Zj Ajaltouni, et al. Measurements of the $B^{+}$, $B^{0}$, $B_{s}^{0}$ meson and $\Lambda b^{0}$ baryon lifetimes. Journal of High Energy Physics, 2014, 04, pp.114. 10.1007/JHEP04(2014)114. in2p3-00945236

HAL Id: in2p3-00945236

https://hal.in2p3.fr/in2p3-00945236

Submitted on 18 Jul 2023

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
# Measurements of the $B^{+}$, $B^{0}$, $B_{s}^{0}$ meson and $\Lambda_{b}^{0}$ baryon lifetimes

The LHCb collaboration†

## Abstract

Measurements of b-hadron lifetimes are reported using pp collision data, corresponding to an integrated luminosity of 1.0 fb⁻¹, collected by the LHCb detector at a centre-of-mass energy of 7 TeV. Using the exclusive decays $B^{+} \to J/\psi K^{+}$, $B^{0} \to J/\psi K^{*}(892)^{0}$, $B^{0} \to J/\psi K_{S}^{0}$, $\Lambda_{b}^{0} \to J/\psi \Lambda$ and $B_{s}^{0} \to J/\psi \phi$ the average decay times in these modes are measured to be

$$ \tau_{B^{+} \to J/\psi K^{+}} = 1.637 \pm 0.004 \pm 0.003 \text{ ps,} $$

$$ \tau_{B^{0} \to J/\psi K^{*0}} = 1.524 \pm 0.006 \pm 0.004 \text{ ps,} $$

$$ \tau_{B^{0} \to J/\psi K_{S}^{0}} = 1.499 \pm 0.013 \pm 0.005 \text{ ps,} $$

$$ \tau_{\Lambda_{b}^{0} \to J/\psi \Lambda} = 1.415 \pm 0.027 \pm 0.006 \text{ ps,} $$

$$ \tau_{B_{s}^{0} \to J/\psi \phi} = 1.480 \pm 0.011 \pm 0.005 \text{ ps,} $$

where the first uncertainty is statistical and the second is systematic. These represent the most precise lifetime measurements in these decay modes. In addition, ratios of these lifetimes, and the ratio of the decay-width difference, $\Delta\Gamma_{d}$, to the average width, $\Gamma_{d}$, in the $B^{0}$ system, $\Delta\Gamma_{d}/\Gamma_{d} = -0.044 \pm 0.025 \pm 0.011$, are reported. All quantities are found to be consistent with Standard Model expectations.

Submitted to JHEP

© CERN on behalf of the LHCb collaboration, license CC-BY-3.0.

<footnote>†Authors are listed on the following pages.</footnote>
LHCb collaboration

R. Aaij$^{40}$, B. Adeva$^{36}$, M. Adinolfi$^{45}$, A. Affolder$^{51}$, Z. Ajaltouni$^{5}$, J. Albrecht$^{9}$, F. Alessio$^{37}$,
M. Alexander$^{50}$, S. Ali$^{40}$, G. Alkhazov$^{29}$, P. Alvarez Cartelle$^{36}$, A.A. Alves Jr$^{24}$, S. Amato$^{2}$,
S. Amerio$^{21}$, Y. Amhis$^{7}$, L. Anderlini$^{17,9}$, J. Anderson$^{39}$, R. Andreassen$^{56}$, M. Andreotti$^{16,7}$,
J.E. Andrews$^{57}$, R.B. Appleby$^{53}$, O. Aquines Gutierrez$^{10}$, F. Archilli$^{37}$, A. Artamonov$^{34}$,
M. Artuso$^{58}$, E. Aslanides$^{6}$, G. Auriemma$^{24,n}$, M. Baalouch$^{5}$, S. Bachmann$^{11}$, J.J. Back$^{47}$,
A. Badalov$^{35}$, V. Balagura$^{30}$, W. Baldini$^{16}$, R.J. Barlow$^{53}$, C. Barschel$^{38}$, S. Barsuk$^{7}$,
W. Barter$^{46}$, V. Batozskaya$^{27}$, Th. Bauer$^{40}$, A. Bay$^{38}$, J. Beddow$^{50}$, F. Bedeschi$^{22}$, I. Bediaga$^{1}$,
S. Belogurov$^{30}$, K. Belous$^{34}$, I. Belyaev$^{30}$, E. Ben-Haim$^{8}$, G. Bencivenni$^{18}$, S. Benson$^{49}$,
J. Benton$^{45}$, A. Berezhnoy$^{31}$, R. Bernet$^{39}$, M.-O. Bettler$^{46}$, M. van Beuzekom$^{40}$, A. Bien$^{11}$,
S. Bifani$^{44}$, T. Bird$^{53}$, A. Bizzeti$^{17,i}$, P.M. Bjørnstad$^{53}$, T. Blake$^{47}$, F. Blanc$^{38}$, J. Blouw$^{10}$,
S. Blusk$^{58}$, V. Bocci$^{24}$, A. Bondar$^{33}$, N. Bondar$^{29}$, W. Bonivento$^{15,37}$, S. Borghi$^{53}$, A. Borgia$^{58}$,
M. Borsato$^{7}$, T.J.V. Bowcock$^{51}$, E. Bowen$^{39}$, C. Bozzi$^{16}$, T. Brambach$^{9}$, J. van den Brand$^{41}$,
J. Bressieux$^{38}$, D. Brett$^{53}$, M. Britsch$^{10}$, T. Britton$^{58}$, N.H. Brook$^{45}$, H. Brown$^{51}$, A. Bursche$^{39}$,
G. Busetto$^{21,r}$, J. Buytaert$^{37}$, S. Cadeddu$^{15}$, R. Calabrese$^{16,f}$, O. Callot$^{7}$, M. Calvi$^{20,k}$,
M. Calvo Gomez$^{35,p}$, A. Camboni$^{35}$, P. Campana$^{18,37}$, D. Campora Perez$^{37}$, A. Carbone$^{14,d}$,
G. Carboni$^{23,l}$, R. Cardinale$^{19,j}$, A. Cardini$^{15}$, H. Carranza-Mejia$^{49}$, L. Carson$^{49}$,
K. Carvalho Akiba$^{2}$, G. Casse$^{51}$, L. Castillo Garcia$^{37}$, M. Cattaneo$^{37}$, Ch. Cauet$^{9}$, R. Cenci$^{57}$,
M. Charles$^{8}$, Ph. Charpentier$^{37}$, S.-F. Cheung$^{54}$, N. Chiapolini$^{39}$, M. Chrzaszcz$^{39,25}$, K. Ciba$^{37}$,
X. Cid Vidal$^{37}$, G. Ciezarek$^{52}$, P.E.L. Clarke$^{49}$, M. Clemencic$^{37}$, H.V. Cliff$^{46}$, J. Closier$^{37}$,
C. Coca$^{28}$, V. Coco$^{37}$, J. Cogan$^{6}$, E. Cogneras$^{5}$, P. Collins$^{37}$, A. Comerma-Montells$^{35}$,
A. Contu$^{15,37}$, A. Cook$^{45}$, M. Coombes$^{45}$, S. Coquereau$^{8}$, G. Corti$^{37}$, I. Counts$^{55}$,
B. Couturier$^{37}$, G.A. Cowan$^{49}$, D.C. Craik$^{47}$, M. Cruz Torres$^{59}$, S. Cunliffe$^{52}$, R. Currie$^{49}$,
C. D'Ambrosio$^{37}$, J. Dalseno$^{45}$, P. David$^{8}$, P.N.Y. David$^{40}$, A. Davis$^{56}$, I. De Bonis$^{4}$,
K. De Bruyn$^{40}$, S. De Capua$^{53}$, M. De Cian$^{11}$, J.M. De Miranda$^{1}$, L. De Paula$^{2}$, W. De Silva$^{56}$,
P. De Simone$^{18}$, D. Decamp$^{4}$, M. Deckenhoff$^{9}$, L. Del Buono$^{8}$, N. Déléage$^{4}$, D. Derkach$^{54}$,
O. Deschamps$^{5}$, F. Dettori$^{41}$, A. Di Canto$^{11}$, H. Dijkstra$^{37}$, S. Donleavy$^{51}$, F. Dordei$^{11}$,
M. Dorigo$^{38}$, P. Dorosz$^{25,o}$, A. Dosil Suárez$^{36}$, D. Dossett$^{47}$, A. Dovbnya$^{42}$, F. Dupertuis$^{38}$,
P. Durante$^{37}$, R. Dzhelyadin$^{34}$, A. Dziurda$^{25}$, A. Dzyuba$^{29}$, S. Easo$^{48}$, U. Egede$^{52}$,
V. Egorychev$^{30}$, S. Eidelman$^{33}$, S. Eisenhardt$^{49}$, U. Eitschberger$^{9}$, R. Ekelhof$^{9}$, L. Eklund$^{50,37}$,
I. El Rifai$^{5}$, Ch. Elsasser$^{39}$, S. Esen$^{11}$, A. Falabella$^{16,f}$, C. Färber$^{11}$, C. Farinelli$^{40}$, S. Farry$^{51}$,
D. Ferguson$^{49}$, V. Fernandez Albor$^{36}$, F. Ferreira Rodrigues$^{1}$, M. Ferro-Luzzi$^{37}$, S. Filippov$^{32}$,
M. Fiore$^{16,f}$, M. Fiorini$^{16,f}$, C. Fitzpatrick$^{37}$, M. Fontana$^{10}$, F. Fontanelli$^{19,j}$, R. Forty$^{37}$,
O. Francisco$^{2}$, M. Frank$^{37}$, C. Frei$^{37}$, M. Frosini$^{17,37,g}$, E. Furfaro$^{23,l}$, A. Gallas Torreira$^{36}$,
D. Galli$^{14,d}$, M. Gandelman$^{2}$, P. Gandini$^{58}$, Y. Gao$^{3}$, J. Garofoli$^{58}$, J. Garra Tico$^{46}$,
L. Garrido$^{35}$, C. Gaspar$^{37}$, R. Gauld$^{54}$, E. Gersabeck$^{11}$, M. Gersabeck$^{53}$, T. Gershon$^{47}$,
Ph. Ghez$^{4}$, A. Gianelle$^{21}$, S. Giani$^{38}$, V. Gibson$^{46}$, L. Giubega$^{28}$, V.V. Gligorov$^{37}$, C. Göbel$^{59}$,
D. Golubkov$^{30}$, A. Golutvin$^{52,30,37}$, A. Gomes$^{1,a}$, H. Gordon$^{37}$, M. Grabalosa Gándara$^{5}$,
R. Graciani Diaz$^{35}$, L.A. Granado Cardoso$^{37}$, E. Graugés$^{35}$, G. Graziani$^{17}$, A. Grecu$^{28}$,
E. Greening$^{54}$, S. Gregson$^{46}$, P. Griffith$^{44}$, L. Grillo$^{11}$, O. Grünberg$^{60}$, B. Gui$^{58}$, E. Gushchin$^{32}$,
Yu. Guz$^{34,37}$, T. Gys$^{37}$, C. Hadjivasiliou$^{58}$, G. Haefeli$^{38}$, C. Haen$^{37}$, T.W. Hafkenscheid$^{62}$,
S.C. Haines$^{46}$, S. Hall$^{52}$, B. Hamilton$^{57}$, T. Hampson$^{45}$, S. Hansmann-Menzemer$^{11}$, N. Harnew$^{54}$,
S.T. Harnew$^{45}$, J. Harrison$^{53}$, T. Hartmann$^{60}$, J. He$^{37}$, T. Head$^{37}$, V. Heijne$^{40}$, K. Hennessy$^{51}$,
P. Henrard$^{5}$, J.A. Hernando Morata$^{36}$, E. van Herwijnen$^{37}$, M. Heß$^{60}$, A. Hicheur$^{1}$, D. Hill$^{54}$,
M. Hoballah$^{5}$, C. Hombach$^{53}$, W. Hulsbergen$^{40}$, P. Hunt$^{54}$, N. Hussain$^{54}$, D. Hutchcroft$^{51}$,
D. Hynds⁵⁰, V. Iakovenko⁴³, M. Idzik²⁶, P. Ilten⁵⁵, R. Jacobsson³⁷, A. Jaeger¹¹, E. Jans⁴⁰,
P. Jaton³⁸, A. Jawahery⁵⁷, F. Jing³, M. John⁵⁴, D. Johnson⁵⁴, C.R. Jones⁴⁶, C. Joram³⁷,
B. Jost³⁷, N. Jurik⁵⁸, M. Kaballo⁹, S. Kandybei⁴², W. Kanso⁶, M. Karacson³⁷, T.M. Karbach³⁷,
I.R. Kenyon⁴⁴, T. Ketel⁴¹, B. Khanji²⁰, C. Khurewathanakul³⁸, S. Klaver⁵³, O. Kochebina⁷,
I. Komarov³⁸, R.F. Koopman⁴¹, P. Koppenburg⁴⁰, M. Korolev³¹, A. Kozlinskiy⁴⁰,
L. Kravchuk³², K. Kreplin¹¹, M. Kreps⁴⁷, G. Krocker¹¹, P. Krokovny³³, F. Kruse⁹,
M. Kucharczyk²⁰,²⁵,³⁷,⁴, V. Kudryavtsev³³, K. Kurek²⁷, T. Kvaratskheliya³⁰,³⁷, V.N. La Thi³⁸,
D. Lacarrere³⁷, G. Lafferty⁵³, A. Lai¹⁵, D. Lambert⁴⁹, R.W. Lambert⁴¹, E. Lanciotti³⁷,
G. Lanfranchi¹⁸, C. Langenbruch³⁷, T. Latham⁴⁷, C. Lazzeroni⁴⁴, R. Le Gac⁶, J. van Leerdam⁴⁰,
J.-P. Lees⁴, R. Lefèvre⁵, A. Leflat³¹, J. Lefrançois⁷, S. Leo²², O. Leroy⁶, T. Lesiak²⁵,
B. Leverington¹¹, Y. Li³, M. Liles⁵¹, R. Lindner³⁷, C. Linn¹¹, F. Lionetto³⁹, B. Liu¹⁵, G. Liu³⁷,
S. Lohn³⁷, I. Longstaff⁵⁰, J.H. Lopes², N. Lopez-March³⁸, P. Lowdon³⁹, H. Lu³, D. Lucchesi²¹,⁴,
J. Luisier³⁸, H. Luo⁴⁹, E. Luppi¹⁶, f, O. Lupton⁵⁴, F. Machefert⁷, I.V. Machikhiliyan³⁰,
F. Maciuc²⁸, O. Maev²⁹,³⁷, S. Malde⁵⁴, G. Manca¹⁵, e, G. Mancinelli⁶, M. Manzali¹⁶, f,
J. Maratas⁵, U. Marconi¹⁴, P. Marino²², t, R. Märki³⁸, J. Marks¹¹, G. Martellotti²⁴, A. Martens⁸,
A. Martín Sánchez⁷, M. Martinelli⁴⁰, D. Martinez Santos⁴¹, D. Martins Tostes², A. Massafferri¹,
R. Matev³⁷, Z. Mathe³⁷, C. Matteuzzi²⁰, A. Mazurov¹⁶,³⁷, f, M. McCann⁵², J. McCarthy⁴⁴,
A. McNab⁵³, R. McNulty¹², B. McSkelly⁵¹, B. Meadows⁵⁶,⁵⁴, F. Meier⁹, M. Meissner¹¹,
M. Merk⁴⁰, D.A. Milanes⁸, M.-N. Minard⁴, J. Molina Rodriguez⁵⁹, S. Monteil⁵, D. Moran⁵³,
M. Morandin²¹, P. Morawski²⁵, A. Mordà⁶, M.J. Morello²², t, R. Mountain⁵⁸, F. Muheim⁴⁹,
K. Müller³⁹, R. Muresan²⁸, B. Muryn²⁶, B. Muster³⁸, P. Naik⁴⁵, T. Nakada³⁸, R. Nandakumar⁴⁸,
I. Nasteva¹, M. Needham⁴⁹, S. Neubert³⁷, N. Neufeld³⁷, A.D. Nguyen³⁸, T.D. Nguyen³⁸,
C. Nguyen-Mau³⁸, q, M. Nicol⁷, V. Niess⁵, R. Niet⁹, N. Nikitin³¹, T. Nikodem¹¹, A. Novoselov³⁴,
A. Oblakowska-Mucha²⁶, V. Obraztsov³⁴, S. Oggero⁴⁰, S. Ogilvy⁵⁰, O. Okhrimenko⁴³,
R. Oldeman¹⁵, e, G. Onderwater⁶², M. Orlandea²⁸, J.M. Otalora Goicochea², P. Owen⁵²,
A. Oyanguren³⁵, B.K. Pal⁵⁸, A. Palano¹³, c, M. Palutan¹⁸, J. Panman³⁷, A. Papanestis⁴⁸,³⁷,
M. Pappagallo⁵⁰, L. Pappalardo¹⁶, C. Parkes⁵³, C.J. Parkinson⁹, G. Passaleva¹⁷, G.D. Patel⁵¹,
M. Patel⁵², C. Patrignani¹⁹, j, C. Pavel-Nicorescu²⁸, A. Pazos Alvarez³⁶, A. Pearce⁵³,
A. Pellegrino⁴⁰, G. Penso²⁴, m, M. Pepe Altarelli³⁷, S. Perazzini¹⁴, d, E. Perez Trigo³⁶, P. Perret⁵,
M. Perrin-Terrin⁶, L. Pescatore⁴⁴, E. Pesen⁶³, G. Pessina²⁰, K. Petridis⁵², A. Petrolini¹⁹, j,
E. Picatoste Olloqui³⁵, B. Pietrzyk⁴, T. Pilař⁴⁷, D. Pinci²⁴, A. Pistone¹⁹, S. Playfer⁴⁹,
M. Plo Casasus³⁶, F. Polci⁸, G. Polok²⁵, A. Poluektov⁴⁷,³³, E. Polycarpo², A. Popov³⁴,
D. Popov¹⁰, B. Popovici²⁸, C. Potterat³⁵, A. Powell⁵⁴, J. Prisciandaro³⁸, A. Pritchard⁵¹,
C. Prouve⁴⁵, V. Pugatch³⁴, A. Puig Navarro³⁸, G. Punzi²², s, W. Qian⁴, B. Rachwal²⁵,
J.H. Rademacker⁴⁵, B. Rakotomiaramanana³⁸, M. Rama¹⁸, M.S. Rangel², I. Raniuk⁴²,
N. Rauschmayr³⁷, G. Raven⁴¹, S. Redford⁵⁴, S. Reichert⁵³, M.M. Reid⁴⁷, A.C. dos Reis¹,
S. Ricciardi⁴⁸, A. Richards⁵², K. Rinnert⁵¹, V. Rives Molina³⁵, D.A. Roa Romero⁵, P. Robbe⁷,
D.A. Roberts⁵⁷, A.B. Rodrigues¹, E. Rodrigues⁵³, P. Rodriguez Perez³⁶, S. Roiser³⁷,
V. Romanovsky³⁴, A. Romero Vidal³⁶, M. Rotondo²¹, J. Rouvinet³⁸, T. Ruf³⁷, F. Ruffini²²,
H. Ruiz³⁵, P. Ruiz Valls³⁵, G. Sabatino²⁴, l, J.J. Saborido Silva³⁶, N. Sagidova²⁹, P. Sail⁵⁰,
B. Saitta¹⁵, e, V. Salustino Guimaraes², B. Sanmartin Sedes³⁶, R. Santacesaria²⁴,
C. Santamarina Rios³⁶, E. Santovetti²³, l, M. Sapunov⁶, A. Sarti¹⁸, C. Satriano²⁴, n, A. Satta²³,
M. Savrie¹⁶, f, D. Savrina³⁰,³¹, M. Schiller⁴¹, H. Schindler³⁷, M. Schlupp⁹, M. Schmelling¹⁰,
B. Schmidt³⁷, O. Schneider³⁸, A. Schopper³⁷, M.-H. Schune⁷, R. Schwemmer³⁷, B. Sciascia¹⁸,
A. Sciubba²⁴, M. Seco³⁶, A. Semennikov³⁰, K. Senderowska²⁶, I. Sepp⁵², N. Serra³⁹, J. Serrano⁶,
P. Seyfert¹¹, M. Shapkin³⁴, I. Shapoval¹⁶,⁴², f, Y. Shcheglov²⁹, T. Shears⁵¹, L. Shekhtman³³,
O. Shevchenko⁴², V. Shevchenko⁶¹, A. Shires<fnref n="9" />, R. Silva Coutinho⁴⁷, G. Simi<fnref n="21" />, M. Sirendi⁴⁶, N. Skidmore⁴⁵, T. Skwarnicki⁵⁸, N.A. Smith⁵¹, E. Smith⁵⁴,⁴⁸, E. Smith⁵², J. Smith⁴⁶, M. Smith⁵³, H. Snoek⁴⁰, M.D. Sokoloff⁵⁶, F.J.P. Soler⁵⁰, F. Soomro³⁸, D. Souza⁴⁵, B. Souza De Paula<fnref n="2" />, B. Spaan<fnref n="9" />, A. Sparkes⁴⁹, F. Spinella<fnref n="22" />, P. Spradlin⁵⁰, F. Stagni³⁷, S. Stahl<fnref n="11" />, O. Steinkamp³⁹, S. Stevenson⁵⁴, S. Stoica²⁸, S. Stone⁵⁸, B. Storaci³⁹, S. Stracka<fnref n="22" />,³⁷, M. Straticiuc²⁸, U. Straumann³⁹, R. Stroili<fnref n="21" />, V.K. Subbiah³⁷, L. Sun⁵⁶, W. Sutcliffe<fnref n="25" />, S. Swientek<fnref n="9" />, V. Syropoulos⁴¹, M. Szczekowski²⁷, P. Szczypka³⁸,³⁷, D. Szilard<fnref n="2" />, T. Szumlak<fnref n="26" />, S. T'Jumpens<fnref n="4" />, M. Teklishyn<fnref n="7" />, G. Tellarini<fnref n="16" />,f, E. Teodorescu²⁸, F. Teubert³⁷, C. Thomas⁵⁴, E. Thomas³⁷, J. van Tilburg<fnref n="11" />, V. Tisserand<fnref n="4" />, M. Tobin³⁸, S. Tolk⁴¹, L. Tomassetti<fnref n="16" />, f, D. Tonelli³⁷, S. Topp-Joergensen⁵⁴, N. Torr⁵⁴, E. Tournefier<fnref n="4" />,⁵², S. Tourneur³⁸, M.T. Tran³⁸, M. Tresch³⁹, A. Tsaregorodtsev<fnref n="6" />, P. Tsopelas⁴⁰, N. Tuning⁴⁰, M. Ubeda Garcia³⁷, A. Ukleja²⁷, A. Ustyuzhanin⁶¹, U. Uwer<fnref n="11" />, V. Vagnoni<fnref n="14" />, G. Valenti<fnref n="14" />, A. Vallier<fnref n="7" />, R. Vazquez Gomez<fnref n="18" />, P. Vazquez Regueiro³⁶, C. Vázquez Sierra³⁶, S. Vecchi<fnref n="16" />, J.J. Velthuis⁴⁵, M. Veltri<fnref n="17" />, h, G. Veneziano³⁸, M. Vesterinen<fnref n="11" />, B. Viaud<fnref n="7" />, D. Vieira<fnref n="2" />, X. Vilasis-Cardona³⁵, p, A. Vollhardt³⁹, D. Volyanskyy<fnref n="10" />, D. Voong⁴⁵, A. Vorobeyev²⁹, V. Vorobeyev³³, C. Voß⁶⁰, H. Voss<fnref n="10" />, J.A. de Vries⁴⁰, R. Waldi⁶⁰, C. Wallace⁴⁷, R. Wallace<fnref n="12" />, S. Wandernoth<fnref n="11" />, J. Wang⁵⁸, D.R. Ward⁴⁶, N.K. Watson⁴⁴, A.D. Webber⁵³, D. Websdale⁵², M. Whitehead⁴⁷, J. Wicht³⁷, J. Wiechczynski<fnref n="25" />, D. Wiedner<fnref n="11" />, L. Wiggers⁴⁰, G. Wilkinson⁵⁴, M.P. Williams⁴⁷,⁴⁸, M. Williams⁵⁵, F.F. Wilson⁴⁸, J. Wimberley⁵⁷, J. Wishahi<fnref n="9" />, W. Wislicki²⁷, M. Witek<fnref n="25" />, G. Wormser<fnref n="7" />, S.A. Wotton⁴⁶, S. Wright⁴⁶, S. Wu<fnref n="3" />, K. Wyllie³⁷, Y. Xie⁴⁹,³⁷, Z. Xing⁵⁸, Z. Yang<fnref n="3" />, X. Yuan<fnref n="3" />, O. Yushchenko³⁴, M. Zangoli<fnref n="14" />, M. Zavertyaev<fnref n="10" />, b, F. Zhang<fnref n="3" />, L. Zhang⁵⁸, W.C. Zhang<fnref n="12" />, Y. Zhang<fnref n="3" />, A. Zhelezov<fnref n="11" />, A. Zhokhov³⁰, L. Zhong<fnref n="3" />, A. Zvyagin³⁷.

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

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

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

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

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

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

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

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

<footnote n="9">⁹Fakultät Physik, Technische Universität Dortmund, Dortmund, Germany</footnote>

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

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

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

<footnote n="13">¹³Sezione INFN di Bari, Bari, Italy</footnote>

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

<footnote n="15">¹⁵Sezione INFN di Cagliari, Cagliari, Italy</footnote>

<footnote n="16">¹⁶Sezione INFN di Ferrara, Ferrara, Italy</footnote>

<footnote n="17">¹⁷Sezione INFN di Firenze, Firenze, Italy</footnote>

<footnote n="18">¹⁸Laboratori Nazionali dell'INFN di Frascati, Frascati, Italy</footnote>

<footnote n="19">¹⁹Sezione INFN di Genova, Genova, Italy</footnote>

<footnote n="20">²⁰Sezione INFN di Milano Bicocca, Milano, Italy</footnote>

<footnote n="21">²¹Sezione INFN di Padova, Padova, Italy</footnote>

<footnote n="22">²²Sezione INFN di Pisa, Pisa, Italy</footnote>

<footnote n="23">²³Sezione INFN di Roma Tor Vergata, Roma, Italy</footnote>

<footnote n="24">²⁴Sezione INFN di Roma La Sapienza, Roma, Italy</footnote>

<footnote n="25">²⁵Henryk Niewodniczanski Institute of Nuclear Physics Polish Academy of Sciences, Kraków, Poland</footnote>

<footnote n="26">²⁶AGH - University of Science and Technology, Faculty of Physics and Applied Computer Science,</footnote>
Kraków, Poland

<sup>27</sup>National Center for Nuclear Research (NCBJ), Warsaw, Poland

<sup>28</sup>Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest-Magurele, Romania

<sup>29</sup>Petersburg Nuclear Physics Institute (PNPI), Gatchina, Russia

<sup>30</sup>Institute of Theoretical and Experimental Physics (ITEP), Moscow, Russia

<sup>31</sup>Institute of Nuclear Physics, Moscow State University (SINP MSU), Moscow, Russia

<sup>32</sup>Institute for Nuclear Research of the Russian Academy of Sciences (INR RAN), Moscow, Russia

<sup>33</sup>Budker Institute of Nuclear Physics (SB RAS) and Novosibirsk State University, Novosibirsk, Russia

<sup>34</sup>Institute for High Energy Physics (IHEP), Protvino, Russia

<sup>35</sup>Universitat de Barcelona, Barcelona, Spain

<sup>36</sup>Universidad de Santiago de Compostela, Santiago de Compostela, Spain

<sup>37</sup>European Organization for Nuclear Research (CERN), Geneva, Switzerland

<sup>38</sup>Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland

<sup>39</sup>Physik-Institut, Universität Zürich, Zürich, Switzerland

<sup>40</sup>Nikhef National Institute for Subatomic Physics, Amsterdam, The Netherlands

<sup>41</sup>Nikhef National Institute for Subatomic Physics and VU University Amsterdam, Amsterdam, The Netherlands

<sup>42</sup>NSC Kharkiv Institute of Physics and Technology (NSC KIPT), Kharkiv, Ukraine

<sup>43</sup>Institute for Nuclear Research of the National Academy of Sciences (KINR), Kyiv, Ukraine

<sup>44</sup>University of Birmingham, Birmingham, United Kingdom

<sup>45</sup>H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom

<sup>46</sup>Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom

<sup>47</sup>Department of Physics, University of Warwick, Coventry, United Kingdom

<sup>48</sup>STFC Rutherford Appleton Laboratory, Didcot, United Kingdom

<sup>49</sup>School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom

<sup>50</sup>School of Physics and Astronomy, University of Glasgow, Glasgow, United Kingdom

<sup>51</sup>Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom

<sup>52</sup>Imperial College London, London, United Kingdom

<sup>53</sup>School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom

<sup>54</sup>Department of Physics, University of Oxford, Oxford, United Kingdom

<sup>55</sup>Massachusetts Institute of Technology, Cambridge, MA, United States

<sup>56</sup>University of Cincinnati, Cincinnati, OH, United States

<sup>57</sup>University of Maryland, College Park, MD, United States

<sup>58</sup>Syracuse University, Syracuse, NY, United States

<sup>59</sup>Pontifícia Universidade Católica do Rio de Janeiro (PUC-Rio), Rio de Janeiro, Brazil, associated to <sup>2</sup>

<sup>60</sup>Institut für Physik, Universität Rostock, Rostock, Germany, associated to <sup>11</sup>

<sup>61</sup>National Research Centre Kurchatov Institute, Moscow, Russia, associated to <sup>30</sup>

<sup>62</sup>KVI - University of Groningen, Groningen, The Netherlands, associated to <sup>40</sup>

<sup>63</sup>Celal Bayar University, Manisa, Turkey, associated to <sup>37</sup>

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

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

<sup>c</sup>Università di Bari, Bari, Italy

<sup>d</sup>Università di Bologna, Bologna, Italy

<sup>e</sup>Università di Cagliari, Cagliari, Italy

<sup>f</sup>Università di Ferrara, Ferrara, Italy

<sup>g</sup>Università di Firenze, Firenze, Italy

<sup>h</sup>Università di Urbino, Urbino, Italy

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

<sup>j</sup>Università di Genova, Genova, Italy

<sup>k</sup>Università di Milano Bicocca, Milano, Italy

<sup>l</sup>Università di Roma Tor Vergata, Roma, Italy
<sup>m</sup>Università di Roma La Sapienza, Roma, Italy

<sup>n</sup>Università della Basilicata, Potenza, Italy

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

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

<sup>q</sup>Hanoi University of Science, Hanoi, Viet Nam

<sup>r</sup>Università di Padova, Padova, Italy

<sup>s</sup>Università di Pisa, Pisa, Italy

<sup>t</sup>Scuola Normale Superiore, Pisa, Italy
# 1 Introduction

Within the framework of heavy quark expansion (HQE) theory [1–7], b-hadron observables are calculated as a perturbative expansion in inverse powers of the b-quark mass, $m_b$. At zeroth order the lifetimes of all weakly decaying b hadrons are equal, with corrections appearing at order $1/m_b^2$. Ratios of b-hadron lifetimes can be theoretically predicted with higher accuracy than absolute lifetimes since many terms in the HQE cancel. The latest theoretical predictions and world-average values for the b-hadron lifetimes and lifetime ratios are reported in Table 1. A measurement of the ratio of the $\Lambda_b^0$ baryon lifetime, using the $\Lambda_b^0 \to J/\psi pK^-$ decay mode<fnref n="1" />, to that of the $B^0$ meson lifetime has recently been made by the LHCb collaboration [8] and is not yet included in the world average.

In this paper, a measurement of the lifetimes of the $B^+$, $B^0$ and $B_s^0$ mesons and $\Lambda_b^0$ baryon is reported using pp collision data, corresponding to an integrated luminosity of 1.0 fb⁻<fnref n="1" />, collected in 2011 with the LHCb detector at a centre-of-mass energy of 7 TeV. The lifetimes are measured from the reconstructed b-hadron decay time distributions of the exclusive decay modes $B^+ \to J/\psi K^+$, $B^0 \to J/\psi K^*(892)^0$, $B^0 \to J/\psi K_S^0$, $B_s^0 \to J/\psi \phi$ and $\Lambda_b^0 \to J/\psi \Lambda$. Collectively, these are referred to as $H_b \to J/\psi X$ decays. In addition, measurements of lifetime ratios are reported.

As a result of neutral meson mixing the decay time distribution of neutral $B_q^0$ mesons ($q \in \{s, d\}$) is characterised by two parameters, namely the average decay width $\Gamma_q$ and the decay width difference $\Delta\Gamma_q$ between the light (L) and heavy (H) $B_q^0$ mass eigenstates. The summed decay rate of $B_q^0$ and $\bar{B}_q^0$ mesons to a final state $f$ is given by [9–11]

$$ \langle \Gamma(B_q^0(t) \to f) \rangle \equiv \Gamma(B_q^0(t) \to f) + \Gamma(\bar{B}_q^0(t) \to f) = R_{q,L}^f e^{-\Gamma_{q,L}t} + R_{q,H}^f e^{-\Gamma_{q,H}t}, \quad (1) $$

where terms proportional to the small flavour specific asymmetry, $a_{fs}^q$, are ignored [12]. Therefore, for non-zero $\Delta\Gamma_q$ the decay time distribution of neutral $B_q^0$ decays is not purely exponential. In the case of an equal admixture of $B_q^0$ and $\bar{B}_q^0$ at $t=0$, the observed average decay time is given by [11]

$$ \tau_{B_q^0 \to f} = \frac{1}{\Gamma_q} \frac{1}{1 - y_q^2} \left( \frac{1 + 2A_{\Delta\Gamma_q}^f y_q + y_q^2}{1 + A_{\Delta\Gamma_q}^f y_q} \right), \quad (2) $$

where $y_q \equiv \Delta\Gamma_q/(2\Gamma_q)$ and $A_{\Delta\Gamma_q}^f \equiv (R_{q,H}^f - R_{q,L}^f)/(R_{q,H}^f + R_{q,L}^f)$ is an observable that depends on the final state, $f$. As such, the lifetimes measured are usually referred to as effective lifetimes. In the $B_s^0$ system, where $\Delta\Gamma_s/\Gamma_s = 0.159 \pm 0.023$ [13], the deviation from an exponential decay time distribution is non-negligible. In contrast, in the $B^0$ system this effect is expected to be small as $\Delta\Gamma_d/\Gamma_d$ is predicted to be $(42 \pm 8) \times 10^{-4}$ in the Standard Model (SM) [14, 15]. Both the BaBar [16, 17] and Belle [18] collaborations have measured $|\Delta\Gamma_d/\Gamma_d|$ and the current world average is $|\Delta\Gamma_d/\Gamma_d| = 0.015 \pm 0.018$ [13]. A deviation in the value of $\Delta\Gamma_d$ from the SM prediction has recently been proposed [19] as

<footnote n="1">¹Charge conjugation is implied throughout this paper, unless otherwise stated.</footnote>
<table><caption>Table 1: Theoretical predictions and current world-average values [13] for b-hadron lifetimes and lifetime ratios.</caption><thead><tr><th>Observable</th><th>Prediction</th><th>World average</th></tr></thead><tbody><tr><td>τ<sub>B<sup>+</sup></sub> [ps]</td><td>–</td><td>1.641 ± 0.008</td></tr><tr><td>τ<sub>B<sup>0</sup></sub> [ps]</td><td>–</td><td>1.519 ± 0.007</td></tr><tr><td>τ<sub>B<sup>0</sup></sub> [ps]</td><td>–</td><td>1.516 ± 0.011</td></tr><tr><td>τ<sub>Λ<sub>b</sub><sup>0</sup></sub> [ps]</td><td>–</td><td>1.429 ± 0.024</td></tr><tr><td>τ<sub>B<sup>+</sup>/τ<sub>B<sup>0</sup></sub></sub></td><td>1.063 ± 0.027 [15, 22, 23]</td><td>1.079 ± 0.007</td></tr><tr><td>τ<sub>B<sup>0</sup>/τ<sub>B<sup>0</sup></sub></sub></td><td>1.00 ± 0.01 [15, 23–25]</td><td>0.998 ± 0.009</td></tr><tr><td>τ<sub>Λ<sub>b</sub><sup>0</sup>/τ<sub>B<sup>0</sup></sub></sub></td><td>0.86–0.95 [3, 23, 26–32]</td><td>0.941 ± 0.016</td></tr></tbody></table>

a potential explanation for the anomalous like-sign dimuon charge asymmetry measured by the D0 collaboration [20]. In this paper, $\Delta\Gamma_d/\Gamma_d$ is measured from the effective lifetimes of $B^0 \to J/\psi K^*(892)^0$ and $B^0 \to J/\psi K_s^0$ decays, as proposed in Ref. [21].

The main challenge in the measurements reported is understanding and controlling the detector acceptance, reconstruction and selection efficiencies that depend upon the b-hadron decay time. This paper is organised as follows. Section 2 describes the LHCb detector and software. The selection criteria for the b-hadron candidates are described in Sec. 3. Section 4 describes the reconstruction efficiencies and the techniques used to correct the decay time distributions. Section 5 describes how the efficiency corrections are incorporated into the maximum likelihood fit that is used to measure the signal yields and lifetimes. The systematic uncertainties on the measurements are described in Sec. 6. The final results and conclusions are presented in Sec. 7.

## 2 Detector and software

The LHCb detector [33] is a single-arm forward spectrometer covering the pseudorapidity range $2 < \eta < 5$, designed for the study of particles containing b or c quarks. The detector includes a high-precision tracking system consisting of a silicon-strip vertex detector (VELO) surrounding the pp interaction region, a large-area silicon-strip detector (TT) 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. The combined tracking system provides a momentum, $p$, measurement with relative uncertainty that varies from 0.4% at 5 GeV/c to 0.6% at 100 GeV/c, and impact parameter resolution of 20 µm for charged particles with high transverse momentum, $p_T$. Charged hadrons are identified using two ring-imaging Cherenkov detectors [34]. Photon, electron and hadron candidates 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 [35]. The right-handed coordinate system adopted has the z-axis along the beam line and the y-axis along the vertical. The trigger [36] 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.

Two distinct classes of tracks are reconstructed using hits in the tracking stations on both sides of the magnet, either with hits in the VELO (long track) or without (downstream track). The vertex resolution of b-hadron candidates reconstructed using long tracks is better than that for candidates reconstructed using downstream tracks. However, the use of long tracks introduces a dependence of the reconstruction efficiency on the b-hadron decay time.

In the simulation, *pp* collisions are generated using PYTHIA 6.4 [37] with a specific LHCb configuration [38]. Decays of hadronic particles are described by EvtGen [39], in which final state radiation is generated using PHOTOS [40]. The interaction of the generated particles with the detector and its response are implemented using the GEANT4 toolkit [41] as described in Ref. [42].

## 3 Candidate selection

The reconstruction of each of the $H_b \to J/\psi X$ decays is similar and commences by selecting $J/\psi \to \mu^+\mu^-$ decays. Events passing the hardware trigger contain dimuon candidates with high transverse momentum. The subsequent software trigger is composed of two stages. The first stage 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^2$. The selection at this stage has a uniform efficiency as a function of decay time. The second stage performs a full event reconstruction, calculating the position of each *pp* interaction vertex (PV) using all available charged particles in the event. The average number of PVs in each event is approximately 2.0. Their longitudinal ($z$) position is known to a precision of approximately 0.05 mm. If multiple PVs are reconstructed in the event, the one with the minimum value of $\chi_{\text{IP}}^2$ is associated with the $J/\psi$ candidate, where $\chi_{\text{IP}}^2$ is the increase in the $\chi^2$ of the PV fit if the candidate trajectory is included. Events are retained for further processing if they contain a $J/\psi \to \mu^+\mu^-$ pair that forms a vertex that is significantly displaced from the PV. This introduces a non-uniform efficiency as function of decay time.

The offline sample of $J/\psi$ meson candidates is selected by requiring each muon to have $p_T$ larger than 500 MeV/$c$ and the $J/\psi$ candidate to be displaced from the PV by more than three times its decay length uncertainty. The invariant mass of the two muons, $m(\mu^+\mu^-)$, must be in the range [3030, 3150] MeV/$c^2$.

The b-hadron candidate selection is performed by applying kinematic and particle identification criteria to the final-state tracks, the details of which are reported in Sec. 3.1 to 3.5. No requirements are placed on variables that are highly correlated to the b-hadron decay time, thereby avoiding the introduction of additional biases. All final-state particles are required to have a pseudorapidity in the range $2.0 < \eta < 4.5$. In addition, the $z$-position of the PV ($z_{\text{PV}}$) is required to be within 100 mm of the nominal interaction point, where the standard deviation of the $z_{\text{PV}}$ distribution is approximately 47 mm. These criteria cause a reduction of approximately 10% in signal yield but define a fiducial region where
the reconstruction efficiency is largely uniform.

The maximum likelihood fit uses the invariant mass, $m(J/\psi X)$, and proper decay time, $t$, of each b-hadron candidate. The decay time of the b-hadron candidate in its rest frame is derived from the relation $t = ml/q$, where $m$ is its invariant mass and the decay length, $l$, and the momentum, $q$, are measured in the experimental frame. In this paper, $t$ is computed using a kinematic decay-tree fit (DTF) [43] involving all final-state tracks from the b-hadron candidate with a constraint on the position of the associated PV. Unlike in the trigger, the position of each PV is calculated using all available charged particles in the event after the removal of the b-hadron candidate final-state tracks. This is necessary to prevent the final-state tracks from biasing the PV position towards the b-hadron decay vertex and helps to reduce the tails of the decay-time resolution function. This prescription does not bias the measured lifetime using simulated events. The $\chi^2$ of the fit, $\chi^2_{DTF}$, is useful to discriminate between signal and background. In cases where there are multiple b-hadron candidates per event, the candidate with the smallest $\chi^2_{DTF}$ is chosen. The $z$-position of the displaced b-hadron vertices are known to a precision of approximately 0.15 mm.

Studies of simulated events show that in the case of $B^0 \to J/\psi K^{*0}$ ($B_s^0 \to J/\psi \phi$) decays, imposing requirements on $\chi^2_{DTF}$ introduces a dependence of the selection efficiency on the decay time if the $K^+$ and $\pi^-$ ($K^+$ and $K^-$) tracks are included in the DTF. If no correction is applied to the decay time distribution, the measured lifetime would be biased by approximately -2 fs relative to the generated value. Using simulated events it is found that this effect is correlated to the opening angle between the $K^+$ and $\pi^-$ ($K^+$ and $K^-$) from the $K^{*0}$ ($\phi$) decay. No effect is observed for the muons coming from the $J/\psi$ decay due to the larger opening angle in this case. To remove the effect, the calculation of $\chi^2_{DTF}$ for the $B^0 \to J/\psi K^{*0}$ and $B_s^0 \to J/\psi \phi$ channels is performed with an alternative DTF in which the assigned track parameter uncertainties of the kaon and pion are increased in such a way that their contribution to the b-hadron vertex position is negligible.

Candidates are required to have $t$ in the range [0.3, 14.0] ps. The lower bound on the decay time suppresses a large fraction of the prompt combinatorial background that is composed of tracks from the same PV, while the upper bound is introduced to reduce the sensitivity to long-lived background candidates. In the case of the $B^0 \to J/\psi K_s^0$ and $\Lambda_b^0 \to J/\psi \Lambda$ decays, the lower bound is increased to 0.45 ps to compensate for the worse decay time resolution in these modes.

In events with multiple PVs, b-hadron candidates are removed if they have a $\chi^2_{\text{IP}}$ with respect to the next best PV smaller than 50. This requirement is found to distort the decay time distribution, but reduces a source of background due to the incorrect association of the b hadron to its production PV.

The invariant mass is computed using another kinematic fit without any constraint on the PV position but with the invariant mass of the $\mu^+\mu^-$ pair, $m(\mu^+\mu^-)$, constrained to the known $J/\psi$ mass [44]. Figures 1 and 2 show the $m(J/\psi X)$ distributions for the selected candidates in each final state and Table 2 gives the corresponding signal yields.
<figure><img src="image_2.png" /><figcaption>Figure 1: Distributions of the (left) mass and (right) decay time of $B^+ \to J/\psi K^+$, $B^0 \to J/\psi K^{*0}$ and $B^0 \to J/\psi K_S^0$ candidates and their associated residual uncertainties (pulls). The data are shown by the black points; the total fit function by the black solid line; the signal contribution by the red dashed line and the background contribution by the blue dotted line.</figcaption></figure>
<figure><img src="image_3.png" /><figcaption>Figure 2: Distributions of the (left) mass and (right) decay time of $B_s^0 \rightarrow J/\psi \phi$ and $\Lambda_b^0 \rightarrow J/\psi \Lambda$ candidates and their associated residual uncertainties (pulls). The data are shown by the black points; the total fit function by the black solid line; the signal contribution by the red dashed line and the background contribution by the blue dotted line.</figcaption></figure>
<table><caption>Table 2: Estimated event yields for the five $b \to J/\psi X$ channels selected using the criteria described in Sec. 3.1 to 3.5.</caption><thead><tr><th>Channel</th><th>Yield</th></tr></thead><tbody><tr><td>$B^{+} \to J/\psi K^{+}$</td><td>229 434 ± 503</td></tr><tr><td>$B^{0} \to J/\psi K^{*0}$</td><td>70 534 ± 312</td></tr><tr><td>$B^{0} \to J/\psi K_{s}^{0}$</td><td>17 045 ± 175</td></tr><tr><td>$B_{s}^{0} \to J/\psi \phi$</td><td>18 662 ± 152</td></tr><tr><td>$\Lambda_{b}^{0} \to J/\psi \Lambda$</td><td>3 960 ± 89</td></tr></tbody></table>

## 3.1 Selection of $B^{+} \to J/\psi K^{+}$ decays

The $B^{+}$ candidates are reconstructed by combining the $J/\psi$ candidates with a charged particle that is identified as a kaon with $p_{T}$ larger than 1 GeV/c and $p$ larger than 10 GeV/c. The invariant mass, $m(J/\psi K^{+})$, must be in the range [5170, 5400] MeV/c$^{2}$, where the lower bound is chosen to remove feed-down from incompletely reconstructed $B^{0} \to J/\psi K^{*0}$ decays. The $\chi_{\text{DTF}}^{2}$ of the fit, which has 5 degrees of freedom, is required to be less than 25. Multiple $B^{+}$ candidates are found in less than 0.02% of selected events.

## 3.2 Selection of $B^{0} \to J/\psi K^{*0}$ decays

The $K^{*0}$ candidates are reconstructed by combining two oppositely charged particles that are identified as a kaon and a pion. The pion and $K^{*0}$ must have $p_{T}$ greater than 0.3 GeV/c and 1.5 GeV/c, respectively. The invariant mass, $m(K^{+}\pi^{-})$, must be in the range [826, 966] MeV/c$^{2}$.

The $B^{0}$ candidates are reconstructed by combining the $J/\psi$ and $K^{*0}$ candidates. The invariant mass, $m(J/\psi K^{+}\pi^{-})$, must be in the range [5150, 5340] MeV/c$^{2}$, where the upper bound is chosen to remove the contribution from $B_{s}^{0} \to J/\psi K^{*0}$ decays. The $\chi_{\text{DTF}}^{2}$ of the fit, which has 3 degrees of freedom, is required to be less than 15. Multiple $B^{0}$ candidates are found in 2.2% of selected events.

## 3.3 Selection of $B^{0} \to J/\psi K_{s}^{0}$ decays

The $K_{s}^{0}$ candidates are formed from the combination of two oppositely charged particles that are identified as pions and reconstructed as downstream tracks. This is necessary since studies of simulated signal decays demonstrate that an inefficiency depending on the b-hadron decay time is introduced by the reconstruction of the long-lived $K_{s}^{0}$ and $\Lambda$ particles using long tracks. Even so, it is found that the acceptance of the TT still depends on the origin of the tracks. This effect is removed by further tightening of the requirement on the position of the PV to $z_{\text{PV}} > -50$ mm.

For particles produced close to the interaction region, this effect is suppressed by the requirements on the fiducial region for the PV, which is further tightened by requiring
that, to account for the additional acceptance introduced by the TT.

The downstream pions are required to have $p_T$ greater than 0.1 GeV/c and $p$ greater than 2 GeV/c. The $K_S^0$ candidate must have $p_T$ greater than 1 GeV/c and be well separated from the $B^0$ decay vertex, to suppress potential background from $B^0 \to J/\psi K^{*0}$ decays where the kaon has been misidentified as a pion. The $\chi^2$ of the $K_S^0$ vertex fit must be less than 25 and the invariant mass of the dipion system, $m(\pi^+\pi^-)$, must be within 15 MeV/c$^2$ of the known $K_S^0$ mass [44]. For subsequent stages of the selection, $m(\pi^+\pi^-)$ is constrained to the known $K_S^0$ mass.

The invariant mass, $m(J/\psi K_S^0)$, of the $J/\psi$ and $K_S^0$ candidate combination must be in the range [5150, 5340] MeV/c$^2$, where the upper bound is chosen to remove the contribution from $B_S^0 \to J/\psi K_S^0$ decays. The $\chi^2_{DTF}$ of the fit, which has 6 degrees of freedom, is required to be less than 30. Multiple $B_S^0$ candidates are found in less than 0.4% of selected events.

### 3.4 Selection of $B_S^0 \to J/\psi \phi$ decays

The $\phi$ candidates are formed from two oppositely charged particles that have been identified as kaons and originate from a common vertex. The $K^+K^-$ pair is required to have $p_T$ larger than 1 GeV/c. The invariant mass of the $K^+K^-$ pair, $m(K^+K^-)$, must be in the range [990, 1050] MeV/c$^2$.

The $B_S^0$ candidates are reconstructed by combining the $J/\psi$ candidate with the $K^+K^-$ pair, requiring the invariant mass, $m(J/\psi K^+K^-)$, to be in the range [5200, 5550] MeV/c$^2$. The $\chi^2_{DTF}$ of the fit, which has 3 degrees of freedom, is required to be less than 15. Multiple $B_S^0$ candidates are found in less than 2.0% of selected events.

### 3.5 Selection of $\Lambda_b^0 \to J/\psi \Lambda$ decays

The selection of $\Lambda_b^0 \to J/\psi \Lambda$ candidates follows a similar approach to that adopted for $B^0 \to J/\psi K_S^0$ decays. Only downstream protons and pions are used to reconstruct the $\Lambda$ candidates. The pions are required to have $p_T$ larger than 0.1 GeV/c, while pions and protons must have $p$ larger than 2 GeV/c. The $\Lambda$ candidate must be well separated from the $\Lambda_b^0$ decay vertex and have $p_T$ larger than 1 GeV/c. The $\chi^2$ of the $\Lambda$ vertex fit must be less than 25 and $m(p\pi^-)$ must be within 6 MeV/c$^2$ of the known $\Lambda$ mass [44]. For subsequent stages of the selection, $m(p\pi^-)$ is constrained to the known $\Lambda$ mass.

The invariant mass, $m(J/\psi \Lambda)$, of the $J/\psi$ and $\Lambda$ candidate combination must be in the range [5470, 5770] MeV/c$^2$. The $\chi^2_{DTF}$ of the fit, which has 6 degrees of freedom, is required to be less than 30. Multiple $\Lambda_b^0$ candidates are found in less than 0.5% of selected events.

## 4 Dependence of efficiencies on decay time

Section 3 described the reconstruction and selection criteria of the $H_b \to J/\psi X$ decays and various techniques that have been used to minimise the dependence of selection efficiencies upon the decay time. After these steps, there remain two effects that distort the $b$-hadron
<figure><img src="image_4.png" /><figcaption>Figure 3: VELO-track reconstruction efficiency for kaon tracks reconstructed using the (a) online and (b) offline algorithms as a function of the kaon $\rho$, as defined in Eq. (3). The red solid lines show the result of an unbinned maximum likelihood fit using the parameterisation in Eq. (4) to the background subtracted data (black points).</figcaption></figure>

decay time distribution. These are caused by the VELO-track reconstruction efficiency, $\varepsilon_{VELO}$, and the combination of the trigger efficiency, $\varepsilon_{trigger}$, and offline selection efficiency, $\varepsilon_{selection|trigger}$. This section will describe these effects and the techniques that are used to evaluate the efficiencies from data control samples.

## 4.1 VELO-track reconstruction efficiency

The largest variation of the efficiency with the decay time is introduced by the track reconstruction in the VELO. The track finding procedure in the VELO assumes that tracks originate approximately from the interaction region [33,45]. In the case of long-lived b-hadron candidates this assumption is not well justified, leading to a loss of reconstruction efficiency for charged particle tracks from the b-hadron decay.

The distance of closest approach of the track to the z-axis is defined as

$$ \rho \equiv \frac{|(\mathbf{d} - \mathbf{v}) \cdot (\mathbf{p} \times \hat{\mathbf{z}})|}{|\mathbf{p} \times \hat{\mathbf{z}}|}, \qquad (3) $$

where $\mathbf{p}$ is the momentum of the final-state track from a b-hadron candidate decaying at point $\mathbf{d}$, $\hat{\mathbf{z}}$ is a unit vector along the z-axis and $\mathbf{v}$ is the origin of the VELO coordinate system. During data taking the position of the LHCb VELO is monitored as a function of time and is centred around the LHC beam line. Using a control sample of $B^+ \to J/\psi K^+$ candidates where the $K^+$ is reconstructed as a downstream track, the VELO-track reconstruction efficiency, $\varepsilon_{VELO}(\rho)$, is computed as the fraction of these tracks that are also reconstructed as long tracks. From samples of simulated b-hadron decays, it is observed that $\varepsilon_{VELO}(\rho)$ can be empirically parameterised by

$$ \varepsilon_{VELO}(\rho) = a(1 + c\rho^2), \qquad (4) $$
<table><caption>Table 3: VELO reconstruction efficiency in data for kaon tracks reconstructed with the online and offline algorithms. In both cases, the correlation coefficient between a and c is 0.2.</caption><thead><tr><th></th><th>a</th><th>c [mm<sup>-2</sup>]</th></tr></thead><tbody><tr><td>Online</td><td>0.9759 &pm; 0.0005</td><td>-0.0093 &pm; 0.0007</td></tr><tr><td>Offline</td><td>0.9831 &pm; 0.0004</td><td>-0.0041 &pm; 0.0005</td></tr></tbody></table>

where the parameters a and c are determined from a fit to the unbinned efficiency distribution.

Figure 3 shows the VELO-track reconstruction efficiency obtained using this method and Table 3 shows the corresponding fit results. Since different configurations of the VELO reconstruction algorithms are used within the LHCb software trigger (online) and during the subsequent processing (offline), it is necessary to evaluate two different efficiencies. The stronger dependence of the online efficiency as a function of $\rho$ is due to the additional requirements used in the first stage of the software trigger such that it satisfies the required processing time.

Applying the same technique to a simulated sample of $B^+ \to J/\psi K^+$ decays yields qualitatively similar behaviour for $\varepsilon_{VELO}(\rho)$. Studies on simulated data show that the efficiency for kaons and pions from the decay of $\phi$ and $K^{*0}$ mesons is smaller than for the kaon in $B^+ \to J/\psi K^+$ decays, due to the small opening between the particles in the $\phi$ and $K^{*0}$ decays, as discussed in Sec. 3. In addition, there are kinematic differences between the calibration $B^+$ sample and the signal samples. Scaling factors on the efficiency parameters are derived from simulation to account for these effects, and have typical sizes in the range [1.04, 1.65], depending on the decay mode and final-state particle being considered.

The distortion to the b-hadron candidate decay time distribution caused by the VELO-track reconstruction is corrected for by weighting each b-hadron candidate by the inverse of the product of the per-track efficiencies. The systematic effect introduced by this weighting is tested using simulated samples of each channel. The chosen efficiency depends on whether the particle is reconstructed with the online or offline variant of the algorithm. Studies on simulated data show that tracks found by the online tracking algorithm are also found by the offline tracking efficiency. For example, the efficiency weight for each $B^0 \to J/\psi K^{*0}$ candidate takes the form

$$w_{B^0 \to J/\psi K^{*0}} = 1 / \left( \varepsilon_{VELO,online}^{\mu^+} \varepsilon_{VELO,online}^{\mu^-} \varepsilon_{VELO,offline}^{K^+} \varepsilon_{VELO,offline}^{\pi^-} \right), \quad (5)$$

since the two muons are required to be reconstructed online, while the kaons and the pions are reconstructed offline.

In the case of the $B^0 \to J/\psi K_S^0$ and $\Lambda_b^0 \to J/\psi \Lambda$ channels, since no VELO information is used when reconstructing the $K_S^0$ and $\Lambda$ particles, the candidate weighting functions take the form $w = 1 / (\varepsilon_{VELO,online}^{\mu^+} \varepsilon_{VELO,online}^{\mu^-})$.
<figure><img src="image_5.png" /><figcaption>Figure 4: Combined trigger and selection efficiency, $\varepsilon_{\text{selection}}(t)$, for $B^{+} \to J/\psi K^{+}$ candidates.</figcaption></figure>

## 4.2 Trigger and selection efficiency

The efficiency of the second stage of the software trigger depends on the b-hadron decay time as it requires that the $J/\psi$ meson is significantly displaced from the PV. A parameterisation of this efficiency, $\varepsilon_{\text{trigger}}(t)$, is obtained for each $b \to J/\psi X$ decay mode by exploiting a corresponding sample of $b \to J/\psi X$ candidates that are selected without any displacement requirement. For each channel, the control sample corresponds to approximately 40% of the total number of signal candidates. A maximum likelihood fit to the unbinned invariant mass distribution $m(J/\psi X)$ is performed to determine the fraction of signal decays that survive the decay-time biasing trigger requirements as a function of decay time.

The same technique is used to determine the decay time efficiency of the triggered candidates caused by the offline selection, $\varepsilon_{\text{selection|trigger}}(t)$, which is introduced by the requirement on the detachment of the $J/\psi$ mesons in the sample used to reconstruct the b-hadron decays. The combined selection efficiency, $\varepsilon_{\text{selection}}(t)$, is given by the product of $\varepsilon_{\text{trigger}}(t)$ and $\varepsilon_{\text{selection|trigger}}(t)$.

Figure 4 shows $\varepsilon_{\text{selection}}(t)$ obtained for the $B^{+} \to J/\psi K^{+}$ channel as a function of decay time. The efficiencies obtained for the other $H_b \to J/\psi X$ channels are qualitatively similar. Studies using simulated events show that the efficiency drop below 0.5 ps is caused by the $J/\psi$ displacement requirement. The dip near 1.5 ps appears because the PV reconstruction in the software trigger is such that some final-state tracks of short-lived b-hadron decays may be used to reconstruct an additional fake PV close to the true b-hadron decay vertex. As a result the reconstructed $J/\psi$ meson does not satisfy the displacement requirement, leading to a decrease in efficiency.

The efficiency parameterisation for each channel is used in the fit to measure the corresponding b-hadron lifetime. An exception is made for the $\Lambda_b^0 \to J/\psi \Lambda$ channel where, owing to its smaller event yield, $\varepsilon_{\text{selection}}(t)$ measured with $B^0 \to J/\psi K_S^0$ decays is used instead. The validity of this approach is checked using large samples of simulated events.
# 5 Maximum likelihood fit

For each channel, the lifetime is determined from a two-dimensional maximum likelihood fit to the unbinned $m(J/\psi X)$ and $t$ distributions. The full probability density function (PDF) is constructed as $\mathcal{P} = f_s(S_m \times S_t) + (1 - f_s)(B_m \times B_t)$, where $f_s$ is the signal fraction, determined in the fit, and $S_m \times S_t$ and $B_m \times B_t$ are the $(m(J/\psi X), t)$ PDFs for the signal and background components, respectively. A systematic uncertainty is assigned to the assumption that the PDFs factorise.

The signal mass PDF, $S_m$, is modelled by the sum of two Gaussian functions. The free parameters in the fit are the common mean, the width of the narrower Gaussian function, the ratio of the second to the first Gaussian width and the fraction of the first Gaussian function. The background mass distribution, $B_m$, is modelled by an exponential function with a single free parameter.

The signal b-hadron decay time distribution is described by an exponential function with decay constant given by the b-hadron lifetime, $\tau_{H_b \to J/\psi X}$. The signal decay time PDF, $S_t$, is obtained by multiplying the exponential function by the combined $t$-dependent trigger and selection efficiency described in Sec. 4.2. From inspection of events in the sidebands of the b-hadron signal peak, the background decay time PDF, $B_t$, is well modelled by a sum of three exponential functions with different decay constants that are free in the fit. These components originate from a combination of prompt candidates, where all tracks originate from the same PV, and long-lived candidates where tracks from the associated PV are combined with other tracks of long-lived particles. For each channel the exponential functions are convolved with a Gaussian resolution function with width $\sigma$ and mean $\Delta$, an offset of the order of a few femtoseconds that is fixed in the fit. Using a sample of prompt $J/\psi$ background events, the decay time resolution for $H_b \to J/\psi X$ channels reconstructed using long tracks has been measured to be approximately 45 fs [46]. For $B^0 \to J/\psi K_S^0$ and $\Lambda_b^0 \to J/\psi \Lambda$ decays, which use downstream tracks to reconstruct the $K_S^0$ and $\Lambda$ particles, a similar study of an event sample composed of prompt $J/\psi$ mesons combined with two downstream tracks, reconstructed as either a $K_S^0$ or $\Lambda$, has determined the resolution to be 65 fs. The systematic uncertainties related to the choice of resolution model are discussed in Sec. 6.

The negative log-likelihood, constructed as

$$-\ln \mathcal{L} = -\alpha \sum_{\text{events } i} w_i \ln \mathcal{P}, \qquad (6)$$

is minimised in the fit, where the weights $w_i$ correspond to the per-candidate correction for the VELO reconstruction efficiency described in Sec. 4.1. The factor $\alpha = \sum_i w_i / \sum_i w_i^2$ is used to include the effect of the weights in the determination of the uncertainties [47]. Figures 1 and 2 show the result of fitting this model to the selected candidates for each channel, projected onto the corresponding $m(J/\psi X)$ and $t$ distributions.

As a consistency check, an alternative fit procedure is developed where each event is given a signal weight, $W_i$, determined using the sPlot [48] method with $m(J/\psi X)$ as the discriminating variable and using the mass model described above. A weighted fit to the
decay time distribution using the signal PDF is then used to measure the b-hadron lifetime. In this case, the negative log-likelihood is given by Eq. (6) where $w_i$ is replaced with $W_i w_i$ and $\alpha = \frac{\sum_i (W_i w_i)}{\sum_i (W_i w_i)^2}$. The difference between the results of the two fitting procedures is used to estimate the systematic uncertainty on the background description.

# 6 Systematic uncertainties

The systematic effects affecting the measurements reported here are discussed in the following and summarised in Tables 4 and 5.

The systematic uncertainty related to the VELO-track reconstruction efficiency can be split into two components. The first uncertainty is due to the finite size of the $B^+ \to J/\psi K^+$ sample, reconstructed using downstream kaon tracks, which is used to determine the per-candidate efficiency weights and leads to a statistical uncertainty on the $\varepsilon_{VELO}(\rho)$ parameterisation. The lifetime fits are repeated after varying the parameters by $\pm 1\sigma$ and the largest difference between the lifetimes is assigned as the uncertainty. The second uncertainty is due to the scaling factors, which are used to correct the efficiency for phase-space effects, obtained from simulated events. The fit is repeated using the unscaled efficiency and half of the variation in fit results is assigned as a systematic uncertainty. These contributions, of roughly the same size, are added in quadrature in Table 4.

A number of additional consistency checks are performed to investigate possible mismodelling of the VELO-track reconstruction efficiency. First, $\varepsilon_{VELO}(\rho)$ is evaluated in two track momentum and two track multiplicity bins and the event weights recalculated. Using both data and simulated events, no significant change in the lifetimes is observed after repeating the fit with the updated weights and, therefore, no systematic uncertainty is assigned. Secondly, to assess the sensitivity to the choice of parameterisation for $\varepsilon_{VELO}(\rho)$ (Eq. 4), the results are compared to those with linear model for the efficiency. The effect is found to be negligible and no systematic uncertainty is applied. Thirdly, the dependence of the VELO-track reconstruction efficiency on the azimuthal angle, $\phi$, of each track is studied by independently evaluating the efficiency in four $\phi$ quadrants for both data and simulation. No dependence is observed. Finally, the efficiency is determined separately for both positive and negative kaons and found to be compatible.

The techniques described in Sec. 4 to correct the efficiency as a function of the decay time are validated on simulated data. The lifetime is fit in each simulated signal mode and the departure from the generated lifetime, $\Delta\tau$, is found to be compatible with zero within the statistical precision of each simulated sample. The measured lifetimes in the data sample are corrected by each $\Delta\tau$ and a corresponding systematic uncertainty is assigned, given by the size of the statistical uncertainty on the fitted lifetime for each simulated signal mode.

The assumption that $m(J/\psi X)$ is independent of the decay time is central to the validity of the likelihood fits used in this study. It is tested by re-evaluating the signal weights of the alternative fit in bins of decay time and then refitting the entire sample using the modified weights. The systematic uncertainty for each decay mode is evaluated
as the sum in quadrature of the lifetime variations, each weighted by the fraction of signal events in the corresponding bin.

For each signal decay mode, the effect of the limited size of the control sample used to estimate the combined trigger and selection efficiency is evaluated by repeating the fits with $\epsilon_{selection}(t)$ randomly fluctuated within its statistical uncertainty. The standard deviation of the distribution of lifetimes obtained is assigned as the systematic uncertainty.

The alternative likelihood fit does not assume any model for the decay time distribution associated with the combinatorial background. Therefore, the systematic uncertainty associated to the modelling of this background is evaluated by taking the difference in lifetimes measured by the nominal and alternative fit methods.

The fit uses a double Gaussian function to describe the $m(J/\psi X)$ distribution of signal candidates. This assumption is tested by repeating the fit using a double-sided Apollonios function [49] where the mean and width parameters are varied in the fit and the remaining parameters are fixed to those from simulation. The differences in lifetime with respect to the default results are taken as systematic uncertainties.

As described in Sec. 5 the dominant background in each channel is combinatorial in nature. It is also possible for backgrounds to arise due to misreconstruction of b-hadron decays where the particle identification has failed. The presence of such backgrounds is checked by inspecting events in the sidebands of the signal and re-assigning the mass hypotheses of at least one of the final-state hadrons. The only contributions that have an impact are $\Lambda_b^0 \to J/\psi pK^-$ decays in the $B_s^0 \to J/\psi \phi$ channel where a proton is misidentified as a kaon and the cross-feed component between $B^0 \to J/\psi K_s^0$ and $\Lambda_b^0 \to J/\psi \Lambda$ decays where pion and protons are misidentified. In the case of $B_s^0 \to J/\psi \phi$ decays, the fit is repeated including a contribution of $\Lambda_b^0 \to J/\psi pK^-$ decays. The two-dimensional PDF is determined from simulation, while the yield is found to be 6% from the sidebands of the $B_s^0$ sample. This leads to the effective lifetime changing by 0.4 fs, which is assigned as a systematic uncertainty. A similar procedure is used to determine the invariant mass shape of the cross-feed background between $B^0 \to J/\psi K_s^0$ and $\Lambda_b^0 \to J/\psi \Lambda$ decays, while the decay time is modelled with the exponential distribution of the corresponding signal mode. A simultaneous fit to both data samples is performed in order to constrain the yield of the cross-feed and the resulting change in lifetime of -0.3 fs and +1.1 fs for $B^0$ and $\Lambda_b^0$, respectively, is assigned as a systematic uncertainty.

Another potential source of background is the incorrect association of signal b hadrons to their PV, which results in an erroneous reconstruction of the decay time. Since the fitting procedure ignores this contribution, a systematic uncertainty is evaluated by repeating the fit after including in the background model a component describing the incorrectly associated candidates. The background distribution is determined by calculating the decay time for each $B^+ \to J/\psi K^+$ decay with respect to a randomly chosen PV from the previous selected event. In studies of simulated events the fraction of this background is less than 0.1%. Repeating the fit with a 1% contribution results in the lifetime changing by 0.1 fs and, therefore, no systematic uncertainty is assigned.

The measurement of the effective lifetime in the $B_s^0 \to J/\psi \phi$ channel is integrated over the angular distributions of the final-state particles and is, in the case of uniform angular
efficiency, insensitive to the different polarisations of the final state [46]. To check if the
angular acceptance introduced by the detector geometry and event selection can affect
the measured lifetime, the events are weighted by the inverse of the angular efficiency
determined in Ref. [46]. Repeating the fit with the weighted dataset leads to a shift of
the lifetime of -1.0 fs, the same as is observed in simulation. The final result is corrected
by this shift, which is also assigned as a systematic uncertainty. The $B_s^0$ effective lifetime
could also be biased due to a small CP-odd S-wave component from $B_s^0 \rightarrow J/\psi K^+ K^-$
decays that is ignored in the fit. For the $m(K^+K^-)$ mass range used here (Sec. 3), Ref. [50]
indicates that the S-wave contribution is 1.1%. The effect of ignoring such a component is
evaluated by repeating the fit on simulated experiments with an additional 1% CP-odd
component. A change in the lifetime of -1.2 fs is observed, which is used to correct the
final lifetime and is also taken as a systematic uncertainty. Finally, as described in Sec. 3,
only events with a decay time larger than 0.3 ps are considered in the nominal fit. This
offset leads to a different relative contribution of the heavy and light mass eigenstates such
that the lifetime extracted from the exponential fit does not correspond to the effective
lifetime defined in Eq. (2). A correction of -0.3 fs is applied to account for this effect.

The presence of a production asymmetry between $B^0$ and $\bar{B}^0$ mesons could bias the measured $B^0 \rightarrow J/\psi K_S^0$ effective lifetime, and therefore $\Delta\Gamma_d/\Gamma_d$, by adding additional terms in Eq. (2). The production asymmetry is measured to be $A_P(B^0) = (0.6 \pm 0.9)\%$ [51], the uncertainty of which is used to estimate a corresponding systematic uncertainty on the $B^0 \rightarrow J/\psi K_S^0$ lifetime of 1.1 fs. No uncertainty is assigned to the $B^0 \rightarrow J/\psi K^{*0}$ lifetime since this decay mode is flavour-specific<fnref n="2" /> and the production asymmetry cancels in the untagged decay rate. For the $B_S^0$ system, the rapid oscillations, due to the large value of $\Delta m_S = 17.768 \pm 0.024$ ps⁻¹ [52], reduce the effect of a production asymmetry, reported as $A_P(B_S^0) = (7 \pm 5)\%$ in Ref. [51], to a negligible level. Hence, no corresponding systematic uncertainty is assigned.

There is a 0.02% relative uncertainty on the lifetime measurements due to the uncer-
tainty on the length scale of LHCb [52], which is mainly determined by the VELO
modules $z$ positions. These are evaluated by a survey, having an accuracy of 0.1 mm over
the full length of the VELO (1000 mm), and refined through a track-based alignment.
The alignment procedure is more precise for the first track hits, that are less affected by
multiple scattering and whose distribution of $z$ positions have an RMS of 100 mm. In
this region, the differences between the module positions obtained from the survey and
track-based alignment are within 0.02 mm, which is taken as systematic uncertainty. The
systematic uncertainty related to the momentum scale calibration affects both the $b$ hadron
candidate mass and momentum and, therefore, cancels when computing the decay time.

The systematic uncertainty related to the choice of 45 fs for the width of the decay-time
resolution function (65 fs in the case of $B^0 \rightarrow J/\psi K_S^0$ and $\Lambda_b^0 \rightarrow J/\psi \Lambda$) is evaluated by
changing the width by $\pm 15$ fs and repeating the fit. This change in width is larger than
the estimated uncertainty on the resolution and leads to a negligible change in the fit

<footnote n="2">²Flavour-specific means that the final state is only accessible via the decay of a $B_{(s)}^0$ meson and accessible by a meson originally produced as a $\bar{B}_{(s)}^0$ only via oscillation.</footnote>
<table><caption>Table 4: Statistical and systematic uncertainties (in femtoseconds) for the values of the b-hadron lifetimes. The total systematic uncertainty is obtained by combining the individual contributions in quadrature.</caption><thead><tr><th>Source</th><th>τ<sub>B<sup>+</sup>→J/ψK<sup>+</sup></sub></th><th>τ<sub>B<sup>0</sup>→J/ψK<sup>*0</sup></sub></th><th>τ<sub>B<sup>0</sup>→J/ψK<sup>0</sup><sub>S</sub></sub></th><th>τ<sub>Λ<sup>0</sup><sub>b</sub>→J/ψΛ</sub></th><th>τ<sub>B<sup>0</sup><sub>S</sub>→J/ψφ</sub></th></tr></thead><tbody><tr><td>Statistical uncertainty</td><td>3.5</td><td>6.1</td><td>12.8</td><td>26.5</td><td>11.4</td></tr><tr><td>VELO reconstruction</td><td>2.0</td><td>2.3</td><td>0.9</td><td>0.5</td><td>2.3</td></tr><tr><td>Simulation sample size</td><td>1.7</td><td>2.3</td><td>2.9</td><td>3.7</td><td>2.4</td></tr><tr><td>Mass-time correlation</td><td>1.4</td><td>1.8</td><td>2.1</td><td>3.0</td><td>0.7</td></tr><tr><td>Trigger and selection eff.</td><td>1.1</td><td>1.2</td><td>2.0</td><td>2.0</td><td>2.5</td></tr><tr><td>Background modelling</td><td>0.1</td><td>0.2</td><td>2.2</td><td>2.1</td><td>0.4</td></tr><tr><td>Mass modelling</td><td>0.1</td><td>0.2</td><td>0.4</td><td>0.2</td><td>0.5</td></tr><tr><td>Peaking background</td><td>—</td><td>—</td><td>0.3</td><td>1.1</td><td>0.4</td></tr><tr><td>Effective lifetime bias</td><td>—</td><td>—</td><td>—</td><td>—</td><td>1.6</td></tr><tr><td>B<sup>0</sup> production asym.</td><td>—</td><td>—</td><td>1.1</td><td>—</td><td>—</td></tr><tr><td>LHCb length scale</td><td>0.4</td><td>0.3</td><td>0.3</td><td>0.3</td><td>0.3</td></tr><tr><td>Total systematic</td><td>3.2</td><td>3.9</td><td>4.9</td><td>5.7</td><td>4.6</td></tr></tbody></table>

results. Consequently, no systematic uncertainty is assigned. Furthermore, to test the sensitivity of the lifetimes to potential mismodelling of the long tails in the resolution, the resolution model is changed from a single Gaussian function to a sum of two or three Gaussian functions with parameters fixed from simulation. Repeating the fit with the new resolution model causes no significant change to the lifetimes and no systematic uncertainty is assigned. The lifetimes are insensitive to the offset, $\Delta$, in the resolution model.

Several consistency checks are performed to study the stability of the lifetimes, by comparing the results obtained using different subsets of the data in terms of magnet polarity, data taking period, b-hadron and track kinematic variables, number of PVs in the event and track multiplicity. In all cases, no trend is observed and all lifetimes are compatible with the nominal results.

The majority of the systematic uncertainties described above can be propagated to the lifetime ratio measurements in Table 7. However, some of the uncertainties are correlated between the individual lifetimes and cancel in the ratio. For the first set of ratios and for $\Delta\Gamma_d/\Gamma_d$, the systematic uncertainty from the VELO-reconstruction efficiency weights and the LHCb length scale are considered as fully correlated. For the second set of ratios, other systematic uncertainties, as indicated in Table 5, cancel, since the ratio is formed from lifetimes measured using the same decay mode. In contrast to the situation for the measurement of the $B^0$ lifetime in the $B^0 \to J/\psi K^{*0}$ mode, the $B^0$ production asymmetry does lead to a systematic uncertainty on the measurement of $\tau_{B^0 \to J/\psi K^{*0}}/\tau_{\bar{B}^0 \to J/\psi K^{*0}}$ since terms like $A_P \cos(\Delta m_d t)$ do not cancel in the decay rates describing the decays of $B^0$ and $\bar{B}^0$ mesons to $J/\psi K^{*0}$ and $J/\psi K^{*0}$ final states. The effect of candidates where the flavour, via the particle identification of the decay products, has not been correctly assigned is investigated and found to be negligible.
<table><caption>Table 5: Statistical and systematic uncertainties (in units of $10^{-3}$) for the lifetime ratios and $\Delta\Gamma_d/\Gamma_d$. For brevity, $\tau_{B^0}$ ($\tau_{\bar{B}^0}$) corresponds to the measurement of $\tau_{B^0 \to J/\psi K^{*0}}$ ($\tau_{\bar{B}^0 \to J/\psi K^{*0}}$). The total systematic uncertainty is obtained by combining the individual contributions in quadrature.</caption><thead><tr><th>Source</th><th>τ<sub>B<sup>+</sup></sub>/τ<sub>B<sup>0</sup></sub></th><th>τ<sub>B<sub>s</sub><sup>0</sup></sub>/τ<sub>B<sup>0</sup></sub></th><th>τ<sub>Λ<sub>b</sub><sup>0</sup></sub>/τ<sub>B<sup>0</sup></sub></th><th>τ<sub>B<sup>+</sup></sub>/τ<sub>B<sup>-</sup></sub></th><th>τ<sub>Λ<sub>b</sub><sup>0</sup></sub>/τ<sub>Λ<sub>b</sub><sup>0</sup></sub></th><th>τ<sub>B<sup>0</sup></sub>/τ<sub>B<sup>0</sup></sub></th><th>ΔΓ<sub>d</sub>/Γ<sub>d</sub></th></tr></thead><tbody><tr><td>Statistical uncertainty</td><td>5.0</td><td>8.5</td><td>18.0</td><td>4.0</td><td>35.0</td><td>8.0</td><td>25.0</td></tr><tr><td>VELO reconstruction</td><td>1.6</td><td>1.7</td><td>1.1</td><td>—</td><td>—</td><td>—</td><td>4.1</td></tr><tr><td>Simulation sample size</td><td>2.0</td><td>2.2</td><td>2.8</td><td>2.1</td><td>5.3</td><td>3.0</td><td>6.3</td></tr><tr><td>Mass-time correlation</td><td>1.6</td><td>1.2</td><td>2.3</td><td>—</td><td>—</td><td>—</td><td>4.7</td></tr><tr><td>Trigger and selection eff.</td><td>1.1</td><td>1.8</td><td>1.5</td><td>—</td><td>—</td><td>—</td><td>4.0</td></tr><tr><td>Background modelling</td><td>0.3</td><td>0.1</td><td>1.5</td><td>0.2</td><td>3.0</td><td>1.4</td><td>3.8</td></tr><tr><td>Mass modelling</td><td>0.2</td><td>0.4</td><td>0.2</td><td>0.1</td><td>0.2</td><td>0.2</td><td>0.8</td></tr><tr><td>Peaking background</td><td>—</td><td>0.3</td><td>0.7</td><td>—</td><td>—</td><td>—</td><td>0.5</td></tr><tr><td>Effective lifetime bias</td><td>—</td><td>1.0</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>B<sup>0</sup> production asym.</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>8.5</td><td>1.9</td></tr><tr><td>Total systematic</td><td>3.2</td><td>3.7</td><td>4.4</td><td>2.1</td><td>6.1</td><td>9.1</td><td>10.7</td></tr></tbody></table>

<table><caption>Table 6: Fit results for the $B^+$, $B^0$, $B_s^0$ mesons and $\Lambda_b^0$ baryon lifetimes. The first uncertainty is statistical and the second is systematic.</caption><thead><tr><th>Lifetime</th><th>Value [ps]</th></tr></thead><tbody><tr><td>τ<sub>B<sup>+</sup>→J/ψK<sup>+</sup></sub></td><td>1.637 ± 0.004 ± 0.003</td></tr><tr><td>τ<sub>B<sup>0</sup>→J/ψK<sup>*0</sup></sub></td><td>1.524 ± 0.006 ± 0.004</td></tr><tr><td>τ<sub>B<sup>0</sup>→J/ψK<sub>S</sub><sup>0</sup></sub></td><td>1.499 ± 0.013 ± 0.005</td></tr><tr><td>τ<sub>Λ<sub>b</sub><sup>0</sup>→J/ψΛ</sub></td><td>1.415 ± 0.027 ± 0.006</td></tr><tr><td>τ<sub>B<sub>s</sub><sup>0</sup>→J/ψφ</sub></td><td>1.480 ± 0.011 ± 0.005</td></tr></tbody></table>

# 7 Results and conclusions

The measured b-hadron lifetimes are reported in Table 6. All results are compatible with existing world averages [13]. The reported $\tau_{\Lambda_b^0 \to J/\psi \Lambda}$ is smaller by approximately $2\sigma$ than a previous measurements from LHCb [8]. With the exception of the $\Lambda_b^0 \to J/\psi \Lambda$ channel, these are the single most precise measurements of the b-hadron lifetimes. The $B_s^0$ meson effective lifetime is measured using the same data set as used in Ref. [46] for the measurement of the $B_s^0$ meson mixing parameters and polarisation amplitudes in $B_s^0 \to J/\psi \phi$ decays. The $B_s^0$ meson effective lifetime computed from these quantities is compatible with the lifetime reported in this paper and a combination of the two results is, therefore, inappropriate.

Table 7 reports the ratios of the $B^+$, $B_s^0$ and $\Lambda_b^0$ lifetimes to the $B^0$ lifetime measured in the flavour-specific $B^0 \to J/\psi K^{*0}$ channel. This decay mode provides a better normalisation than the $B^0 \to J/\psi K_S^0$ channel due to the lower statistical uncertainty on the $B^0$ meson lifetime and due to the fact that the $B^0 \to J/\psi K^{*0}$ lifetime only depends quadratically on $\Delta\Gamma_d/\Gamma_d$, as shown in Eq. (7). The statistical and systematic uncertainties from
<table><caption>Table 7: Lifetime ratios for the $B^+$, $B^0$, $B_s^0$ mesons and $\Lambda_b^0$ baryon. The first uncertainty is statistical and the second is systematic.</caption><thead><tr><th>Ratio</th><th>Value</th></tr></thead><tbody><tr><td>τ<sub>B<sup>+</sup></sub>/τ<sub>B<sup>0</sup>→J/ψK<sup>*0</sup></sub></td><td>1.074 ± 0.005 ± 0.003</td></tr><tr><td>τ<sub>B<sub>s</sub><sup>0</sup></sub>/τ<sub>B<sup>0</sup>→J/ψK<sup>*0</sup></sub></td><td>0.971 ± 0.009 ± 0.004</td></tr><tr><td>τ<sub>Λ<sub>b</sub><sup>0</sup></sub>/τ<sub>B<sup>0</sup>→J/ψK<sup>*0</sup></sub></td><td>0.929 ± 0.018 ± 0.004</td></tr><tr><td>τ<sub>B<sup>+</sup></sub>/τ<sub>B<sup>-</sup></sub></td><td>1.002 ± 0.004 ± 0.002</td></tr><tr><td>τ<sub>Λ<sub>b</sub><sup>0</sup></sub>/τ<sub>Λ<sub>b</sub><sup>0</sup></sub></td><td>0.940 ± 0.035 ± 0.006</td></tr><tr><td>τ<sub>B<sup>0</sup>→J/ψK<sup>*0</sup></sub>/τ<sub>B<sup>0</sup>→J/ψK<sup>*0</sup></sub></td><td>1.000 ± 0.008 ± 0.009</td></tr></tbody></table>

the absolute lifetime measurements are propagated to the ratios, taking into account the correlations between the systematic uncertainties. All ratios are consistent with SM predictions [15, 22–25, 30–32] and with previous measurements [13]. Furthermore, the ratios $\tau_{B^+}/\tau_{B^-}$, $\tau_{\Lambda_b^0}/\tau_{\Lambda_b^0}$ and $\tau_{B^0\to J/\psi K^{*0}}/\tau_{B^0\to J/\psi K^{*0}}$ are reported. Measuring any of these different from unity would indicate a violation of CPT invariance or, for $B^0 \to J/\psi K^{*0}$ decays, could also indicate that $\Delta\Gamma_d$ is non-zero and $B^0 \to J/\psi K^{*0}$ is not 100% flavour-specific. No deviation from unity of these ratios is observed.

The effective lifetimes of $B^0 \to J/\psi K^{*0}$ and $B^0 \to J/\psi K_S^0$ decays are used to measure $\Delta\Gamma_d/\Gamma_d$. Flavour-specific final states such as $B^0 \to J/\psi K^{*0}$ have $A_{\Delta\Gamma_d}^{B^0\to J/\psi K^{*0}} = 0$, while $A_{\Delta\Gamma_d}^{B^0\to J/\psi K_S^0} = \cos(2\beta)$ to a good approximation in the SM, where $\beta = \arg[-(V_{cd}V_{cb}^*)/(V_{td}V_{tb}^*)]$ is one of the CKM unitarity triangle angles. Hence, the two effective lifetimes can be expressed as

$$ \tau_{B^0 \to J/\psi K^{*0}} = \frac{1}{\Gamma_d} \frac{1}{1 - y_d^2} (1 + y_d^2), \qquad (7) $$

$$ \tau_{B^0 \to J/\psi K_S^0} = \frac{1}{\Gamma_d} \frac{1}{1 - y_d^2} \left( \frac{1 + 2 \cos(2\beta) y_d + y_d^2}{1 + \cos(2\beta) y_d} \right). \qquad (8) $$

Using the effective lifetimes reported in Table 6 and $\beta = (21.5_{-0.7}^{+0.8})^\circ$ [13], a fit of $\Delta\Gamma_d$ and $\Gamma_d$ to the expressions in Eq. (7) and Eq. (8) leads to

$$ \Gamma_d = 0.656 \pm 0.003 \pm 0.002 \text{ ps}^{-1}, \qquad (9) $$

$$ \Delta\Gamma_d = -0.029 \pm 0.016 \pm 0.007 \text{ ps}^{-1}, \qquad (10) $$

where the first uncertainty is statistical and the second is systematic. The correlation coefficient between $\Delta\Gamma_d$ and $\Gamma_d$ is 0.43 when including statistical and systematic uncertainties. The combination gives

$$ \frac{\Delta\Gamma_d}{\Gamma_d} = -0.044 \pm 0.025 \pm 0.011, \qquad (11) $$

consistent with the SM expectation [14, 15] and the current world-average value [13].
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
and Region Auvergne (France); BMBF, DFG, HGF and MPG (Germany); SFI (Ireland);
INFN (Italy); FOM and NWO (The Netherlands); SCSR (Poland); MEN/IFA (Romania);
MinES, Rosatom, RFBR and NRC "Kurchatov Institute" (Russia); MinECo, XuntaGal
and GENCAT (Spain); SNSF and SER (Switzerland); NAS Ukraine (Ukraine); STFC
(United Kingdom); NSF (USA). We also acknowledge the support received from the ERC
under FP7. The Tier1 computing centres are supported by IN2P3 (France), KIT and
BMBF (Germany), INFN (Italy), NWO and SURF (The Netherlands), PIC (Spain),
GridPP (United Kingdom). We are indebted to the communities behind the multiple open
source software packages we depend on. We are also thankful for the computing resources
and the access to software R&D tools provided by Yandex LLC (Russia).

References

[1] V. A. Khoze and M. A. Shifman, *Heavy quarks*, Sov. Phys. Usp. **26** (1983) 387.

[2] M. A. Shifman and M. Voloshin, *Preasymptotic effects in inclusive weak decays of charmed particles*, Sov. J. Nucl. Phys. **41** (1985) 120.

[3] M. A. Shifman and M. Voloshin, *Hierarchy of lifetimes of charmed and beautiful hadrons*, Sov. Phys. JETP **64** (1986) 698.

[4] I. I. Bigi, N. Uraltsev, and A. Vainshtein, *Nonperturbative corrections to inclusive beauty and charm decays: QCD versus phenomenological models*, Phys. Lett. **B293** (1992) 430, arXiv:hep-ph/9207214.

[5] I. I. Bigi, *The QCD perspective on lifetimes of heavy-flavor hadrons*, arXiv:hep-ph/9508408.

[6] N. Uraltsev, *Heavy quark expansion in beauty and its decays*, arXiv:hep-ph/9804275.

[7] M. Neubert, *B decays and the heavy-quark expansion*, Adv. Ser. Direct. High Energy Phys. **15** (1998) 239, arXiv:hep-ph/9702375.

[8] LHCb collaboration, R. Aaij et al., *Precision measurement of the $\Lambda_b^0$ baryon lifetime*, Phys. Rev. Lett. **111** (2013) 102003, arXiv:1307.2476.

[9] K. Hartkorn and H. Moser, *A new method of measuring $\Delta\Gamma/\Gamma$ in the $B_s^0$ anti-$B_s^0$ system*, Eur. Phys. J. **C8** (1999) 381.
[10] I. Dunietz, R. Fleischer, and U. Nierste, *In pursuit of new physics with $B_s$ decays*, Phys. Rev. **D63** (2001) 114015, [arXiv:hep-ph/0012219](https://arxiv.org/abs/0012219).

[11] R. Fleischer and R. Knegjens, *Effective lifetimes of $B_s$ decays and their constraints on the $B_s^0\bar{B}_s^0$ mixing parameters*, Eur. Phys. J. **C71** (2011) 1789, [arXiv:1109.5115](https://arxiv.org/abs/1109.5115).

[12] U. Nierste, *CP asymmetry in flavor-specific B decays*, [arXiv:hep-ph/0406300](https://arxiv.org/abs/0406300).

[13] Heavy Flavor Averaging Group, Y. Amhis et al., *Averages of b-hadron, c-hadron, and $\tau$-lepton properties as of early 2012*, [arXiv:1207.1158](https://arxiv.org/abs/1207.1158), updated results and plots available at: [http://www.slac.stanford.edu/xorg/hfag/](http://www.slac.stanford.edu/xorg/hfag/).

[14] A. Lenz and U. Nierste, *Theoretical update of $B_s - \bar{B}_s$ mixing*, JHEP **06** (2007) 072, [arXiv:hep-ph/0612167](https://arxiv.org/abs/0612167).

[15] A. Lenz and U. Nierste, *Numerical updates of lifetimes and mixing parameters of B mesons*, [arXiv:1102.4274](https://arxiv.org/abs/1102.4274).

[16] BaBar collaboration, B. Aubert et al., *Limits on the decay-rate difference of neutral B mesons and on CP, T, and CPT violation in $B^0\bar{B}^0$ oscillations*, Phys. Rev. Lett. **92** (2004) 181801, [arXiv:hep-ex/0311037](https://arxiv.org/abs/0311037).

[17] BaBar collaboration, B. Aubert et al., *Limits on the decay-rate difference of neutral B mesons and on CP, T, and CPT violation in $B^0\bar{B}^0$ oscillations*, Phys. Rev. **D70** (2004) 012007, [arXiv:hep-ex/0403002](https://arxiv.org/abs/0403002).

[18] Belle collaboration, T. Higuchi et al., *Search for time-dependent CPT violation in hadronic and semileptonic B decays*, Phys. Rev. **D85** (2012) 071105, [arXiv:1203.0930](https://arxiv.org/abs/1203.0930).

[19] G. Borissov and B. Hoeneisen, *Understanding the like-sign dimuon charge asymmetry in $p\bar{p}$ collisions*, Phys. Rev. **D87** (2013) 074020, [arXiv:1303.0175](https://arxiv.org/abs/1303.0175).

[20] D0 collaboration, V. M. Abazov et al., *Study of CP-violating charge asymmetries of single muons and like-sign dimuons in $p\bar{p}$ collisions*, Phys. Rev. **D89** (2014) 012002, [arXiv:1310.0447](https://arxiv.org/abs/1310.0447).

[21] T. Gershon, $\Delta\Gamma_d$: a forgotten null test of the standard model, J. Phys. **G38** (2011) 015007, [arXiv:1007.5135](https://arxiv.org/abs/1007.5135).

[22] M. Beneke et al., *The $B^+ - B_d^0$ lifetime difference beyond leading logarithms*, Nucl. Phys. **B639** (2002) 389, [arXiv:hep-ph/0202106](https://arxiv.org/abs/0202106).

[23] E. Franco, V. Lubicz, F. Mescia, and C. Tarantino, *Lifetime ratios of beauty hadrons at the next-to-leading order in QCD*, Nucl. Phys. **B633** (2002) 212, [arXiv:hep-ph/0203089](https://arxiv.org/abs/0203089).
[24] M. Beneke, G. Buchalla, and I. Dunietz, *Width Difference in the $B_s - \bar{B}_s$ System*, Phys. Rev. **D54** (1996) 4419, arXiv:hep-ph/9605259.

[25] Y.-Y. Keum and U. Nierste, *Probing penguin coefficients with the lifetime ratio $\tau(B_s)/\tau(B_d)$*, Phys. Rev. **D57** (1998) 4282, arXiv:hep-ph/9710512.

[26] N. Uraltsev, *On the problem of boosting nonleptonic b baryon decays*, Phys. Lett. **B376** (1996) 303, arXiv:hep-ph/9602324.

[27] I. I. Bigi, M. A. Shifman, and N. Uraltsev, *Aspects of heavy quark theory*, Ann. Rev. Nucl. Part. Sci. **47** (1997) 591, arXiv:hep-ph/9703290.

[28] D. Pirjol and N. Uraltsev, *Four fermion heavy quark operators and light current amplitudes in heavy flavor hadrons*, Phys. Rev. **D59** (1999) 034012, arXiv:hep-ph/9805488.

[29] M. Voloshin, *Reducing model dependence of spectator effects in inclusive decays of heavy baryons*, Phys. Rev. **D61** (2000) 074026, arXiv:hep-ph/9908455.

[30] C. Tarantino, *Beauty hadron lifetimes and B meson CP violation parameters from lattice QCD*, Eur. Phys. J. **C33** (2004) S895, arXiv:hep-ph/0310241.

[31] F. Gabbiani, A. I. Onishchenko, and A. A. Petrov, $\Lambda_b^0$ lifetime puzzle in heavy quark expansion, Phys. Rev. **D68** (2003) 114006, arXiv:hep-ph/0303235.

[32] F. Gabbiani, A. I. Onishchenko, and A. A. Petrov, *Spectator effects and lifetimes of heavy hadrons*, Phys. Rev. **D70** (2004) 094031, arXiv:hep-ph/0407004.

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

[34] M. Adinolfi et al., *Performance of the LHCb RICH detector at the LHC*, Eur. Phys. J. **C73** (2013) 2431, arXiv:1211.6759.

[35] A. A. Alves Jr et al., *Performance of the LHCb muon system*, JINST **8** (2013) P02022, arXiv:1211.1346.

[36] R. Aaij et al., *The LHCb trigger and its performance in 2011*, JINST **8** (2013) P04022, arXiv:1211.3055.

[37] T. Sjöstrand, S. Mrenna, and P. Skands, *PYTHIA 6.4 physics and manual*, JHEP **05** (2006) 026, arXiv:hep-ph/0603175.

[38] I. Belyaev et al., *Handling of the generation of primary events in GAUSS*, the LHCb simulation framework, Nuclear Science Symposium Conference Record (NSS/MIC) **IEEE** (2010) 1155.
[39] D. J. Lange, *The EVTGEN particle decay simulation package*, Nucl. Instrum. Meth. **A462** (2001) 152.

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

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

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

[43] W. D. Hulsbergen, Decay chain fitting with a Kalman filter, Nucl. Instrum. Meth. **A552** (2005) 566, arXiv:physics/0503191.

[44] Particle Data Group, J. Beringer et al., Review of particle physics, Phys. Rev. **D86** (2012) 010001, and 2013 partial update for the 2014 edition.

[45] O. Callot, FastVelo, a fast and efficient pattern recognition package for the Velo, LHCb-PUB-2011-001.

[46] LHCb collaboration, R. Aaij et al., Measurement of CP-violation and the $B_s^0$-meson decay width difference with $B_s^0 \to J/\psi K^+ K^-$ and $B_s^0 \to J/\psi \pi^+ \pi^-$ decays, Phys. Rev. **D87** (2013) 112010, arXiv:1304.2600.

[47] Y. Xie, sFit: a method for background subtraction in maximum likelihood fit, arXiv:0905.0724.

[48] M. Pivk and F. R. Le Diberder, sPlot: a statistical tool to unfold data distributions, Nucl. Instrum. Meth. **A555** (2005) 356, arXiv:physics/0402083.

[49] D. M. Santos and F. Dupertuis, Mass distributions marginalized over per-event errors, submitted to Nucl. Instrum. Meth. A (2013) arXiv:1312.5000.

[50] LHCb collaboration, R. Aaij et al., Amplitude analysis and branching fraction measurement of $\bar{B}_s^0 \to J/\psi K^+ K^-$, Phys. Rev. **D87** (2013) 072004, arXiv:1302.1213.

[51] LHCb collaboration, R. Aaij et al., First measurement of time-dependent CP violation in $B_s^0 \to K^+ K^-$ decays, JHEP **10** (2013) 183, arXiv:1308.1428.

[52] LHCb collaboration, R. Aaij et al., Precision measurement of the $B_s^0 - \bar{B}_s^0$ oscillation frequency $\Delta m_s$ in the decay $B_s^0 \to D_s^+ \pi^-$, New J. Phys. **15** (2013) 053021, arXiv:1304.4741.