We also include in our analysis LHC constraints from searches for
heavy long-lived charged particles (HLCP)
that are, in general, relevant to coannihilation regions where the
mass difference between the lightest SUSY particle (LSP) and the
next-to-lightest SUSY particle (NLSP)
may be small and the NLSP may therefore be long-lived. As we discuss
below, important roles
are played in our analysis by \(\tilde{\tau }_1\) , \(\tilde{\chi }^\pm _{1}\)  and \(\tilde{u}_R/\tilde{c}_R\)  coannihilation,
but only in the \(\tilde{\tau }_1\)  case is the NLSP - LSP mass difference small enough to
offer the possibility of a long-lived charged particle.
We implement in our global analysis the preliminary CMS 13-TeV result [1]}
using tracking and time-of-flight measurements,
based on the recipe and the efficiency map as a function of the pseudo-rapidity
and velocity of the HLCP given in [2]}.
We use Pythia 8 [3]} and Atom [4]} to generate and analyse the events,
and assume that the efficiencies for detecting slow-moving \(\tilde{\tau }_1\) s are similar at 8 and 13 TeV.
A similar recasting method was used in [5]}. See
also [6]} for another
approach using simplified model topologies.
The efficiency contains a lifetime-dependent factor \(\propto {\rm exp}(- d m/ p \tau )\) , where \(d\)  is a distance \(d \simeq 10\)  m
that depends on the pseudorapidity, and \(m, p\)  and \(\tau \)  are the mass, momentum and
lifetime of the long-lived particle. This factor drops rapidly for particles with
lifetimes \(\lesssim 10\)  ps, corresponding to \(m_{\tilde{\tau }_1} - m_{\tilde{\chi }^0_{1}} \gtrsim 1.6 \,\, \mathrm {GeV}\) .
