Search for heavy charged long-lived particles in the ATLAS detector in 36.1 fb⁻¹ of proton-proton collision data at = 13 TeV

A search for heavy charged long-lived particles is performed using a data sample of 36 . 1 fb − 1 of proton-proton collisions at ﬃﬃﬃ s p ¼ 13 TeV collected by the ATLAS experiment at the Large Hadron Collider. The search is based on observables related to ionization energy loss and time of flight, which are sensitive to the velocity of heavy charged particles traveling significantly slower than the speed of light. Multiple search strategies for a wide range of lifetimes, corresponding to path lengths of a few meters, are defined as model independently as possible, by referencing several representative physics cases that yield long-lived particles within supersymmetric models, such as gluinos/squarks ( R -hadrons), charginos and staus. No significant deviations from the expected Standard Model background are observed. Upper limits at 95% confidence level are provided on the production cross sections of long-lived R -hadrons as well as directly pair-produced staus and charginos. These results translate into lower limits on the masses of long-lived gluino, sbottom and stop R -hadrons, as well as staus and charginos of 2000, 1250, 1340, 430, and 1090 GeV, respectively.


I. INTRODUCTION
The search for heavy charged long-lived particles presented in this paper is based on a data sample of 36.1 fb −1 of proton-proton (pp) collisions at ffiffi ffi s p ¼ 13 TeV collected in 2015 and 2016.It utilizes observables related to large ionization energy loss (dE=dx) and time of flight (ToF), which are signatures of heavy charged particles traveling significantly slower than the speed of light.The mass of the particles is estimated using the dE=dx (m dE=dx ) or ToF (m ToF ) measurements together with the reconstructed momentum.The background is estimated in a purely data-driven manner, and multiple signal regions are defined to address the different possible signatures of heavy longlived particles (LLPs) that reach at least the ATLAS hadronic calorimeter, which corresponds to decay lengths of a few meters.Previous searches for LLPs that are stable within the detector were performed by the CMS Collaboration using 2.5 fb −1 of data at 13 TeV [1], as well as by the ATLAS Collaboration using 3.2 fb −1 of data at 13 TeV [2] and 19.8 fb −1 of data at 8 TeV [3].
LLPs are predicted in a variety of theories that extend the Standard Model (SM) [4].Theories with supersymmetry (SUSY) [5][6][7][8][9][10], which either violate [11][12][13] or conserve [4,[14][15][16][17][18][19]] R-parity, allow for the existence of charged long-lived sleptons ( l), squarks ( q), gluinos (g) and charginos ( χAE 1 ).Colored LLPs (e.g., q and g) would hadronize forming so-called R-hadrons [14], which are bound states composed of the LLP and light SM quarks or gluons, and may emerge from the collision as charged or neutral states.Through hadronic interactions of the light-quark constituents with the detector material, especially inside the calorimeters, R-hadrons can change to states with a different electric charge.Thus they might not be reconstructed as a consistently charged track in the inner tracking detector (ID) and in the muon spectrometer (MS), even if the lifetime is long enough to traverse the entire detector.Searches for Rhadrons are performed following two different approaches: using all available detector information ("full-detector Rhadron search"), or disregarding all information from the muon spectrometer ("MS-agnostic R-hadron search") to minimize the dependence on the modeling of R-hadron interactions with the material of the detector.Long-lived gluinos are motivated for example by split-SUSY models [18,19], in which high-mass squarks can lead to very long gluino lifetimes.Long-lived squarks, in particular a light top squark (stop) as the next-to-lighest SUSY particle, is motivated for example by electroweak baryogenesis [20,21], where nonuniversal squark mass terms can lead to a small mass difference between the stop and the neutralino as the lightest supersymmetric particle (LSP), and the lightest chargino is heavy, leading to suppressed radiative decays and long stop lifetimes.
If the charged LLP does not interact hadronically, it would predominantly lose energy via ionization as it passes through the ATLAS detector.Searches for long-lived charginos and sleptons (focusing on staus, as they are expected in most models to be the lightest) identified in both the ID and MS are therefore performed.The searches for staus are motivated by gauge-mediated SUSY breaking (GMSB) [15,[22][23][24][25][26][27] assuming the LSP to be a gravitino and with the light stau (τ) as an LLP, decaying to a τ-lepton and a gravitino with an unconstrained lifetime.While essentially all events in the GMSB models that include long-lived staus involve cascade decays that end in two LLPs, in this paper direct di-stau production through a Drell-Yan process is taken as a benchmark.Results for pair-produced charginos are motivated by a minimal anomaly-mediated supersymmetry breaking (mAMSB) model, where often the supersymmetric partners of the SM W-boson fields, the wino fermions, are the lightest gaugino states.In this particular case, the lightest of the charged mass eigenstates, a chargino, and the lightest of the neutral mass eigenstates, a neutralino, are both almost pure wino and nearly mass-degenerate, resulting in longlived charginos (see Refs. [28,29] for details).
This paper is organized as follows.A brief description of the ATLAS detector is given in Sec.II with emphasis on the parts relevant for this analysis, followed by details of the calibration of key observables in Sec.III.The dataset and simulated event samples and the subsequent event selection are described in Secs.IV and V, respectively.Section VI explains the method of background estimation and compares these estimates with data.Section VII details the origin and estimation of systematic uncertainties.The results including upper cross-section limits are shown in Sec.VIII, and finally the conclusions are summarized in Sec.IX.

II. ATLAS DETECTOR
The ATLAS detector [30] is a multipurpose particle detector consisting of the ID immersed in a 2 T solenoidal magnetic field, electromagnetic as well as hadronic calorimeters and a MS based on three large air-core toroid superconducting magnets with eight coils each.The ID comprises a silicon pixel detector, a silicon microstrip detector (SCT) and a transition-radiation tracker.With almost 4π coverage in solid angle, 1 the ATLAS detector is sensitive to the missing transverse momentum associated with each event.Several components are used to determine either dE=dx or ToF in this search and are discussed in more detail below.
The innermost component of the ATLAS detector is the pixel detector, consisting of four radial layers of pixelsensors in the barrel region and three disks on each side in the end cap region.The pixel detector measures the ionization energy loss of charged particles traversing it via a time-over-threshold (ToT) technique [31].
The ATLAS calorimeter in the central detector region consists of a liquid-argon electromagnetic calorimeter followed by a steel-absorber scintillator-tile sampling calorimeter.The latter serves as hadronic calorimeter covering the region up to jηj ¼ 1.7.In the ϕ-direction the tile calorimeter is segmented into 64 wedges with a size of Δϕ ¼ 2π=64 ≈ 0.1.In the r-z plane it is formed by three radial layers in the central barrel region and three in an extended barrel on each side, where each layer is further segmented into cells.The cells in the first and third layer in the barrel and in all layers in the extended barrel have a rectangular shape, while in the second layer of the barrel they are composed of two shifted rectangles (see Fig. 1).Overall, the tile calorimeter consists of 73 cells in 64 ϕ-segments resulting in 4672 individual cells.The optimal filtering algorithm (OFA) [32] used for the readout provides, besides a precise measurement of the energy deposits of particles and jets, a timing measurement.The resolution of the single-cell timing in the tile calorimeter is 1.3-1.7 ns.
The outermost part of the ATLAS detector is the MS.Immersed in toroidal magnetic fields, the MS provides particle tracking and momentum reconstruction for charged particles with jηj < 2.7, as well as triggering information in the range jηj < 2.4.Three layers of muon detectors are arranged in concentric shells at distances between 5 and 10 m from the interaction point (IP) in the barrel region, and in wheels perpendicular to the beam axis at distances between 7.4 and 21.5 m in the end cap regions where jηj > 1.05.Four different detector technologies are used in the MS.In the barrel region the muon trigger relies on resistive-plate chambers (RPCs), while thin-gap chambers are employed in the end cap wheels.High-precision tracking for the momentum measurements is performed by monitored drift tubes (MDTs), except for the innermost layer in the forward region of 2.0 < jηj < 2.7, where cathode-strip chambers are installed.A muon traversing the detector in the barrel region typically has around 20 hits in the MDTs and 14 hits in the RPCs.Both systems exhibit a sufficiently precise single-hit timing resolution (MDTs 3.2 ns, RPCs 1.8 ns) to distinguish relativistic muons traveling at almost the speed of light from slowly propagating stable massive particles.For the reconstruction of slow particles traversing the full detector, a dedicated tracking algorithm that treats the velocity β as a free parameter is used [33,34].
1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the center of the detector and the z axis coinciding with the axis of the beam pipe.The x axis points from the IP to the center of the LHC ring, and the y axis points upward.Cylindrical coordinates (r, ϕ) are used in the transverse plane, ϕ being the azimuthal angle around the beam pipe.The pseudorapidity is defined in terms of the polar angle θ as η ¼ − ln tanðθ=2Þ.Object distances in the η-ϕ plane are given by ΔR ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi ðΔηÞ 2 þ ðΔϕÞ 2 p .
The ATLAS trigger system [35] consists of a hardwarebased level-1 trigger followed by a software-based highlevel trigger, which runs reconstruction and calibration software similar to the offline reconstruction, and reduces the event rate to about 1 kHz.

III. CALIBRATION OF MAIN OBSERVABLES
To achieve optimal identification performance for heavy, charged, LLPs, measurements of the specific ionization energy loss and ToF have to be calibrated.
The pixel detector provides measurements of specific ionization energy losses.The dE=dx is corrected for ηdependence and on a run-by-run basis, to reduce the effect of the degradation of charge collection in the silicon sensor due to accumulated irradiation.Individual measurements resulting from a particle traversing the detector are combined in a truncated mean to reduce the effects of the Landau tails on the estimate of the most probable value (MPV).Finally, a value of the βγ of the particle is estimated from the MPV and the momentum, using a three-parameter empirical function, calibrated using protons, kaons and pions with low transverse momentum (p T ).The expected resolution for βγ from the dE=dx and momentum measurement is about 14%, which is taken as the uncertainty in the measured values.A detailed description of the βγ estimation using the pixel detector can be found in Ref. [31].
The particle velocity, β, is determined via ToF measurements in the tile calorimeter, the MDTs and the RPCs.To ensure an optimal β resolution, a series of custom calibrations is performed.High-p T muons are used for the calibration as, given the timing resolutions, they travel effectively at the speed of light.
For the tile calorimeter, only cells with a minimum energy deposit of 500 MeV are used.The time difference (t 0 ) relative to a particle traveling at the speed of light from the IP to the cell center (cell distance) is not allowed to be larger than 25 ns to reduce the effect of out-of-time pileup.As the first calibration step, a bias introduced by the OFA is corrected.The signal simulation reproduces the OFA bias reliably.The correction is very small for an in-time signal, while for late-arriving particles it ranges up to 10 ns.This is followed by a correction for the particle path in the cell.An effective spatial position corresponding to the timing measurement is estimated using an extrapolation of the particle track from the production vertex to the respective cell in the tile calorimeter (effective distance).The difference between the effective distance and the cell distance is used as a correction for the timing calibrations, while the effective distance is used for the estimation of the individual measurements of the velocity (β HIT TILE ).The corresponding correction is largest (3 ns) for the edge η regions of the largest cells in the tile calorimeter, which are located in the outermost layer of the extended barrel.This is followed by separate η-dependent corrections for data and simulation, which accounts for small remaining η-dependences of the timing for the outer parts of the cells in data, and ensures an η-independent timing measurement in simulation.Furthermore, a calibration of the timing as a function of the energy deposit per ϕ-projected cell (averaging cells of identical geometry over ϕ) is applied, which does not exceed 0.2 ns.Additionally, calibration constants in data are estimated for each cell followed by a run-by-run correction of the overall timing.These calibrations account for differences between single-tile calorimeter cells and for misalignments of the ATLAS and LHC clocks.The cellwise calibration factors are mostly well below 1 ns, while for some cells they range up to 2 ns.The run-by-run correction constants are between −0.6 and 0.2 ns.A smearing of the timing obtained in simulation is applied to achieve the same resolution as observed in data.The width of a Gaussian parametrization of the t 0 distribution serves as the uncertainty.The uncertainties are estimated as a function of the energy deposit per ϕ-projected cell.Finally, the uncertainties are adjusted to achieve a unit Gaussian width for the pull distribution, which gives a correction of 1% to the uncertainty.The individual measurements of β HIT TILE , with the ϕ-projected resolutions shown in Fig. 1, are combined in a weighted average β TILE using the inverse squared uncertainties as weights.The final resolution achieved for β TILE is σ β TILE ¼ 0.068. 2  For the MDTs and RPCs, each of the 323799 drift tubes and 362262 RPC readout strips is calibrated individually by performing a Gaussian parametrization of the timing information, and correcting for the offset to the expected value for a particle traveling at the speed of light.The width of the fitted distribution is taken as the uncertainty for the measurements in the respective tube/strip.The uncertainties are, in a manner similar to the tile calorimeter, adjusted to give a unit Gaussian width for the pull distributions.The corrections to the uncertainties are 16% for the MDTs, less than 1% for the RPC η-strips and 2% for the ϕ-strips.The squared inverse of the uncertainty is used as a weight for the calculation of a weighted-average β.Time-dependent phase-shift variations between the ATLAS and LHC clocks are addressed similarly to the tile calorimeter in a separate correction of timing information for each LHC run, and separately for MDTs, RPC readout strips measuring η and RPC strips measuring ϕ.The calibration results in a β resolution of σ β MDT ¼ 0.026 for the MDTs and σ β RPC ¼ 0.022 for the RPCs (with single-hit t 0 resolutions of about 3.2 and 1.8 ns, respectively).
Finally, the ToF-based β measurements in the different subsystems are combined into an overall β ToF , which is estimated as a weighted average of the β measurements from the different subsystems using the inverse squared uncertainties as weights.Furthermore, the combined 2 The resolutions (σ x ) are determined by performing a Gaussian fit to the core of the respective distributions.
uncertainty and the probability of compatibility between the measurements from the different subsystems are calculated for each candidate.The final distributions of the β measurements are shown in Fig. 2 for muons from Z → μμ selected events in data and simulation.The final resolution achieved for β ToF is σ β ToF ¼ 0.021.

IV. DATA AND SIMULATED EVENTS
The analysis presented in this paper is based on a data sample of 36.1 fb −1 of pp collisions at ffiffi ffi s p ¼ 13 TeV collected in 2015 and 2016, with a bunch spacing of 25 ns.
Reconstructed Z → μμ events in data and simulation are used for timing resolution studies.Simulated signal events are used to study the expected signal behavior.Pair production of gluinos (squarks), with masses between 400 and 3000 GeV (600 and 1400 GeV), was simulated in PYTHIA6 [36] (version 6.427) with the AUET2B [37] tuned set of underlying event and hadronization parameters (tune) and the CTEQ6L1 [38] PDF set, incorporating specialized hadronization routines [39,40] to produce final states containing R-hadrons.The other SUSY particle masses, except that of the lightest neutralino, were set to very high values to ensure negligible effects on gluino or squark pair-production.The fraction of gluino-balls, i.e., bound states of a gluino and gluon, was conservatively set to 10% [4,39], in order to account for the possibility of other (neutral and hence invisible) final states than those searched for.While the search is optimized for R-hadrons long-lived enough to reach at least the hadronic calorimeter, samples with gluino lifetimes of 10, 30 and 50 ns, where a significant fraction of the LLPs will decay before the calorimeter, are also investigated.
PYTHIA6 relies on a parton shower to add additional high-p T partons to the event.To achieve a more accurate description of QCD radiative effects, the PYTHIA6 events were reweighted to match the transverse-momentum distribution of the gluino-gluino or squark-squark system to the distribution obtained in dedicated leading-order MG5_AMC@NLO (version 2.2.3) [41] simulations with one additional parton in the matrix-element calculation.Pair production of staus in a GMSB scenario [m Messenger ¼ 500 TeV, C grav ¼ 100000, tan β ¼ 10, signðμÞ ¼ 1] was simulated with two additional partons at leading order, using MG5_AMC@NLO (version 2.3.3) in combination with PYTHIA8 [42] (version 8.212) and EvtGen [43] (version 1.2.0) with the A14 [44] tune and NNPDF23LO [45] PDF set.The pair-produced lightest stau mass eigenstate is a mixture of the left-handed and right-handed superpartners of the τ-lepton, although it is the partner of the right-handed lepton at the 99% level in the model considered here.
Pair production of charginos in an mAMSB scenario (m 0 ¼ 5 TeV, tan β ¼ 5) was simulated with two additional partons at leading order, using MG5_AMC@NLO (version 2.3.3) in combination with PYTHIA8 (version 8.212) and EvtGen (version 1.2.0) with the A14 tune and NNPDF23LO PDF set.
Samples of Z → μμ events, which are used for general testing, calorimeter and MS timing calibrations as well as the evaluation of systematic uncertainties, were simulated using POWHEG-BOX [46] (r2856) in combination with PYTHIA8 (version 8.186) and EvtGen (version 1.2.0) with the AZNLO [47] tune and CTEQ6L1 PDF set.
All events were passed through a full detector simulation [48] based on the GEANT4 framework [49].For the R-hadron simulations, hadronic interactions with matter were handled by dedicated GEANT4 [49] routines based on different scattering models: the model used to describe gluino (squark) R-hadron interactions is referred to as the generic (Regge) model [50].The R-hadrons interact only moderately with the detector material, as most of the R-hadron momentum is carried by the heavy gluino or squark, which has a small interaction cross section.Typically, the energy deposit in the calorimeters is less than 10 GeV.
All simulated events included a modeling of contributions from pileup by overlaying minimum-bias pp interactions from the same (in-time pileup) and nearby (out-of-time pileup) bunch crossings simulated in PYTHIA8 (version 8.186) and EvtGen (version 1.2.0) with the A2 [51] tune and MSTW2008LO [52] PDF set.The simulated events were reconstructed using the same software used for collision data, and were reweighted so that the distribution of the number of collisions per bunch crossing matched that of the data.

V. EVENT SELECTION
Five dedicated signal regions (SRs), imposing requirements on the entire event and on the individual candidate tracks, are defined in this section, addressing differences in topology and expected interactions with the detector for three different benchmark scenarios: staus, charginos and R-hadrons.
Events are selected by online triggers based on large missing transverse momentum ( ⃗p miss T , with magnitude denoted by E miss T ) or signatures of single (isolated) highmomentum muons.Large E miss T values are produced mainly when QCD initial-state radiation (ISR) boosts the R-hadron system, resulting in an imbalance between ISR and Rhadrons whose momenta are not fully accounted for in the E miss T calculation.The adopted triggers impose thresholds from 70 to 110 GeV on E miss T and 20 to 26 GeV on single muons, depending on the data-taking period.
The offline event selection requires all relevant detector components to be fully operational, a primary vertex (PV) built from at least two well-reconstructed charged-particle tracks, each with p T > 400 MeV, and at least one candidate track that meets the criteria specified below.
A "common track selection" is implemented for candidates in all SRs and is described below.ID tracks (denoted by the superscript "trk") are required to have a minimum p trk T of 50 GeV and a momentum measurement p trk < 6.5 TeV.The candidate track is required to be matched to the PV using 'loose' requirements on the transverse (d 0 ) and longitudinal (z 0 ) impact parameters (jd 0 j < 2 mm, jz 0 sin θj < 3 mm). 3To ensure good track reconstruction, the candidate track must have at least seven silicon clusters 4 (N clusters silicon > 6), no shared or split clusters in the pixel detector (N shared pix þ N split pix ¼ 0) [55], and at least three clusters5 in the SCT (N clustersþdead SCT > 2).Mainly to ensure a reliable timing measurement in the calorimeters, it is required that the sum of the track-p T in a cone of ΔR ¼ 0.2 around the candidate track is below 5 GeV.Jets reconstructed in the calorimeter are used to veto electrons and SM hadrons.Jets are reconstructed with the anti-k t clustering algorithm [56] with radius parameter R ¼ 0.4 and using as inputs clusters of energy deposits in the calorimeter, calibrated such that the average response of an electron is unity.The jets are then calibrated using the method described in Ref. [57].An electron veto is imposed by rejecting any candidate track for which the nearest jet with p T > 20 GeV and within a ΔR ¼ 0.05 cone around the track has at least 95% of its energy deposited in the electromagnetic calorimeter.A veto against SM hadrons is imposed by rejecting any candidate track for which any associated jet within a ΔR ¼ 0.05 cone of the track has an energy larger than the track momentum.Candidate tracks are required to have a cluster in the innermost 3 The transverse impact parameter is defined as the distance of closest approach between a track and the beam-line in the transverse plane.The longitudinal impact parameter corresponds to the z-coordinate distance between the primary vertex and the point along the track at which the transverse impact parameter is defined.The charge released by a moving charged particle is rarely contained within just one pixel; neighboring pixels/strips registering hits are joined together using a connected component analysis [53,54] to form clusters. pixel-detector layer (if expected) or a cluster in the second layer (if the innermost is not expected, but the second layer is).A "cosmics veto" rejects candidate tracks that have an opposite-sign track/muon on the other side of the detector, satisfying ΔR cosmics < 0.04. 6A "Z veto" rejects each candidate track that together with the highest-p T muon in the event forms an invariant mass within 10 GeVof the Z-boson mass.
Three different selection criteria for each track in an event, hereafter referred to as ID+CALO, LOOSE and TIGHT, are defined.The ID+CALO selection, which does not use MS information, is designed for the R-hadron searches.The LOOSE and TIGHT criteria are applied to ID-MS combined tracks in the search using the full detector information, with TIGHT providing stricter requirements.
An ID+CALO candidate selection starts from an ID track that fulfills the requirements of the above-mentioned common track selection, and which has at least one timing measurement in the tile calorimeter.In addition, the candidate tracks need to have jηj < 1.65, to ensure a minimal sensitivity in the tile calorimeter.To obtain a reliable dE=dx measurement in the pixel detector, the tracks are required to have at least two clusters used in the respective estimation (N clusters good dE=dx > 1), a dE=dx value between 0 and 20 MeV g −1 cm 2 , as well as an estimate of βγ (0.2 < ðβγÞ dE=dx < 10).To ensure a reliable ToF measurement in the calorimeter, individual β HIT TILE estimates and their individual uncertainties from calorimeter cells associated with the track are combined using a χ 2 , and the resulting weighted average β TILE is required to be between 0.2 and 2, with an uncertainty σ β TILE < 0.06 and a χ 2 probability Pð χ 2 ; NDFÞ > 0.01.
A LOOSE candidate selection starts from a combined ID-MS track with at least one timing measurement in the tile calorimeter or MS, and fulfilling the requirements of the above-mentioned common track selection for its ID track.In addition, the candidate tracks need to have a large combined transverse momentum (p cand T > 70 GeV), a combined momentum p cand < 6.5 TeV, jηj < 2, and hits in at least two MS stations.All ToF -based β measurements must be consistent, i.e., for candidates with β measurements in more than one ToF system (N systems > 1), the weighted means of all systems have to be consistent to within 5σ; and for candidates with only one system (N systems ¼ 1) the respective internal measurements have to be consistent to within 5σ.In cases where a dE=dx measurement from the pixel detector exists, the β estimates based on dE=dx and ToF have to be consistent to within 5σ.The uncertainty in the final β, whether from a single system or a weighted average of several systems, has to be below 0.025 (σ β ToF < 0.025) for the candidate to be accepted.At least one system has to yield a β measurement and a final value for β ToF between 0.2 and 2 is required.
A TIGHT candidate selection is identical to the LOOSE selection, except for a tighter pseudorapidity requirement (jηj < 1.65), an additional requirement on the pixel dE=dx measurement (1 < dE=dx < 20 MeVg −1 cm 2 ), and requiring at least two systems to yield a β measurement.
To target the three different benchmark scenarios this analysis uses five distinct selections, as shown in Table I.Signal regions are defined by imposing requirements on estimated masses, in addition to the criteria in Table I.The signal regions that are combined are designed to be orthogonal to allow for a combination in the statistical interpretation of the results.
The search for stable R-hadrons is performed in both an MS-agnostic SR (SR-Rhad-MSagno) and a full-detector SR (SR-Rhad-FullDet) approach.The former is much less dependent on the hadronic-interaction model for Rhadrons.It is based on events solely selected through E miss T triggers and candidates stemming from ID tracks fulfilling the ID+CALO selection plus the final selection requirements  I.If no LOOSE candidate is found in the event, then it is checked whether there is any candidate satisfying the ID+CALO requirements with the additional selections of Table I.The two selections are therefore mutually exclusive.In the rare case of events with more than one candidate, the candidate with the highest p T is chosen to improve background suppression.
The masses are derived from m ¼ p=βγ using momentum and ToF measurements to give m ToF , and, where applicable, ionization measurements to give m dE=dx .The final selection requirements on the masses are then obtained by fitting the reconstructed mass distribution in signal events for each simulated mass hypothesis with a Gaussian function, and taking the mean minus twice the width as the lower requirement on the respective mass hypothesis.The upper requirement on the mass is left open, and hence is constrained only by kinematics to be < 6.5 TeV.For the R-hadron SRs the lower requirements are evaluated in the two-dimensional m ToF -m dE=dx -plane, while for SRs SR-2Cand-FullDet and SR-1Cand-FullDet only the m ToF distribution is used.This choice is based on the fact that especially for low-mass candidates a significant fraction of candidates have dE=dx < 0.945 MeV g −1 cm 2 , at which point there is essentially no separation power between the various mass hypotheses.The final lower requirements on the masses are shown in Fig. 3.
The searches for pair-produced stable staus and charginos are performed using two orthogonal SRs.The region SR-2Cand-FullDet contains events with exactly two candidates fulfilling the LOOSE selection and the corresponding final selection requirements stated in Table I, while SR-1Cand-FullDet contains events with exactly one candidate, which in this case must satisfy the TIGHT selection and the corresponding final selection requirements stated in Table I.For the region SR-2Cand-FullDet, which has precedence over SR-1Cand-FullDet in the event categorization, the candidate with the lower m ToF is used to derive limits.
A set of 16 discovery regions (DRs) is defined for setting model-independent upper limits on cross sections and stating p 0 values.These DRs are indicated by larger markers in Fig. 3.The resulting DRs are: four for SR-Rhad-MSagno, four for SR-Rhad-FullDet (combined ID+CALO þ LOOSE), as well as four each for exclusive SR-2Cand-FullDet and SR-1Cand-FullDet regions.

VI. BACKGROUND ESTIMATION
The background is estimated with a fully data-driven method.First, the probability density functions (pdfs) of the key variables are determined from data, using sideband regions where possible.The key variables are momentum, β ToF and ðβγÞ dE=dx for the R-hadron SRs, and momentum and β ToF for the chargino and stau SRs.Distributions of expected background in m ToF (and m dE=dx ) are obtained by randomly sampling the pdfs and using the equation m ¼ p=βγ.
For this procedure to be valid, β ToF and ðβγÞ dE=dx must not be correlated with momentum.In principle this is true, as the primary background contribution is from highmomentum muons with mismeasured β ToF and ðβγÞ dE=dx , but there can be an implicit correlation via η, in particular as the p T requirement for reconstructing candidates translates into a momentum requirement p req > 200 GeV even for jηj < 2.0.This means that momentum is correlated with η for jηj > 1.75 due to the lack of low-momentum tracks.
As β ToF and ðβγÞ dE=dx are correlated with η due to the different resolutions in different detector regions, the result is some correlation of β ToF and ðβγÞ dE=dx with momentum.To remove the effect of these correlations, the pdfs are estimated in five (SR-Rhad-MSagno, SR-Rhad-FullDet, SR-1Cand-FullDet) or six (SR-2Cand-FullDet) jηj bins.A variable binning is used to account for the different regions of the subsystems used.For the sampling of the background, η from the candidate is used to get the corresponding pdfs.This is a safe procedure as effects from signal contamination are negligible.
To derive the SR-Rhad-MSagno and SR-Rhad-FullDet momentum pdfs, the final requirements on β ToF and ðβγÞ dE=dx are inverted, while for the β ToF and ðβγÞ dE=dx pdfs the final requirement on the momentum is inverted and a minimum momentum of 50 GeV is required.For the SR-2Cand-FullDet, the momentum pdf can be estimated from a sideband in β ToF , but not the β ToF pdf, because in the high-jηj region no candidates with p < 100 GeV are left.SR-1Cand-FullDet instead uses sidebands for the momentum pdfs, since for the tight candidates jηj < 1.65 is required and hence tracks with p < 200 GeV for the full range are available.For the β pdf in SR-1Cand-FullDet and SR-2Cand-FullDet, a minimum momentum of 70 GeV is required.The final requirements are summarized in Table I.
The background yield is normalized to data using lowmass control regions (CRs).The number of estimated events is scaled, so that the number of events in the CR matches the number observed in data.In the R-hadron searches, the CR is identical to the SR except that the minimum mass requirement is replaced by an upper limit on m dE=dx and m ToF , which both have to be below 300 GeV.Similarly, in the chargino/stau search, the minimum mass requirement on m ToF is replaced by an upper limit of 150 GeV and 200 GeV for SR-2Cand-FullDet and SR-1Cand-FullDet, respectively.
The background estimates overlaid with MC signal events at the expected mass limit are shown in Figs.4-7.The statistical uncertainty from the pdfs is propagated to the background estimate and shown as gray bands.For all SRs, agreement between data and estimated background in the low-mass regions is found.

VII. SYSTEMATIC UNCERTAINTIES
The two major uncertainties in the signal yields are from the theoretical cross section and the modeling of ISR, as well as the dedicated full-detector track reconstruction and the ToF-based β measurement in the MS in some cases.All individual contributions are outlined below and summarized in Table II.

A. Theoretical cross sections
R-hadron production cross sections are calculated to next-to-leading order (NLO) in the strong coupling constant, adding the resummation of soft-gluon emission at next-to-leading-logarithm accuracy (NLO þ NLL) [58][59][60][61][62]. Stau and chargino signal cross sections are calculated to NLO in the strong coupling constant (NLO) using PROSPINO2 [63].The nominal cross section and the uncertainty is taken from an envelope of cross-section predictions using different PDF sets and factorization and renormalization scales, as described in Ref. [64].This prescription results in an uncertainty in the cross section of between 14% (at 600 GeV) and 57% (at 3000 GeV) for gluino R-hadrons, and between 14% (at 600 GeV) and 23% (at 1400 GeV) for squark R-hadrons.For direct pairproduction of staus and chargino pair-production the uncertainty is between 6% (at 290 GeV) and 10% (at 910 GeV) and between 4% (at 200 GeV) and 10% (at 1500 GeV), respectively.

B. Signal efficiency
Missing-transverse-momentum triggers used in these searches rely solely on calorimeter-energy deposits to calculate the transverse energy, and are thus largely blind to muons, which can therefore be used for calibration and systematic uncertainties.To evaluate the trigger efficiency, the trigger turn-on curve is obtained by fitting the measured efficiency as a function of E miss T in Z → μμ events, in both data and simulation.These efficiency turn-on curves are then applied to the E miss T spectrum from simulated events.The total uncertainty is estimated from four contributions: the relative difference between the efficiencies obtained using the fitted threshold curves from Z → μμ in data and  simulation, the differences in efficiency obtained from independent AE1σ variations in fit parameters relative to the unchanged turn-on-curve fit for both Z → μμ data and simulation, and a 10% variation of the E miss T to assess the scale uncertainty.The E miss T trigger is estimated to contribute a total uncertainty of 1.6% and 2% to the signal efficiency for R-hadrons and staus/charginos, respectively.
To account for a possible mismodeling of the singlemuon trigger timing in simulation, all simulated events were reweighted, depending on the β and jηj of the candidate, to match the data and to compute systematic uncertainties.Systematic uncertainties of 1%, 4% and 3% are assigned to R-hadrons, staus and charginos, respectively, by taking the difference in signal trigger efficiency between unweighted and reweighted events.
To address a possible mismodeling of ISR, and hence E miss T in the signal events, half of the difference between the selection efficiency for the nominal PYTHIA6 events and those reweighted to match MG5_AMC@NLO predictions is taken as an uncertainty in the expected signal yield, and found to be below 15% in all cases.For staus and charginos, the uncertainties in the amount of ISR (and FSR) are evaluated by varying generator parameters in the simulation.The uncertainty from the choice of renormalizaton/factorization scale is evaluated by changing the default scale by a factor of two in MG5_AMC@NLO.The uncertainty from the choice of CKKW-L merging [65,66] scale is evaluated by taking the maximum deviation from the nominal MC sample when varying it by a factor of two.The uncertainty from the parton-shower generator tuning is evaluated using PYTHIA8 tune variations.The systematic uncertainty on the ISR/FSR is calculated by adding all three components in quadrature.The uncertainties are evaluated for three mass points of staus and charginos and are estimated to be 4% and 5%, respectively.The latter procedure was also implemented for the MG5_AMC@NLO gluino-gluino and squark-squark samples, and the resulting uncertainties are similar to the nominal-vs-reweighted ones.
The uncertainty in the pileup modeling in simulation is found to affect the signal efficiency by between 0.1% and 5.5%, typically decreasing as a function of the simulated LLP mass and varying with benchmark model.
An additional systematic uncertainty in the signal efficiency is estimated for the dedicated tracking algorithm to cover all discrepancies between data and simulation, by randomly rejecting 10% of the reconstructed objects.The effect on the final signal efficiency is found to be between 0.2% and 14.8%.TABLE III.Expected signal yield (N exp ) and acceptance (a) × efficiency (ε), estimated background (N est ) and observed number of events in data (N obs ) for the full range of simulated masses in the MS-agnostic R-hadron search.

R-hadron
Simulated   Uncertainties in the measurement of dE=dx in the pixel detector result in uncertainties in the signal yield ranging from 6% to 13% for R-hadrons and less than 0.5% for staus and charginos.They account for both the shape difference between the ionization distribution in data and simulation, and the scale shift in data due to radiation damage.
The systematic uncertainty in the calorimeter-based β estimation is assessed by scaling the calorimeter-cell-time smearing of simulated events by AE5% and by varying the cell-time correction introduced to correct for the bias due to the OFA by AE50%, and is found to be below 2% in all cases.
The systematic uncertainty in the MS-based β estimation is derived by varying the MDT smearing constants by AE10% to bracket the distribution seen in data.
By comparing the signal efficiency with and without the correction for incorrectly modeled timing behavior in the RPCs, the overall uncertainty is found to be between 1.0% and 21.7%.

C. Integrated luminosity
The uncertainty in the combined 2015 þ 2016 integrated luminosity is 2.1%.It is derived, following a methodology similar to that detailed in Ref. [67], from a calibration of the luminosity scale using x-y beam-separation scans performed in August 2015 and May 2016.

D. Background estimation
To estimate the systematic uncertainty of the background estimate, three main contributions are considered: the  systematic uncertainty of the normalization as well as the influence of the jηj binning and the definition of the sidebands.The normalization uncertainty, derived by applying the relative statistical uncertainty in the number of data events in the CR to the estimated number of background events in the SR, is the dominant contribution for most SRs.Only at very high masses do the other contributions become significant and eventually take over.
To test the stability of the choice of jηj binning, the number of jηj bins in the pdfs is varied from nominal five (six) to three (four) and six (nine) for SR-Rhad-MSagno, SR-Rhad-FullDet and SR-1Cand-FullDet (SR-2Cand-FullDet).For each pdf, the background is estimated and half the maximal difference in the number of background events in the SR is taken as a systematic uncertainty.To check for an influence of the sideband definition, the selection requirements are varied, the pdfs re-evaluated correspondingly, and half the maximal difference of background counts in the SR, using the nominal and varied pdfs, is taken as the systematic uncertainty.FIG. 9. Expected (dashed lines) and observed (marked solid lines) upper cross-section limits in the gluino (top), sbottom (bottom left) and stop (bottom right) R-hadron searches, respectively, using two independent and not to be combined approaches based on SRs SR-Rhad-MSagno (light blue) and SR-Rhad-FullDet (red).The shaded light-blue/light-red bands represent the AE1σ uncertainties in the expected limits.The result obtained using the former SR has a much reduced dependence on the modeling of R-hadron interaction with matter with respect to the other, therefore the two results must not combined.The theory prediction along with its AE1σ uncertainty is shown as a white line and a dark-blue band, respectively.

Mass
the MS-agnostic and full-detector R-hadron, stau and chargino searches, respectively.
As can be seen in Tables III-VIII, no significant excess of observed data events above the expected background is found in the examined mass ranges and signal regions.The yields are summarized for a subset of discovery regions (see Sec. V) in Fig. 8 and Table V, the latter also showing the p 0 values.
Upper limits at 95% confidence level (C.L.) are placed on the production cross sections for various benchmark models, as shown in Figs. 9 and 10.These limits are obtained from the expected signal and the estimated background in signal region SR-Rhad-MSagno (SR-Rhad-FullDet or SR-1Cand-FullDet/SR-2Cand-FullDet) using a one-bin (two-bin) counting experiment applying the CL s prescription [68].Model-independent upper limits defined as cross section × acceptance × efficiency for the abovementioned discovery regions are shown in Table V.Given the predicted theoretical cross sections, also shown in Figs. 9 and 10, the cross-section limits are translated into lower limits on masses for the various benchmark models.
The MS-agnostic search yields expected lower limits at 95% C.L. on the R-hadron masses of 2060, 1270 and 1345 GeV for the production of long-lived gluino, sbottom and stop R-hadrons, respectively.The corresponding observed lower limits on the masses are 1950, 1190 and 1265 GeV.The expected signal yield (N exp ) and efficiency, estimated background (N est ) and observed number of events in data (N obs ) for the full range of simulated masses can be found in Table III.The sensitivity first increases and then decreases with increasing R-hadron mass.The same effect is visible in the total efficiency of the E miss T trigger, and is due to the change in production channel for gluino R-hadrons from gluon-initiated to quark-initiated with increasing mass.
For metastable gluino R-hadrons, the MS-agnostic search yields expected lower limits on mass of 1980, 2080 and 2090 GeV for lifetimes of 10, 30 and 50 ns, respectively.The corresponding observed lower limits are 1860, 1960 and 1980 GeV.The expected signal yield (N exp ) and efficiency, estimated background (N est ) and observed number of events in data (N obs ) for the full range of simulated masses can be found in Table VI.As additional E miss T can arise from one of the gluino R-hadrons decaying before the calorimeters, the sensitivity increases for mass regions that have limited trigger efficiency (low and high masses).With decreasing lifetime, and thereby increasing probability for both gluino R-hadrons to decay before the calorimeters, the signal efficiency drops, as does the sensitivity of the search optimized for R-hadrons longlived enough to exit the detector.However, the expected lower limits on mass at lifetimes of 50 and 30 ns remain more stringent than those of the search targeting lowlifetime metastable gluino R-hadrons [31].Expected and observed lower limits on mass as a function of lifetime for both the MS-agnostic and full-detector searches for gluino R-hadrons are shown in Fig. 11.
Using the full-detector search, the expected (observed) lower limits on the mass are 2050 GeV (2000 GeV), 1280 GeV (1250 GeV) and 1370 GeV (1340 GeV) for the production of long-lived gluino, sbottom and stop Rhadrons, respectively.The expected signal yield (N exp ) and efficiency, estimated background (N est ) and observed number of events in data (N obs ) for the full range of simulated masses can be found in Table IV.The sensitivity of the R-hadron full-detector search in mass regions with noticeable background yields is slightly better than that of the MS-agnostic search, due mainly to the increased signal efficiency when including the single-muon trigger and an improved β-resolution when using full-detector candidates.However, this specific search is consequently slightly more model-dependent, especially on the modeling of hadronic interactions.
Expected (observed) lower limits on mass for direct production of staus and charginos are set at 420 GeV (430 GeV) and 1070 GeV (1090 GeV), respectively.The expected signal yield (N exp ) and efficiency, estimated background (N est ) and observed number of events in data (N obs ) for the full range of simulated masses can be found in Tables VII and VIII.

IX. CONCLUSION
A search for heavy, charged, long-lived particles is performed using a data sample of 36.1 fb −1 of protonproton collisions at ffiffi ffi s p ¼ 13 TeV collected by the ATLAS experiment at the Large Hadron Collider at CERN.The search is based on observables related to large ionization losses, measured in the innermost tracking detector, and slow propagation velocities, measured in the tile calorimeter and muon spectrometer.Both observables are signatures of heavy charged particles traveling significantly slower than the speed of light.No significant deviations from the expected background are observed.Upper limits at 95% confidence level are provided on the production cross sections of long-lived R-hadrons, as well as directly pair-produced staus and charginos.These results translate into lower limits on the masses of long-lived gluino, sbottom and stop R-hadrons, as well as staus and charginos of 2000, 1250, 1340, 430 and 1090 GeV, respectively.America.In addition, individual groups and members have received support from BCKDF, CANARIE, CRC and Compute Canada, Canada; COST, ERC, ERDF, Horizon 2020, and Marie Skłodowska-Curie Actions, European Union; Investissements d' Avenir Labex and Idex, ANR, France; DFG and AvH Foundation, Germany; Herakleitos, Thales and Aristeia programmes co-financed by EU-ESF and the Greek NSRF, Greece; BSF-NSF and GIF, Israel; CERCA Programme Generalitat de Catalunya, Spain; The Royal Society and Leverhulme Trust, United Kingdom.The crucial computing support from all WLCG partners is acknowledged gratefully, in particular from CERN, the ATLAS Tier-1 facilities at TRIUMF (Canada), NDGF (Denmark, Norway, Sweden), CC-IN2P3 (France), KIT/GridKA (Germany), INFN-CNAF (Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Taiwan), RAL (UK) and BNL (USA), the Tier-2 facilities worldwide and large non-WLCG resource providers.Major contributors of computing resources are listed in Ref. [69].FIG.11.Expected (dashed blue line) and observed (solid blue line) lower mass limits at the 95% C.L. level for gluino R-hadrons with different mean lifetimes derived using the MS-agnostic analysis (SR-Rhad-MSagno).Expected (empty circles) and observed (filled circles) limits on stable gluino R-hadrons are shown for both the MS-agnostic (SR-Rhad-MSagno) and the full-detector analysis (SR-Rhad-FullDet).

FIG
FIG. 1. Resolution (σ β ) of the β HIT TILE measurement for the different tile-calorimeter cells, which are shown with their actual shape.The gray circles indicate, the distance, d, from the interaction point and the gray straight lines the direction in η.

4
FIG.3.The lower mass requirement (m min ) defining the final counting regions in the m ToF -m dE=dx plane used in the R-hadron searches (left) and the m ToF -distribution used in the chargino and stau searches (right) for the respective simulated mass of the signal at particle level (m truth ).The larger markers indicate the lower mass requirements for the "discovery regions."

FIG. 4 .
FIG. 4. Background estimate for the MS-agnostic analysis targeting gluino R-hadrons [SR-Rhad-MSagno (ID+CALO)] in the m ToF -m dE=dx plane (top), the m ToF -projection (bottom left) and the m dE=dx -projection (bottom right).The last bin of each distribution includes the overflow.The dashed red lines in the upper figure indicate the lower bounds of the signal region for a representative signal choice, while the dotted gray lines illustrate the upper bound for the control region.The signal (2200 GeV gluino R-hadron) is indicated by green markers in the upper and green dash-dotted lines in the lower plots.The lower panels show the ratio of observed data to estimated background.The shaded gray area shows the statistical uncertainty of the background estimate.

FIG. 5 .
FIG.5.Background estimate for the full-detector analysis (LOOSE part) targeting gluino R-hadrons (SR-Rhad-FullDet) in the m ToF -m dE=dx plane (top), the m ToF -projection (bottom left) and the m dE=dx -projection (bottom right).The last bin of each distribution includes the overflow.The dashed red lines in the upper figure indicate the lower bounds of the signal region for a representative signal choice, while the dotted gray lines illustrate the upper bound for the control region.The signal (2200 GeV gluino R-hadron) is indicated by green markers in the upper and green dash-dotted lines in the lower plots.The lower panels show the ratio of observed data to estimated background.The shaded gray area shows the statistical uncertainty of the background estimate.

FIG. 6 .
FIG.6.Background estimate for the full-detector analysis (ID+CALO part) targeting gluino R-hadrons (SR-Rhad-FullDet) in the m ToF -m dE=dx plane (top), the m ToF -projection (bottom left) and the m dE=dx -projection (bottom right).The last bin of each distribution includes the overflow.The dashed red lines in the upper figure indicate the lower bounds of the signal region for a representative signal choice, while the dotted gray lines illustrate the upper bound for the control region.The signal (2200 GeV gluino R-hadron) is indicated by green markers in the upper and green dash-dotted lines in the lower plots.The lower panels show the ratio of observed data to estimated background.The shaded gray area shows the statistical uncertainty of the background estimate.

FIG. 7 .
FIG.7.Background estimate for the analysis targeting pair-produced staus and charginos in the one-TIGHT -candidate SR (left) and the two-LOOSE -candidates SR (right).The last bin of each distribution includes the overflow.The dashed red lines indicate the lower bound of the signal region for a representative signal choice, while the dotted gray lines illustrate the upper bound for the control region.The signal (1200 GeV chargino) is indicated by green dash-dotted lines.The lower panels show the ratio of observed data to estimated background.The shaded gray area shows the statistical uncertainty of the background estimate.

FIG. 8 .
FIG.8.Expected and observed events in the 16 discovery regions introduced in Sec.V and for the respective CRs, including the systematic uncertainties.For data the Poisson error is shown.The lower plot shows the significance for SRs with a surplus of events between observed and estimated.The bin contents of SRs with inclusive mass ranges are correlated.
distributions observed in data together with the background estimate, its statistical uncertainty and a representative expected signal are shown in Figs.4 FIG.10.Expected (dashed red line) and observed (marked solid red line) upper cross-section limits using combined two-and onecandidate SRs for stau pair production (left) and chargino pair production (right).The shaded dark-red (light-red) bands represent the AE1σ (AE2σ) uncertainties in the expected limits.The theory prediction along with its AE1σ uncertainty is shown as a white line and a blue band, respectively.

TABLE I .
Summary of the five sets of SRs.The trigger as well as the track candidate selection and the number of candidates per event required for the respective SR are given.Also the final selection requirements together with the mass window (one-or two-dimensional) for the final counting are stated.SRs in one block (delimited by horizontal lines) are combined for the statistical interpretation of the results.For SR-Rhad-FullDet, the ID+CALO is used as a fallback only if no LOOSE candidates are found in the event; hence the two SRs are mutually exclusive.
stated in TableI.The full-detector search, on the other hand, takes full advantage of the MS, both in terms of triggering events and in using additional ToF measurements in the MDTs and RPCs.Candidates originate from combined ID-MS tracks fulfilling the LOOSE selection and passing the final requirements states in Table

TABLE II .
Summary of systematic uncertainties.Ranges indicate a dependence on the mass hypothesis.
mass [GeV] N exp AE σ N exp a × ε AE σ a×ε N est AE σ N est N obs

TABLE IV .
Expected signal yield (N exp ) and acceptance (a) × efficiency (ε), estimated background (N est ) and observed number of events in data (N obs ) for the full range of simulated masses in the full-detector R-hadron search.N exp AE σ N exp a × ε AE σ a×ε N est AE σ N est N obs N exp AE σ N exp a × ε AE σ a×ε N est AE σ N est N obs

TABLE V .
p 0 values and model-independent upper limits on cross section (σ) × acceptance (a) × efficiency (ε) for the 16 discovery regions.

TABLE VI .
Expected signal yield (N exp ) and acceptance (a) × efficiency (ε), estimated background (N est ) and observed number of events in data (N obs ) for the full range of simulated masses in the MS-agnostic search for metastable gluino R-hadrons.N exp AE σ N exp a × ε AE σ a×ε N exp AE σ N exp a × ε AE σ a×ε N exp AE σ N exp a × ε AE σ a×ε N est AE σ N est N obs

TABLE VII .
Expected signal yield (N exp ) and acceptance (a) × efficiency (ε), estimated background (N est ) and observed number of events in data (N obs ) for the full range of simulated masses in the full-detector direct-stau search.N exp AE σ N exp a × ε AE σ a×ε N exp AE σ N exp N obs N exp AE σ N exp a × ε AE σ a×ε N est AE σ N est