PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of springeropenLink to Publisher's site
The European Physical Journal. C, Particles and Fields
 
Eur Phys J C Part Fields. 2015; 75: 595.
Published online 2015 December 15. doi:  10.1140/epjc/s10052-015-3809-7
PMCID: PMC4684630

Search for long-lived heavy charged particles using a ring imaging Cherenkov technique at LHCb

R. Aaij, B. Adeva, M. Adinolfi, A. Affolder, Z. Ajaltouni, S. Akar, J. Albrecht, F. Alessio, M. Alexander, S. Ali, G. Alkhazov, P. Alvarez Cartelle, A. A. Alves Jr, S. Amato, S. Amerio, Y. Amhis, L. An, L. Anderlini, J. Anderson, M. Andreotti, J. E. Andrews, R. B. Appleby, O. Aquines Gutierrez, F. Archilli, P. d’Argent, A. Artamonov, M. Artuso, E. Aslanides, G. Auriemma, M. Baalouch, S. Bachmann, J. J. Back, A. Badalov, C. Baesso, W. Baldini, R. J. Barlow, C. Barschel, S. Barsuk, W. Barter, V. Batozskaya, V. Battista, A. Bay,corresponding author L. Beaucourt, J. Beddow, F. Bedeschi, I. Bediaga, L. J. Bel, I. Belyaev, E. Ben-Haim, G. Bencivenni, S. Benson, J. Benton, A. Berezhnoy, R. Bernet, A. Bertolin, M.-O. Bettler, M. van Beuzekom, A. Bien, S. Bifani, T. Bird, A. Birnkraut, A. Bizzeti, T. Blake, F. Blanc, J. Blouw, S. Blusk, V. Bocci, A. Bondar, N. Bondar, W. Bonivento, S. Borghi, A. Borgia, M. Borsato, T. J. V. Bowcock, E. Bowen, C. Bozzi, D. Brett, M. Britsch, T. Britton, J. Brodzicka, N. H. Brook, A. Bursche, J. Buytaert, S. Cadeddu, R. Calabrese, M. Calvi, M. Calvo Gomez, P. Campana, D. Campora Perez, L. Capriotti, A. Carbone, G. Carboni, R. Cardinale, A. Cardini, P. Carniti, L. Carson, K. Carvalho Akiba, R. Casanova Mohr, G. Casse, L. Cassina, L. Castillo Garcia, M. Cattaneo, Ch. Cauet, G. Cavallero, R. Cenci, M. Charles, Ph. Charpentier, M. Chefdeville, S. Chen, S. F. Cheung, N. Chiapolini, M. Chrzaszcz, X. Cid Vidal, G. Ciezarek, P. E. L. Clarke, M. Clemencic, H. V. Cliff, J. Closier, V. Coco, J. Cogan, E. Cogneras, V. Cogoni, L. Cojocariu, G. Collazuol, P. Collins, A. Comerma-Montells, A. Contu, A. Cook, M. Coombes, S. Coquereau, G. Corti, M. Corvo, I. Counts, B. Couturier, G. A. Cowan, D. C. Craik, A. Crocombe, M. Cruz Torres, S. Cunliffe, R. Currie, C. D’Ambrosio, J. Dalseno, P. N. Y. David, A. Davis, K. De Bruyn, S. De Capua, M. De Cian, J. M. De Miranda, L. De Paula, W. De Silva, P. De Simone, C. T. Dean, D. Decamp, M. Deckenhoff, L. Del Buono, N. Déléage, D. Derkach, O. Deschamps, F. Dettori, B. Dey, A. Di Canto, F. Di Ruscio, H. Dijkstra, S. Donleavy, F. Dordei, M. Dorigo, A. Dosil Suárez, D. Dossett, A. Dovbnya, K. Dreimanis, G. Dujany, F. Dupertuis, P. Durante, R. Dzhelyadin, A. Dziurda, A. Dzyuba, S. Easo, U. Egede, V. Egorychev, S. Eidelman, S. Eisenhardt, U. Eitschberger, R. Ekelhof, L. Eklund, I. El Rifai, Ch. Elsasser, S. Ely, S. Esen, H. M. Evans, T. Evans, A. Falabella, C. Färber, C. Farinelli, N. Farley, S. Farry, R. Fay, D. Ferguson, V. Fernandez Albor, F. Ferrari, F. Ferreira Rodrigues, M. Ferro-Luzzi, S. Filippov, M. Fiore, M. Fiorini, M. Firlej, C. Fitzpatrick, T. Fiutowski, P. Fol, M. Fontana, F. Fontanelli, R. Forty, O. Francisco, M. Frank, C. Frei, M. Frosini, J. Fu, E. Furfaro, A. Gallas Torreira, D. Galli, S. Gallorini, S. Gambetta, M. Gandelman, P. Gandini, Y. Gao, J. García Pardiñas, J. Garofoli, J. Garra Tico, L. Garrido, D. Gascon, C. Gaspar, R. Gauld, L. Gavardi, G. Gazzoni, A. Geraci, D. Gerick, E. Gersabeck, M. Gersabeck, T. Gershon, Ph. Ghez, A. Gianelle, S. Gianì, V. Gibson, L. Giubega, V. V. Gligorov, C. Göbel, D. Golubkov, A. Golutvin, A. Gomes, C. Gotti, M. Grabalosa Gándara, R. Graciani Diaz, L. A. Granado Cardoso, E. Graugés, E. Graverini, G. Graziani, A. Grecu, E. Greening, S. Gregson, P. Griffith, L. Grillo, O. Grünberg, B. Gui, E. Gushchin, Yu. Guz, T. Gys, C. Hadjivasiliou, G. Haefeli, C. Haen, S. C. Haines, S. Hall, B. Hamilton, T. Hampson, X. Han, S. Hansmann-Menzemer, N. Harnew, S. T. Harnew, J. Harrison, J. He, T. Head, V. Heijne, K. Hennessy, P. Henrard, L. Henry, J. A. Hernando Morata, E. van Herwijnen, M. Heß, A. Hicheur, D. Hill, M. Hoballah, C. Hombach, W. Hulsbergen, T. Humair, N. Hussain, D. Hutchcroft, D. Hynds, M. Idzik, P. Ilten, R. Jacobsson, A. Jaeger, J. Jalocha, E. Jans, A. Jawahery, F. Jing, M. John, D. Johnson, C. R. Jones, C. Joram, B. Jost, N. Jurik, S. Kandybei, W. Kanso, M. Karacson, T. M. Karbach, S. Karodia, M. Kelsey, I. R. Kenyon, M. Kenzie, T. Ketel, B. Khanji, C. Khurewathanakul, S. Klaver, K. Klimaszewski, O. Kochebina, M. Kolpin, I. Komarov, R. F. Koopman, P. Koppenburg, L. Kravchuk, K. Kreplin, M. Kreps, G. Krocker, P. Krokovny, F. Kruse, W. Kucewicz, M. Kucharczyk, V. Kudryavtsev, K. Kurek, T. Kvaratskheliya, V. N. La Thi, D. Lacarrere, G. Lafferty, A. Lai, D. Lambert, R. W. Lambert, G. Lanfranchi, C. Langenbruch, B. Langhans, T. Latham, C. Lazzeroni, R. Le Gac, J. van Leerdam, J. P. Lees, R. Lefèvre, A. Leflat, J. Lefrançois, O. Leroy, T. Lesiak, B. Leverington, Y. Li, T. Likhomanenko, M. Liles, R. Lindner, C. Linn, F. Lionetto, B. Liu, S. Lohn, I. Longstaff, J. H. Lopes, D. Lucchesi, H. Luo, A. Lupato, E. Luppi, O. Lupton, F. Machefert, I. V. Machikhiliyan, F. Maciuc, O. Maev, S. Malde, A. Malinin, G. Manca, G. Mancinelli, P. Manning, A. Mapelli, J. Maratas, J. F. Marchand, U. Marconi, C. Marin Benito, P. Marino, R. Märki, J. Marks, G. Martellotti, M. Martinelli, D. Martinez Santos, F. Martinez Vidal, D. Martins Tostes, A. Massafferri, R. Matev, Z. Mathe, C. Matteuzzi, A. Mauri, B. Maurin, A. Mazurov, M. McCann, J. McCarthy, A. McNab, R. McNulty, B. McSkelly, B. Meadows, F. Meier, M. Meissner, M. Merk, D. A. Milanes, M. N. Minard, D. S. Mitzel, J. Molina Rodriguez, S. Monteil, M. Morandin, P. Morawski, A. Mordà, M. J. Morello, J. Moron, A. B. Morris, R. Mountain, F. Muheim, J. Müller, K. Müller, V. Müller, M. Mussini, B. Muster, P. Naik, T. Nakada, R. Nandakumar, I. Nasteva, M. Needham, N. Neri, S. Neubert, N. Neufeld, M. Neuner, A. D. Nguyen, T. D. Nguyen, C. Nguyen-Mau, V. Niess, R. Niet, N. Nikitin, T. Nikodem, A. Novoselov, D. P. O’Hanlon, A. Oblakowska-Mucha, V. Obraztsov, S. Ogilvy, O. Okhrimenko, R. Oldeman, C. J. G. Onderwater, B. Osorio Rodrigues, J. M. Otalora Goicochea, A. Otto, P. Owen, A. Oyanguren, A. Palano, F. Palombo, M. Palutan, J. Panman, A. Papanestis, M. Pappagallo, L. L. Pappalardo, C. Parkes, G. Passaleva, G. D. Patel, M. Patel, C. Patrignani, A. Pearce, A. Pellegrino, G. Penso, M. Pepe Altarelli, S. Perazzini, P. Perret, L. Pescatore, K. Petridis, A. Petrolini, E. Picatoste Olloqui, B. Pietrzyk, T. Pilař, D. Pinci, A. Pistone, S. Playfer, M. Plo Casasus, T. Poikela, F. Polci, A. Poluektov, I. Polyakov, E. Polycarpo, A. Popov, D. Popov, B. Popovici, C. Potterat, E. Price, J. D. Price, J. Prisciandaro, A. Pritchard, C. Prouve, V. Pugatch, A. Puig Navarro, G. Punzi, W. Qian, R. Quagliani, B. Rachwal, J. H. Rademacker, B. Rakotomiaramanana, M. Rama, M. S. Rangel, I. Raniuk, N. Rauschmayr, G. Raven, F. Redi, S. Reichert, M. M. Reid, A. C. dos Reis, S. Ricciardi, S. Richards, M. Rihl, K. Rinnert, V. Rives Molina, P. Robbe, A. B. Rodrigues, E. Rodrigues, P. Rodriguez Perez, S. Roiser, V. Romanovsky, A. Romero Vidal, M. Rotondo, J. Rouvinet, T. Ruf, H. Ruiz, P. Ruiz Valls, J. J. Saborido Silva, N. Sagidova, P. Sail, B. Saitta, V. Salustino Guimaraes, C. Sanchez Mayordomo, B. Sanmartin Sedes, R. Santacesaria, C. Santamarina Rios, E. Santovetti, A. Sarti, C. Satriano, A. Satta, D. M. Saunders, D. Savrina, M. Schiller, H. Schindler, M. Schlupp, M. Schmelling, T. Schmelzer, B. Schmidt, O. Schneider, A. Schopper, M. H. Schune, R. Schwemmer, B. Sciascia, A. Sciubba, A. Semennikov, I. Sepp, N. Serra, J. Serrano, L. Sestini, P. Seyfert, M. Shapkin, I. Shapoval, Y. Shcheglov, T. Shears, L. Shekhtman, V. Shevchenko, A. Shires, R. Silva Coutinho, G. Simi, M. Sirendi, N. Skidmore, I. Skillicorn, T. Skwarnicki, 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, B. Spaan, P. Spradlin, S. Sridharan, F. Stagni, M. Stahl, S. Stahl, O. Steinkamp, O. Stenyakin, F. Sterpka, S. Stevenson, S. Stoica, S. Stone, B. Storaci, S. Stracka, M. Straticiuc, U. Straumann, R. Stroili, L. Sun, W. Sutcliffe, K. Swientek, S. Swientek, V. Syropoulos, M. Szczekowski, P. Szczypka, T. Szumlak, S. T’Jampens, T. Tekampe, M. Teklishyn, G. Tellarini, F. Teubert, C. Thomas, E. Thomas, J. van Tilburg, V. Tisserand, M. Tobin, J. Todd, S. Tolk, L. Tomassetti, D. Tonelli, S. Topp-Joergensen, N. Torr, E. Tournefier, S. Tourneur, K. Trabelsi, M. T. Tran, M. Tresch, A. Trisovic, A. Tsaregorodtsev, P. Tsopelas, N. Tuning, M. Ubeda Garcia, A. Ukleja, A. Ustyuzhanin, U. Uwer, C. Vacca, V. Vagnoni, G. Valenti, A. Vallier, R. Vazquez Gomez, P. Vazquez Regueiro, C. Vázquez Sierra, S. Vecchi, J. J. Velthuis, M. Veltri, G. Veneziano, M. Vesterinen, B. Viaud, D. Vieira, M. Vieites Diaz, X. Vilasis-Cardona, A. Vollhardt, D. Volyanskyy, D. Voong, A. Vorobyev, V. Vorobyev, C. Voß, J. A. de Vries, R. Waldi, C. Wallace, R. Wallace, J. Walsh, S. Wandernoth, J. Wang, D. R. Ward, N. K. Watson, D. Websdale, A. Weiden, M. Whitehead, D. Wiedner, G. Wilkinson, M. Wilkinson, M. Williams, M. P. Williams, M. Williams, F. F. Wilson, J. Wimberley, J. Wishahi, W. Wislicki, M. Witek, G. Wormser, S. A. Wotton, S. Wright, K. Wyllie, Y. Xie, Z. Xu, Z. Yang, X. Yuan, O. Yushchenko, M. Zangoli, M. Zavertyaev, L. Zhang, Y. Zhang, A. Zhelezov, A. Zhokhov, and L. Zhong

Abstract

A search is performed for heavy long-lived charged particles using 3.0 fb-1 of proton–proton collisions collected at s 7 and 8  TeV with the LHCb detector. The search is mainly based on the response of the ring imaging Cherenkov detectors to distinguish the heavy, slow-moving particles from muons. No evidence is found for the production of such long-lived states. The results are expressed as limits on the Drell–Yan production of pairs of long-lived particles, with both particles in the LHCb pseudorapidity acceptance, 1.8 < η < 4.9. The mass-dependent cross-section upper limits are in the range 2–4 fb (at 95 % CL) for masses between 14 and 309  GeV /c2.

Introduction

Several extensions of the Standard Model (SM) propose the existence of charged massive stable particles (CMSP). Stable particles, in this context, are long-lived particles that can travel through a detector without decaying. These particles can have long lifetimes for a variety of reasons, e.g. a new (approximately) conserved quantum number, a weak coupling or a limited phase space in any allowed decay. In supersymmetric (SUSY) models, CMSPs can be sleptons (~), charginos, or R-hadrons. R-hadrons are colourless states combining squarks (q~) or gluinos (g~) and SM quarks or gluons. In the gauge-mediated supersymmetry breaking (GMSB) model [13] the breakdown of SUSY is mediated by gauge interactions and can occur at a relatively low energy scale. For a particular range of parameter space in the minimal model (mGMSB) the next-to-lightest supersymmetric particle can be a long-lived stau (τ1~), with a mass of the order of 100 GeV /c2 or higher. The τ1~ is the lightest mass eigenstate, resulting from the mixture of right-handed and left-handed superpartners of the τ, dominated by the right-handed component.

A CMSP loses energy mainly via ionisation; strongly interacting CMSPs are not considered here. In a detector such as LHCb a CMSP with a kinetic energy above about 5  GeV should be able to traverse the muon chambers. Those particles would often be produced with a relatively low velocity and could be identified by their time-of-flight, and by their specific energy loss, dE/dx, in the detectors; Cherenkov radiation would be absent in Cherenkov counters tuned for ultra-relativistic particles.

Several experiments have searched for CMSPs [412]. With the exception of DELPHI [5], which had Cherenkov counters, the analyses are based on dE/dx and time-of-flight measurements. The primary interest here is to show the potential of the identification technique based on ring imaging Cherenkov (RICH) detectors, in addition to the exploration of the forward pseudorapidity region only partially covered by the central detectors at the Tevatron and the LHC.

The analysis described in this study is mainly based on the absence of Cherenkov radiation in the RICH detectors. This technique is used to search for pairs of CMSPs in LHCb, produced by a Drell–Yan mechanism.

The LHCb detector and the detection of slow particles

The LHCb detector [13, 14] is a single-arm forward spectrometer covering the approximate pseudorapidity range 1.8 < η < 4.9, 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 (the vertex locator, VELO) surrounding the proton–proton interaction region [15], a large-area silicon-strip detector located upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip detectors and straw drift tubes [16] placed downstream of the magnet. The tracking system provides a measurement of momentum, p , of charged particles with a relative uncertainty that varies from 0.5 % at low momentum to 1.0 % at 200  GeV /c. The minimum distance of a track to a primary vertex (PV), the impact parameter (IP), is measured with a resolution of (15 + 29/pT ) μm, where pT is the component of the momentum transverse to the beam, in  GeV /c. Photons, electrons and hadrons are identified by a calorimeter system consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers [17].

Different types of charged particles are distinguished using information from two RICH detectors [18]. The RICH system, which plays a crucial role in this analysis, consists of an upstream detector with silica aerogel and C4F10 gas radiators, positioned directly after the VELO, and a downstream detector with a CF4 gas radiator, located just after the tracking system.

The online event selection is performed by a trigger [19], which consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, which applies a full event reconstruction.

The analysis presented here is based on two data sets collected in 2011 and 2012 corresponding to integrated luminosities of 1.0  and 2.0 badhbox-1 from proton–proton collisions recorded at centre-of-mass energies of 7 and 8  TeV, respectively.

In the production process considered, CMSPs can have velocities β ≡ v/c as low as 0.7, and their arrival time at the subdetectors can differ by several nanoseconds with respect to lighter particles with β ≃ 1. For illustration, the β spectrum is shown in Fig. 1, for two values of the CMSP mass, at centre-of-mass energy of 7  TeV. The effects of such delayed detection on the efficiencies of the subdetectors are determined from simulation in which the timing information is modelled according to dedicated electronic measurements and tests in beam. The muon chambers have the largest inefficiency for slow-particle reconstruction. The maximal delay for a particle to be accepted by the front-end electronics is 12.5 ns [17]. In the most downstream muon chamber, this delay corresponds to the arrival of a particle with β = 0.83. To be identified as a muon, the charged particle must be associated with hits in the last four muon chambers, a requirement that particles with β < 0.8 fail to meet. The large time-of-flight can also bias the reconstructed position of the particle passing through the tracker straw tubes, which accept a maximal drift-time of about 35 ns for tracks passing close to the straw radius of 2.5 mm [16]. These combined effects result in a vanishingly small reconstruction efficiency for particles with β < 0.8 but an efficiency above 95 % if β > 0.95, as shown in Fig. 1.

Fig. 1
CMSP velocity spectrum for the CMSP masses of 124 and 309  GeV /c2. The proton–proton centre-of-mass energy is 7 TeV. The dots with error bars show the efficiency to detect tracks as a function of the β of the particle (right scale ...

Simulation

CMSP signal

The adopted framework is stau pair production, τ1~+τ1~-, in mGMSB via a Drell–Yan process. Pairs of CMSPs originating from cascade decays of heavier particles are explicitly not considered. In the following the symbol τ1~ is used when the context is explicitly the mGMSB model, while CMSP is kept for the more general context.

The mGMSB model has six parameters [2, 3]: the SUSY breaking scale (Λ), the mass scale of the SUSY loop messengers (Mm), the number of messenger supermultiplets (N5), the ratio of the vacuum expectation values of the two neutral Higgs fields (tanβ), the sign of the Higgs mass parameter (μ), and the parameter Cgrav, which affects the gravitino mass. The Spheno3.0 SUSY spectrum generator [20] is used to compute the masses of the τ1~ as a function of the above six parameters. The SPS7 benchmark scenario [21] is used to determine the parameter space, where N5 = 3, tanβ = 15, μ > 0, Mm = 2Λ, and the parameter Cgrav = 4000 are fixed. Variation of Λ then uniquely determines the τ1~ mass and lifetime, which is of the order of 100 ns. In this study the τ1~ is considered stable.

The predictions for τ1~ pair production are based on next-to-leading order (NLO) cross-section calculations by the Prospino2.1 program [22] using the CTEQ6.6M parton distribution function (PDF) set [23]. These predictions at s 7 and 8  TeV are presented in Table 1. The relative theoretical uncertainties vary between 5 and 8 %, and are determined following Ref. [24].

Table 1
Values of the mGMSB Λ parameters in the SPS7 scenario used in this study, the corresponding masses of the τ1~, mτ~, and the cross-section of the pair production at next-to-leading order. The last two columns give the detector acceptance ...

Fully simulated signal samples, with masses varying from 124 to 309 GeV/c2, have been produced for proton–proton collisions at s 7 and 8  TeV. The τ1~ pairs generated by Pythia 6.423 [25], with both τ1~ particles in the fiducial range 1.8 < η < 4.9 are passed to Geant4  [26, 27] for detector simulation. The fraction of τ1~ pairs within the fiducial range is defined as the acceptance, A. The acceptance factor obtained from Pythia with the MSTW2008 PDF set [28] is also shown in Table 1, with model uncertainties ranging from 5 to 9 % for τ1~ mass from 124 to 309  GeV /c2, mainly associated to the choice of PDF.

For larger τ1~ masses, the Drell–Yan process results in a lower forward boost of the τ1~ pair, with a subsequent increase in the pair opening angle in the detector frame. The decrease of A for an increasing τ1~ mass is due to a higher probability for one of the particles to escape the LHCb geometrical acceptance.

Background

The main background is from the Drell–Yan production of muon pairs, Z/γ → μ+μ-. Samples of Z/γ → μ+μ- events have been produced with Pythia and fully simulated with Geant4. The cross-section for this process has been calculated with DYNNLO [29] at next-to-next-to-leading order with the MSTW2008 PDF set. The preselection requirements (see Sect. 4.1) lead to values of the predicted cross-section in LHCb for s 7 and 8  TeV of 1.08 ± 0.10 and 1.36 ± 0.12 pb, respectively. These values are nearly two orders of magnitude larger than the predicted τ1~ pair cross-section in the most favourable case, corresponding to Λ = 40 Te V.

Other background sources include muons produced by top quark pairs, and from τ pairs. To study the background contributions from these processes, samples of Z/γ → τ+τ- and top quark pair decays have been simulated.

Data selection

The event selection is performed in two steps: a preselection aimed at suppressing the most prominent backgrounds, followed by a multivariate analysis, based on an artificial neural network that is trained using calibrated simulation.

Preselection

CMSP candidates are identified as high-momentum charged particles with hits in the VELO, all the tracking stations and the four last muon detectors.

Events are selected that contain two or more such particles where one of the particles passes the high-pT single muon trigger with a threshold of 15  GeV /c. The trigger efficiency is estimated from simulation to be 92 % for a mass of 124  GeV /c2, and 89 % for 309  GeV /c2. The two candidates must have opposite charge and each have pT > 50  GeV /c. To reject background from Z/γ → μ+μ- the pair must have a dimuon mass larger than 100  GeV /c2. A mass-dependent lower threshold on momentum is applied to select particles with β > 0.8.

Several criteria are used to reject muons, electrons and hadrons. Pions and kaons in jets may be identified as muons if they decay in flight or if shower fragments escape from the calorimeters to the muon stations. As hadrons and electrons deposit more energy in the calorimeters than that expected for CMSPs, an efficient rejection of these backgrounds is achieved by requiring the sum of the ECAL and HCAL energies associated with the extrapolation of the charged particle to the calorimeters to be less than 1 % of the momentum of that particle. The background from misidentified muons contributes approximately equally to same- and opposite-charge pairs [30]. No same-charge event is found in the preselected data, showing that this contribution is negligible.

CMSPs, as well as muons from Z/γ decays, would be produced at the PV and should have a smaller IP with respect to the PV than muons from heavy quark or tau decays. Requiring an IP of less than 50 μm selects efficiently CMSP candidates. After preselection, the contribution from the Z/γ → τ+τ- process where both taus decay leptonically to muons is estimated from simulation to contribute less than 0.1 events in total. Pairs of muons produced from top quark decays into b quarks and W± bosons, with the W± bosons decaying leptonically into muons, contribute less than one event, as determined from simulation.

In summary, after preselection the only significant source of background is from Z/γ → μ+μ- decays. The predicted number of dimuon events in the 7  Te V(8  TeV) data set is 249 ± 49 (570 ± 110) which is in good agreement with the 239 (713) observed candidate events. The uncertainties comprise contributions from the preselection cuts (Sect. 6), and the uncertainty on the Z/γ → μ+μ- cross-section (Sect. 3.2). The expected number of events with τ1~ pairs is 2.7 events in the full data set of 3.0 badhbox-1, according to the cross-section calculated with Prospino2.1, with SPS7 parameters and a τ1~ mass of 124  GeV /c2.

Selection

An artificial neural network (ANN) is used to distinguish CMSPs from muons by exploiting the difference in interactions that these particles have in matter. To reduce model dependence, the ANN is applied to the individual CMSP candidates, rather than to CMSP-pairs, and a minimum requirement is placed on the product of the two ANN responses. Four variables of the CMSP candidates are used as ANN inputs, computed from the energy deposited in the VELO sensors (ΔE VELO), in the ECAL (ΔE ECAL), in the HCAL (ΔE HCAL), and a likelihood variable associated with the RICH information (DLLx). Model dependence is reduced as much as possible by the absence of kinematical observables in the ANN. The energy loss of a charged particle traversing a VELO sensor follows a Landau distribution. The most probable energy deposition in a sensor is estimated using a truncated mean where only the 60 % lowest depositions are averaged.

Particle identification for a CMSP candidate, using RICH information, is provided by the DLLx variable. A particle identification hypothesis is assigned to a track using a likelihood method. The information from the three radiators is combined and a “delta log-likelihood” (DLL) value computed. The DLL gives, for each track, the change in the overall event log-likelihood when the particle ID hypothesis is changed from π to μ, e, K, p. The DLLx classification has been added to account for high momentum particles which do not radiate, or have a Cherenkov angle which is too small to fit one of the five particle hypotheses. A positive DLLx indicates a high probability that the candidate has a relatively low velocity. More details are given in Sect. 5.

Events with both candidate CMSPs with DLLx >  - 5 are used in the analysis, with no loss of signal, as deduced from simulation. The numbers of selected events with CMSP-pairs are given in Table 2.

Table 2
Number of events with both CMSPs candidates satisfying DLLx >  - 5, and the final efficiency, ϵ, after the multivariate analysis selection, given for each mass hypothesis

Simulated events are used to train the ANN. The first three variables defined above are calibrated using muons from Z/γ → μ+μ- to ensure that simulation agrees with data. A total of 25k Z events from the 2011 data set and 65k from the 2012 data set are used for the calibration. In the Z mass region the expected amount of signal is smaller than one event and cannot bias the procedure. The DLLx variable is by far the most discriminating, and its calibration procedure is presented in detail in Sect. 5.

The ANN training is carried out independently for the 7 and 8  Te Vdata sets, and for all the CMSP mass hypotheses considered. Figure 2 shows the distribution of the four ANN input variables for the 8  Te Vdata set, compared to the background and signal predictions; good agreement can be seen between data and simulated background.

Fig. 2
Number of CMSP candidates, as a function of the four variables used as inputs to the ANN. There are two CMSP candidates per event. The black dots with error bars show the 2012 data. The dashed red histogram is the expected shape for 124  GeV / ...

The discriminating variable is the product of the ANN outputs obtained for the two CMSP candidates. This “pair significance” is shown in Fig. 3 for the 124 and 309  GeV /c2 CMSP mass hypotheses for both 7 and 8  Te Vdata sets. The requirement placed on the pair significance is determined by the value needed to achieve a 95 % signal efficiency. After the pair significance selection, the signal efficiency for candidate events in the LHCb acceptance is 50 % for CMSPs with a mass of 124  GeV /c2, decreasing for increasing mass to a minimum of 27 % at 309  GeV /c2. The signal efficiency values for CMSPs in the acceptance, after the ANN selection, are given in Table 2. After the full selection is applied, the dimuon background is suppressed by a factor of 10-5.

Fig. 3
Number of CMSP pairs, N, as a function of the pair significance. The left and right figures correspond to the CMSP masses of 124 and 309  GeV /c2, respectively. The black points with their statistical uncertainty show the 7 TeV (top) and ...

CMSP identification with Cherenkov detectors

The present study uses the Cherenkov radiation produced in the RICH detectors to identify CMSPs. The Cherenkov momentum thresholds for muons, protons, and CMSPs with masses of 124 and 309  GeV /c2, are given in Table 3 for the three radiators in the LHCb detectors. Only CMSP candidates with momenta above 200  GeV /c are considered. For this momentum range, particles with masses of the order of MeV /c2 to GeV /c2, have Cherenkov angles very close to the saturation value arccos(1/(nβ)), where n is the refractive index of the medium. The fraction of CMSPs with momentum above 2  TeV /c is negligible, and the CMSPs are therefore expected not to produce Cherenkov radiation in the gaseous radiators. Around half of the 124  GeV /c2 CMSPs have a momentum above the Cherenkov threshold for aerogel, and Cherenkov angles smaller than the saturation value. This allows them to be separated from the background. Only a few percent are expected to be in the momentum range 1.4–2.0  TeV /c, corresponding to Cherenkov angles from 0.225 to 0.234  rad. It is possible to distinguish these angles from the saturation value of 0.242  rad in the aerogel as the angular resolution is about 5.6  mrad.

Table 3
Refractive indices and Cherenkov β thresholds for the three radiators. The momentum threshold is given for muons, protons, and 124 and 309  GeV /c2 CMSPs

As previously said, the variable DLLx has been introduced to identify high momentum particles which do not radiate, or have a Cherenkov angle which is too small to fit one of the five particle hypotheses, π, μ, e, K, p. The DLLx value is positive for the momentum distributions of the CMSPs, for all of the masses considered.1

Simulated dimuon background events and CMSP signal samples used to train the ANN are first validated with data.

The study of the background samples is performed on a set of muons above the Cherenkov threshold and selected from Z decays. Such events have an event topology and kinematics that are very close to those of the dimuon background expected in the CMSP analysis. The DLLx distribution is shown in Fig. 4a for muons from data and simulated Z decays. For illustration, the expected signal shapes for CMSPs with masses of 124 and 309  GeV /c2 are superimposed. The small difference is due to a change in the underlying event and some light from the aerogel for the 124  GeV /c2 case. A clear separation between the signal and background muon DLLx distributions is observed. The difference in the data and simulated muon distributions is mainly due to the lack of precision in the mapping of the photon detection efficiency in the RICH system. In particular, the peak at DLLx >  - 5 is produced by the decrease of the photon detection efficiency when approaching boundaries in the RICH modules. The simulation only partially reproduces this behaviour and the number of candidates above DLLx =  - 5 is too low by around a factor of two. To compensate for this, 15 % of simulated muon events with DLLx falling close to zero have been shifted by anad hoc value to obtain the best agreement between data and simulation in the DLLx >  - 5 region, resulting in the distribution shown in Fig. 4b. It is expected that a correct efficiency map should produce such a shift, moving above zero the slightly negative DLLx values. This set of simulated background events is used to train the ANN. Note that only candidates with DLLx >  - 5 are used in the ANN. In order to assess the systematic uncertainty associated to this correction method, two other procedures are considered. In the first procedure, the ANN training is performed on the original background simulated data set. In the second, the DLLx values for each muon are randomly chosen following a set of templates inferred from the DLLx data distributions as function of p and η. Despite the fact that this operation is done in bins of p and η, it is obvious that most of the correlation is lost in the randomisation process. The three methods are found to provide the same final discrimination power and their contributions to the systematic uncertainties are small. This is due to the strong separation between signal and background that is provided intrinsically by the DLLx variable.

Fig. 4
a DLLx distributions from muons selected from Z decays in 7 TeV data (black points) and simulation (blue histogram). The expected shape from 124 and 309  GeV /c2 CMSPs are also shown. b The simulated DLLx values are shifted in such a ...

The validation of the signal sample is more complex due to the absence of a SM process that can be used for calibration. The quality of the simulation was studied using protons from Λ → pπ decays with a velocity below the Cherenkov threshold. The differences between the data and simulation for these protons are extrapolated to the CMSP kinematics using a fast simulation method, and the contribution to the systematic uncertainty estimated.

The proton is below the Cherenkov threshold in all of the RICH radiators for p < 3.8 GeV /c, and above the threshold for p ≳ 30 GeV /c. Pairs of opposite-charge tracks are selected from data and from simulated events passing a minimum bias trigger. The pair must combine to form a particle with a mass compatible with the known mass of the Λ baryon, and the reconstructed vertex must be more than 3  mm from the beam axis.

Samples of protons below and above Cherenkov threshold are obtained by choosing the momentum regions below 3.8  GeV /c and above 30  GeV /c, respectively. Figure 5 shows the corresponding DLLx distributions, indicating a reasonable agreement between the DLLx distributions from data collected at 7  TeV and simulation.

Fig. 5
The DLLx variable for: a protons selected above the Cherenkov threshold, p  > 30  GeV /c, in 7 TeV data (black line) and simulation (blue points), b and c below-threshold protons, p  < 3.8 GeV /c, from data and simulation, respectively. ...

A Gaussian function plus a polynomial, to account for the tail at low DLLx, is fitted to the DLLx distributions for below-Cherenkov-threshold particles. The mean and width of the Gaussian functions are 6.5 ± 0.3 and 3.4 ± 0.3 for data, and 6.0 ± 0.3 and 3.6 ± 0.4 for simulation. The DLLx value in data is 0.5 ± 0.4 units higher, which may indicate that there is a lower photon detection efficiency in data compared to simulation. A maximal deviation of ±1 DLLx units is considered in the following to assess the systematic effects. The extrapolation from the low momentum proton result to the CMSP regime is made using a fast simulation.

In addition to the Geant4-based full simulation, a fast simulation describing the main features of the RICH measurement process is also used. This allows the impact of varying parameters and the algorithms to be studied in a more efficient way. The fast simulation generates a target particle (a proton from Λ decays, a muon or a CMSP) with a momentum distribution representing the phenomenon under study. The underlying event is represented by a number of pions with a momentum distribution obtained from minimum bias events. The simulation of Cherenkov emission in the radiators is then performed for each particle. The number of Cherenkov photons generated by a particle of velocity β follows a Poisson distribution of average N0(β2n2 - 1)/(β2(n2 - 1)), where N0 is the maximal number of photons for a saturated ring and n is the refractive index. The ring has an average radius corresponding to the expected Cherenkov angle and a Gaussian profile of width σc representing the angular resolution of the detector. Finally, random noise is added using the probability for a pixel to fire, probnoise. The nominal values of the parameters used in the fast simulation are given in Table 4. The event log-likelihood for each target particle hypothesis is

LL=-pixeliall pixelsνi+pixeliactive pixelsln(eνi-1)
1

where νi is the probability to have photons in the pixel i, including the random background. Note that the formula is valid for the binary readout implemented in the RICH electronics. The centre of each Cherenkov ring is defined by the true particle direction. The DLL values are subsequently computed.

Table 4
Nominal values of the parameters used in the fast simulation

Simulated distributions for protons from Λ decays with p < 3.8 GeV /c are shown in Fig. 6. The average DLLx is 6.2, with an RMS of 4.0, for a fast simulation made using the nominal parameters. Figure 6a also shows the distributions after varying the detection efficiency by ±20 %. The corresponding distributions are shifted by ∓1.2 units. The random noise probability was changed by ±40 % from its nominal value which produces negligible variation as seen in Fig. 6b. This study shows that a variation of ±1 DLLx units is obtained by changing the photon detection efficiency by 15 %. A variation of the same size can be obtained by changing the angular resolution σc by 50 %.

Fig. 6
DLLx for protons with p < 3.8 GeV /c. In black the result with nominal simulation parameters are shown. a The red, dashed, and green, dash-dot, plots are for a change by + 20 and -20 % of the photon yield, b is for a change by + 40 % ...

The DLLx distributions for CMSPs obtained from the fast simulation with nominal parameters are consistent with those obtained from full simulation. Changing the photon detection efficiency by ±15 %, as inferred from the Λ study, a variation of ±2 DLLx units is obtained. An identical result can be obtained by changing the angular resolution.

In summary, an uncertainty of two DLLx units is inferred from the comparison of below-threshold protons in data and simulation when using Λ decays. The extrapolation to the CMSP regime is obtained by varying the simulation parameters and leads to an uncertainty of four DLLx units on the average DLLx value for such particles.

Uncertainties and results

After the ANN selection the signal prediction for the chosen model is 2.5, 0.9 and 0.3 events for the 124, 154 and 185  GeV /c2τ1~ masses, and below 0.1 for the other mass values. The expected background is negligible, less than 0.02 events for all the mass hypotheses.

A summary of the systematic uncertainties is given in Table 5. The total systematic uncertainties are approximately 5 % for the signal yield and 50 % for the background yield.

Table 5
Systematic uncertainties (in %) for the selection of the signal of CMSP pairs and on the background retention. Where relevant, the lower value corresponds to the CMSP mass of 124  GeV /c2 and the higher one to 309  GeV /c2

Two methods are used to determine the luminosity: a Van der Meer scan and a beam-gas imaging method [31]. The uncertainties on the integrated luminosities are 1.7 % for the 7  Te Vdata set and 1.2 % for the 8  TeV data set.

The efficiency for triggering, reconstructing and identifying high-pT muons has been studied in detail for the LHCb Z and W boson cross-section measurements [32, 33], and the agreement between data and simulation was found to be better than 2 %. This percentage is taken as the corresponding uncertainty for this analysis.

A further efficiency uncertainty is considered due to the delayed signals in the tracking and muon systems. The timing precision affects the amplitude recorded by the front-end electronic boards and the measurement of the drift time in the straw tubes. The effect on the efficiency due to a timing uncertainty of ±1 ns is determined from simulation as a function of the β of the particle. Subsequently, a weighted average of the uncertainty is obtained from the β distributions for each mass hypothesis, providing values varying from 2.2 to 3.7 % for CMSP masses from 124 to 309  GeV /c2.

The comparison of the IP distributions in data and simulated Z/γ → μ+μ- events indicates a maximal discrepancy of ±5 μm. By changing the requirement on the IP parameter by this amount, the corresponding efficiency variation is ±0.7 % for the background and negligible for the signal.

The hadron and electron rejection is affected by the calibration of the calorimeters. From the comparison of Z/γ → μ+μ- decays in data and simulation a relative uncertainty of 10 % is inferred. This translates into a 15 % change on the background yield, while the signal is almost unaffected.

The training of the ANN is affected by the uncertainty on the background and signal models. A 2.7 % contribution to the signal efficiency uncertainty is associated with the calibration procedures, determined by the comparison of data and simulation for Z/γ → μ+μ- and Λ → pπ decays. Error propagation is performed by modifying the ANN training sets, while keeping the test sets and the pair significance selection fixed. Adding the uncertainties in quadrature with the statistical uncertainty of 1 %, gives a total of 2.9 %.

The ANN selection leaves a very small amount of simulated background. The binomial uncertainty on the background retention is large, at approximately 50 %. This value is assigned as the uncertainty on the background selection efficiency.

As already stated, the acceptance A is affected by model uncertainties in the range from 5 to 9 % for τ1~ mass from 124 to 309  GeV /c2. In addition, the choice of the PDF affects the efficiency by modifying the momentum of the CMSP. By scanning various PDFs, we have found that this effect is small, not larger that 0.4 %, for all the models.

The cross-section upper limits are computed using the Feldman–Cousins method [34] for zero observed candidates, taking into account the expected number of background events and the uncertainties [35]. The predicted amount of background is so small that it has no sizeable effect on the result. The upper limits at a 95 % confidence level (CL) for CMSP pair production in the LHCb geometrical acceptance at s 7 and 8  TeV are listed in Table 6 and shown in Fig. 7 together with the theoretical cross-sections calculated for the particular model described in Sect. 3.1.

Table 6
Cross-section upper limits at 95 % CL for CMSP pair production in the LHCb acceptance in the 7 and 8 TeV
Fig. 7
Upper limits at 95 % CL on the cross-sections for the pair production of CMSPs in the LHCb acceptance (points) and the corresponding predictions assuming the Drell–Yan production of τ1~ (bands representing ±1σ uncertainty) with ...

Conclusions

A search for pairs of long-lived charged particles, with masses in the range 124–309  GeV /c2, using τ1~ pairs predicted by the mGMSB model as a benchmark scenario, is performed using data from proton–proton collisions at 7 and 8  Te V, corresponding to an integrated luminosity of 3.0 badhbox-1, collected with the LHCb detector in the forward pseudorapidity range 1.8 < η < 4.9. The candidates are assumed to interact only through the electroweak interaction in the detector. Hence, they behave like heavy muons and their main signature is the absence of a signal in the RICH detectors. The detection efficiency is limited to particles with β > 0.8 due to the acceptance in time after beam crossing. The main background contribution comes from Z/γ → μ+μ- and is reduced to less than  ∼ 0.02 events. No events have been observed. Upper limits are set on the Drell–Yan CMSP pair production cross-section. For proton–proton collisions at s 7  Te V, the 95 % CL upper limits for the production cross-section of a pair of CMSPs in the LHCb acceptance vary from 6.1 fb for a mass of 124  GeV /c2 up to 10.7 fb for a mass of 309  GeV /c2. At s 8  Te V, they vary from 3.4 to 5.7 fb for the same masses.

In LHCb the identification of CMSPs relies on the energy deposited in the subdetectors, the main discrimination power being provided by the RICH system. Together with the forward pseudorapidity coverage, this unique feature allows LHCb to complement the searches undertaken by the central detectors at the Tevatron and LHC.

Acknowledgments

We express our gratitude to our colleagues in the CERN accelerator departments for the excellent performance of the LHC. We thank the technical and administrative staff at the LHCb institutes. We acknowledge support from CERN and from the national agencies: CAPES, CNPq, FAPERJ and FINEP (Brazil); NSFC (China); CNRS/IN2P3 (France); BMBF, DFG, HGF and MPG (Germany); INFN (Italy); FOM and NWO (The Netherlands); MNiSW and NCN (Poland); MEN/IFA (Romania); MinES and FANO (Russia); MinECo (Spain); SNSF and SER (Switzerland); NASU (Ukraine); STFC (UK); NSF (USA). The Tier1 computing centres are supported by IN2P3 (France), KIT and BMBF (Germany), INFN (Italy), NWO and SURF (The Netherlands), PIC (Spain), GridPP (UK). We are indebted to the communities behind the multiple open source software packages on which we depend. We are also thankful for the computing resources and the access to software R&D tools provided by Yandex LLC (Russia). Individual groups or members have received support from EPLANET, Marie Skłodowska-Curie Actions and ERC (European Union), Conseil général de Haute-Savoie, Labex ENIGMASS and OCEVU, Région Auvergne (France), RFBR (Russia), XuntaGal and GENCAT (Spain), Royal Society and Royal Commission for the Exhibition of 1851 (UK).

Footnotes

1An anomalous signature in the RICH detectors could also be produced by tracks with wrongly assigned momenta. This can happen if the particle has an absolute electric charge that is different from the proton charge. For instance, a 1  GeV /c2 proton-like particle is expected to produce Cherenkov light in the two gaseous radiators when the measured momentum is above 30  GeV /c. If the particle has one third of the proton charge, the measured momentum is overestimated by a factor of three and this will lead to an incorrect calculation of the Cherenkov emission.

Authors are listed at the end of this paper.

References

1. Dimopoulos S, Thomas SD, Wells JD. Sparticle spectroscopy and electroweak symmetry breaking with gauge mediated supersymmetry breaking. Nucl. Phys. 1997;488:B39. doi: 10.1016/S0550-3213(97)00030-8. [Cross Ref]
2. Giudice GF, Rattazzi R. Theories with gauge-mediated supersymmetry breaking. Phys. Rep. 2011;332:419.
3. S. P. Martin, A supersymmetry primer. arXiv:hep-ph/9709356
4. ALEPH Collaboration, V. Barate et al., Search for pair-production of long-lived heavy charged particles in e+e-e+e- annihilation. Phys. Lett. B 405, 379 (1997). arXiv:hep-ex/9706013
5. DELPHI Collaboration, P. Abreu et al., Search for heavy stable and long-lived particles in e+e-e+e-=189 GeV. Phys. Lett. B 478, 65 (2000). arXiv:hep-ex/0103038
6. L3 Collaboration, P. Achard et al., Search for heavy neutral and charged leptons in e+e-e+e- annihilation at LEP. Phys. Lett. B 517, 75 (2001). arXiv:hep-ex/0107015
7. OPAL Collaboration, G. Abbiendi et al., Search for stable and long-lived massive charged particles in e+e-e+e-=130–209 GeV. Phys. Lett. B 572, 8 (2003). arXiv:hep-ex/0305031
8. H1 Collaboration, G. Aktas et al., Measurement of anti-deuteron photoproduction and a search for heavy stable charged particles at HERA. Eur. Phys. J. C 36, 413 (2004). arXiv:hep-ex/0403056
9. D0 Collaboration, V.M. Abazov et al., Search for long-lived charged massive particles with the D0 detector. Phys. Rev. Lett. 102, 161802 (2009). arXiv:0809.4472 [PubMed]
10. CDF Collaboration, T. Aaltonen et al., Search for long-lived massive charged particles in 1.96 TeV pp¯ collisions. Phys. Rev. Lett. 103, 021802 (2009). arXiv:0902.1266 [PubMed]
11. ATLAS Collaboration, G. Aad et al., Searches for heavy long-lived charged particles with the ATLAS detector in proton–proton collisions at s=8 TeV. JHEP 01, 68 (2015). arXiv:1411.6795
12. CMS Collaboration, S. Chatrchyan et al., Searches for long-lived particles in pppp and 8 TeV. JHEP 07, 122 (2013). arXiv:1305.0491
13. LHCb Collaboration, A.A. Alves Jr. et al., The LHCb detector at the LHC. JINST 3, S08005 (2008)
14. LHCb Collaboration, R. Aaij et al., LHCb detector performance. Int. J. Mod. Phys. A 30, 1530022 (2015). arXiv:1412.6352
15. Aaij R, et al. Performance of the LHCb vertex locator. JINST. 2014;9:P09007. doi: 10.1088/1748-0221/9/09/P09007. [Cross Ref]
16. Arink R, et al. Performance of the LHCb outer tracker. JINST. 2014;9:P01002. doi: 10.1088/1748-0221/9/01/P01002. [Cross Ref]
17. Alves AA, Jr, et al. Performance of the LHCb muon system. JINST. 2013;8:P02022.
18. Adinolfi M, et al. Performance of the LHCb RICH detector at the LHC. Eur. Phys. J. C. 2013;73:2431. doi: 10.1140/epjc/s10052-013-2431-9. [PMC free article] [PubMed] [Cross Ref]
19. Aaij R, et al. The LHCb trigger and its performance in 2011. JINST. 2013;8:P04022. doi: 10.1088/1748-0221/8/04/P04022. [Cross Ref]
20. Porod W. SPheno, a program for calculating supersymmetric spectra, SUSY particle decays and SUSY particle production at e+e-e+e- colliders. Comput. Phys. Commun. 2003;153:275. doi: 10.1016/S0010-4655(03)00222-4. [Cross Ref]
21. Allanach BC. The Snowmass points and slopes: benchmarks for SUSY searches. Eur. Phys. J. C. 2002;25:113. doi: 10.1007/s10052-002-0949-3. [Cross Ref]
22. W. Beenakker, R. Hoepker, M. Spira, PROSPINO: a program for the PROduction of Supersymmetric Particles In Next-to-leading Order QCD. arXiv:hep-ph/9611232
23. Nadolsky PM, et al. Implications of CTEQ global analysis for collider observables. Phys. Rev. D. 2008;78:013004. doi: 10.1103/PhysRevD.78.013004. [Cross Ref]
24. M. Krämer et al., Supersymmetry production cross sections in pp collisions at s=7. arXiv:1206.2892
25. Sjöstrand T, Mrenna S, Skands P. PYTHIA 6.4 physics and manual. JHEP. 2006;05:026. doi: 10.1088/1126-6708/2006/05/026. [Cross Ref]
26. Geant4 collaboration, S. Agostinelli et al., Geant4: a simulation toolkit. Nucl. Instrum. Methods A 506, 250 (2003)
27. Clemencic M, et al. The LHCb simulation application, Gauss: design, evolution and experience. J. Phys. Conf. Ser. 2011;331:032023. doi: 10.1088/1742-6596/331/3/032023. [Cross Ref]
28. Martin AD, et al. Parton distributions for the LHC. Eur. Phys. J. C. 2009;63:189. doi: 10.1140/epjc/s10052-009-1072-5. [Cross Ref]
29. Catani S, et al. Vector boson production at hadron colliders: a fully exclusive QCD calculation at next-to-next-to-leading order. Phys. Rev. Lett. 2009;103:082001. doi: 10.1103/PhysRevLett.103.082001. [PubMed] [Cross Ref]
30. LHCb Collaboration, R. Aaij et al., Inclusive WW TeV. JHEP 06, 058 (2012). arXiv:1204.1620
31. LHCb Collaboration, R. Aaij et al., Precision luminosity measurements at LHCb. JINST 9, P12005 (2014). arXiv:1410.0149
32. LHCb Collaboration, R. Aaij et al., Measurement of the forward ZZ TeV. arXiv:1505.07024 (submitted to JHEP)
33. LHCb Collaboration, R. Aaij et al., Measurement of the forward WW TeV. JHEP 12, 079 (2014). arXiv:1408.4354
34. Feldman GJ, Cousins RD. A unified approach to the classical statistical analysis of small signals. Phys. Rev. D. 1998;57:3873. doi: 10.1103/PhysRevD.57.3873. [Cross Ref]
35. Conrad J. Including systematic uncertainties in confidence interval construction for Poisson statistics. Phys. Rev. D. 2003;67:012002. doi: 10.1103/PhysRevD.67.012002. [Cross Ref]

Articles from Springer Open Choice are provided here courtesy of Springer