Constraints on the pMSSM, AMSB model and on other models from the search for long-lived charged particles in proton-proton collisions at sqrt(s) = 8 TeV

Stringent limits are set on the long-lived lepton-like sector of the phenomenological minimal supersymmetric standard model (pMSSM) and the anomaly-mediated supersymmetry breaking (AMSB) model. The limits are derived from the results presented in a recent search for long-lived charged particles in proton-proton collisions, based on data collected by the CMS detector at a centre-of-mass energy of 8 TeV at the Large Hadron Collider. In the pMSSM parameter sub-space considered, 95.9% of the points predicting charginos with a lifetime of at least 10 ns are excluded. These constraints on the pMSSM are the first obtained at the LHC. Charginos with a lifetime greater than 100 ns and masses up to about 800 GeV in the AMSB model are also excluded. The method described can also be used to set constraints on other models.


Introduction
We present new constraints on long-lived chargino production in the phenomenological minimal supersymmetric standard model (pMSSM) [1] and on the anomaly-mediated supersymmetry breaking (AMSB) model [2][3][4].In the pMSSM, we consider the parameter sub-space for particle masses up to 3 TeV and charginos ( χ ± ) with a mean proper decay length (cτ) greater than 50 cm.In the AMSB model, the small mass difference between the lightest chargino and neutralino ( χ 0 1 ) often leads to a long chargino lifetime.Long-lived charged particles are predicted by various extensions of the standard model (SM) [5][6][7], such as supersymmetry (SUSY) [8] and theories with extra dimensions [9,10].If such particles have a mass lighter than a few TeV they could be produced by the CERN LHC.The energy available in the proton-proton (pp) collisions at the LHC is such that particles with mass 100 GeV and lifetime greater than O(1) ns could be observed with the CMS detector [11] as high-momentum tracks with anomalously large rates of energy loss through ionization (dE/dx).These particles could also be highly penetrating such that the fraction reaching the CMS muon system would be sizable.The muon system could therefore be used to help in identification and in the measurement of the time-of-flight (TOF) of the particles.The signature described above is exploited in a previous CMS search [12], which sets the most stringent limits to date on a number of representative models predicting massive long-lived charged particles such as tau sleptons (staus), top squarks, gluinos, and leptons with an electric charge between e/3 and 8e.
The main thrust of this paper is to present constraints on the pMSSM and the AMSB model, obtained using the results from the search for heavy stable charged particles (HSCP) in protonproton collisions at √ s = 7 and 8 TeV [12].The method applied relies on the factorization of the acceptance in terms of two probabilities, calculated using the standard CMS simulation and reconstruction tools at 8 TeV.Moreover, the results of the acceptance calculations have been tabulated and made publicly available [13].Thus this technique, which allows the signal acceptance for a model with long-lived particles to be computed using the kinematic properties of the particles at their production point, may in the future be used by others to evaluate constraints on other extensions of the SM without use of the CMS software.

The CMS detector
The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter providing a field of 3.8 T. Within the field volume are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter, and a brass and scintillator hadron calorimeter.Muons are measured in gas-ionization detectors embedded in the steel flux-return yoke of the magnet.The tracker measures charged particles within the pseudorapidity range |η| < 2.5 and provides a transverse momentum p T resolution of about 2.8% for 100 GeV particles.The analog readout of the tracker also enables the particle ionization energy loss to be measured with a resolution of about 5%.Muons are measured in the range |η| < 2.4, with detection planes based on three technologies: drift tubes, cathode strip chambers, and resistive-plate chambers.The muon system extends out to 11 m from the interaction point in the z direction and to 7 m radially.Matching tracks in the muon system to tracks measured in the silicon tracker results in a transverse momentum resolution between 1 and 5%, for p T values up to 1 TeV.The time resolution of the muon system is of the order of 1 ns.This provides a time-of-flight measurement that can be used to determine the inverse of the long-lived particle velocity as a fraction of the velocity of light (1/β) with a resolution of 0.065 over the full η range [12].A more detailed description of the CMS detector, of the coordinate system, and of the kinematic variables can be found in Ref. [11].
The first level of the CMS trigger system, composed of custom hardware processors, uses information from the calorimeters and muon detectors to select the events of interest.The high-level trigger processor farm further decreases the event rate from around 100 kHz to around 400 Hz for data storage.

Estimation of signal acceptance
Transcribing the results presented in Ref. [12], in terms of limits on models other than those considered in the reference, requires the tabulation of the relevant signal acceptance.This acceptance can then be used to constrain the beyond the standard model (BSM) scenarios, as described in the next sections.There can be significant differences in signal acceptance between the models investigated in Ref. [12] and the model to be tested.These differences arise because the dE/dx and TOF measurements are affected by the distribution and orientation of the material encountered by particles travelling within the CMS detector.The combination of these effects with the differences in the kinematic properties between models can result in large differences in signal acceptance.
The signal acceptance can be accurately computed by fully simulating and reconstructing signal events in the CMS detector and by applying the same selection criteria as adopted in Ref. [12].We refer to this method of determining the signal acceptance as the "full simulation".This procedure requires extensive knowledge of the CMS detector and in particular use of CMS software that accurately models detector response to simulated signal events and employs the full CMS reconstruction routines such that identical selection criteria can be used on simulated signal events as on data collected from collisions.
An alternate method, which only requires information on the kinematic properties of the longlived particles at their production point, is presented here.Such a method can be used if the efficiency for triggering and selecting events can be expressed in terms of probabilities associated with each individual long-lived particle.This is the case for models with lepton-like massive long-lived particles since the event selection specified in Ref. [12] imposes only requirements on measurements performed on individual particles.The adjective "lepton-like", defined in Ref. [12], indicates particles that do not interact strongly and are therefore not subject to hadronization.
The event selection requirements of Ref. [12] are expressed in terms of measured p T , dE/dx, TOF, and mass values of individual particles in an event.The probability that a long-lived particle in an event passes the online or offline selection requirements in Ref. [12] can be expressed as a function of the true (generator-level) kinematic properties (k) of the particle: β, η, and p T .Any other equivalent set of kinematic and directional variables could also be used to express this dependence.The offline selection of Ref. [12] has fixed values for p T , dE/dx, and TOF thresholds but uses different reconstructed mass thresholds (m thresh ) depending on the mass (m) of the long-lived particle in the model being tested.
With these individual particle probabilities, the acceptance A for a model to pass the online and offline selections can be computed with a Monte Carlo technique by generating a large number of events N such that where is the probability that the event with index i containing M longlived particles with true kinematic properties k 1 i , k 2 i , . . ., k M i passes the online selection, and is the probability that the event with the same kinematic properties passes the offline selection with mass threshold m thresh , after having passed the online selection.Given the mass resolution of the detector, a mass threshold of m thresh 0.6m has an efficiency of about 95% for the benchmark models considered in Ref. [12].Consequently, throughout this paper, mass thresholds of 0, 100, 200, and 300 GeV are used for true long-lived particle masses of m ≤ 166, ≤330, ≤500, and ≥500 GeV, respectively.The choice of 100 GeV steps for the mass thresholds is made so that the information provided in Table 3 of Ref. [12] regarding the background expectation and the observed count in the signal region may be used, thus allowing a more general application of the factorization method.In the case where only one long-lived particle is present in each event, the probabilities have the simplest form P on (k) and P off (m thresh , k).If each event contains two long-lived particles, the probabilities can be expressed using the probabilities for events with either single long-lived particle: ( The expression for events with more than two long-lived particles per event can also be expressed in terms of the probabilities P on (k) and P off (m thresh , k) associated with each individual long-lived particle.
The probabilities P on (k) and P off (m thresh , k), derived using the full simulation as described in the Appendix, are evaluated using the selection requirements adopted by the "tracker+TOF" analysis described in Ref. [12].This selection, where tracks are required to be reconstructed in both the tracker and the muon system, has been shown to be the most sensitive to signatures with lepton-like long-lived particles.The probabilities thus obtained are stored in the form of look-up tables.
The factorization method described above is validated by comparing the estimated signal acceptance values with those obtained with full simulation for a few benchmark models predicting long-lived leptons.In the rest of this paper we refer to the former as "the fast technique".
In the full simulation case, pileup due to multiple interactions per bunch crossing is also simulated.Agreement better than 10% between full simulation and the fast technique presented in this paper is generally observed for the considered values of m thresh .
More details on the determination of the probabilities P on (k) and P off (m thresh , k) can be found in the Appendix, which also contains the details of the validation of the fast technique based on Eq. ( 1).The Appendix also explains how this technique can be used in the future to estimate the CMS exclusion limits for extensions of the SM other than the ones considered in this paper using publicly available look-up tables for P on (k) and P off (m thresh , k).
The probabilities P on (k) and P off (m thresh , k) are computed with simulated stable particles, but the method is easily extended to particles with finite lifetimes by correcting the P on (k) probability for the probability that the long-lived particle, with mass m, lifetime τ, and momentum p, travels at least the distance x required to produce the minimum number of track measurements in the CMS muon system, as required in Ref. [12].The correction consists of an exponential factor, exp[−mx/(τ p)], to be applied to P on (k).This correction is not applied also to P off (m thresh , k), since it needs to be made just once for each long-lived particle.The distance x only depends on the pseudorapidity of the particle: ( These values of x ensure that the particle traverses the entire muon system before decaying.This choice results in a conservative estimate of the signal acceptance since it ignores the contribution to the acceptance from particles that decay before the end of the muon detector but still pass the selection.In Fig. 1, the acceptance obtained with the fast technique is compared with the acceptance obtained with the full simulation of the detector, as a function of the lifetime of the particle.It can be seen from the ratio panels that in most cases the agreement between the two methods is within 10%, corresponding to the systematic uncertainty in the fast estimation technique used to compute the acceptance.For lifetimes less than 10 ns, the spread is somewhat larger, but the tendency to underestimate the acceptance is less than 15%.The offline event selection in Ref. [12] includes two isolation requirements.The first is defined by Σp T < 50 GeV, where the sum is over all tracks (except the candidate track) within a radius ∆R = √ (∆η) 2 + (∆φ) 2 = 0.3 around the candidate track.The second requirement is that E/p < 0.3, where E is the sum of energy deposited in the calorimeters within a radius ∆R = 0.3 around the candidate track (including the candidate energy deposit) and p is the candidate track momentum reconstructed in the tracker.The probabilities P on (k) and P off (m thresh , k) are estimated with single-particle events and thus do not account for the possibility that a longlived particle might fail the isolation requirements.In order to accurately model the isolation requirements, the following procedure, which uses generator-level information from a simulation of the BSM model under test, should be used.The isolation requirements must be determined for each long-lived particle at the generator-level.The following conditions are imposed: The equations ( 4) represent the generator-level equivalents of the tracker and calorimeter isolation requirements, respectively.The sum of the tracks p T around the long-lived candidate is replaced by a sum over transverse momenta of all charged particles around the direction of the long-lived particle.Similarly, the sum of the calorimeter energy deposits around the long-lived candidate is replaced by a sum over the energies of all visible particles, except the long-lived one, around the direction of the long-lived particle of momentum magnitude p. Muons are not considered in the sum over visible particles since they deposit very little energy in the calorimeter.If a long-lived particle does not satisfy both isolation requirements expressed by Eqs.
(4), the corresponding factor P off (m thresh , k) in Eq. ( 2) is set to zero.While pileup effects are not taken into account in the above prescriptions, it was checked that these omissions do not have a significant effect on the results.The robustness of the fast estimate of the acceptance against the hadronic activity surrounding the long-lived particle is assessed by testing different chargino production mechanisms: direct pair production of charginos, pair production of gluinos that each decay to a heavy quark (b,t) and a chargino, pair production of gluinos that each decay to a light quark (u,d,s,c) and a chargino, and pair production of squarks that each decay to a quark and a chargino.A ∼10% agreement between the fast estimate of the model acceptance and the full simulation predictions is observed as shown in Fig. 2. The estimate of the acceptance for chargino pair production is not modified by the isolation requirement since the charginos are always isolated for that production mechanism.
In the two following sections, the fast technique is used to estimate the signal acceptance of the CMS detector for two extensions of the SM.Given that no significant excess of events over the predicted backgrounds is observed in Ref. [12], the signal acceptance estimated with the technique described above is used in calculating cross section limits of these models, at 95% confidence level (CL).The limits in this paper are established using the CL S approach [14,15] where p-values are computed with a hybrid Bayesian-frequentist technique [16].A log-normal model [17,18] is adopted for the nuisance parameters.The nuisance parameters are: the expected background in the signal region, the integrated luminosity, and the signal acceptance.
The expected background and the integrated luminosity, as well as their uncertainties, are provided in Ref. [12].The uncertainty in the signal acceptance is assumed to be 25% for all signal models.This value results from adding 10% to the ≤15% uncertainty reported in Section 6 of Ref. [12].The additional 10% accounts for the systematic uncertainty incurred by the use of the fast estimation technique to compute the acceptance.Table 1 summarizes the information needed, in addition to the signal acceptance evaluated with the fast technique, to set limits on a signal model predicting lepton-like charged long-lived particles.
To demonstrate the validity of this procedure, we have used it to translate the estimate of the signal acceptance to the 95% CL limit on the production cross section for two of the models considered in Ref. [12].The limits obtained on the pair production and inclusive production of staus in the context of the gauge mediated supersymmetry breaking (GMSB) model are shown in Table 2 for signal acceptance values estimated with full simulation and with the procedure described above.The limits obtained with the two techniques agree within 8%.The differences are due only to a small difference in the signal acceptance computed with the two different techniques and to the larger uncertainty assigned to the acceptance from the fast technique.
The results in Table 1 may be compared with those shown in Table 7 and Fig. 8 of Ref. [12].In making the comparisons, it should be noted that there are important differences in the way that  m thresh is chosen in the two cases.In Ref. [12] the value was varied in steps of 10 GeV in order to optimise the resultant limit.In this paper a simpler approach has been followed with m thresh varied in 100 GeV steps and chosen to be the largest value satisfying the condition m thresh ≤ 0.6m.The latter approach generally results in a somewhat higher estimation of the background and therefore a more conservative limit.In the most extreme case, the pair production of staus, with m = 308 GeV, the background is estimated to be 5.6 ± 1.1 events with m thresh = 100 GeV, compared with the estimate of 0.7 ± 0.1 events obtained with m thresh = 190 GeV in Ref. [12], resulting in a cross section limit that is about three times higher in the former case.Nonetheless, the limits agree within ∼15% in almost all cases, allowing restrictive limits to be set on a general class of models.) trigger in combination with a single-muon trigger, there is one caveat to the proposed factorization method: the efficiency of the E miss T trigger cannot be modeled accurately in terms of single long-lived particle kinematic properties.Accounting for the presence of other undetectable particles using a Monte Carlo method would not help because E miss T often depends significantly on detector effects due to the other particles.The assumption that the E miss T trigger adds negligibly to the event selection performed by the muon trigger must therefore be satisfied in order to apply the method to a given signal.Deviations from this assumption would result in an underestimation of the signal acceptance.The assumption is satisfied by models with lepton-like long-lived particles.Models with long-lived colored particles, such as top squarks or gluinos, do not satisfy this condition and thus cannot currently be tested with the technique presented in this paper.Long-lived colored particles hadronize in color singlet bound states [19] that could interact with the detector material leading to complex situations, detailed in Refs.[12,20], and inducing significant instrumental E miss T .For instance, pair-produced colored long-lived particles may hadronize to a charged and a neutral hadron.In this case the E miss T is strongly modified by the presence of the neutral hadron since it is not visible in the tracker and deposits only O(1) GeV in the calorimeter.Moreover, the interactions of the color singlet bound states with matter may lead to a modification of the electric charge of the bound states [19].In this case, the hadron containing a long-lived colored particle can be electrically charged at production, but neutral in the muon system and therefore fail the muon reconstruction.The probabilities P on (k) and P off (m thresh , k) do not account for the possible modification of the electric charge experienced by the hadron-like particles.

Constraints on the pMSSM
The pMSSM is a 19-parameter realization of the minimal supersymmetric model (MSSM) [1] that captures most of the phenomenological features of the R-parity conserving MSSM.The free parameters of the pMSSM, in addition to the SM parameters, are: (1) the gaugino mass parameters M 1 , M 2 , and M 3 ; (2) the ratio of the vacuum expectation values of the two Higgs doublets tan β = v 2 /v 1 ; (3) the higgsino mass parameter µ and the pseudoscalar Higgs mass m A ; (4) the 10 sfermion mass parameters m F , where . . .); and 5) the trilinear couplings A t , A b , and A τ .To minimize the theoretical uncertainties in the Higgs sector, these parameters are defined at the electroweak symmetry breaking scale.
In the pMSSM, all MSSM parameters are allowed to vary freely, subject to the requirement that the model is consistent with some basic constraints: first, the sparticle spectrum must be free of tachyons and cannot lead to color or charge breaking minima in the scalar potential.We also require that the model is consistent with electroweak symmetry breaking and that the Higgs potential is bounded from below.Finally, in this study, we also require the lightest SUSY particle to be the lightest neutralino ( χ 0 1 ).Furthermore, for practical reasons, we limit our study to the pMSSM sub-space chosen to cover sparticle masses up to about 3 TeV.Table 3 summarizes the boundaries of the considered sub-space.
Table 3: The pMSSM parameter space used in the scan.
In Ref. [21], some previous CMS published results were reinterpreted in the context of this pMSSM sub-space.This analysis, however, did not consider the region of parameter space in which long-lived charginos are predicted.Using the technique described in Section 3, we can extend the results of Ref. [21] to regions of the parameter space leading to long-lived particles.
We sample 20 million points in a pMSSM parameter space from a prior probability density function that encodes results from indirect SUSY searches and pre-LHC searches as done in Ref. [21].From this set we select 7205 points with a Higgs boson mass in the range 120 ≤ m h ≤ 130 GeV and predicting long-lived (cτ > 50 cm) charginos.Tightening further the mass window to 123 ≤ m h ≤ 128 GeV to reflect more recent constraints from the Higgs boson mass measurements [22][23][24] reduces the size of the subset by 45% but does not further constrain the chargino mass or lifetime.We therefore use the 120 ≤ m h ≤ 130 GeV window in order to minimize the statistical uncertainty.
For each of the 7205 points in this subset, we have generated 10 000 events using PYTHIA v6.426 [25].Both direct pair production of charginos and indirect chargino production through the decay of heavier SUSY particles were considered for this study.The generated events have been used to evaluate the signal acceptance of the HSCP search, given the chargino kinematic properties predicted by PYTHIA for the considered pMSSM sub-space.
The fast technique is used to obtain acceptance values expected to be in good agreement with the full simulation prediction.The predicted signal acceptance is then used to compute 95% CL limits on the 7205 analysed pMSSM parameter points.A parameter point is excluded if the observed limit obtained on the cross section is less than the theoretical prediction at leading-order as calculated by PYTHIA.The use of leading-order instead of next-to-leading-order theoretical cross section is driven by practical considerations given the large number of parameter points considered.
Figure 3 shows the fraction of parameter points excluded as a function of the chargino lifetime.The fraction of excluded model points with a chargino lifetime longer than 1000 ns (10 ns) is 100.0%(95.9%).Although these values depend on the random point sampling in the 19-dimensional pMSSM parameter space, it is remarkable that a high fraction of the points predicting long-lived charginos are excluded.
Figure 4 shows the number of parameter points predicted and excluded by the analysis of the results obtained in Ref. [12] as a function of the chargino mass, chargino lifetime, and the mass difference between the chargino and the neutralino.

Constraints on the AMSB model
In the AMSB model [2][3][4] the lightest chargino and neutralino are almost degenerate (m , where the neutralino is the lightest SUSY particle.In this model, the chargino lifetime, expected to be of the order of a nanosecond or larger, is determined by the mass splitting with the neutralino.Previous searches for AMSB charginos [26,27] looked for a chargino decaying within the tracking volume into a neutralino and a soft charged pion.The pion has a momentum of ∼ 100 MeV and is generally not reconstructed.The experimental chargino signature therefore takes the form of a disappearing track inside the tracking system.The main limitation of that search technique is that the sensitivity drops quickly as the lifetime increases because the chargino is required to decay within the tracking region.
In contrast, the sensitivity of the search for HSCP [12] is maximal when the chargino decay occurs outside the detector.These two searches are therefore complementary.In this context, the fast technique discussed in the previous sections can be used to assess the limits set by Ref. [12] on long-lived charginos in the AMSB model.
The minimal version of the AMSB model is fully characterized by four parameters: the ratio of Higgs doublet vacuum expectation values at the electroweak scale, the sign of the higgsino mass term, the universal scalar mass (m 0 ), and the gravitino mass (m 3/2 ), which dictates the value of the chargino mass.The values of the first two parameters are set to tan β = 5 and µ > 0. The scalar mass is set to a large value (1 TeV) in order to prevent the appearance of tachyonic sleptons.Gravitino masses ranging from 3.5 to 32 TeV are used in order to scan chargino masses from 100 to 900 GeV.
Samples of simulated charginos with masses from 100 to 900 GeV and lifetimes from 1 ns to 10 µs are produced with PYTHIA v6.426.The SUSY mass spectrum and the decay tables are calculated using ISASUSY with ISAJET v7.80 [28].For each sample, 10 000 events are generated using an inclusive SUSY production and the acceptance of the search for long-lived particles is estimated using the technique described in Section 3. The estimated signal acceptance is shown in Fig. 5 (left).The acceptance is reduced for short lifetimes because the probability that a particle reaches the muon system before decaying, is exponentially smaller.As explained in Section 3, a systematic uncertainty of 25% is assigned to the signal acceptance.A point in the mass-lifetime parameter space of the AMSB model is considered to be excluded when the 95% CL observed limit on the cross section is lower than the leading-order theoretical cross section.The excluded region in this plane is shown in Fig. 5 (right).These results extend those from previous searches at LHC experiments [26,27] by excluding charginos with lifetimes 100 ns up to masses of about 800 GeV.While the signal acceptance remains nearly constant over a wide mass range, heavier charginos cannot be excluded because of their smaller production cross section.The sensitivity of the search for HSCP [12] is limited to charginos with an average lifetime ≥ 3 ns while previous searches based on short track signatures are sensitive to lifetimes of ∼ 0.1 ns [26,27].

Summary
The results of the CMS search for long-lived charged particles have been analysed to set constraints on the phenomenological minimal supersymmetric standard model (pMSSM) and on the anomaly-mediated SUSY breaking model (AMSB), both of which predict the existence of long-lived massive particles in certain regions of their parameter space.A novel technique for estimating the signal acceptance with an accuracy of 10% is presented.This technique only uses generator-level information, while the integrated luminosity, the expected standard model background yields, and the corresponding uncertainties are taken from a previous CMS search [12].The technique and the tabulated probabilities, available as supplementary material to this paper [13], can be used by others to estimate the CMS exclusion limits for different models predicting long-lived lepton-like particles.In the context of the AMSB model, charginos with lifetimes 100 ns (3 ns) and masses up to about 800 GeV (100 GeV) are excluded at 95% confidence level.The most stringent limits to date are set on the long-lived sector of the pMSSM sub-space that covers SUSY particle masses up to about 3 TeV.In this sub-space, 95.9% (100%) of the points with a chargino lifetime τ ≥ 10 ns (1000 ns) are excluded by the present analysis of the results from the CMS search in Ref. [12].These are the first constraints on the pMSSM obtained at the LHC.

A Details of the fast technique
The probabilities P on (k) and P off (m thresh , k), introduced in Section 3, are computed using samples of single long-lived particles, uniformly distributed in η and β, and produced in PYTHIA v6.426 [25].Twenty samples, each containing one million stau particles, were produced for the following long-lived particle masses: 100, 126, 156, 200, 247, 308, 370, 432, 494, 557, 651, 745, 871, 1029, 1200, 1400, 1600, 1800, 2000, and 2500 GeV.These generated events were then processed by the full CMS simulation and reconstruction software.The full simulation includes pileup effects.The selection requirements adopted by the "tracker+TOF" analysis described in Ref. [12] were then applied to the reconstructed events in order to evaluate the probability that a particle with kinematics k can pass the selection criteria.
The probabilities P on (k) and P off (m thresh , k) were evaluated in 3D bins of |η|, β, and p T .For |η| we consider the following bin boundaries: 0.0, 0.25, 0.50, 0.75, 1.00, 1.10, 1.125, 1.50, 1.75, 2.00, 2.10.The granularity is finer around the transition between the tracker barrel and the tracker endcap to better model the efficiency in this region.Since the detector is symmetric in η, the probabilities for the positive and negative η regions were averaged in order to reduce the statistical uncertainties.A constant bin width of 0.05 is considered between 0.0 and 1.0 for the binning in β.For the binning along p T the following bin boundaries are considered 5,50,60,70,80,90,100,110,120,130,140,150,160,170,180,190,200,220,240,260,280,300,325,350,375,400,425,450,475,500,550,600,650,700,750,800,850,900,950,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2250, 2500, ∞ GeV.Figures A.1, A.2, and A.3 show the estimated values of P on (k), P off (0 GeV, k), and P off (300 GeV, k), respectively.The complete evaluation of P on (k) and P off (m thresh , k) for the bin ranges previously described is available as supplemental material to this paper [13].
The fast technique is validated by comparing the estimated signal acceptance values with those obtained from the full simulation for three benchmark models predicting long-lived leptons: pair production of staus, inclusive production of staus in the context of a gauge mediated symmetry breaking (GMSB) model, and pair production of long-lived leptons with a null weak isospin [29].These three models, having significantly different kinematic properties, are studied in Ref. [12].
The signal acceptance obtained with the two methods is shown in Figs.A.4 and A.5 for the benchmark models considered for two thresholds on the reconstructed mass, m thresh .The acceptance computed with the full simulation and with the technique presented in this paper agrees within 10%.The agreement between the two techniques is worse when the mass threshold is close to the true mass of the long-lived particle.The larger disagreement close to the mass threshold is a consequence of the coarse binning of probability P off (m thresh , k), which only partially reflects the sharp variation in probability with the reconstructed values of p T , η, and mass.The reconstructed particle mass is generally lower than the actual particle mass because of a cutoff applied in the data acquisition electronics of the tracker detector.In Ref. [12], the threshold on the reconstructed mass is optimized for each signal model considered.In this paper, for the sake of simplicity, the threshold on the reconstructed mass used in the fast technique is set at 60% of the true particle mass.The hatched area in Figs.A.4 and A.5 indicates the range not satisfying the requirement on the reconstructed mass threshold.Outside of this region the estimated signal acceptance is always compatible within 10% with the results from full simulation.
Several points of the pMSSM subset have also been fully simulated, as described in Ref. [12], in order to further validate the fast estimation of the acceptance in the context of pMSSM.The acceptances estimated using the full simulation prediction and the proposed technique are compared in Fig. A.6 for points leading to quasi-stable (cτ ≥100 m) charginos.An agreement at the level of 10% is observed.

Figure 1 :
Figure 1: Signal acceptance as a function of the chargino lifetime for a benchmark model having a chargino of mass 100 GeV (left) and 700 GeV (right), with a mass threshold of 0 GeV and 300 GeV, respectively.The panel below each figure shows the ratio of acceptance from the fast technique to the acceptance obtained from a full simulation of the detector.

Figure 2 :
Figure2: Signal acceptance as a function of the chargino production mechanism for a benchmark model having a chargino of mass 200 GeV (left) and 700 GeV (right), with a mass threshold of 100 GeV and 300 GeV, respectively.From left to right, the production mechanisms considered are: direct pair production of charginos; pair production of gluinos that each decay to a heavy quark (b,t) and a chargino; pair production of gluinos that each decay to a light quark (u,d,s,c) and a chargino; and pair production of squarks that each decay to a quark and a chargino.The panel below each figure shows the ratio of acceptance from the fast technique with the isolation requirements to the acceptance obtained from a full simulation of the detector.The estimated acceptance is given with and without the generator-level isolation.Pileup is present only in the full simulation samples.

Figure 3 :
Figure 3: (left) Number of pMSSM points, in the sub-space covering sparticle masses up to about 3 TeV, that are excluded at a 95% CL (hatched red) or allowed (white) as a function of the chargino lifetime.(right) Enlargement of the long-lived region.The bottom panel shows the fraction of pMSSM points excluded by the analysis based on the results from the HSCP search [12].

Figure 4 :
Figure 4: Number of pMSSM parameter points in the sub-space covering sparticle masses up to about 3 TeV shown as a function of the chargino mass and (upper row) of the mass difference between the chargino and the neutralino, and (lower row) chargino lifetime.The left panels show the entire set of points considered while the right panels show the set of points excluded by the analysis based on the results from the HSCP search [12].

Figure 5 :
Figure 5: (left) Signal acceptance as a function of chargino mass for the AMSB model as predicted by the fast technique.(right) Observed and expected excluded region on the chargino mass and lifetime parameter space in the context of the AMSB model with tan β = 5 and µ ≥ 0. The excluded region is indicated by the hatched area.

search
Council and EPLANET (European Union); the Leventis Foundation; the A. P. Sloan Foundation; the Alexander von Humboldt Foundation; the Belgian Federal Science Policy Office; the Fonds pour la Formation à la Recherche dans l'Industrie et dans l'Agriculture (FRIA-Belgium); the Agentschap voor Innovatie door Wetenschap en Technologie (IWT-Belgium); the Ministry of Education, Youth and Sports (MEYS) of the Czech Republic; the Council of Science and Industrial Research, India; the HOMING PLUS programme of Foundation for Polish Science, cofinanced from European Union, Regional Development Fund; the Compagnia di San Paolo (Torino); the Consorzio per la Fisica (Trieste); MIUR project 20108T4XTM (Italy); the Thalis and Aristeia programmes cofinanced by EU-ESF and the Greek NSRF; and the National Priorities Research Program by Qatar National Research Fund.

Figure A. 1 :Figure A. 2 :
Figure A.1: Values taken by the probability P on (k) as a function of the true particle-variables p T , β, and |η|.The binning in p T in these figures is coarser than the one used in the study.

Figure A. 3 :
Figure A.3: Values taken by the probability P off (300 GeV, k) as a function of the true particlevariables p T , β, and |η|.The binning in p T in these figures is coarser than the one used in the study.

Figure A. 4 :
Figure A.4: Signal acceptance for a mass threshold of 0 (left) and 300 GeV (right).The upper and lower sets of distributions show the pair production and inclusive production of staus as predicted by the GMSB model, respectively.The panel below each figure shows the ratio of acceptance from the fast technique to the acceptance obtained from a full simulation of the detector.The hatched area indicates the range not satisfying the requirement on the reconstructed mass threshold.

Figure A. 5 :
Figure A.5: Signal acceptance for a mass threshold of 0 (left) and 300 GeV (right).The tested model is the pair production of leptons with no weak isospin.The panel below each figure shows the ratio of acceptance from the fast technique to the acceptance obtained from a full simulation of the detector.The hatched area indicates the range not satisfying the requirement on the reconstructed mass threshold.

Figure A. 6 :
Figure A.6: Signal acceptance for a mass threshold of 0 (left) and 300 GeV (right) on a few representative pMSSM points predicting quasi-stable (cτ ≥ 100 m) charginos.The panel below each figure shows the ratio of acceptance from the fast technique to the acceptance obtained from a full simulation of the detector.The hatched area indicates the range not satisfying the requirement on the reconstructed mass threshold.

Table 1 :
Summary of the information needed to set limits on a signal model predicting leptonlike charged long-lived particles.The mass threshold, the corresponding expected background, and the observed numbers of events, as well as the uncertainty in the signal acceptance evaluated with the fast technique, are provided as a function of the long-lived particle mass.

Table 2 :
Signal acceptance estimated from the fast technique and with the full simulation of the detector, as well as the corresponding expected and observed cross section limits.Results are provided for both the pair production and the inclusive production of staus as predicted by the GMSB model.The mass threshold, the corresponding expected background and the observed numbers of events is also shown.