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

A search is performed for heavy long-lived charged particles using 3.0 fb-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rm fb}^{-1}$$\end{document} of proton–proton collisions collected at s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}$$\end{document}=\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$=$$\end{document} 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\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.8 < \eta < 4.9$$\end{document}. 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\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm {\,GeV\!/}c^2}$$\end{document}.


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 [1][2][3] 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/c 2 or higher. The τ 1 is the lightest mass eigenstate, resulting from the mixture of right-handed and e-mail: aurelio.bay@epfl.ch left-handed superpartners of the τ , dominated by the righthanded 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 [4][5][6][7][8][9][10][11][12]. With the exception of DELPHI [5], which had Cherenkov counters, the analyses are based on dE/dx and time-offlight 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/ p T ) µm, where p T 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 C 4 F 10 gas radiators, positioned directly after the VELO, and a downstream detector with a CF 4 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 fb −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 inef-ficiency 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-offlight 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.

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 (M m ), the number of messenger supermultiplets (N 5 ), the ratio of the vacuum expectation values of the two neutral Higgs fields (tan β), the sign of the Higgs mass parameter (μ), and the parameter C grav , 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 N 5 = 3, tan β = 15, μ > 0, M m = 2 , and the parameter C grav = 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 nextto-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].
Fully simulated signal samples, with masses varying from 124 to 309 GeV/c 2 , 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 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 A (TeV) m τ (GeV/c 2 ) σ (fb) A (%) 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/c 2 , 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 crosssection 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 TeV. 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, fol-lowed 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 highp T 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/c 2 , and 89 % for 309 GeV/c 2 . The two candidates must have opposite charge and each have p T > 50 GeV/c. To reject background from Z /γ → μ + μ − the pair must have a dimuon mass larger than 100 GeV/c 2 . 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 TeV (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 L = 3.0 fb −1 , according to the cross-section calculated with Prospino2.1, with SPS7 parameters and a τ 1 mass of 124 GeV/c 2 .

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. 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 TeV data sets, and for all the CMSP mass hypotheses considered. Figure 2 shows the distribution of the four ANN input variables for the 8 TeV data set, compared to the background and signal predictions; good agreement can be seen between data and simulated background.
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/c 2 CMSP mass hypotheses for both 7 and 8 TeV data 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/c 2 , decreasing for increasing mass to a minimum of 27 % at 309 GeV/c 2 . 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 .

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/c 2 , 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/c 2 to GeV/c 2 , 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/c 2 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. The energy in the VELO is given in units of minimum ionising particle (MIP) deposition. The first bin of the histogram for E in the ECAL has been multiplied by a factor 0.25 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. 1 An 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/c 2 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.
For illustration, the expected signal shapes for CMSPs with masses of 124 and 309 GeV/c 2 are superimposed. The small difference is due to a change in the underlying event and some light from the aerogel for the 124 GeV/c 2 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 an ad 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.
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. 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 N 0 (β 2 n 2 − 1)/(β 2 (n 2 −1)), where N 0 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, prob noise . 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 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. 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 %.
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/c 2 τ 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.
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 TeV data set and 1.2 % for the 8 TeV data set.  The efficiency for triggering, reconstructing and identifying highp T 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/c 2 .
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 train- ing 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/c 2 . 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.

Conclusions
A search for pairs of long-lived charged particles, with masses in the range 124-309 GeV/c 2 , using τ 1 pairs predicted by the mGMSB model as a benchmark scenario, is performed using data from proton-proton collisions at 7 and 8 TeV, corresponding to an integrated luminosity of 3.0 fb −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 TeV, the 95 % CL upper limits for the production crosssection of a pair of CMSPs in the LHCb acceptance vary from 6.1 fb for a mass of 124 GeV/c 2 up to 10.7 fb for a mass of 309 GeV/c 2 . At √ s = 8 TeV, 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.