EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) CERN-PH-EP-2011-186 LHCb-PAPER-2011-025 2 1 0 Search for the rare decays 2 n 0 + − 0 + − B → µ µ and B → µ µ a J s 5 1 ] x e p- The LHCb Collaboration1 e h [ 3 v 0 Abstract 0 6 1 A search for the decays B0 → µ+µ− and B0 → µ+µ− is performed with 0.37 fb−1 of . √ s 2 pp collisions at s = 7 TeV collected by the LHCb experiment in 2011. The upper 1 limits on the branching fractions are B(B0 → µ+µ−) < 1.6 × 10−8 and B(B0 → µ+µ−) 1 s 1 < 3.6×10−9 at95%confidencelevel. AcombinationoftheseresultswiththeLHCblimits : v obtained with the 2010 dataset leads to B(B0 → µ+µ−) < 1.4×10−8 and B(B0 → µ+µ−) s Xi < 3.2×10−9 at 95% confidence level. r a Keywords: LHC, b-hadron, FCNC, rare decays, leptonic decays. 1Authors are listed on the following pages. R. Aaij23, C. Abellan Beteta35,n, B. Adeva36, M. Adinolfi42, C. Adrover6, A. Affolder48, Z. Ajaltouni5, J. Albrecht37, F. Alessio37, M. Alexander47, G. Alkhazov29, P. Al- varez Cartelle36, A.A. Alves Jr22, S. Amato2, Y. Amhis38, J. Anderson39, R.B. Appleby50, O. Aquines Gutierrez10, F. Archilli18,37, L. Arrabito53, A. Artamonov 34, M. Artuso52,37, E. Aslanides6, G. Auriemma22,m, S. Bachmann11, J.J. Back44, D.S. Bailey50, V. Balagura30,37, W. Baldini16, R.J. Barlow50, C. Barschel37, S. Barsuk7, W. Barter43, A. Bates47, C. Bauer10, Th. Bauer23, A. Bay38, I. Bediaga1, S. Belogurov30, K. Belous34, I. Belyaev30,37, E. Ben-Haim8, M. Benayoun8, G. Bencivenni18, S. Benson46, J. Benton42, R. Bernet39, M.-O. Bettler17, M. van Beuzekom23, A. Bien11, S. Bifani12, T. Bird50, A. Bizzeti17,h, P.M. Bjørnstad50, T. Blake37, F. Blanc38, C. Blanks49, J. Blouw11, S. Blusk52, A. Bobrov33, V. Bocci22, A. Bondar33, N. Bondar29, W. Bonivento15, S. Borghi47,50, A. Borgia52, T.J.V. Bowcock48, C. Bozzi16, T. Brambach9, J. van den Brand24, J. Bressieux38, D. Brett50, M. Britsch10, T. Britton52, N.H. Brook42, H. Brown48, A. Bu¨chler-Germann39, I. Burducea28, A. Bursche39, J. Buytaert37, S. Cadeddu15, O. Callot7, M. Calvi20,j, M. Calvo Gomez35,n, A. Camboni35, P. Campana18,37, A. Carbone14, G. Carboni21,k, R. Cardinale19,i,37, A. Cardini15, L. Carson49, K. Carvalho Akiba2, G. Casse48, M. Cattaneo37, Ch. Cauet9, M. Charles51, Ph. Charpentier37, N. Chiapolini39, K. Ciba37, X. Cid Vidal36, G. Ciezarek49, P.E.L. Clarke46,37, M. Clemencic37, H.V. Cliff43, J. Closier37, C. Coca28, V. Coco23, J. Cogan6, P. Collins37, A. Comerma-Montells35, F. Constantin28, G. Conti38, A. Contu51, A. Cook42, M. Coombes42, G. Corti37, G.A. Cowan38, R. Currie46, B. D’Almagne7, C. D’Ambrosio37, P. David8, P.N.Y. David23, I. De Bonis4, S. De Capua21,k, M. De Cian39, F. De Lorenzi12, J.M. De Miranda1, L. De Paula2, P. De Simone18, D. Decamp4, M. Deckenhoff9, H. Degaudenzi38,37, M. Deissenroth11, L. Del Buono8, C. Deplano15, D. Derkach14,37, O. Deschamps5, F. Dettori24, J. Dickens43, H. Dijkstra37, P. Diniz Batista1, F. Domingo Bonal35,n, S. Donleavy48, F. Dordei11, P. Dornan49, A. Dosil Su´arez36, D. Dossett44, A. Dovbnya40, F. Dupertuis38, R. Dzhelyadin34, A. Dziurda25, S. Easo45, U. Egede49, V. Egorychev30, S. Eidelman33, D. van Eijk23, F. Eisele11, S. Eisenhardt46, R. Ekelhof9, L. Eklund47, Ch. Elsasser39, D. Elsby55, D. Esperante Pereira36, L. Est`eve43, A. Falabella16,14,e, E. Fanchini20,j, C. F¨arber11, G. Fardell46, C. Farinelli23, S. Farry12, V. Fave38, V. Fernandez Albor36, M. Ferro-Luzzi37, S. Filippov32, C. Fitzpatrick46, M. Fontana10, F. Fontanelli19,i, R. Forty37, M. Frank37, C. Frei37, M. Frosini17,f,37, S. Furcas20, A. Gallas Torreira36, D. Galli14,c, M. Gandelman2, P. Gandini51, Y. Gao3, J-C. Garnier37, J. Garofoli52, J. Garra Tico43, L. Garrido35, D. Gascon35, C. Gaspar37, N. Gauvin38, M. Gersabeck37, T. Gershon44,37, Ph. Ghez4, V. Gibson43, V.V. Gligorov37, C. G¨obel54, D. Golubkov30, A. Golutvin49,30,37, A. Gomes2, H. Gordon51, M. Grabalosa Ga´ndara35, R. Graciani Diaz35, L.A. Granado Cardoso37, E. Graug´es35, G. Graziani17, A. Grecu28, E. Greening51, S. Gregson43, B. Gui52, E.Gushchin32, Yu.Guz34, T.Gys37, G.Haefeli38, C.Haen37, S.C.Haines43, T.Hampson42, S.Hansmann-Menzemer11, R.Harji49, N.Harnew51, J.Harrison50, P.F.Harrison44, J.He7, V. Heijne23, K. Hennessy48, P. Henrard5, J.A. Hernando Morata36, E. van Herwijnen37, E. Hicks48, K. Holubyev11, P. Hopchev4, W. Hulsbergen23, P. Hunt51, T. Huse48, R.S. Huston12, D. Hutchcroft48, D. Hynds47, V. Iakovenko41, P. Ilten12, J. Imong42, ii R. Jacobsson37, A. Jaeger11, M. Jahjah Hussein5, E. Jans23, F. Jansen23, P. Jaton38, B. Jean-Marie7, F. Jing3, M. John51, D. Johnson51, C.R. Jones43, B. Jost37, M. Kaballo9, S. Kandybei40, M. Karacson37, T.M. Karbach9, J. Keaveney12, I.R. Kenyon55, U. Kerzel37, T. Ketel24, A. Keune38, B. Khanji6, Y.M. Kim46, M. Knecht38, P. Koppenburg23, A. Kozlinskiy23, L. Kravchuk32, K. Kreplin11, M. Kreps44, G. Krocker11, P. Krokovny11, F. Kruse9, K. Kruzelecki37, M. Kucharczyk20,25,37,j, T. Kvaratskheliya30,37, V.N. La Thi38, D. Lacarrere37, G. Lafferty50, A. Lai15, D. Lambert46, R.W. Lambert24, E. Lanciotti37, G. Lanfranchi18, C. Langenbruch11, T. Latham44, C. Lazzeroni55, R. Le Gac6, J. van Leerdam23, J.-P. Lees4, R. Lef`evre5, A. Leflat31,37, J. Lefran¸cois7, O. Leroy6, T. Lesiak25, L. Li3, L. Li Gioi5, M. Lieng9, M. Liles48, R. Lindner37, C. Linn11, B. Liu3, G. Liu37, J.H. Lopes2, E. Lopez Asamar35, N. Lopez-March38, H. Lu38,3, J. Luisier38, A. Mac Raighne47, F. Machefert7, I.V. Machikhiliyan4,30, F. Maciuc10, O. Maev29,37, J. Magnin1, S. Malde51, R.M.D. Mamunur37, G. Manca15,d, G. Mancinelli6, N. Mangiafave43, U. Marconi14, R. Ma¨rki38, J. Marks11, G. Martellotti22, A. Martens8, L. Martin51, A. Mart´ın S´anchez7, D. Martinez Santos37, A. Massafferri1, Z. Mathe12, C. Matteuzzi20, M. Matveev29, E. Maurice6, B. Maynard52, A. Mazurov16,32,37, G. McGregor50, R. McNulty12, C. Mclean14, M. Meissner11, M. Merk23, J. Merkel9, R. Messi21,k, S. Miglioranzi37, D.A. Milanes13,37, M.-N. Minard4, J. Molina Rodriguez54, S. Monteil5, D. Moran12, P. Morawski25, R. Mountain52, I. Mous23, F. Muheim46, K. Mu¨ller39, R. Muresan28,38, B. Muryn26, B. Muster38, M. Musy35, J. Mylroie- Smith48, P. Naik42, T. Nakada38, R. Nandakumar45, I. Nasteva1, M. Nedos9, M. Needham46, N. Neufeld37, C. Nguyen-Mau38,o, M. Nicol7, V. Niess5, N. Nikitin31, A. Nomerotski51, A. Novoselov34, A. Oblakowska-Mucha26, V. Obraztsov34, S. Oggero23, S. Ogilvy47, O. Okhrimenko41, R. Oldeman15,d, M. Orlandea28, J.M. Otalora Goicochea2, P. Owen49, K. Pal52, J. Palacios39, A. Palano13,b, M. Palutan18, J. Panman37, A. Papanestis45, M. Pappagallo47, C. Parkes47,37, C.J. Parkinson49, G. Passaleva17, G.D. Patel48, M. Patel49, S.K. Paterson49, G.N. Patrick45, C. Patrignani19,i, C. Pavel- Nicorescu28, A. Pazos Alvarez36, A. Pellegrino23, G. Penso22,l, M. Pepe Altarelli37, S. Perazzini14,c, D.L. Perego20,j, E. Perez Trigo36, A. P´erez-Calero Yzquierdo35, P. Perret5, M. Perrin-Terrin6, G. Pessina20, A. Petrella16,37, A. Petrolini19,i, A. Phan52, E. Pi- catoste Olloqui35, B. Pie Valls35, B. Pietrzyk4, T. Pilaˇr44, D. Pinci22, R. Plackett47, S. Playfer46, M. Plo Casasus36, G. Polok25, A. Poluektov44,33, E. Polycarpo2, D. Popov10, B. Popovici28, C. Potterat35, A. Powell51, T. du Pree23, J. Prisciandaro38, V. Pugatch41, A. Puig Navarro35, W. Qian52, J.H. Rademacker42, B. Rakotomiaramanana38, M.S. Rangel2, I. Raniuk40, G. Raven24, S. Redford51, M.M. Reid44, A.C. dos Reis1, S. Ricciardi45, K. Rinnert48, D.A. Roa Romero5, P. Robbe7, E. Rodrigues47,50, F. Rodrigues2, P. Rodriguez Perez36, G.J. Rogers43, S. Roiser37, V. Romanovsky34, M. Rosello35,n, J. Rouvinet38, T. Ruf37, H. Ruiz35, G. Sabatino21,k, J.J. Sa- borido Silva36, N. Sagidova29, P. Sail47, B. Saitta15,d, C. Salzmann39, M. Sannino19,i, R. Santacesaria22, C. Santamarina Rios36, R. Santinelli37, E. Santovetti21,k, M. Sapunov6, A. Sarti18,l, C. Satriano22,m, A. Satta21, M. Savrie16,e, D. Savrina30, P. Schaack49, M. Schiller24, S. Schleich9, M. Schlupp9, M. Schmelling10, B. Schmidt37, O. Schneider38, A. Schopper37, M.-H. Schune7, R. Schwemmer37, B. Sciascia18, A. Sciubba18,l, iii M. Seco36, A. Semennikov30, K. Senderowska26, I. Sepp49, N. Serra39, J. Serrano6, P. Seyfert11, B. Shao3, M. Shapkin34, I. Shapoval40,37, P. Shatalov30, Y. Shcheglov29, T. Shears48, L. Shekhtman33, O. Shevchenko40, V. Shevchenko30, A. Shires49, R. Silva Coutinho44, T. Skwarnicki52, A.C. Smith37, N.A. Smith48, E. Smith51,45, K. Sobczak5, F.J.P. Soler47, A. Solomin42, F. Soomro18, B. Souza De Paula2, B. Spaan9, A. Sparkes46, P. Spradlin47, F. Stagni37, S. Stahl11, O. Steinkamp39, S. Stoica28, S. Stone52,37, B. Storaci23, M. Straticiuc28, U. Straumann39, V.K. Subbiah37, S. Swientek9, M. Szczekowski27, P. Szczypka38, T. Szumlak26, S. T’Jampens4, E. Teodorescu28, F. Teubert37, C. Thomas51, E. Thomas37, J. van Tilburg11, V. Tisserand4, M. Tobin39, S. Topp-Joergensen51, N. Torr51, E. Tournefier4,49, M.T. Tran38, A. Tsaregorodtsev6, N. Tuning23, M. Ubeda Garcia37, A. Ukleja27, P. Urquijo52, U. Uwer11, V. Vagnoni14, G. Valenti14, R. Vazquez Gomez35, P. Vazquez Regueiro36, S. Vecchi16, J.J. Velthuis42, M.Veltri17,g, B.Viaud7, I.Videau7, X.Vilasis-Cardona35,n, J.Visniakov36, A.Vollhardt39, D. Volyanskyy10, D. Voong42, A. Vorobyev29, H. Voss10, S. Wandernoth11, J. Wang52, D.R. Ward43, N.K. Watson55, A.D. Webber50, D. Websdale49, M. Whitehead44, D. Wiedner11, L. Wiggers23, G. Wilkinson51, M.P. Williams44,45, M. Williams49, F.F. Wilson45, J. Wishahi9, M. Witek25, W. Witzeling37, S.A. Wotton43, K. Wyllie37, Y. Xie46, F. Xing51, Z. Xing52, Z. Yang3, R. Young46, O. Yushchenko34, M. Zavertyaev10,a, F. Zhang3, L. Zhang52, W.C. Zhang12, Y. Zhang3, A. Zhelezov11, L. Zhong3, E. Zverev31, A. Zvyagin37. 1Centro Brasileiro de Pesquisas F´ısicas (CBPF), Rio de Janeiro, Brazil 2Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil 3Center for High Energy Physics, Tsinghua University, Beijing, China 4LAPP, Universit´e de Savoie, CNRS/IN2P3, Annecy-Le-Vieux, France 5Clermont Universit´e, Universit´e Blaise Pascal, CNRS/IN2P3, LPC, Clermont-Ferrand, France 6CPPM, Aix-Marseille Universit´e, CNRS/IN2P3, Marseille, France 7LAL, Universit´e Paris-Sud, CNRS/IN2P3, Orsay, France 8LPNHE, Universit´e Pierre et Marie Curie, Universit´e Paris Diderot, CNRS/IN2P3, Paris, France 9Fakult¨at Physik, Technische Universita¨t Dortmund, Dortmund, Germany 10Max-Planck-Institut fu¨r Kernphysik (MPIK), Heidelberg, Germany 11Physikalisches Institut, Ruprecht-Karls-Universit¨at Heidelberg, Heidelberg, Germany 12School of Physics, University College Dublin, Dublin, Ireland 13Sezione INFN di Bari, Bari, Italy 14Sezione INFN di Bologna, Bologna, Italy 15Sezione INFN di Cagliari, Cagliari, Italy 16Sezione INFN di Ferrara, Ferrara, Italy 17Sezione INFN di Firenze, Firenze, Italy 18Laboratori Nazionali dell’INFN di Frascati, Frascati, Italy 19Sezione INFN di Genova, Genova, Italy 20Sezione INFN di Milano Bicocca, Milano, Italy 21Sezione INFN di Roma Tor Vergata, Roma, Italy 22Sezione INFN di Roma La Sapienza, Roma, Italy iv 23Nikhef National Institute for Subatomic Physics, Amsterdam, The Netherlands 24Nikhef National Institute for Subatomic Physics and Vrije Universiteit, Amsterdam, The Netherlands 25Henryk Niewodniczanski Institute of Nuclear Physics Polish Academy of Sciences, Krac´ow, Poland 26AGH University of Science and Technology, Krac´ow, Poland 27Soltan Institute for Nuclear Studies, Warsaw, Poland 28Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest-Magurele, Romania 29Petersburg Nuclear Physics Institute (PNPI), Gatchina, Russia 30Institute of Theoretical and Experimental Physics (ITEP), Moscow, Russia 31Institute of Nuclear Physics, Moscow State University (SINP MSU), Moscow, Russia 32Institute for Nuclear Research of the Russian Academy of Sciences (INR RAN), Moscow, Russia 33Budker Institute of Nuclear Physics (SB RAS) and Novosibirsk State University, Novosibirsk, Russia 34Institute for High Energy Physics (IHEP), Protvino, Russia 35Universitat de Barcelona, Barcelona, Spain 36Universidad de Santiago de Compostela, Santiago de Compostela, Spain 37European Organization for Nuclear Research (CERN), Geneva, Switzerland 38Ecole Polytechnique F´ed´erale de Lausanne (EPFL), Lausanne, Switzerland 39Physik-Institut, Universita¨t Zu¨rich, Zu¨rich, Switzerland 40NSC Kharkiv Institute of Physics and Technology (NSC KIPT), Kharkiv, Ukraine 41Institute for Nuclear Research of the National Academy of Sciences (KINR), Kyiv, Ukraine 42H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom 43Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom 44Department of Physics, University of Warwick, Coventry, United Kingdom 45STFC Rutherford Appleton Laboratory, Didcot, United Kingdom 46School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom 47School of Physics and Astronomy, University of Glasgow, Glasgow, United Kingdom 48Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom 49Imperial College London, London, United Kingdom 50School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom 51Department of Physics, University of Oxford, Oxford, United Kingdom 52Syracuse University, Syracuse, NY, United States 53CC-IN2P3, CNRS/IN2P3, Lyon-Villeurbanne, France, associated member 54Pontif´ıcia Universidade Cat´olica do Rio de Janeiro (PUC-Rio), Rio de Janeiro, Brazil, associated to 2 55University of Birmingham, Birmingham, United Kingdom aP.N. Lebedev Physical Institute, Russian Academy of Science (LPI RAS), Moscow, Russia bUniversit`a di Bari, Bari, Italy cUniversit`a di Bologna, Bologna, Italy dUniversit`a di Cagliari, Cagliari, Italy eUniversit`a di Ferrara, Ferrara, Italy fUniversit`a di Firenze, Firenze, Italy gUniversit`a di Urbino, Urbino, Italy hUniversit`a di Modena e Reggio Emilia, Modena, Italy iUniversit`a di Genova, Genova, Italy v jUniversit`a di Milano Bicocca, Milano, Italy kUniversit`a di Roma Tor Vergata, Roma, Italy lUniversit`a di Roma La Sapienza, Roma, Italy mUniversit`a della Basilicata, Potenza, Italy nLIFAELS, La Salle, Universitat Ramon Llull, Barcelona, Spain oHanoi University of Science, Hanoi, Viet Nam vi 1 Introduction Measurements of low-energy processes can provide indirect constraints on particles that are too heavy to be produced directly. This is particularly true for Flavour Changing Neutral Current (FCNC) processes which are highly suppressed in the Standard Model (SM) and can only occur through higher-order diagrams. The SM predictions for the branchingfractionsoftheFCNCdecays2 B0 → µ+µ− andB0 → µ+µ− areB(B0 → µ+µ−) s s = (3.2±0.2)×10−9 and B(B0 → µ+µ−) = (0.10±0.01)×10−9 [1]. However, contributions from new processes or new heavy particles can significantly enhance these values. For example, within Minimal Supersymmetric extensions of the SM (MSSM), in the large tanβ regime, B(B0 → µ+µ−) is found to be approximately proportional to tan6β [2], s where tanβ is the ratio of the vacuum expectation values of the two neutral CP-even Higgs fields. The branching fractions could therefore be enhanced by orders of magnitude for large values of tanβ. The best published limits from the Tevatron are B(B0 → µ+µ−) < 5.1 × 10−8 s at 95% confidence level (CL) by the D0 collaboration using 6.1 fb−1 of data [3], and B(B0 → µ+µ−) < 6.0 × 10−9 at 95% CL by the CDF collaboration using 6.9 fb−1 of data [4]. In the same dataset the CDF collaboration observes an excess of B0 → µ+µ− s candidates compatible with B(B0 → µ+µ−) = (1.8+1.1)×10−8 and with an upper limit of s −0.9 B(B0 → µ+µ−) < 4.0×10−8 at 95% CL. The CMS collaboration has recently published s B(B0 → µ+µ−) < 1.9×10−8 at 95% CL and B(B0 → µ+µ−) < 4.6×10−9 at 95% CL using s 1.14 fb−1 of data [5]. The LHCb collaboration has published the limits [6] B(B0 → µ+µ−) s < 5.4 × 10−8 and B(B0 → µ+µ−) < 1.5 × 10−8 at 95% CL based on about 37 pb−1 of integrated luminosity collected in the 2010 run. This Letter presents an analysis of the data recorded by LHCb in the first half of 2011 which correspond to an integrated luminosity of ∼ 0.37 fb−1. The results of this analysis are then combined with those published from the 2010 dataset. 2 The LHCb detector The LHCb detector [7] is a single-arm forward spectrometer designed to study production and decays of hadrons containing b or c quarks. The detector consists of a vertex loca- tor (VELO) providing precise locations of primary pp interaction vertices and detached vertices of long lived hadrons. The momenta of charged particles are determined using information from the VELO together with the rest of the tracking system, composed of a large area silicon tracker located before a warm dipole magnet with a bending power of ∼ 4 Tm, and a combination of silicon strip detectors and straw drift chambers located after the magnet. Two Ring Imaging Cherenkov (RICH) detectors are used for charged hadron identification in the momentum range 2–100 GeV/c. Photon, electron and hadron candidates are identified by electromagnetic and hadronic calorimeters. A muon system of alternating layers of iron 2Inclusion of charged conjugated processes is implied throughout. 1 and drift chambers provides muon identification. The two calorimeters and the muon system provide the energy and momentum information to implement a first level (L0) hardware trigger. An additional trigger level (HLT) is software based, and its algorithms are tuned to the experimental operating condition. Events with a muon final states are triggered using two L0 trigger decisions: the single-muon decision, which requires one muon candidate with a transverse momentum p larger than 1.5 GeV/c, and the di-muon decision, which requires two muon candidates T √ with transverse momenta p and p satisfying the relation p ·p > 1.3 GeV/c. T,1 T,2 T,1 T,2 The single muon trigger decision in the second trigger level (HLT) includes a cut on the impact parameter (IP) with respect to the primary vertex, which allows for a lower p requirement (p > 1.0 GeV/c, IP > 0.1 mm). The di-muon trigger decision requires T T muon pairs of opposite charge with p > 500 MeV/c, forming a common vertex and T with an invariant mass m > 4.7 GeV/c2. A second trigger decision, primarily to select µµ J/ψ events, requires 2.97 < m < 3.21 GeV/c2. The remaining region of the di-muon µµ invariant mass range is also covered by trigger decisions that in addition require the di-muon secondary vertex to be well separated from the primary vertex. Events with purely hadronic final states are triggered by the L0 trigger if there is a calorimeter cluster with transverse energy E > 3.6 GeV. Other HLT trigger decisions T select generic displaced vertices, providing high efficiency for purely hadronic decays. 3 Analysis strategy ¯ Assuming the branching fractions predicted by the SM, and using the bb cross-section measured by LHCb in the pseudorapidity interval 2 < η < 6 and integrated over all transverse momenta of σ = 75 ± 14µb [8], approximately 3.9 B0 → µ+µ− and 0.4 bb s B0 → µ+µ− eventsareexpectedtobetriggered,reconstructedandselectedintheanalysed sample embedded in a large background. The general structure of the analysis is based upon the one described in Ref. [6]. First a very efficient selection removes the biggest amount of background while keeping most of the signal within the LHCb acceptance. The number of observed events is compared to thenumberofexpectedsignalandbackgroundeventsinbinsoftwoindependentvariables, the invariant mass and the output of a multi-variate discriminant. The discriminant is a Boosted Decision Tree (BDT) constructed using the TMVA package [9]. It supersedes the Geometrical Likelihood (GL) used in the previous analysis [6] as it has been found more performant in discriminating between signal and background events in simulated samples. No data were used in the choice of the multivariate discriminant in order not to bias the result. The combination of variables entering the BDT discriminant is optimized using sim- ulated events. The probability for a signal or background event to have a given value of the BDT output is obtained from data using B0 → h+h(cid:48)− candidates (where h((cid:48)) can be (s) a pion or a kaon) as signal and sideband B0 → µ+µ− candidates as background. (s) The invariant mass line shape of the signals is described by a Crystal Ball function [10] 2 withparametersextractedfromdatacontrolsamples. Thecentralvaluesofthemassesare obtained from B0 → K+π− and B0 → K+K− samples. The B0 and B0 mass resolutions s s are estimated by interpolating those obtained with di-muon resonances (J/ψ,ψ(2S) and Υ(1S,2S,3S)) and cross-checked with a fit to the invariant mass distributions of both inclusive B0 → h+h(cid:48)− decays and exclusive B0 → K+π− decays. The central values of (s) the masses and the mass resolution are used to define the signal regions. The number of expected signal events, for a given branching fraction hypothesis, is obtained by normalizing to channels of known branching fractions: B+→ J/ψK+, B0→ s J/ψφ and B0 → K+π−. These channels are selected in a way as similar as possible to the signals in order to minimize the systematic uncertainty related to the different phase space accessible to each final state. TheBDToutputandinvariantmassdistributionsforcombinatorialbackgroundevents in the signal regions are obtained using fits of the mass distribution of events in the mass sidebands in bins of the BDT output. The two-dimensional space formed by the invariant mass and the BDT output is binned. For each bin we count the number of candidates observed in the data, and compute the expected number of signal events and the expected number of background events. The binning is unchanged with respect to the 2010 analysis [6]. The compatibility of the observed distribution of events in all bins with the distribution expected for a given branching fraction hypothesis is computed using the CL method [11], which allows a s given hypothesis to be excluded at a given confidence level. 4 Selection The B0 → µ+µ− selections require two muon candidates of opposite charge. Tracks are (s) required to be of good quality and to be displaced with respect to any primary vertex. The secondary vertex is required to be well fitted (χ2/nDoF < 9) and must be sepa- rated from the primary vertex in the forward direction by a distance of flight significance (L/σ(L)) greater than 15. When more than one primary vertex is reconstructed, the one that gives the minimum impact parameter significance for the candidate is chosen. The reconstructed candidate has to point to this primary vertex (IP/σ(IP) < 5). Improvements have been made to the selection developed for 2010 data [6]. The RICH is used to identify kaons in the B0 → J/ψφ normalization channel and the Kullback- s Leibler (KL) distance [12] is used to suppress duplicated tracks created by the recon- struction. This procedure compares the parameters and correlation matrices of the recon- structed tracks and where two are found to be similar, in this case with a symmetrized KL divergence less than 5000, only the one with the higher track fit quality is considered. The inclusive B0 → h+h(cid:48)− sample is the main control sample for the determina- (s) tion from data of the probability distribution function (PDF) of the BDT output. This sample is selected in exactly the same way as the B0 → µ+µ− signals apart from the (s) muon identification requirement. The same selection is also applied to the B0 → K+π− normalization channel. 3 The muon identification efficiency is uniform within ∼ 1% in the considered phase space therefore no correction is added to the BDT PDF extracted from the B0 → h+h(cid:48)− (s) sample. The remaining phase space dependence of the muon identification efficiency is instead taken into account in the computation of the normalization factor when the B0 → K+π− channel is considered. The J/ψ → µµ decay in the B+ → J/ψK+ and B0 → J/ψφ normalization channels s is selected in a very similar way to the B0 → µ+µ− channels, apart from the pointing (s) requirement. K± candidates are required to be identified by the RICH detector and to pass track quality and impact parameter cuts. To avoid pathological events, all tracks from selected candidates are required to have a momentum less than 1TeV/c. Only B candidates with decay times less than 5τ , B0 (s) where τ is the B lifetime [13], are accepted for further analysis. Di-muon candidates B0 (s) coming from elastic di-photon production are removed by requiring a minimum transverse momentum of the B candidate of 500MeV/c. 5 Determination of the mass and BDT distributions The variables entering the BDT discriminant are the six variables used as input to the GL in the 2010 analysis plus three new variables. The six variables used in the 2010 anal- ysis are the B lifetime, impact parameter, transverse momentum, the minimum impact parameter significance (IP/σ(IP)) of the muons, the distance of closest approach between the two muons and the isolation of the two muons with respect to any other track in the event. The three new variables are: 1. the minimum p of the two muons; T 2. the cosine of the angle between the muon momentum in the B rest frame and the vector perpendicular to the B momentum and the beam axis: p p −p p y,µ1 x,B x,µ1 y,B cosP = (1) p (m /2) T,B µµ where µ labels one of the muons and m is the reconstructed B candidate mass3; 1 µµ 3. the B isolation [14] p (B) T I = , (2) B (cid:80) p (B)+ p T i T,i where p (B) is the B transverse momentum with respect to the beam line T and the sum is over all the tracks, excluding the muon candidates, that satisfy 3As the B is a (pseudo)-scalar particle, this variable is uniformely distributed for signal candidates whileispeakedatzeroforb¯b→µ+µ−X backgroundcandidates. Infact,muonsfromsemi-leptonicdecays are mostly emitted in the direction of the b’s and, therefore, lie in a plane formed by the B momentum and the beam axis. 4