Search for long-lived charged particles in proton-proton collisions at sqrt(s) = 13 TeV

Results are presented of a search for heavy stable charged particles produced in proton-proton collisions at sqrt(s) = 13 TeV using a data sample corresponding to an integrated luminosity of 2.5 inverse femtobarns collected in 2015 with the CMS detector at the CERN LHC. The search is conducted using signatures of anomalously high energy deposits in the silicon tracker and long time of flight measurements by the muon system. The data are consistent with the expected background, and upper limits are set on the cross sections for production of long-lived gluinos, top squarks, tau sleptons, and leptonlike long-lived fermions. These upper limits are equivalently expressed as lower limits on the masses of new states; the limits for gluinos, ranging up to 1610 GeV, are the most stringent to date. Limits on the cross sections for direct pair production of long-lived tau sleptons are also determined.


Introduction
Many extensions of the standard model (SM) include heavy long-lived charged particles that might have high momentum, but speed significantly smaller than the speed of light [1][2][3] and/or charge, Q, not equal to the elementary charge ±1e [4][5][6][7].Those particles with lifetimes greater than a few nanoseconds can travel distances larger than the size of a typical collider detector and appear quasi-stable like the pion or kaon.These particles are generally referred to as heavy stable charged particles (HSCPs) and can be singly (|Q| = 1e), fractionally (|Q| < 1e), or multiply (|Q| > 1e) charged.Without dedicated searches, HSCPs may be misidentified or unobserved, since charged particle identification algorithms at hadron collider experiments generally assume that particles have speeds close to the speed of light and charges of ±1e.Additionally, HSCPs may be charged during only a part of their passage through detectors [8] further limiting the ability of standard algorithms to identify them.
For HSCP masses greater than about 100 GeV, a significant fraction of particles produced at the LHC will have a relative velocity β ≡ v/c < 0.9.It is possible to distinguish |Q| ≥ 1e particles with β < 0.9 from light SM particles traveling close to the speed of light through their higher rate of energy loss via ionization (dE/dx) or through their longer time of flight (TOF) to the outer detectors.This paper describes a search for HSCPs using the CMS detector in two ways: (i) requiring tracks to be reconstructed only in the silicon detectors, the tracker-only analysis; (ii) requiring tracks to be reconstructed in both the silicon detectors and the muon system, referred to as the tracker+TOF analysis.

CMS
Figure 1: Distribution of the dE/dx estimator, I h (see Section 3.1), versus particle momentum for tracks in the 13 TeV data, and for simulation of HSCP for singly or multiply charged particles with masses of 400 and 1000 GeV.The vertical scale shows the density of entries for data only.
The dependence of dE/dx on the particle momentum is described by the Bethe-Bloch formula [9].This dependence can be seen in Fig. 1, which shows the dE/dx estimator versus momentum for good quality (Section 4) high transverse momentum (p T > 55 GeV) tracks from data and the generated Monte Carlo (MC) samples for HSCP signals with various charges.In the momentum range of interest at the LHC (10-1000 GeV), SM particles have nearly uniform ionization energy loss (≈3 MeV/cm).Searching for candidates with larger dE/dx gives sensitivity to massive particles with |Q| ≥ 1e.
In the present paper, results of searches for singly and multiply charged HSCPs in 2.5 fb −1 of data collected with the CMS detector at √ s = 13 TeV in 2015 are presented.Similar limits on HSCPs were recently obtained by the ATLAS experiment [30,31]

Signal benchmarks
The analyses described in this paper employ several HSCP models as benchmarks, to account for a range of signatures that are experimentally accessible.
The first type of signal consists of HSCPs that interact via the strong force and hadronize with SM quarks to form R-hadrons [2,3].As in Ref. [27], events involving direct pair production of gluinos ( g) and top squarks ( t 1 ), with mass values in the range 300-2600 GeV, are generated according to the Split Supersymmetry (Split SUSY) scenarios [32][33][34][35].Gluinos are generated assuming the squark mass is 10 TeV [32,36].In the region of parameter space where squarks are too heavy to be produced at the LHC, the gluino-gluino production cross section and kinematic distributions depend only on the gluino mass, thus the cross section limits are modelindependent.PYTHIA 8.153 [37], with the underlying event tune CUETP8M1 [38], is used to generate the 13 TeV MC samples.The fraction, f , of produced g hadronizing into a g-gluon state is an arbitrary value of the hadronization model.It determines the fraction of R-hadrons that are neutral at production.For this search, results are obtained for two different values of f , 0.1 and 0.5.As in Ref. [27], two scenarios of R-hadron strong interactions with nuclear matter are considered.The first scenario follows the cloud model in Refs.[8,39], which assumes that the R-hadron is surrounded by a cloud of colored, light constituents that interact during scattering.Therefore, the R-hadron interacting inside the detector may change its charge sign.The second scenario adopts a model of complete charge suppression [40] where the R-hadron becomes a neutral particle before it enters the muon system.Both the tracker-only and tracker+TOF analyses are used to search for these signals, although only the tracker-only analysis is expected to have sensitivity in the charge-suppressed scenario.In the case of a discovery, a comparison of the numbers of events found in the two analyses could give a hint about the nature of the new long-lived particle.
The second type of signal consists of HSCPs that behave like leptons.The minimal gauge mediated supersymmetry breaking (mGMSB) model [41] is selected as a benchmark for leptonlike HSCPs.Production of quasi-stable sleptons at the LHC can proceed either directly or via production of heavier supersymmetric particles (mainly squarks and gluinos) that decay and lead to two sleptons at the end of the decay chain.This latter process is dominant because the direct production process is electroweak.Direct production is relevant only if squarks and gluinos are too heavy to be produced at the LHC.The mGMSB model is explored using the SPS7 slope [42], which has the tau slepton (stau τ 1 ) as the next-to-lightest supersymmetric particle (NLSP).The particle mass spectrum and the decay table are produced with the program ISASUGRA 7.69 [43].The mGMSB model is characterized by six fundamental parameters.The mGMSB parameter Λ, which corresponds to the effective supersymmetry breaking scale, is varied from 31 to 510 TeV.It is proportional to the sparticle masses.The range of its values gives a tau slepton mass of 100 to 1600 GeV.Other parameters are fixed to the following values.The number of the mes-senger SU(5) multiplets N mes = 3 and their mass scale M is set as M mes /Λ = 2.The ratio of the vacuum expectation values of the Higgs doublets is tan β = 10 and a positive sign of the higgsino mass term, µ > 0, is assumed.The large value of the scale factor of the gravitino coupling, C grav = 10000, results in a long-lived τ 1 .The SUSY mass spectrum produced is input to PYTHIA 6.4 [37] with the Z2* tune [44] as the generator for a MC simulation at 13 TeV.Two tau slepton samples are generated for each SUSY point: one with all processes (labeled "GMSB stau") and one with only direct pair production (labeled "Pair prod.stau").The pair-produced stau includes only τ 1 , which is predominantly τ R for these model parameters.The direct production of long-lived stau is model independent.Both cross section and kinematics depend only on the stau mass and the scan over the stau mass parameter shows the effect of variations in center-of-mass energy and integrated luminosity.The tracker-only and tracker+TOF analyses are both used to search for these signals.
The last type of signal is based on modified Drell-Yan (DY) production of long-lived leptonlike fermions.In this scenario, new massive spin-1/2 particles have arbitrary electric charge but are neutral under SU(3) Colour and SU(2) Left , and therefore couple only to the photon and the Z boson.PYTHIA v6.4 [37] with the Z2 * tune [44] is used to generate these 13 TeV MC signal samples.Simulations of events with leptonlike fermions are generated with masses ranging from 100 to 2600 GeV and for electric charges |Q| = 1e and 2e.
Different PYTHIA tunes were studied and the effects on the kinematic distribution were negligible for the HSCPs considered.The tracker-only and tracker+TOF analyses are both expected to have sensitivity to |Q| = 2e HSCPs.
In all signal samples, simulated minimum bias events are overlaid with the primary collision to produce the effect of additional interactions in the same LHC beam crossing (pileup).

The CMS detector
The central feature of the CMS [45] apparatus is a 3.8 T superconducting solenoid of 6 m internal diameter.Within the solenoid volume are a silicon tracker, a lead tungstate crystal electromagnetic calorimeter, and a brass and scintillator hadron calorimeter, each composed of a barrel and two endcap sections.Outside the solenoid, forward calorimeters extend the pseudorapidity (η) coverage provided by the barrel and endcap detectors.Muons are measured in gas-ionization detectors embedded in the steel flux-return yoke outside of the solenoid.The missing transverse momentum vector p miss T is defined as the projection on the plane perpendicular to the beam axis of the negative vector sum of the momenta of all reconstructed particles in an event.Its magnitude is referred to as E miss T .
The silicon tracker, consisting of 1440 silicon pixel and 15 148 silicon strip detector modules , measures charged particles within the pseudorapidity range |η| < 2.5.Isolated particles of transverse momentum p T = 100 GeV and with |η| < 1.4 have track resolutions of 2.8% in p T and 10 (30) µm in the transverse (longitudinal) impact parameter [46].Muons are measured in the pseudorapidity range |η| < 2.4, using three technologies: drift tubes (DTs), cathode strip chambers (CSCs), and resistive-plate chambers (RPCs).Matching muons to tracks measured in the silicon tracker results in a relative transverse momentum resolution for muons with 20 < p T < 100 GeV of 1.3-2.0% in the barrel and better than 6% in the endcaps.The p T resolution in the barrel is better than 10% for muons with p T up to 1 TeV [47].
The first level (L1) of the CMS trigger system, composed of custom hardware processors, uses information from the calorimeters and muon detectors to select events of interest within a fixed time interval of less than 4 µs.The high-level trigger (HLT) processor farm further decreases the event rate from around 100 kHz to less than 1 kHz, before data storage.A more detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [45].

dE/dx measurements
For the reconstructed track, information about dE/dx can be gained from measurements of ionization deposited in layers of the pixel and silicon tracker.The ionization charge measured is compared with that expected from a Minimum-Ionizing Particle (MIP), and its level of compatibility can provide a probability, using a dE/dx discriminator.As in Ref. [24], to distinguish SM particles from HSCP candidates the I as discriminator is used and is given by where N is the number of measurements in the silicon-tracker detectors, P i is the probability for a MIP to produce a charge smaller or equal to that of the ith measurement for the observed path length in the detector, and the sum is over the track measurements ordered in terms of increasing P i .
In addition, the dE/dx of a track is estimated using a harmonic-2 estimator: where c i is the charge per unit path length in the sensitive part of the silicon detector of the ith track measurement.The harmonic-2 estimator has units MeV/cm and the summation includes just the top 85% of the charge measurements.Ignoring the low charge measurements increases the resilience of the estimator against instrumental biases.This procedure is not necessary for I as which is, by construction, robust against that type of bias.
The mass of a candidate particle can be calculated [27] from its momentum and its I h dE/dx estimate, based on the relation: where the empirical parameters K = 2.684 ± 0.001 MeV cm −1 and C = 3.375 ± 0.001 MeV cm −1 are determined from data using a sample of low-momentum protons.As the momentum reconstruction is done assuming |Q| = 1e particles, Eq. ( 3) leads to an accurate mass reconstruction only for singly charged particles.
The HSCP candidates are preselected using the I as discriminator because it has a better signalto-background discriminating power compared to the I h estimator or the mass.Nonetheless, the mass is used at the last stage of the analysis, after the I as selection, to further discriminate between signal and backgrounds since the latter tend to have a low reconstructed mass.

Time of flight measurements
The time of flight to the muon system can be used to discriminate between particles travelling at near the speed of light and slower candidates.Both the DT and the CSC muon systems measure the time of each hit.In the DT, the precision position is obtained from this time measurement.The synchronization works in a such a way that a relativistic muon produced at the interaction point gives an aligned pattern of hits in consecutive DT layers.For a slower HSCP particle, hits in each DT layer will be reconstructed as shifted with respect to its true position and will form a zigzag pattern with an offset proportional to the particle delay, δ t .In the CSC the delay is measured for each hit separately.Each δ t measurement can be used to determine the track β via the equation: where L is the flight distance.The track β −1 value is calculated as the weighted average of the β −1 measurements from the DT and CSC systems associated with the track.The weight for the i th DT measurement is given by where n is the number of φ projection measurements found in the muon chamber producing the measurement and σ DT is the time resolution of the DT measurements, for which the measured value of 3 ns is used.The factor (n − 2)/n accounts for residuals computed using the two parameters of a straight line determined from the same n measurements.The minimum number of hits in a given DT chamber that allows for at least one residual calculation is n = 3.
The weight for the ith CSC measurement is given by where σ i , the measured time resolution, is 7.0 ns for cathode strip measurements and 8.6 ns for anode wire measurements.
The resolution on the weighted average β −1 measurement is approximately 0.065 in both the DT and CSC subsystems.

Data selection
All events pass a trigger requiring either a reconstructed muon with high transverse momentum or large E miss T , calculated using an online particle-flow algorithm [48][49][50].
The muon trigger is more efficient than the E miss T trigger for all HSCP models considered with the exception of the charge-suppressed R-hadron model, but it is not efficient for particles that are slow (β < 0.6).
The E miss T trigger can recover some events in which the HSCP is charged in the tracker and neutral in the muon subsystem.The particle-flow algorithm rejects tracks reconstructed only in the tracker and having a track p T significantly greater than the matched energy deposited in the calorimeter [49], as would be the case for HSCPs that become neutral in the calorimeter.Thus only an HSCP's energy deposit in the calorimeter, roughly 10-20 GeV, will be included in the E miss T calculation.Where one or more HSCPs fail to be reconstructed as muon candidates, the events may appear to have significant E miss T .
For both the tracker-only and the tracker+TOF analyses, the muon high-level trigger requires a muon candidate with p T > 50 GeV and the E miss T trigger requires E miss T > 170 GeV.Using these two triggers for both analyses allows for increased sensitivity to HSCP candidates that arrive in the muon system very late, as well as for hadronlike HSCPs, which may be charged only in the tracker.
Offline, for the tracker-only analysis, all events are required to have a candidate track with p T > 55 GeV as measured in the tracker, relative uncertainty in p T (σ p T /p T ) less than 0.25, |η| < 2.1, and the track fit χ 2 /dof < 5.The magnitudes of the impact parameters d z and d xy must both be less than 0.5 cm, where d z and d xy are the longitudinal and transverse impact parameters with respect to the vertex with the smallest d z .The requirements on the impact parameters are very loose compared to the resolutions for tracks in the tracker.Candidates must pass isolation requirements in the tracker and calorimeter.The tracker isolation criterion is ∑ p T < 50 GeV, where the sum is over all tracks (except the candidate) within ∆R = √ (∆η) 2 + (∆φ) 2 = 0.3 of the candidate track.The calorimeter isolation criterion is E/p = 0.3, where E is the sum of energy deposited in the calorimeter towers within ∆R = 0.3 and p is the track momentum reconstructed from the tracker.Candidate tracks must have at least two measurements in the silicon pixel detector and at least six measurements in the strip detectors.In addition, there must be measurements in at least 80% of the silicon layers between the first and last measurements of the track.To reduce the contamination from clusters with a large energy deposition due to overlapping tracks, a filtering procedure is applied to remove clusters in the silicon strip tracker that are not consistent with the passage of a singly charged particle (i.e., a narrow cluster with most of the energy deposited in one or two strips).After cluster filtering, there must be at least six measurements in the silicon tracker that are used for the dE/dx calculation.
The tracker+TOF analysis applies the same criteria, but additionally requires a reconstructed muon matched to the track in the inner detectors.At least eight independent time measurements are needed for the TOF computation.Finally, 1/β > 1 and σ 1/β < 0.15 are required.

Background estimation
For background estimation we follow the procedure described in our previous work [27].Candidates passing the preselection (Section 4) are subject to either two or three additional criteria to improve the signal-to-background discrimination.By choosing two uncorrelated criteria it is possible to predict the background using the ABCD (matrix) method.In this approach, the expected background in the signal region, D, is estimated by BC/A, where B and C are the number of candidates that fail the first or second criterion, respectively, while A is the number of candidates that fail both criteria.
Results are based upon a comparison of the number of candidates passing the selection criteria defining the signal region with the number of predicted background events in that region.Fixed selections on the appropriate set of I as , p T , and 1/β are used to define the final signal region (and the regions for the background prediction).The values are chosen to give the best discovery potential over the signal mass regions of interest.
For the tracker-only analysis, the two criteria are p T > 65 GeV and I as >0.3.The candidates passing only the I as requirement fall into the B region and those passing only the p T requirement fall into the C region.The B and C candidates are then used to form binned probability density functions in I h and p, respectively, such that, using the mass value (Eq.( 3)), the full mass spectrum of the background in the signal region D can be predicted.However, the η distribution of candidates with low dE/dx differs from the distribution of candidates with high dE/dx.To correct for this, events in the C region are weighted such that the η distribution matches that in the B region.For the tracker+TOF analysis, a three-dimensional matrix method is used with p T > 65 GeV, I as > 0.175, and 1/β > 1.25, creating eight regions labeled A-H.Region D represents the signal region, with events passing all three criteria.The candidates in the A, F, and G regions pass only the 1/β, I as , and p T criteria, respectively, while the candidates in the B, C, and H regions fail only the p T , I as , and 1/β criteria, respectively.The E region contains events that fail all three criteria.Background estimates can be made from several different combinations of these regions.The combination D = AGF/E 2 is used because it yields the smallest statistical uncertainty.As in the tracker-only analysis, events in the G region are reweighted to match the η distribution in the B region.The spread in background estimated from the other combinations is less than 20%, which is taken as the systematic uncertainty in the collision background estimate.The same 20% systematic uncertainty is used for the tracker-only analysis.
In order to check the background prediction, samples with a loose selection, which would be dominated by background tracks, are used for the tracker-only and tracker+TOF analyses.The loose selection sample for the tracker-only analysis is defined as p T > 60 GeV and I as > 0.10.The loose selection sample for the tracker+TOF analysis is defined by p T > 60 GeV, I as > 0.05, and 1/β > 1.05.Figure 2 shows the observed and estimated mass spectra for these samples.
For both analyses, an additional requirement on the reconstructed mass is applied.The specific requirement is adapted to each HSCP model.For a given signal mass and model, the mass requirement is M ≥ M reco − 2σ, where M reco is the average reconstructed mass for the given mass M HSCP and σ is the expected resolution.Simulation is used to determine M reco and σ.
Table 1 lists the final selection criteria, the predicted number of background events, and the number of events observed in the signal region.Agreement between prediction and observation is seen for both the tracker-only and the tracker+TOF analyses.Figure 3 shows the predicted and observed mass distributions for the tracker-only and the tracker+TOF analyses with the final selection.
Table 1: Selection criteria for the two analyses with the number of predicted and observed events.In the background prediction, the statistical and systematic uncertainties are added in quadrature.

Systematic uncertainties
The sources of systematic uncertainty considered are those related to the background prediction, the signal acceptance, and the integrated luminosity.The uncertainty in the integrated luminosity is 2.7% at √ s = 13 TeV [51].The uncertainties in the collision background predictions are estimated to be at the level of 20% for the tracker-only and the tracker+TOF analyses, as described in Section 5.
The signal acceptance is obtained from MC samples of the various signals processed through the full detector simulation (Section 2).Systematic uncertainties are derived by comparing the response of the detector in the data and simulation.The relevant sources of uncertainty are discussed below.
The signal trigger efficiency is dominated by the muon triggers efficiency, for all the models except the charge-suppressed ones.The uncertainty in the muon trigger efficiency has many contributions.It is estimated from the difference between the trigger efficiency in data and that seen in simulation, using Z(µµ) data.For genuine muons, the trigger efficiency uncertainty is 3%.
For slow moving particles, the effect of the timing synchronization of the muon system is tested by shifting the arrival times in simulation by the synchronization accuracy observed in data, resulting in an efficiency change of less than 4% for most samples but up to 8% for the 2.4 TeV gluino sample.The uncertainty in the E miss T trigger efficiency is found by varying the jet energy scale in the simulation of the high-level trigger by its uncertainty in data.The E miss T uncertainty is found to be less than 12% for all samples.The total trigger uncertainty is found to be less than 13% for all the samples, since the muon trigger inefficiencies are often compensated by the E miss T trigger and vice versa.
Low-momentum protons are used to compare the observed and simulated distributions of I h and I as that reflect the energy loss in the silicon tracker.The dE/dx distributions of signal samples are varied by the observed differences in order to estimate the systematic uncertainty.The uncertainty in the signal acceptance is usually less than 10%, and is at most 15%.
Bias in the energy loss measurement due to highly ionizing particles (HIP), such as lowmomentum protons produced in pp collisions earlier than the triggering collision, was also considered as a source of uncertainty in the I h estimate.In 2015, the LHC collision frequency was doubled, with bunches colliding every 25 ns compared to collisions every 50 ns in 2012, causing an increase of the HIP rate.The contribution of HIPs was included in simulations with the rate observed during the 2015 data taking.The uncertainty in this rate is found to be 25% and 80% for pixel and strip sensors, respectively.Varying the HIP rate in the simulation by these amounts leads to a change in signal acceptance of at most 4% for both analyses.
Dimuon events are used to test the MC simulation of 1/β by comparing with data.An offset of at most 1.5% is found for the muon system.The resulting uncertainty (labeled "Time of flight" in Table 2) in the signal acceptance is found to be less than 5% by shifting 1/β by this amount.
As in Ref. [26], the uncertainties in the efficiencies for muon [47] and track [52] reconstruction are each less than 2%.The track momentum uncertainty is estimated by shifting the momentum of the inner track, as in Ref. [26].This uncertainty is found to be less than 5% for most of the samples, increasing to 20% for masses above 2 TeV.
The uncertainty in the number of pileup events is evaluated by varying ±5% the minimum bias cross section used to calculate the weights applied to signal events in order to reproduce the pileup observed in data.The uncertainties due to pileup estimated with this procedure are The total systematic uncertainty in the signal acceptance is the sum in quadrature of the uncertainties due to the sources discussed above.For almost all signal models, it is less than 20% for both analyses.Only for the tracker+TOF analysis of the gluino ( f = 0.5) sample it is larger, but does not exceed 25%.
Table 2 summarizes the systematic uncertainties for the two analyses.As the uncertainty often depends on the model and HSCP mass, the largest systematic uncertainty is reported for each source.

Results
No significant excess of events is observed above the predicted background.Cross section limits are placed at 95% CL using a CL s approach [53][54][55] where a profile likelihood technique [56] is used.It utilizes a log-normal model [57,58] for the nuisance parameters, which are the integrated luminosity, the signal acceptance, and the expected background in the signal region.
The observed limits are shown in Fig. 4 for both the tracker-only and the tracker+TOF analyses along with the theoretical predictions.The theoretical cross sections are computed at NLO or NLO+NLL [59][60][61][62] using PROSPINO [63] with CTEQ6.6MPDFs [64].The uncertainty bands of the theoretical cross sections include the PDF uncertainty, the renormalization and factorization scale uncertainties, and the uncertainty in α s .The 95% CL limits on the production cross sections are shown in Tables 3, 4, 5, and 6 for long-lived gluino, top squark, tau slepton, and modified Drell-Yan signals, respectively.The limits were determined from the numbers of events passing all final criteria (including the mass criteria).
Mass limits are obtained from the intersection of the observed limit and the central value of the theoretical cross section.The tracker-only analysis excludes f = 0.1 gluino masses below 1610 (1580) GeV for the cloud interaction model (charge-suppressed model).Top squark masses below 1040 (1000) GeV are excluded for the cloud (charge-suppressed) models.In addition, the tracker+TOF analysis excludes τ 1 masses below 490 (240) GeV for the GMSB (direct pair production) model.Drell-Yan signals with |Q| = 1e (2e) are excluded below 550 (680) GeV.
The mass limits obtained at √ s = 13 TeV for various HSCP signal models are summarized in Table 7 and compared with earlier results at √ s = 7 and 8 TeV [27].A significant increase  Figure 4: Results of the HSCP search as the cross section upper limits at 95% CL for various signal models for the tracker-only analysis (left) and tracker+TOF analysis (right) at √ s = 13 TeV.In the legend, "CS" stands for charge-suppressed interaction model. in mass limit is obtained for all models with large QCD production cross section (gluinos, top squarks, and inclusive production of GMSB tau sleptons), arising from the higher centerof-mass energy pp collisions delivered by the LHC.For scenarios with much smaller crosssections, directly pair-produced tau sleptons and Drell-Yan signals with |Q| = 1e, the results do not improve, because the larger integrated luminosity at 7 and 8 TeV with respect to that at 13 TeV prevails over the effect of the increase of the centre-of-mass energy.For the |Q| = 2e analysis, results from the previous analysis optimized for multiply charged signals [27] are also provided.

Summary
A search for heavy stable charged particles produced in proton-proton collisions at √ s = 13 TeV using the CMS detector is presented.Two complementary analyses were performed: using only the tracker and using both the tracker and the muon system.The data are found to be compatible with the expected background.Mass limits for long-lived gluinos, top squarks, tau sleptons, and multiply charged particles are calculated.The models for R-hadronlike HSCPs include a varying fraction of g-gluon hadronization and two different interaction models leading to a variety of exotic experimental signatures.The limits are significantly improved over those from Run 1 of the LHC, and the limits on long-lived gluinos, ranging up to 1610 GeV, are the most stringent to date.Table 3: Summary of the search for long-lived gluinos: the p T (GeV), I as , 1/β, and mass thresholds M (GeV) requirements, the predicted and observed yields passing these criteria, and the resulting expected (exp.) and observed (obs.)cross section limits.The signal efficiencies and theoretical (theo.)cross sections are also listed.Table 6: Summary of the search for long-lived particles from modified Drell-Yan models of various charge: the p T (GeV), I as , 1/β, and mass thresholds M (GeV) requirements, the predicted and observed yields passing these criteria, and the resulting expected (exp.) and observed (obs.)cross section limits.The signal efficiencies and theoretical (theo.)cross sections are also listed.[27].In the model name, "CS" stands for charged suppressed interaction model and "DY" for Drell-Yan.The limits for doubly charged particles are also compared to the earlier CMS results obtained with the 'multiply charged' analysis, which was specifically designed to search for multiply charged particles.

Mass
using 3.2 fb −1 of 13 TeV data collected in 2015.

Figure 2 :
Figure 2: Observed and predicted mass spectra for loose selection candidates in the tracker-only (left) and tracker+TOF (right) analyses.The expected distributions for representative signals are shown as histograms.

Figure 3 :
Figure 3: Observed and predicted mass spectra for candidates passing the final selection in the tracker-only (left) and tracker+TOF (right) analyses.The expected distributions for representative signals are shown as histograms.

Table 2 :
Systematic uncertainties for the two HSCP searches.All values are relative uncertainties in the signal acceptance for the tracker-only and tracker+TOF analyses.

Table 4 :
Summary of the search for long-lived top squarks: the p T (GeV), I as , 1/β, and mass thresholds M (GeV) requirements, the predicted and observed yields passing these criteria, and the resulting expected (exp.) and observed (obs.)cross section limits.The signal efficiencies and theoretical (theo.)cross sections are also listed.

Table 5 :
Summary of the search for long-lived tau sleptons: the p T (GeV), I as , 1/β, and mass thresholds M (GeV) requirements, the predicted and observed yields passing these criteria, and the resulting expected (exp.) and observed (obs.)cross section limits.The signal efficiencies and theoretical (theo.)cross sections are also listed.

Table 7 :
Mass limits obtained at √ s = 13 TeV for various HSCP candidate models compared with earlier results for √ s = 7 + 8 TeV