Search for the electroweak production of supersymmetric particles in √ s = 8 TeV pp collisions with the ATLAS detector

Disclaimer/Complaints regulations If you believe that digital publication of certain material infringes any of your rights or (privacy) interests, please let the Library know, stating your reasons. In case of a legitimate complaint, the Library will make the material inaccessible and/or remove it from the website. Please Ask the Library: https://uba.uva.nl/en/contact, or a letter to: Library of the University of Amsterdam, Secretariat, Singel 425, 1012 WP Amsterdam, The Netherlands. You will be contacted as soon as possible.

The ATLAS experiment has performed extensive searches for the electroweak production of charginos, neutralinos, and staus. This article summarizes and extends the search for electroweak supersymmetry with new analyses targeting scenarios not covered by previously published searches. New searches use vectorboson fusion production, initial-state radiation jets, and low-momentum lepton final states, as well as multivariate analysis techniques to improve the sensitivity to scenarios with small mass splittings and lowproduction cross sections. Results are based on 20 fb −1 of proton-proton collision data at ffiffi ffi s p ¼ 8 TeV recorded with the ATLAS experiment at the Large Hadron Collider. No significant excess beyond Standard Model expectations is observed. The new and existing searches are combined and interpreted in terms of 95% confidence-level exclusion limits in simplified models, where a single production process and decay mode is assumed, as well as within phenomenological supersymmetric models. DOI: 10.1103/PhysRevD.93.052002
The charginos and neutralinos are mixtures of the bino, winos, and higgsinos, collectively referred to as the electroweakinos, that are superpartners of the U(1), SU(2) gauge bosons and the Higgs bosons, respectively. Their mass eigenstates are referred to asχ AE i (i ¼ 1, 2) andχ 0 j (j ¼ 1, 2, 3, 4) in order of increasing mass. The direct production of charginos, neutralinos, and sleptons (l) through electroweak (EW) interactions may dominate the SUSY production at the Large Hadron Collider (LHC) if the masses of the gluinos and squarks are large. Previous searches for electroweak SUSY production at ATLAS targeted the production ofl þl− ,τ þτ− ,χ þ 1χ − 1 (decaying throughl or W bosons),χ AE 1χ 0 2 (decaying throughl or W and Z=h bosons), andχ 0 2χ 0 3 (decaying throughl or Z bosons) [19][20][21][22][23], and found no significant excess beyond SM expectations. These searches are typically sensitive to scenarios where there is a relatively large Oðm W;Z Þ splitting between the produced sparticles and the LSP, leaving uncovered territory for smaller mass splittings.
This article addresses EW SUSY production based on the 20.3 fb −1 of ffiffi ffi s p ¼ 8 TeV proton-proton collisions collected by the ATLAS experiment in 2012. A series of new analyses targeting regions in parameter space not covered by previous ATLAS analyses [19][20][21][22][23] are presented. The results from new and published searches are combined and reinterpreted to provide the final 8 TeV ATLAS limits on the production of EW SUSY particles in a variety of models. The dependence of the limits on the mass of the intermediate slepton in models of electroweakino production withl-mediated decays is also studied, thus generalizing the results of Refs. [19][20][21]. In cases where the LSP is wino or higgsino dominated, the lighter electroweakino statesχ AE 1 ,χ 0 2 can have mass differences with theχ 0 1 ranging from a few GeV to a few tens of GeV, depending on the values of the other parameters in the mixing matrix [24]. In particular, in naturalness-inspired models [25,26] the higgsino must be light, so theχ 0 1 ,χ 0 2 , andχ AE 1 are usually higgsino-dominated and have a small mass splitting. Therefore, a situation with a lightχ 0 1 approximately mass degenerate with theχ AE 1 and χ 0 2 has a strong theoretical motivation. A relatively low mass splitting between the produced sparticles and the LSP (referred to as compressed scenarios) results in lowmomentum decay products that are difficult to reconstruct efficiently, and probing these signatures is experimentally challenging. The new analyses introduced in this article improve the sensitivity to the compressed spectra. The twoand three-lepton searches forχ þ 1χ − 1 andχ AE 1χ 0 2 production in Refs. [19,20] are extended by lowering the transverse momentum threshold on reconstructed leptons, and by boosting the electroweak SUSY system through the requirement of QCD initial state radiation (ISR). The search for the vector-boson fusion (VBF) production of χ AE 1χ AE 1 uses the signature of a same-sign light lepton (e, μ) pair with two jets to probe compressed spectra.
In many SUSY scenarios with large tan β, the stau (τ) is lighter than the selectron and smuon [27], resulting in taurich final states. Coannihilation processes [28] favor a light τ that has a small mass splitting with a bino LSP, as it can set the relic density to the observed value [29]. An additional new search is presented here, which uses a final state with two hadronically decaying τ leptons and multivariate techniques to improve the sensitivity to directτ production compared to the search presented in Ref. [22].
The article is organized as follows: Sec. II describes the signal models studied in this article; Sec. III provides a brief description of the ATLAS detector; Secs. IV and V outline the Monte Carlo (MC) simulation and event selection, respectively; Sec. VI discusses the analysis strategy common to all analyses studied in this article; Sec. VII presents the direct stau production search; Sec. VIII presents the compressed spectra searches in direct production; Sec. IX presents the search for same-sign charginopair production via VBF; Sec. X provides a global overview of the results of the ATLAS searches for electroweakino production at 8 TeV, integrating the results of the new analyses with published analyses in the framework of several relevant signal models; and finally conclusions are drawn in Sec. XI.

II. SUSY SCENARIOS
The SUSY scenarios considered in this article can be divided into two categories: simplified models and phenomenological models. The simplified models [40] target the production of charginos, neutralinos, and sleptons, where the masses and the decay modes of the relevant particles are the only free parameters. In each of the simplified models, a single production process with a fixed decay chain is considered for optimization of the event selection and interpretation of the results. To illustrate the range of applicability of the searches, several classes of phenomenological models that consider all relevant SUSY production and decay processes are also used to interpret the results. These models include the five-dimensional EW phenomenological minimal supersymmetric standard model (pMSSM) [41], the nonuniversal Higgs masses (NUHM) model [42,43], and a gauge-mediated SUSY breaking (GMSB) model [44][45][46][47][48][49].
R-parity is assumed to be conserved in all SUSY scenarios considered in this article. The LSP is assumed to be the lightest neutralinoχ 0 1 except in the GMSB scenarios, where it is the gravitinoG. The next-to-LSP (NLSP) is usually one or more of the charginos, neutralinos, or sleptons. All SUSY particles are assumed to decay promptly, with the exception of the LSP, which is stable. Finally, SUSY particles that are not considered in a given model are decoupled by setting their masses to values inaccessible at the LHC.
Unless stated otherwise, signal cross sections are calculated to next-to-leading order (NLO) in the strong coupling constant using PROSPINO2 [50] and are shown in Fig. 1 for a number of selected simplified-model production modes. The cross sections for the production of charginos and neutralinos are in agreement with the NLO calculations matched to resummation at next-to-leading logarithmic (NLL) accuracy (NLO þ NLL) within about 2% [51][52][53]. The nominal cross section and the uncertainty are taken from the center and spread, respectively, of the envelope of crosssection predictions using different parton distribution function (PDF) sets and factorization and renormalization scales, as described in Ref. [54].
A. Direct stau-pair production simplified model Two simplified models describing the direct production ofτ þτ− are used in this article: one considers stau partners of the left-handed τ lepton (τ L ), and a second considers stau partners of the right-handed τ lepton (τ R ). In both models, the stau decays with a branching fraction of 100% to the SM tau lepton and the LSP. The diagram for this model can be seen in Fig. 2(a).
B. Direct chargino-pair, chargino-neutralino, and neutralino-pair production simplified models In the simplified models describing the direct production ofχ þ 1χ − 1 andχ AE 1χ 0 2 , both theχ AE 1 andχ 0 2 are assumed to be pure wino and mass degenerate, while theχ 0 1 is assumed to be pure bino. However, it is possible to reinterpret the results from these simplified models by assuming different compositions of theχ 0 1 ,χ 0 2 , andχ AE 1 for the same masses of the states. Two different scenarios for the decays of theχ AE 1 andχ 0 2 are considered, as shown in Figs. 2(b) and 2(c): 2 production withl L -mediated decays: Theχ AE 1 andχ 0 2 decay with a branching fraction of 1=6 viaẽ L ,μ L ,τ L ,ν e ,ν μ , orν τ with masses mν l ¼ ml 25, 0.5, 0.75, or 0.95, (ii)χ AE 1χ 0 2 production withτ L -mediated decay: The firstand second-generation sleptons and sneutrinos are assumed to be very heavy, so that theχ AE 1 andχ 0 2 decay with a branching fraction of 1=2 viaτ L orν τ with masses mν τ ¼ mτ L ¼ 0.5ðm~χAE 1 þ m~χ0 1 Þ. In the simplified models considered here, the slepton mass is assumed to lie between theχ 0 1 andχ AE 1 =χ 0 2 masses, which increases the branching fraction to leptonic final states compared to scenarios without sleptons. The compressed spectra searches in this article are less sensitive to scenarios where theχ AE 1 =χ 0 2 decay through SM W, Z, or Higgs bosons, as the branching fraction to leptonic final states is significantly suppressed. The results of the ATLAS searches forχ þ 1χ − 1 production with WW-mediated decays [19],χ AE 1χ 0 2 production with WZ-mediated decays [20], andχ AE 1χ 0 2 production with Wh-mediated decays [23] are summarized in Sec. X E. In these scenarios with decays mediated by SM bosons, the W, Z, and h bosons are assumed to decay with SM branching fractions.
In the simplified models of the direct production ofχ 0 2χ 0 3 , theχ 0 2 andχ 0 3 are assumed to be pure higgsino and mass degenerate, while theχ 0 1 is assumed to be pure bino. Theχ 0 2 andχ 0 3 are assumed to decay with a branching fraction of one-half viaẽ R ,μ R with mass ml 0.25, 0.5, 0.75, or 0.95 (χ 0 2χ 0 3 production withl R -mediated decay). The associated diagram is shown in Fig. 2(d). In thisχ 0 2χ 0 3 simplified model, the choice of right-handed sleptons in the decay chain ensures high lepton multiplicities in the final state while suppressing the leptonic branching fraction of any associated chargino, thus enhancing the rate of four-lepton events with respect to events with lower lepton multiplicities.

C. Simplified model of same-sign chargino-pair production via vector-boson fusion
A simplified model forχ AE 1χ AE 1 production via VBF [55,56] is also considered. As in the case of direct production, theχ AE 1 is assumed to be pure wino, and mass degenerate with theχ 0 2 , and theχ 0 1 is assumed to be pure bino. Theχ AE 1 decays with a branching fraction of 1=6 viã e L ,μ L ,τ L ,ν e ,ν μ , orν τ with masses mν l ¼ ml 1 production via VBF, where the sparticles are produced along with two jets, is shown in Fig. 2(e). The jets are widely separated in pseudorapidity 1 η and have a relatively high dijet invariant mass m jj . Because of the VBF topology, the charginos are often boosted in the transverse plane, forcing the decay products to be more collinear and energetic, even in highly compressed spectra. This feature of VBF production makes it a good candidate to probe compressed SUSY scenarios that are experimentally difficult to explore via the direct production modes. The signal cross sections are calculated to leading order (LO) in the strong coupling constant using  (more details on the cross-section calculation are given in the Appendix).
The uncertainties on the signal cross sections are calculated by using different PDF sets (2%) and by varying the renormalization and factorization scales between 0.5 and 2 times the nominal values (6%) [58]. For aχ AE 1 with mass of 120 GeV, the cross section forχ AE 1χ AE 1 production in association with two jets satisfying the criteria m jj > 350 GeV and jΔη jj j > 1.6 is 1.1 fb. For the assumed mixings in the chargino-neutralino sector, and the mass values considered in the analysis, the cross section for χ AE 1χ AE 1 VBF production is found to be independent of thẽ χ 0 1 mass.

D. Phenomenological minimal supersymmetric standard model
The analysis results are interpreted in a pMSSM scenario. The masses of the sfermions, the gluino, and the CPodd Higgs boson are set to high values (2 TeV, 2 TeV, and 500 GeV, respectively), thus decoupling the production of these particles and allowing only the direct production of charginos and neutralinos decaying via SM gauge bosons and the lightest Higgs boson. The remaining four parameters, the ratio of the expectation values of the two Higgs doublets (tan β), the gaugino mass parameters M 1 and M 2 , and the higgsino mass parameter μ, determine the phenomenology of direct electroweak SUSY production. For the analysis presented here, μ and M 2 are treated as free parameters. The remaining parameters are fixed to tan β ¼ 10 and M 1 ¼ 50 GeV, so that the relic dark-matter density is below the cosmological bound [29] across most of the μ-M 2 grid. The lightest Higgs boson has a mass close to 125 GeV, which is set by tuning the mixing in the top squark sector, and decays to SUSY as well as SM particles where kinematically allowed.

E. Two-parameter nonuniversal Higgs masses model
Radiatively driven natural SUSY [59] allows the Z and Higgs boson masses to be close to 100 GeV, with gluino and squark masses beyond the TeV scale. In the twoparameter NUHM model (NUHM2) that is considered in this article, the direct production of charginos and neutralinos is dominant in a large area of the parameter space considered. The mass hierarchy, composition, and production cross section of the SUSY particles are governed by the universal soft SUSY-breaking scalar mass m 0 , the soft SUSY-breaking gaugino mass m 1=2 , the trilinear SUSYbreaking parameter A 0 , the pseudoscalar Higgs boson mass m A , tan β, and μ. Both μ and m 1=2 are treated as free parameters, and the other parameters are fixed to and signðμÞ > 0. These conditions ensure a low level of electroweak fine-tuning, while keeping the lightest Higgs boson mass close to 125 GeVand the squark masses to a few TeV. The gluino mass typically satisfies mg ≃ 2.5m 1=2 . For low gluino masses, the production of strongly interacting 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 along 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 z axis. The pseudorapidity is defined in terms of the polar angle θ as η ¼ − ln tanðθ=2Þ. SUSY particles dominates; as the gluino mass increases, the production of electroweakinos becomes more important. The charginos and neutralinos decay via W, Z, and Higgs bosons.

F. Gauge-mediated SUSY breaking model
Minimal GMSB models are described by six parameters: the SUSY-breaking mass scale in the low-energy sector (Λ), the messenger mass (M mess ), the number of SU(5) messenger fields (N 5 ), the scale factor for the gravitino mass (C grav ), tan β, and μ. In the model presented here, Λ and tan β are treated as free parameters, and the remaining parameters are fixed to M mess ¼ 250 TeV, N 5 ¼ 3, C grav ¼ 1, and signðμÞ > 0. For high Λ values, the EW production of SUSY particles dominates over other SUSY processes. In most of the relevant parameter space, the NLSP is theτ for large values of tan β (tan β > 20), and the final states contain two, three, or four tau leptons. In the region where the mass difference between the stau and selectron/smuon is smaller than the sum of the tau and the electron/muon masses, the stau, selectron, and smuon decay directly into the LSP and a lepton, defining the phenomenology. The charginos and neutralinos decay as and the LSP is the gravitinoG.

III. THE ATLAS DETECTOR
The ATLAS detector [60] is a multipurpose particle physics detector with forward-backward symmetric cylindrical geometry. The inner tracking detector (ID) covers jηj < 2.5 and consists of a silicon pixel detector, a semiconductor microstrip detector, and a transition radiation tracker. The ID is surrounded by a thin superconducting solenoid providing a 2 T axial magnetic field. A highgranularity lead/liquid-argon sampling calorimeter measures the energy and the position of electromagnetic showers within jηj < 3.2. Sampling calorimeters with liquid argon as the active medium are also used to measure hadronic showers in the end cap (1.5 < jηj < 3.2) and forward (3.1 < jηj < 4.9) regions, while a steel/scintillator tile calorimeter measures hadronic showers in the central region (jηj < 1.7). The muon spectrometer (MS) surrounds the calorimeters and consists of three large superconducting air-core toroid magnets, each with eight coils, a system of precision tracking chambers (jηj < 2.7), and fast trigger chambers (jηj < 2.4). A three-level trigger system [61] selects events to be recorded for off-line analysis.

IV. MONTE CARLO SIMULATION
Monte Carlo generators are used to simulate SM processes and new physics signals. The SM processes considered are those that can lead to leptonic signatures. Details of the signal and background MC simulation samples used in this article, as well as the order of cross-section calculations in perturbative QCD used for yield normalization, are shown in Table I. For all MC simulation samples, the propagation of particles through the ATLAS detector is modeled with GEANT 4 [96] using the full ATLAS detector simulation [97], or a fast simulation using a parametric response of the electromagnetic and hadronic calorimeters [98] and GEANT 4 elsewhere. The effect of multiple proton-proton collisions in the same or nearby beam bunch crossings (in-time and out-of-time pileup) is incorporated into the simulation by overlaying additional minimum-bias events generated with PYTHIA -8 onto hard-scatter events. Simulated events are weighted to match the distribution of the mean number of interactions per bunch crossing in data and are reconstructed in the same manner as data. The simulated MC samples are corrected to account for differences with respect to the data in the heavy-flavor quark jet selection efficiencies and misidentification probabilities, lepton efficiencies, tau misidentification probabilities, as well as the energy and momentum measurements of leptons and jets. 0 2 ) signal samples simulated with Herwig++ are reweighted to match theχ þ 1χ − 1 (χ AE 1χ 0 2 ) system transverse momentum distribution obtained from the MADGRAPH samples that are generated with an additional parton in the matrix element to give a better description of the ISR.

V. EVENT RECONSTRUCTION
Events recorded during stable data-taking conditions are analyzed if the reconstructed primary vertex has five or more tracks with transverse momentum p T > 400 GeV associated with it. The primary vertex of an event is identified as the vertex with the highest Σp T 2 of associated tracks. After the application of beam, detector, and dataquality requirements, the total luminosity considered in these analyses corresponds to 20.3 fb −1 (20.1 fb −1 for the direct stau production analysis due to a different trigger requirement).
Electron candidates are required to have jηj < 2.47 and p T > 7 GeV, where the p T and η are determined from the calibrated clustered energy deposits in the electromagnetic calorimeter and the matched ID track, respectively. Electrons must satisfy "medium" identification criteria, following Ref. [99]. Muon candidates are reconstructed by combining tracks in the ID and tracks in the MS [100], and are required to have jηj < 2.5 and p T > 5 GeV. Events containing one or more muons that have transverse impact parameter with respect to the primary vertex jd 0 j > 0.2 mm or longitudinal impact parameter with respect to the primary vertex jz 0 j > 1 mm are rejected to suppress cosmic-ray muon background. In the direct stau production analysis, and the two-lepton compressed spectra analyses, electrons and muons are required to have p T > 10 GeV.
Jets are reconstructed with the anti-k t algorithm [101] with a radius parameter of R ¼ 0.4. Three-dimensional calorimeter energy clusters are used as input to the jet reconstruction. The clusters are calibrated using the local hadronic calibration [102], which gives different weights to the energy deposits from the electromagnetic and hadronic components of the showers. The final jet energy calibration corrects the calorimeter response to the particle-level jet energy [102,103], where correction factors are obtained from simulation and then refined and validated using data.
Corrections for in-time and out-of-time pileup are also applied based on the jet area method [102]. Central jets must have jηj < 2.4 and p T > 20 GeV, and a "jet vertex fraction" (JVF) [102] larger than 0.5 if p T < 50 GeV. The JVF is the p T -weighted fraction of the tracks in the jet that are associated with the primary vertex. Requiring large JVF values suppresses jets from pileup. Forward jets are those with 2.4 < jηj < 4.5 and p T > 30 GeV. Events containing jets failing to satisfy the quality criteria described in Ref. [102] are rejected to suppress events with large calorimeter noise and noncollision backgrounds. Central jets are identified as containing b-hadrons (referred to as b-tagged) using a multivariate technique based on quantities related to reconstructed secondary vertices. The chosen working point of the b-tagging algorithm [104] correctly identifies b-hadrons in simulated tt samples with an efficiency of 80%, with a light-flavor jet misidentification probability of about 4% and a c-jet misidentification probability of about 30%.
Hadronically decaying τ leptons (τ had ) are reconstructed using jets described above with jηj < 2.47 and a lower p T threshold of 10 GeV. The τ had reconstruction algorithm uses information about the tracks within ΔR ≡ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi ðΔϕÞ 2 þ ðΔηÞ 2 p ¼ 0.2 of the seed jet, in addition to the electromagnetic and hadronic shower shapes in the calorimeters. The τ had candidates are required to have one or three associated tracks (prongs), as τ leptons predominantly decay to either one or three charged pions together with a neutrino and often additional neutral pions. The τ had candidates are required to have p T > 20 GeV and unit total charge of their constituent tracks. A boosted decision tree algorithm (BDT) uses discriminating track and cluster variables to optimize τ had identification, where "loose," "medium," and "tight" working points are defined [105]. Electrons misidentified as τ had candidates are vetoed using transition radiation and calorimeter information. The τ had candidates are corrected to the τ energy scale [105] using an ηand p T -dependent calibration. Kinematic variables built using taus in this article use only the visible decay products from the hadronically decaying tau.
The missing transverse momentum is the negative vector sum of the transverse momenta of all muons with p T > 10 GeV, electrons with p T > 10 GeV, photons with p T > 10 GeV [99], jets with p T > 20 GeV, and calibrated calorimeter energy clusters with jηj < 4.9 not associated with these objects. Hadronically decaying τ leptons are included in the E miss T calculation as jets. Clusters associated with electrons, photons, and jets are calibrated to the scale of the corresponding objects. Calorimeter energy clusters not associated with these objects are calibrated using both calorimeter and tracker information [106]. For jets, the calibration includes the pileup correction described above, while the JVF requirement is not considered when selecting jet candidates.
To avoid potential ambiguities among objects, "tagged" leptons are candidate leptons separated from each other and from jets in the following order: (1) If two electron candidates are reconstructed with ΔR < 0.1, the lower energy candidate is discarded. (2) Jets within ΔR ¼ 0.2 of an electron candidate, and τ had candidates within ΔR ¼ 0.2 of an electron or muon, are discarded. Tagged leptons satisfying additional identification criteria are called signal leptons. To maximize the search sensitivity, some analyses presented in this article require different additional criteria for signal leptons, and these are highlighted where necessary. Signal τ leptons must satisfy medium identification criteria [105], while for the final signal-region selections, both the medium and tight criteria are used. Unless stated otherwise, signal electrons (muons) are tagged electrons (muons) for which the scalar sum of the transverse momenta of tracks within a cone of ΔR ¼ 0.3 around the lepton candidate is less than 16% (12%) of the lepton p T . Tracks used for the electron (muon) isolation requirement defined above are those that have p T > 0.4 (1.0) GeV and jz 0 j < 2 mm with respect to the primary vertex of the event. Tracks of the leptons themselves as well as tracks closer in z 0 to another vertex (that is not the primary vertex) are not included. The isolation requirements are imposed to reduce the contributions from semileptonic decays of hadrons and jets misidentified as leptons. Signal electrons must also satisfy tight identification criteria [99], and the sum of the extra transverse energy deposits in the calorimeter (corrected for pileup effects) within a cone of ΔR ¼ 0.3 around the electron candidate must be less than 18% of the electron p T . To further suppress electrons and muons originating from secondary vertices, the d 0 normalized to its uncertainty is required to be small, with jd 0 j=σðd 0 Þ < 5ð3Þ, and jz 0 sin θj < 0.4 mm (1 mm) for electrons (muons).
Events must satisfy the relevant trigger for the analysis and satisfy the corresponding p T -threshold requirements shown in Table II.

VI. GENERAL ANALYSIS STRATEGY
The broad range of EW SUSY scenarios considered by the ATLAS experiment is accompanied by a large number of experimental signatures: from the two-tau signature from direct stau production to three-lepton signatures fromχ AE 1χ 0 2 production. As much as possible the individual analyses follow a common approach. Signal regions (SR) are defined to target one or more EW SUSY scenarios, using kinematic variables with good signal-background separation, as described in Sec. VI A. The optimization of key selection variables is performed by maximizing the expected sensitivity to the signal model. A common background estimation strategy is used for the analyses in this article: the main SM backgrounds are estimated by normalizing MC simulation samples to data in dedicated control regions (CRs); backgrounds due to nonprompt and fake leptons are derived from data as outlined in Sec. VI B, while small backgrounds are estimated purely using MC simulation samples. The HISTFITTER [107] software framework is used in all analyses for constraining the background normalizations and the statistical interpretation of the results. The CRs are defined with kinematic properties similar to the SRs, yet are disjoint from the SR, and have high purity for the background process under consideration. The CRs are designed in a way that minimizes the contamination from the signal model, and cross contamination between multiple CRs is taken into account in the normalization to data. To validate the modeling of the SM backgrounds, the yields and shapes of key kinematic variables are compared to data in validation regions (VR). The VRs are defined to be close to, yet disjoint from the SR and CR, and be dominated by the background process under consideration. The VRs are designed such that the contamination from the signal model is low. Three different fit configurations are used. The "background-only fit" is used for estimating the expected background in the SRs and VRs using observations in the CRs, with no assumptions made on any signal model. In the absence of an observed excess of events in one or more signal regions, the "model-dependent signal fit" is used to set exclusion limits in a particular model, where the signal contribution from the particular model that is being tested is taken into account in all CR and SR. Finally, in the "model-independent signal fit," both the CRs and SRs are used in the same manner as for the modeldependent signal fit, but signal contamination is not accounted for in the CRs. A likelihood function is built as the product of Poisson probability functions, describing the observed and expected number of events in the CRs and SRs. The observed number of events in various CRs and SRs are used in a combined profile likelihood fit to determine the expected SM background yields in each of the SRs. The systematic uncertainties on the expected background yields described in Sec. VI C are included as nuisance parameters, constrained to be Gaussian with a width determined by the size of the uncertainty. Correlations between control and signal regions, and background processes, are taken into account with common nuisance parameters. The free parameters and the nuisance parameters are determined by maximizing the product of the Poisson probability functions and the Gaussian constraints on the nuisance parameters.
After the background modeling is understood and validated, the predicted background in the SR is compared to the observed data. In order to quantify the probability for the background-only hypothesis to fluctuate to the observed number of events or higher, the one-sided p 0value is calculated. For this calculation, the profile likelihood ratio is used as a test statistic to exclude the signalplus-background hypothesis if no significant excess is observed. A signal model can be excluded at 95% CL if the CL s [108] of the signal plus background hypothesis is < 0.05. For each signal region, the expected and observed upper limits at 95% CL on the number of beyond-the-SM TABLE II. The triggers used in the analyses and the off-line p T threshold used, ensuring that the lepton(s) or E miss T triggering the event are in the plateau region of the trigger efficiency. Where multiple triggers are listed for an analysis, events are used if any of the triggers is passed. Muons are triggered within a restricted range of jηj < 2.4.

Trigger p T threshold [GeV] Analysis
Single τ 150 Direct stau production Double τ 40,25 Single isolated e 25 Compressed spectra l þ l − , 3l Single isolated μ 25 Double e 14, 14 25,10 Compressed spectra l þ l − , l AE l AE , 3l Double μ 14, 14 18,10 Compressed spectra l þ l − , l AE l AE , 3l Triple e 20,9,9 Compressed spectra 3l Triple μ 7,7, 7 19,5,5 Compressed spectra 3l Combined eμ 14ðeÞ,10ðμÞ 18ðμÞ,10ðeÞ 9ðeÞ, 9ðeÞ, 7ðμÞ 9ðeÞ, 7ðμÞ, 7ðμÞ Compressed spectra 3l E miss events (S 95 exp and S 95 obs ) are calculated using the model-independent signal fit. The 95% CL upper limits on the signal cross section times efficiency (hϵσi 95 obs ) and the CL b value for the background-only hypothesis are also calculated for each analysis in this article. : In the three-lepton analysis, the lepton not forming the SFOS lepton pair with mass closest to the Z boson mass is used. In cases where the second lepton or tau is used, the variable is labeled as m X T , where X is the object used with the E miss T to form the transverse mass.

m SFOS
The invariant mass of the SFOS lepton pair in the event. In the three-lepton analysis, the SFOS pair with mass closest to the Z boson mass is used. The two-tau invariant mass.

m T2
The "stransverse mass" is calculated as where l1=τ1 and l2=τ2 denote the highest-and second-highest-p T leptons or taus in the event, respectively, andq T is a test transverse vector that minimizes the larger of the two transverse masses m T . The m T2 distribution has a kinematic end point for events where two massive pair-produced particles each decay to two particles, one of which is detected and the other escapes undetected [109,110].
m eff The scalar sum of the transverse momenta of the signal leptons, taus, jets, and E miss T in the event In the case of the two-tau analysis, only the sum of the E miss T and two taus is used.
The quantity R 2 is defined as The R 2 distribution is shifted toward unity for signal events compared to the background, due to the existence of the LSPs that results in a larger E miss T .
M R Δ , Δϕ β R The super-razor quantities M R Δ and Δϕ β R are defined in Ref. [111]. These variables are motivated by the generic process of the pair production of two massive particles, each decaying into a set of visible and invisible particles (i.e.χ AE 1 → lν lχ 0 1 ). Similar to m T2 , M R Δ is sensitive to the squared mass difference of the pair-produced massive particle and the invisible particle, via a kinematic end point. These two variables are expected to provide a similar performance for discriminating the signal from the background. For systems where the invisible particle has a mass that is comparable to the pair-produced massive particle (i.e. compressed spectra), the variable Δϕ β R has a pronounced peak near π. The effect is magnified as the spectrum becomes more and more compressed, making this variable a good discriminator for compressed spectra searches.

B. Common reducible background estimation
Electron and muon candidates can be classified into three main types, depending on their origin: "real" leptons are prompt and isolated leptons from a W or Z boson, a prompt tau, or a SUSY particle decay; "fake" leptons can originate from a misidentified light-flavor quark or gluon jet (referred to as "light flavor"); "nonprompt" leptons can originate from a semileptonic decay of a heavy-flavor quark, from the decay of a meson, or from an electron from a photon conversion. The background due to nonprompt and fake electrons and muons, collectively referred to as "reducible," is commonly estimated using the matrix method described in Ref. [112]. The matrix method extracts the number of events with one or two fake or nonprompt leptons from a system of linear equations relating the number of events with two signal or tagged leptons (before signal lepton identification requirements are applied) to the number of events with two candidates that are either real, fake, or nonprompt. The coefficients of the linear equations are functions of the real-lepton identification efficiencies and of the fake and nonprompt lepton misidentification probabilities, both defined as a fraction of the corresponding tagged leptons satisfying the signal lepton requirements.
The real-lepton identification efficiencies are obtained from MC simulation samples in the region under consideration to account for detailed kinematic dependencies and are multiplied by correction factors to account for residual differences with respect to the data. The correction factors are obtained from a control region rich in Z → e þ e − and Z → μ þ μ − decays. The fake and nonprompt lepton misidentification probabilities are calculated as the weighted averages of the corrected type-and process-dependent misidentification probabilities defined below according to their relative contributions in a given signal or validation region. The type-and process-dependent misidentification probabilities for each relevant fake and nonprompt lepton type (heavy-flavor, light-flavor, or conversion) and for each reducible background process are corrected using the ratio ("correction factor") of the misidentification probability in data to that in simulation obtained from dedicated control samples. The correction factors are assumed to be independent of the selected regions and of any potential composition or kinematic differences. For nonprompt electrons and muons from heavy-flavor quark decays, the correction factor is measured in a bb-dominated control sample. The correction factor for the conversion candidates is determined in events with a converted photon radiated from a muon in Z → μμ decays.

C. Common systematic uncertainties
Several sources of systematic uncertainty are considered for the SM background estimates and signal yield predictions. When the MC simulation samples are normalized to data yields in the CR, there is a partial cancellation of both the experimental and theoretical modeling systematic uncertainties.
The experimental systematic uncertainties affecting the simulation-based estimates include the following: the uncertainties due to the jet energy scale and resolution [100,102]; the uncertainties due to the lepton energy scale, energy resolution, and identification efficiency [99,100,105]; the uncertainty due to the hadronic tau misidentification probability [105]; the uncertainty on the E miss T from energy deposits not associated with reconstructed objects (E miss T soft-term resolution) [106]; and the uncertainties due to b-tagging efficiency and mistag probability [104]. The uncertainty on the integrated luminosity is AE2.8% and is derived following the same methodology as that detailed in Ref. [113]. The uncertainty due to the modeling of the pileup in the MC simulation samples is estimated by varying the distribution of the number of interactions per bunch crossing overlaid in the MC samples by AE10%. An uncertainty is applied to MC samples to cover differences in efficiency observed between the trigger in data and the MC trigger simulation.
The systematic uncertainties due to the limitations in theoretical models or calculations affecting the simulationbased background estimates include the cross-section uncertainties that are estimated by varying the renormalization and factorization scales and the PDFs, and the acceptance uncertainties due to PDFs and the choice of MC generator and parton shower. The cross-section uncertainties for the irreducible backgrounds used here are 30% for ttV [76,77], 50% for tZ, 5% for ZZ, 7% for WZ, and 100% for the triboson samples. For the Higgs boson samples, a 20% uncertainty is used for VH and VBF production, while a 100% uncertainty is assigned to ttH and Higgs boson production via gluon fusion [74]. For theχ þ signal simulations that are sensitive to ISR, the impact of the choice of renormalization scales, factorization scales, the scale for the first emission in the so-called MLM matching scheme [114], and MLM matching scale are evaluated by varying these individually between 0.5 and 2 times the nominal values in MadGraph.

VII. DIRECT STAU PRODUCTION
This section presents a search for direct stau-pair production with subsequent decay into final states with two taus and E miss T . The search for direct stau production is very challenging, as the final state is difficult to trigger on and to separate from the SM background. In Ref. [22], the best observed upper limit on the direct stau production cross section was found for a stau mass of 80 GeV and a masslessχ 0 1 , where the theoretical cross section at NLO is 0.07 (0.17) pb for right-handed (left-handed) stau-pair production and the excluded cross section is 0. 22 (0.28) pb. This analysis is an update of Ref. [22], using a multivariate analysis technique instead of a simple cut-based method to improve the sensitivity to direct stau-pair production.

A. Event selection
Events are selected using the basic reconstruction, object, and event selection criteria described in Sec. V. In addition, if taus form an SFOS pair with m SFOS < 12 GeV, the event is rejected. Events with exactly two hadronically decaying tau candidates are selected, where the two tau candidates are required to have opposite-sign (OS) charge. At least one tau must satisfy the tight tau identification BDT requirement, and events with additional tagged light leptons are vetoed. Events must satisfy either the single-tau or ditau trigger criteria, as described in Sec. V.
To suppress events from Z boson decays, events are rejected if the invariant mass of the tau pair lies within AE10 GeV of the peak value of 81 GeV for Z boson candidates. 2 To suppress background from events containing a top quark, events with b-tagged jets are vetoed. To further select SUSY events from direct stau production and suppress WW and tt production, m T2 is calculated using the two taus and the E miss T in the event. The additional requirement of m T2 > 30 GeV is applied to select events for the training and optimization of the multivariate analysis (MVA).
After applying the preselection listed above, both the signal and background MC samples are split in two. Half is used for the BDT training and the other half for testing. Twelve variables with good discriminatory power are considered as input for the BDT training procedure: E miss T , m eff , m T2 , m ττ , Δϕðτ; τÞ, Δηðτ; τÞ, p τ1 T , p τ2 T , m Tτ1 , m Tτ2 , ΔϕðE miss T ; τ1Þ, and ΔϕðE miss T ; τ2Þ. The MC simulation samples are compared to data for these variables and their correlations to ensure that they are modeled well.
A direct stau production scenario with mðτ R ;χ 0 1 Þ ¼ ð109; 0Þ GeV is used for the training and optimization of the BDT, and the BDT response requirement (t cut ) is chosen based on the best expected sensitivity for discovery. The two-tau MVA SR definition is shown in Table III.

B. Background determination
The main SM backgrounds in the two-tau MVA SR are W þ jets and diboson production. Contributions from diboson, tt, and Z þ jets processes are estimated using MC simulation samples and validated using data in WWrich, tt-rich, or Z-rich validation regions, as defined in Ref. [22]. 2 The Z boson mass in ditau decays is reconstructed lower than the Z boson mass value due to the neutrinos from the tau decay.
The W þ jets contribution in the signal region is dominated by events where the W decays to a tau lepton and a jet is misidentified as another tau. The contribution is estimated by normalizing the yields from MC simulation samples to data in a dedicated control region. The W þ jets control region selects events with the W boson decaying to a muon and neutrino to suppress the multijet background, which is larger for the electron channel. Events containing exactly one isolated muon and one tau satisfying the tight identification requirement are selected, where the muon and tau must have opposite electrical charge. To reduce the contribution from Z þ jets production, m τ T þ m μ T > 80 GeV is required, and the reconstructed invariant mass of the muon and tau must be outside the Z mass window (12 GeV < m τμ < 40 GeV or m τμ > 100 GeV). To further suppress multijet and Z þ jets processes, E miss T > 40 GeV is required, and the muon and tau must not be back-to-back [Δϕðτ; μÞ < 2.7 and Δηðτ; μÞ < 2.0]. The contribution from events with top quarks is suppressed by rejecting events containing btagged jets. The multijet background in the W þ jets control region is estimated using a region with the same requirements, but with a same-sign muon and tau. The contribution from other SM processes is subtracted using MC simulation samples, and the ratio of opposite-sign muon and tau events to same-sign events is assumed to be unity for the multijet background.
The contribution from multijet events in the signal region, where both selected taus are misidentified jets, is small and is estimated using the so-called ABCD method. Four exclusive regions (A, B, C, D) are defined in a twodimensional plane as a function of the two uncorrelated TABLE IV. Numbers of events observed in data and expected from SM processes and the SUSY reference point mðτ R ;χ 0 1 Þ ¼ ð109; 0Þ GeV in the two-tau MVA validation and signal regions. The uncertainties shown include both statistical and systematic components. The "top" contribution includes the single top, tt, and ttV processes. The multijet background estimation is taken from data, as described in the text. In the VR, the multijet scale factor from fitting the background is not applied, while the W þ jets scale factor is applied. In the SR, both the multijet and the W þ jets scale factors are applied. Also shown are the model-independent limits calculated from the signal region observations: the one-sided p 0 -value; the expected and observed upper limits at 95% CL on the number of beyond-the-SM events (S 95 exp and S 95 obs ) for each signal region, calculated using pseudoexperiments and the CL s prescription; the observed 95% CL upper limit on the signal cross section times efficiency (hϵσi 95 obs ); and the CL b value for the background-only hypothesis.

SM process
discriminating variables m T2 and the tau identification criterion. Regions A and B are required to have two medium taus where at least one meets the tight tau identification criteria, while regions C and D are required to have two loose taus that fail to satisfy the tight tau identification criteria. In regions A and C (B and D) m T2 > 30 GeV (m T2 < 20 GeV) is also required. The multijet background in signal region A can be estimated from and N D are the numbers of events in regions A, B, C, and D, respectively.
The assumption that the ratios N A =N C and N B =N D are the same is confirmed using MC simulation samples and in validation regions using data.
A simultaneous likelihood fit to the multijet estimation and W þ jets CR is performed to normalize the corresponding background estimates and obtain the expected yields in the SR (as described in Sec. VI). After the simultaneous fit, the multijet and W þ jets normalization factors are found to be 1.4 þ2.5 −1.4 and 0.98 AE 0.30, respectively. Because of the small number of events in some of the ABCD regions, the uncertainty on the multijet normalization factor is large; however, the multijet contribution to the total background is very small and the effect on the total signal region background uncertainty is small.
Two multijet validation regions are defined with the same selection as for the signal region, but with t cut < 0.07 and intermediate m T2 . These multijet validation regions are enriched in events with jets misidentified as hadronic tau decays, and good agreement is seen between the data and expectation across the BDT input kinematic variables. A further two validation regions are defined to check the modeling of the W þ jets background. The intermediate BDT region −0.2 < t cut < 0.07 is used, with a high E miss T selection, where the W þ jets background is seen to be  Table III. Table IV and Figs. 3(a), 3(b), 3(c), and 3(d) show the agreement between data and expectation in the validation regions. The purity of the multijet and W þ jets validation regions is ∼90% and ∼50%, respectively, while the signal contamination from the mðτ R ;χ 0 1 Þ ¼ ð109;0Þ GeV scenario is < 1% and < 10%, respectively.

C. Results
The observed number of events in the signal region is shown in Table IV along with the background expectations, uncertainties, p 0 -value, S 95 exp , S 95 obs , hϵσi 95 obs , and the CL b value. The individual sources of uncertainty on the background estimation in the SR are shown in Table V, where the dominant sources are the statistical uncertainty on the MC simulation samples, the uncertainty on the E miss T from energy deposits not associated with reconstructed objects, and the statistical uncertainty on the normalization factor applied to the W þ jets background. Generator modeling uncertainties for the W þ jets background are estimated by varying the renormalization and factorization scales individually between 0.5 and 2 times the nominal values in Alpgen. Additionally, the impact of the jet p T threshold used for parton-jet matching in Alpgen W þ jets simulation is assessed by changing the jet p T threshold from 15 GeV to 25 GeV. Figures 4(a), 4(b), 4(c), and 4(d) show the distributions of the BDT response prior to the t cut selection, and the E miss T , m eff , and m T2 quantities in the SR, where good agreement between the expected background and the observed data is seen.  In many SUSY scenarios, one or more of the mass differences between the charginos and neutralinos is small, resulting in final states with low-momentum leptons that require dedicated searches. The two-lepton analysis in Ref. [19] excludedχ þ 1χ − 1 scenarios withl L -mediated decays withχ AE 1 -χ 0 1 mass splittings down to approximately 100 GeV, while the three-lepton analysis in Ref. [20] excludedχ AE 1χ 0 2 scenarios withl L -mediated decays down toχ 0 2 -χ 0 1 mass splittings of 20 GeV. The analyses presented in this section focus on event selections based on lowmomentum leptons, and also on the production in association with ISR jets to provide improved sensitivity to the compressed spectra scenarios not covered by previous searches. As discussed in Sec. II, simplified models describingχ þ 1χ − 1 andχ AE 1χ 0 2 production are considered for these compressed spectra searches, where theχ AE 1 =χ 0 2 decay only through sleptons or sneutrinos. The compressed spectra searches are less sensitive to scenarios where thẽ χ AE 1 =χ 0 2 decay through SM W, Z, or Higgs bosons, as the branching fraction to leptonic final states is significantly suppressed. The experimental sensitivity to these scenarios is expected to be recovered with a larger data set.

A. Searches with two opposite-sign light leptons
Previous searches for directχ þ 1χ − 1 production using two opposite-sign light-lepton final states are extended here to increase the sensitivity to compressed SUSY scenarios. The opposite-sign, two-lepton analysis presented here probesχ AE 1 -χ 0 1 mass splittings below 100 GeV using an ISR-jet selection.

Event selection
Events are reconstructed as described in Sec. V, with the signal light-lepton p T threshold raised to p T ¼ 10 GeV. In addition, in events where tagged light leptons form an SFOS pair with m SFOS < 12 GeV, both leptons in the pair are rejected. Events must have exactly two signal light leptons with opposite charge and satisfy the symmetric or asymmetric dilepton trigger criteria, as described in Sec. V.
To suppress the top-quark (tt and Wt) production contribution to the background, events containing central b-tagged jets or forward jets are rejected. To suppress events from Z boson decays, events with invariant mass of the reconstructed SFOS pair within 10 GeV of the Z boson mass (91.2 GeV) are rejected in the same-flavor channel.
Two SRs, collectively referred to as SR2l-1, are defined. Both are designed to provide sensitivity toχ þ 1χ − 1 production withl L -mediated decays and lowχ AE 1 -χ 0 1 mass splittings and rely on a high-p T ISR jet to boost the leptons, which would otherwise have too low momentum to be reconstructed. The super-razor variables that are discussed in Sec. VI A are used to discriminate between signal and backgrounds. Both the same-flavor (SF) and different-flavor (DF) channels are used. The first SR, SR2l-1a, requires R 2 > 0.5ð0.7Þ in the SF (DF) channel, whereas the second SR, SR2l-1b, requires R 2 > 0.65ð0.75Þ. Both SRs require M R Δ > 20 GeV to reduce SM Z þ jets background, and Δϕ β R > 2ð2.5Þ in the SF (DF) to further increase the signal sensitivity. Table VI summarizes the complete definitions of the SRs. SR2l-1a provides sensitivity for moderatẽ χ AE 1 -χ 0 1 mass splittings from 50 GeV to 100 GeV, while SR2l-1b provides sensitivity forχ AE 1 -χ 0 1 mass splittings less than 50 GeV.

Background determination
The SM background is dominated by WW diboson and top-quark production. The MC predictions for these SM sources, in addition to contributions from ZV production, where V ¼ W or Z, are normalized in dedicated control regions for each background. The reducible background is estimated using the matrix method as described in Sec. VI B. Finally, contributions from remaining sources of SM background, which include Higgs boson production and Z þ jets, are small and are estimated from simulation. These are collectively referred to as "Others." ¼1 Veto >20 The top CR is defined using the DF sample in order to suppress events from SM Z boson production. Events are required to have exactly one central light-flavor jet with p T > 80 GeV, no forward jet, and M R Δ > 20 GeV. At least one b-tagged jet is required to enrich the purity in topquark production and ensure orthogonality to the SRs. Figures 5(a) and 5(b) show the M R Δ and Δϕ β R distributions in this CR, respectively. The estimated signal contamination in this CR is less than 1% for the signal models considered.
The WW CR is also defined using the DF sample. Events are required to have exactly one central light jet, no forward jet or b-tagged jet, p ll T > 70 GeV, and M R Δ > 20 GeV. In order to ensure orthogonality to the SRs, Δϕ β R < 2 is required. Figure 5(c) shows the R 2 distribution in this CR. The estimated signal contamination in this CR is less than 20% for the signal models considered.
The ZV CR is defined by using the SF samples and by requiring exactly one central light jet, no forward jet, or btagged jet, p ll T > 70 GeV, Δϕ β R > 2, and M R Δ > 20 GeV. In order to increase the purity in ZV production, events with invariant mass of the reconstructed SFOS pair within 10 GeV of the Z boson mass are used. This requirement also ensures orthogonality to the SRs. Figure 5(d) shows the p ll T distribution in this CR. The estimated signal contamination in this CR is less than 10% for the signal models considered.
A simultaneous likelihood fit to the top, WW, and ZV CRs is performed to normalize the corresponding background estimates to obtain yields in the SR (as described in Sec. VI). Table VI summarizes the definitions of  the CRs, and Table VII summarizes the numbers of observed and predicted events in these CRs, data/MC normalizations, and CR compositions obtained from the simultaneous fit.
Systematic uncertainties affect the estimates of the backgrounds and signal event yields in the control and signal regions. A breakdown of the different sources of systematic uncertainty on the background estimate as described in Sec. VI C is shown in Table VIII. Generator modeling uncertainties are estimated by comparing the results from the POWHEG BOX and MC@NLO event generators for top-quark events, and POWHEG BOX and aMC@NLO for WW events, using HERWIG for parton showering in all cases. Parton showering uncertainties are estimated in top-quark and WW events by comparing POWHEG BOX+HERWIG with POWHEG BOX+PYTHIA. Both generator modeling and parton showering uncertainties are estimated for ZV events by comparing Powheg Box +Pythia to Sherpa. Top-quark samples are generated using AcerMC+Pythia to evaluate the uncertainties related to the amount of initial-and final-state radiation [115]. The impact of the choice of renormalization and factorization scales is evaluated by varying these individually between 0.5 and 2 times the nominal values in POWHEG BOX for topquark events and in aMC@NLO for diboson events. The dominant contributions among the "Theory and modeling" uncertainties come from the generator modeling and parton showering uncertainties.

Results
The observed number of events in each signal region is shown in Table IX    simplified models withl L -mediated decays and differentχ AE 1 andχ 0 1 masses in GeV. The uncertainties shown include both statistical and systematic components. Also shown are the modelindependent limits calculated from the opposite-sign two-lepton signal region observations: the one-sided p 0 values; the expected and observed upper limits at 95% CL on the number of beyondthe-SM events (S 95 exp and S 95 obs ) for each signal region, calculated using pseudoexperiments and the CL s prescription; the observed 95% CL upper limit on the signal cross section times efficiency (hϵσi 95 obs ); and the CL b value for the background-only hypothesis.

B. Searches with two same-sign light leptons
In compressed mass scenarios, one or more of the three leptons fromχ AE 1χ 0 2 production may have momentum too low to be reconstructed. Therefore, the search forχ AE 1χ 0 2 production using two same-sign leptons can complement the three-lepton search documented in Ref. [20] and extend the reach for small mass splittings. The search for same-sign lepton pairs is preferable to opposite-sign pairs, due to the comparatively small SM background. A multivariate analysis technique is used here to discriminate between signal and backgrounds.

Event selection
Events are selected using the basic reconstruction, object, and event selection criteria described in Sec. V. In addition, if tagged light leptons form an SFOS pair with m SFOS < 12 GeV, both leptons in the pair are rejected. Signal electrons with p T < 60 GeV have a tightened track (calorimeter) isolation of 7% (13%) of the electron p T applied, whereas for electrons with p T > 60 GeV, a track isolation requirement of 4.2 GeV (7.8 GeV) is used. For signal muons, the track (calorimeter) isolation requirement is tightened to 6% (14%) of the muon p T for p T < 60 GeV, and 4.2 GeV (8.4 GeV) otherwise. The stricter lepton isolation requirements are optimized to suppress the Events must have exactly two light leptons with the same charge, e AE e AE , μ AE μ AE , or e AE μ AE , and satisfy the symmetric or asymmetric dilepton trigger criteria, as described in Sec. V. Eight BDTs are independently trained to define eight signal regions optimized for four mass splitting scenarios, mðχ 0 2 Þ − mðχ 0 1 Þ ¼ 20, 35, 65, 100 GeV, referred to as ΔM20, ΔM35, ΔM65, and ΔM100, respectively, each with and without the presence of a central light jet with p T > 20 GeV, referred to as ISR and no-ISR. For the BDT training, signal scenarios ofχ AE 1χ 0 2 production withl Lmediated decays are used, where the slepton mass is set at 95% between theχ AE 1 and theχ 0 1 masses. Seven variables are considered as input for the BDT training procedure: m T2 , p T ll , E miss;rel T , H T , m lep1 T , m lep2 T , and Δϕðl; lÞ. Three further variables are also considered for the ISR signal regions: Δϕð; jet1Þ and the ratios E miss;rel T =p T jet1 and p T lep1 =p T jet1 . These variables exploit the kinematic properties of a compressed mass SUSY system, with and without a high-p T ISR jet. The MC simulation samples are compared to data for these variables and their correlations to ensure that they are modeled well.
For the training and testing of the BDT, the signal and background samples are split into two halves, including those backgrounds estimated from data as described in Sec. VIII B 2. The eight signal region definitions are shown in Table X. Since the selection on the BDT output, t cut , is independent for each SR, the overlap between SRs with looser and tighter selections is small.

Background determination
Several SM processes produce events with two samesign signal leptons. The SM background processes are classified as irreducible background if they lead to events with two real, prompt, same-sign leptons, reducible background if the event has at least one fake or nonprompt lepton, or "charge flip" if the event has one lepton with mismeasured charge.
Irreducible processes include diboson (W AE W AE , WZ, ZZ), triboson (VVV), ttV, tZ, and Higgs boson production and are determined using the corresponding MC samples. The reducible Wγ process is estimated with MC simulation samples; other reducible processes are estimated with the matrix method, similar to that described in Sec. VI B.
In this implementation of the matrix method, the fake and nonprompt lepton misidentification probabilities are measured in control regions that are kinematically close and similar in composition to the signal regions. The regions where the misidentification probabilities are measured are required to have large H T (H T > 50 GeV) and large transverse mass using the leading lepton (m T > 50 GeV). The contamination from signal events in these measurement regions is < 1%. The charge-flip, irreducible, and Wγ backgrounds are subtracted from the control regions before calculating lepton misidentification probabilities.
Charge-flip processes include sources of opposite-sign prompt leptons for which the charge of one lepton is mismeasured (Z, tt, W þ W − ). In the relevant momentum range the muon charge-flip background is found to be negligible. Control samples of e þ e − and e AE e AE with invariant mass near the Z boson mass 2 production with small mass splittings between theχ AE 1 =χ 0 2 and LSP. The selection on the BDT output, t cut , is independent for each SR. (75 < m ll < 100 GeV) are used to extract the electron charge-flip rate. A small background due to misidentified jets is subtracted by interpolating the mass sidebands and subtracting them from the observed data events. A likelihood fit is used that takes the numbers of e þ e − and e AE e AE pairs observed in the charge-flip control regions as input.
The charge-flip probability is a free parameter of the fit and is extracted as a function of the electron p T and η. The charge-flip background event yield is found by applying the charge-flip probability to control regions in data with the same kinematic requirements as the signal and validation regions, but with opposite-sign light lepton pairs. The contamination from fake and nonprompt leptons, and from signal events, is negligible in the e þ e − and e AE e AE control regions.
Generator modeling uncertainties for the diboson processes are estimated by comparing the results from the POWHEG BOX and MC@NLO event generators, while parton showering uncertainties are estimated by comparing MC@NLO+HERWIG with MC@NLO+PYTHIA. The impact of the choice of renormalization and factorization scales is evaluated by varying these individually between 0.5 and 2 times the nominal values in aMC@NLO for diboson events.
To test the background prediction methods, two validation regions with looser selection on the BDT output than the SRs are defined; the definitions are shown in Table X. The light-lepton flavor content (ee, μμ, or eμ) is checked separately in each validation region. Table XI

Results
The observed number of events in each signal region is shown in

C. Searches with three light leptons
Previous searches forχ AE 1χ 0 2 production using the threelepton final state are extended here to increase the sensitivity to compressed SUSY scenarios. The three-lepton analysis presented here probesχ 0 2 -χ 0 1 mass splittings below 25 GeV using low-p T leptons and ISR jets.

Event selection
Events are selected as described in Sec. V. In addition, signal muons with p T < 15 GeV have tightened track and calorimeter isolation requirements of 7% of the muon p T . The stricter muon isolation requirements suppress SM backgrounds with semileptonically decaying b=c-hadrons, which are larger for muons rather than electrons due to the lower muon-p T threshold. Events must satisfy a singlelepton, dilepton, or trilepton trigger.
Four signal regions are defined with exactly three light leptons, all with p T < 30 GeV, and at least one SFOS pair present among the leptons. All signal regions veto events with b-tagged jets to reduce the tt SM background and events with 8.4 < m SFOS < 10.4 GeV to suppress backgrounds with leptonic ϒ decays. The three-lepton signal region selections are summarized in Table XIV. The first two signal regions, SR3l-0a and SR3l-0b, closely follow the selection in Ref. [20], using E miss T , m T , and m SFOS selections. SR3l-0a and SR3l-0b are defined with E miss T > 50 GeV and 30 < m lll < 60 GeV to reject diboson processes. Events with a jet with p T > 50 GeV are vetoed to be disjoint from the ISR signal region. The first signal region, SR3l-0a, targets the smallestχ 0 2 -χ 0 1 mass splittings by selecting events with m min SFOS between 4 and 15 GeV. In addition, SR3l-0a requires small m T to reduce the WZ SM background. The second signal region, SR3l-0b, targets the slightly largerχ 0 2 -χ 0 1 mass splittings by selecting events with m min SFOS between 15 and 25 GeV. The third and fourth signal regions, SR3l-1a and SR3l-1b, both require the presence of a p T > 50 GeV jet to target signal production with ISR. The leptons from a 2 production with small mass splittings between theχ AE 1 =χ 0 2 and LSP.
Common Common compressed SUSY decay chain would have too low p T to be reconstructed; however, because of the recoil against the high-p T ISR jet, all three leptons can be boosted enough to meet the selection requirements. The third signal region, SR3l-1a, targets the smallestχ 0 2 -χ 0 1 mass splittings and selects events with 5 < m min SFOS < 15 GeV. Here the leading jet is required to be back-to-back in the transverse plane with the E miss T , ΔϕðE miss T ; jet1Þ > 2.7 rad, and the ratio of leading lepton p T to the jet p T is required to be small, p T lep1 =p T jet1 < 0.2, to suppress the diboson and tt backgrounds. The fourth signal region, SR3l-1b, targets the slightly largerχ 0 2 -χ 0 1 mass splittings by selecting events with 15 < m min SFOS < 25 GeV. To suppress the WZ and tt backgrounds in SR3l-1b, the angle between the E miss T and the three-lepton system is required to be large, ΔϕðE miss T ; 3lÞ > 0.7π rad.

Background determination
Several SM processes produce events with three signal leptons. The SM background processes are classified as irreducible background if they lead to events with three or more real leptons, or as reducible background if the event has at least one fake or nonprompt lepton. The predictions for irreducible and reducible backgrounds are tested in validation regions. For this search, irreducible processes include diboson (WZ and ZZ), VVV, ttV, tZ, and Higgs boson production and are determined from MC simulation samples.
Reducible processes include single and pair production of top quarks, WW production, and a single W or Z boson produced in association with jets or photons. The dominant reducible background component is tt, followed by Z þ jets. The reducible background is estimated using the matrix method, similar to that described in Sec. VI B. In this implementation of the matrix method, the highest-p T signal electron or muon is taken to be real, and only the second and third leptons are used in the matrix method. Simulation studies show that neglecting the case that the leading lepton is nonprompt or fake is valid in more than 95% of the events. The uncertainty on the reducible background includes the MC statistical uncertainty on the weights for the processdependent misidentification probabilities, the uncertainty on the correction factors for the misidentification probability, the statistical uncertainty on the data events to which the matrix equation is applied, and the statistical uncertainty from the misidentification probability measured in simulation.
The systematic uncertainty related to the theoretical modeling of the WZ and ZZ backgrounds is assessed by comparing MC estimates with data in dedicated regions. The WZ region requires three light leptons with p T >30GeV, an SFOS pair among the three leptons, 30 < E miss T < 50 GeV, and one jet with p T > 50 GeV. Events with an SFOS pair or three-lepton invariant mass within 10 GeV of the Z boson mass are vetoed. The ZZ region is defined with four light leptons with p T > 10 GeV, two SFOS pairs with invariant mass within 10 GeV of the Z boson mass, and E miss T < 50 GeV. This approach for estimating the systematic uncertainties is used here instead of the MC-based approach discussed in Sec. VI C. The WZ and ZZ MC simulation samples are both found to agree with observations in the dedicated regions within 15%, which is applied as a systematic uncertainty in the three-lepton validation and signal regions.
The background predictions are tested in validation regions that are defined to be adjacent to, yet disjoint from, the signal regions. Low-E miss T validation regions ("a" regions) and high-E miss T þ b-jet validation regions ("b" regions) are defined to target different background processes. The definition of the regions and the targeted processes are shown in Table XV. In the three-lepton validation regions, the observed data counts and SM expectations are in good agreement within statistical and systematic uncertainties, as shown in Table XVI   XVII. Expected and observed yields in the three-lepton signal regions. The uncertainties shown include both statistical and systematic components. The "Others" background category includes ttV, VVV, and SM Higgs boson production. Also shown are the model-independent limits calculated from the three-lepton signal region observations: the one-sided p 0 -values; the expected and observed upper limits at 95% CL on the number of beyond-the-SM events (S 95 exp and S 95 obs ) for each signal region, calculated using pseudoexperiments and the CL s prescription; the observed 95% CL upper limit on the signal cross section times efficiency (hϵσi 95 obs ); and the CL b value for the background-only hypothesis.

Results
The observed number of events in each signal region is shown in Table XVII Table XVIII, where the dominant sources of uncertainty are the statistical uncertainty on the data for the reducible background estimate, and the uncertainty on the electron and muon misidentification probabilities. Figures 9(a), 9(b), 9(c), and 9(d) show the distributions of the quantities E miss T , m lll , ΔϕðE miss T ; jet1Þ, and p T jet1 in SR3l-0a, SR3l-0b, SR3l-1a, and SR3l-1b regions, respectively, prior to the requirements on these variables. For illustration, the distributions are also shown for aχ AE   This section presents a search for the same-sign charginopair production via VBF with subsequentl L -mediated chargino decays into final states with two same-sign light leptons, at least two jets, and E miss T . Although the cross section for VBF production is significantly lower than that for direct production, the two additional jets in the event provide a means to separate the signal from the background for compressed spectra scenarios and complement the direct production searches that use low-momentum leptons and ISR jets.

A. Event selection
Events are selected using the basic reconstruction, object, and event selection criteria described in Sec. V. In addition, signal muons with p T < 15 GeV have tightened isolation requirements as in the three-lepton analysis described in Sec. VIII C. A tighter isolation is needed for muons rather than electrons due to the lower p T threshold for muons. The stringent lepton isolation suppresses the dominant reducible background processes. Events are required to satisfy an E miss T trigger. One signal region, SR2l-2, is defined with exactly two same-sign light leptons, at least two jets (central light or forward) and large missing transverse momentum E miss T > 120 GeV. In order to select events that originate from VBF production, the highest-p T jet (jet 1) and the second highest-p T jet (jet 2) are required to have large invariant mass, m jj > 350 GeV, to be well separated in pseudorapidity, jΔη jj j > 1.6, and to be in opposite sides of the detector, η jet1 · η jet2 < 0. The last requirement greatly reduces the SM background originating from non-VBF diboson and Higgs boson production. The residual SM background originating from diboson and top-quark production is minimized by requiring the events to have no b-tagged jets, moderate invariant mass of the two leptons (m ll < 100 GeV), small stransverse mass (m T2 < 40 GeV), and a high-p T jet (p T jet1 > 95 GeV). In addition, requirements are made on the ratios of the jet p T , E miss T , p jj T , and p ll T . The SR definition is summarized in Table XIX.

B. Background determination
Several SM processes lead to events with two same-sign signal leptons. The irreducible background is dominated by diboson production, which is estimated using MC simulation samples. The dominant reducible background component is from W þ jets production, followed by tt production, and these are estimated using a data-driven technique called the "fake factor method," similar to that described in Ref. [116]. The production of Wγ is also an important background component and is modeled using MC simulation samples. The charge-flip background is estimated by applying data-driven corrections to the MC simulation samples, following the procedure outlined in Sec. VIII B 2.
The fake factor method estimates the contributions from processes that produce one or two fake or nonprompt   Target process  Dibosons  Nonprompt and fake leptons   TABLE XXI. Observed and expected number of events in the same-sign, two-lepton VBF validation and signal regions. The numbers of signal events are shown for theχ AE 1χ AE 1 VBF simplified model withl L -mediated decays, with theχ AE 1 andχ 0 1 masses in GeV. The uncertainties shown include both statistical and systematic components. The model-independent limits are also shown: the one-sided p 0 value; the expected and observed upper limits at 95% CL on the number of beyond-the-SM events (S 95 exp and S 95 obs ) for the signal region, calculated using pseudoexperiments and the CL s prescription; the observed 95% CL upper limit on the signal cross section times efficiency (hϵσi 95 obs ); and the CL b value for the background-only hypothesis.

VR-VV
VR-Fakes SR2l-2 l flavor/sign l AE l AE , l AE l 0AE l AE l AE , l AE l 0AE l AE l AE , l AE l 0AE leptons using data events that contain one signal lepton and one lepton failing to satisfy the signal lepton requirements. These events are scaled by a "fake factor" to predict the reducible background in the signal region. The fake factor is defined as the ratio of events with two signal leptons to events with one signal lepton and one lepton failing the signal lepton requirements. It is measured in data using a control sample of jets faking leptons in Z → ll events. The SM background process dependence of the fake factor is studied using simulation, and no strong dependence is observed. Residual differences are covered by assigning a 30% uncertainty, independent of the lepton p T , to the fake factor. The uncertainty on the reducible background estimate ranges from 37% to 42%, depending on the channel (ee, μμ, or eμ), and is dominated by the prompt lepton contamination in the control sample and the uncertainty on the extrapolation of fake factors into the signal region.
The contributions from diboson processes are estimated using MC simulation samples. SHERPA is used to produce all diboson samples, taking into account both the strong and the electroweak production of associated jets. The W AE W AE þ 2jets and WZ þ 2jets processes are normalized to NLO cross sections using corrections evaluated in dedicated VBF fiducial regions at the parton level. The corrections are calculated separately for strong and electroweak jet production. For the W AE W AE þ 2jets production, the fiducial cross section is calculated using POWHEG BOX +PYTHIA [62,63,117], and the fiducial region is defined to be identical to the signal region at the parton level, except for the lepton isolation requirement. For the WZ þ 2jets production, the fiducial cross sections are calculated using . Since it is not possible to define a fiducial region that is identical to the signal region using VBFNLO-2.7.0, a looser set of requirements is imposed. The generator modeling uncertainty is estimated by comparing POWHEG BOX+PYTHIA with VBFNLO-2.7.0 for W AE W AE +2jets production, and parton showering uncertainties are estimated by comparing POWHEG BOX+HERWIG with POWHEG BOX+PYTHIA. The impact of the choice of renormalization and factorization scales is evaluated by varying each between 0.5 and 2 times the nominal values.
The uncertainties due to the PDFs are evaluated using 90% CL CT10 PDF eigenvectors. Finally, the interference between the strong and electroweak jet production is studied at LO accuracy using SHERPA and is found to have a negligible effect on the combined fiducial cross section in the signal region.
The background predictions are tested in VRs that are defined to be as kinematically close to the SR as possible. The first VR, VR-Fakes, is defined with two signal light  leptons, large E miss T , and at least two jets to test backgrounds with fake and nonprompt leptons modeled by the fake factor method. The second VR, VR-VV, adopts the same requirements as the VR-Fakes, in addition to higher leptonp T thresholds and a b-jet veto that allow it to test the MC modeling of the diboson background. By definition, the VRs are not disjoint from the SR, but have negligible overlaps. The overlap between the VR-Fakes (VR-VV) and the SR is 2.4% (0.2%), and the largest signal contamination is 1.9% (0.9%) of the total expected background in the VR-Fakes (VR-VV). The definitions of the validation regions are shown in Table XX, along with the targeted processes. The yields in the VRs are shown in Table XXI, where the background expectation is in good agreement with the observed data, within the total uncertainties. Figures 10(a), 10(b), 10(c), and 10(d) show the distributions of p T lep2 and m jj in VR-VV, along with p T lep2 and E miss T in VR-Fakes, with good agreement observed.

C. Results
The observed number of events in the signal region is shown in Table XXI along with the background expectation and uncertainties, p 0 -value, S 95 exp , S 95 obs , hϵσi 95 obs , and the CL b value. No significant excess with respect to the SM expectation is observed. A breakdown of the different sources of systematic uncertainty in the signal region, including those described in Sec. VI C, is shown in Table XXII.  show the distributions of the quantities m jj , jΔη jj j, E miss T , and p T lep2 in the signal region.

X. INTERPRETATION OF RESULTS
Previous ATLAS searches for EW SUSY production [19][20][21][22][23] are combined with the new analyses presented in Secs. VII-IX. The combined results are interpreted in the SUSY models discussed in Sec. II. The analyses combined for each SUSY model are shown in Table XXIII. Limits in the simplified models targeted by the analysis presented in the previous sections are presented in Secs. X A-X D. A summary is provided in Sec. X E, including the limits previously obtained from the ATLAS searches forχ þ 1χ − 1 production with WW-mediated decays [19],χ AE 1χ 0 2 production with WZ-mediated decays [20], andχ AE 1χ 0 2 production with Wh-mediated decays [23]. Finally, limits on phenomenological models are presented in Secs. X F-X H. For these models, the new searches presented in this article are not included, since they target very specific areas of parameter space and their sensitivity is small. Exclusion limits are calculated by statistically combining results from a number of disjoint signal regions. In general, the analyses in Table XXIII are mutually exclusive by design (the exceptions are indicated in the table), using the lepton multiplicity and charge, and are statistically combined. Where overlapping signal regions exist within an analysis, the signal region with the best-expected exclusion is used. During the combinations, all experimental uncertainties are treated as correlated between regions and processes, with the exception of the experimental uncertainties on data-driven backgrounds, which are correlated between regions only. Theoretical uncertainties on the irreducible background and signal are treated as correlated between regions, while statistical uncertainties are treated as TABLE XXIII. Searches used to probe each of the models described in Sec. II.

Model
Wh [23] 2l † [19] 2τ * [22] 3l ⋄ [20] 4l [21] 2τ uncorrelated between regions and processes. For the exclusion limits, the observed and expected 95% CL limits are calculated using asymptotic formulas for each SUSY model point, taking into account the theoretical and experimental uncertainties on the SM background and the experimental uncertainties on the signal. Where the three-lepton [20] analysis is used in the combination, 95% CL limits are calculated using pseudoexperiments as the asymptotic approximation becomes inappropriate where the expected and observed yields are close to zero. The impact of the theoretical uncertainties on the signal cross section is shown for the observed mass limit; where quoted in the text, mass limits refer to the −1σ variation on the observed limit.

A. Direct stau production
The combination of the two-tau MVA results in Sec. VII with the simple cut-based analysis from Ref. [22] is used to set limits on the direct production of stau pairs. For each   1 production withl L -mediated decays, (a) where theχ 0 1 is massless and the intermediate slepton mass is set to 5%, 25%, 50%, 75%, and 95% of theχ AE 1 mass, and (b) as a function of theχ AE 1 andχ 0 1 masses, where the slepton mass is halfway between theχ 0 2 andχ 0 1 masses. The limits in (a) are set using the 2l analysis from Ref. [19], while the limits in (b) use the opposite-sign, two-lepton analysis from this article. The limit from Ref. [19] is also shown in (b). signal point, the signal region with the best expected limit is used. The upper limits on the cross section for direct stau production are shown in Fig. 12 for combinedτ LτL and τ RτR production, where the observed limit is nearly always above the theoretical prediction. One scenario of combined τ LτL andτ RτR production is excluded, where theτ R mass is 109 GeV and theχ 0 1 is massless. For this scenario, cross sections above 0.115 pb are excluded, where the theoretical cross section at NLO is 0.128 pb. No scenarios can be excluded where onlyτ RτR production orτ LτL production is considered. Cross sections above 0.06 (0.21) pb are excluded forτ RτR (τ LτL ) production with aτ R (τ L ) mass of 109 GeV and a masslessχ 0 1 , where the theoretical cross section at NLO is 0.04 (0.09) pb. For this scenario [mðτ R Þ ¼ 109 GeV, mðχ 0 1 Þ ¼ 0 GeV], the expected yields fromτ RτR production are larger than fromτ LτL in the signal region, making the experimental limits stronger forτ RτR production. However, for other mass points the experimental limit is generally weaker forτ RτR production due to   1 production for (a) mðχ AE 1 Þ ¼ 110 GeV and (b) mðχ AE 1 Þ ¼ 120 GeV. The limits are set with respect to the mass difference between theχ AE 1 andχ 0 1 , and use the results from the same-sign, two-lepton VBF analysis. 3 production withl R -mediated decays, (a) where theχ 0 1 is massless and the intermediate slepton mass is set to 5%, 25%, 50%, 75%, and 95% of theχ 0 2 mass, and (b) as a function of theχ 0 2 andχ 0 1 masses, where the slepton mass is halfway between theχ 0 2 andχ 0 1 masses. The limits in (a) and (b) are set using a combination of the 3l analysis from Ref. [20] and the 4l analysis from Ref. [21]. the lower production cross section. These limits on direct production of stau pairs improve upon the previous limits in Ref. [22], particularly for stau masses below ∼150 GeV.

B. Direct chargino production
The opposite-sign, two-lepton analysis in Ref. [19] is used to reinterpret the limits onχ þ 1 production decaying through sleptons, where the slepton mass is varied between theχ AE 1 andχ 0 1 masses. Scenarios where the slepton mass is 5%, 25%, 50%, 75%, and 95% of theχ AE 1 mass are studied for a masslessχ 0 1 , and the limits are shown in Fig. 13(a). For the majority of theχ AE 1 masses considered, the slepton mass does not have a significant effect on the sensitivity, andχ AE 1 masses are excluded up to ∼500 GeV. The sensitivity is reduced for a very small mass splitting between the chargino and the slepton (x ¼ 0.95), as in this case leptons from theχ AE 1 →νl decays have low momentum, making these events difficult to reconstruct in the two lepton final state.
Limits are also set in theχ þ 1χ − 1 scenario withl Lmediated decays, with slepton masses set halfway between theχ AE 1 and theχ 0 1 masses, where both theχ AE 1 and theχ 0 1 masses are varied. Figure 13(b) shows the opposite-sign, two-lepton analysis presented in Sec. VIII A, which provides new sensitivity to compressed scenarios forχ AE 1 masses below ∼220 GeV. The 2l analysis in Ref. [19] continues to dominate the sensitivity to scenarios with large mass splittings, excludingχ AE 1 masses up to ∼465 GeV. The same-sign, two-lepton VBF analysis described in Sec. IX is used to set limits on VBFχ AE 1χ AE 1 production, where theχ AE 1 decays through sleptons. Figures 14(a) and 14(b) show the 95% CL upper limits on the cross section for mðχ AE 1 Þ ¼ 110 GeV and mðχ AE 1 Þ ¼ 120 GeV, as a function of the mass splitting between the chargino and the neutralino. The best observed upper limit on the VBF χ AE 1χ AE 1 production cross section is found for aχ AE 1 mass of 120 GeV and mðχ AE 1 Þ − mðχ 0 1 Þ ¼ 25 GeV, where the theoretical cross section at LO is 4.33 fb and the excluded cross section is 10.9 fb. The sensitivity is slightly stronger for higherχ AE 1 masses, since these scenarios were used for optimizing the signal selection. production withl L -mediated decays, where theχ 0 1 is massless and the intermediate slepton mass is set to 5%, 25%, 50%, 75%, and 95% of theχ AE 1 mass. The limits are set using the 3l analysis from Ref. [20].
[GeV] 2 production withl L -mediated decays, as a function of theχ AE 1 andχ 0 1 masses, where the intermediate slepton mass is set to theχ 0 1 mass plus (a) 50% or (b) 95% of the difference between theχ AE 1 and theχ 0 1 masses. The limits in (a) are set using a combination of the 3l analysis from Ref. [20] and the same-sign, two-lepton analysis from this article, while the limits in (b) use the combination of the three-lepton and same-sign, two-lepton analyses from this article.

C. Direct neutralino production
The combination of the three-lepton analysis in Ref. [20] and four-lepton analysis in Ref. [21] is used to set limits oñ χ 0 2χ 0 3 production withl R -mediated decays, where the slepton mass is varied between theχ 0 2 andχ 0 1 masses. Scenarios where the slepton mass is 5%, 25%, 50%, 75%, and 95% of theχ 0 2 mass are studied for a masslessχ 0 1 , and the limits are shown in Fig. 15(a). For the majority ofχ 0 2 masses considered, the slepton mass does not have a significant effect on the sensitivity, andχ 0 2 masses are excluded up to ∼600 GeV. The sensitivity is reduced for a very small mass splitting between theχ 0 2 and slepton (x ¼ 0.95) as the lepton produced in theχ 0 2 → ll R decay has low momentum. The reduced sensitivity is not seen for a very small mass splitting between the slepton and the LSP (x ¼ 0.05) as the lepton produced in thel R → lχ 0 1 decay can carry some of the momentum of the slepton.
Limits are also set in theχ 0 2χ 0 3 scenario withl R -mediated decays, with slepton masses set halfway between theχ 0 2 and theχ 0 1 masses, where both theχ 0 2 and theχ 0 1 masses are varied. The combination of the three-and four-lepton analysis is again used here and limits are shown in Fig. 15(b), whereχ 0 2 ,χ 0 3 masses up to 670 GeV are excluded, improving the previous limits by 30 GeV for χ 0 1 masses below 200 GeV.

D. Direct neutralino-chargino production
The three-lepton analysis in Ref. [20] is used to reinterpret the limits onχ AE 1χ 0 2 production decaying through sleptons. Scenarios where the slepton mass is 5%, 25%, 50%, 75%, and 95% of theχ AE 1 mass are studied for a masslessχ 0 1 . The limits on these variable slepton mass scenarios are shown in Fig. 16. For the majority ofχ AE 1 masses considered, the slepton mass does not have a significant effect on the sensitivity, andχ AE 1 masses are excluded up to ∼700 GeV. The same reduction in sensitivity is seen for a small mass splitting between theχ 0 2 and slepton (x ¼ 0.95) as in theχ 0 2χ 0 3 interpretation in Sec. X C. Forχ AE 1χ 0 2 production scenarios decaying through SM W, Z, or Higgs bosons [20], the results in Fig. 16 would be degraded due to lower branching fractions into leptonic final states. The pMSSM scenario in Sec. X F shows the sensitivity to SUSY scenarios without sleptons in theχ AE 1χ 0 2 decay chain. Limits are also set in theχ AE 1χ 0 2 scenarios withl Lmediated decays, with slepton masses set halfway and at 95% between theχ AE 1 and theχ 0 1 masses, where both theχ AE 1 and theχ 0 1 masses are varied. Figures 17(a) and 17(b) show that the combination of the published and new analyses gives an improved sensitivity to compressed scenarios up tõ χ AE 1 masses of ∼250 GeV. In scenarios with large mass splittings,χ AE 1 masses are excluded up to ∼700 GeV for slepton masses set to theχ 0 1 mass plus 50% or 95% of the difference between theχ AE 1 and theχ 0 1 masses. In the compressed areas of theχ AE 1χ 0 2 scenario withl L -mediated decays, and slepton masses set halfway (95%) between thẽ χ AE 1 and theχ 0 1 masses, the three-lepton (same-sign, twolepton) analysis has the strongest sensitivity.
Finally, limits are set in theχ AE 1χ 0 2 scenario withτmediated decays, using combined results from the twotau analysis in Ref. [22] and the three-lepton analysis in Ref. [20]. Figure 18 shows that the sensitivity to largeχ AE 1 masses is improved by 20 GeV with the new combination, whereχ AE 1 masses are excluded up to ∼400 GeV for masslessχ 0 1 .

E. Summary of simplified electroweakino production
The ATLAS results for electroweakino searches at 8 TeV in the framework of simplified models are summarized in Figs. 19(a) and 19(b) in the mðχ AE 1 ;χ 0 2 Þ-mðχ 0 1 Þ plane. As explained in Sec. II, each of theχ AE 1 =χ 0 2 =χ 0 3 decays considered in the plot is assumed to have a 100% branching fraction, and the production cross section is for pure winõ χ þ 2 production with decays mediated by SM bosons are summarized in Fig. 19(a). All of the limits are from the two-lepton, three-lepton, and Wh analyses from Refs. [19,20,23]. The new analyses targeting compressed spectra presented in this article have only a small sensitivity to these scenarios and did not significantly improve upon published limits. The limits forχ þ  2 production with τ-mediated decays, as a function of theχ AE 1 andχ 0 1 masses. The limits are set using a combination of the 3l analysis from Ref. [20] and the 2τ analysis from Ref. [22]. summarized in Fig. 19(b). The limits are from the new analyses in Secs. X B-X D and the previously published analyses.

F. pMSSM
The two-lepton, three-lepton, and Wh analyses from Refs. [19,20,23] are combined to improve the sensitivity in the considered pMSSM scenario where the EW SUSY production and the decays through W, Z, or h bosons are dominant. The 95% CL exclusion in the pMSSM μ − M 2 plane for the scenario of heavy sleptons, tan β ¼ 10, and M 1 ¼ 50 GeV is shown in Fig. 20. Including the Wh analysis in the new combination results in a stronger limit at high values of M 2 , in particular in the intermediate μ region.

G. NUHM2
The two-, three-, and four-lepton analyses from Refs. [19][20][21] are combined to set limits in a new interpretation for the NUHM2 model. The 95% CL exclusion in the NUHM2 m 1=2 − μ plane is shown in Fig. 21, where the three-lepton analysis offers the best sensitivity All limits at 95% CL FIG. 20. The 95% CL exclusion limit in the pMSSM scenario, using a combination of the 2l and 3l analyses from Ref. [19] and the Wh analysis from Ref. [23]. The areas excluded by the −1σ expected limit are shown in green. The blue contour corresponds to the limits from the combination of the 2l and 3l analyses from Ref. [19]. The grey dotted contours show the chargino mass isolines. l Observed limit 2 l Observed limit 3 l Observed limit 4 All limits at 95% CL FIG. 21. The 95% CL exclusion limit in the NUHM2 scenario, using a combination of the 2l, 3l, and 4l analyses from Refs. [19][20][21]. The areas excluded by the −1σ expected limit are shown in green. The black, pink, and blue contours correspond to the limits from the 2l, 3l, and 4l analyses, respectively. and drives the combined limit. The results in the three-lepton signal regions lead to a weaker observed exclusion than expected for the compressed scenarios in the high-m 1=2 , low-μ region. In general, m 1=2 values up to 300 GeV are excluded in the NUHM2 model.

H. GMSB
The four-lepton analysis from Ref. [21] is reinterpreted in the GMSB model described in Sec. II. The 95% CL exclusion in the GMSB Λtan β plane is shown in Fig. 22, where Λ values up to 94 TeV are excluded for all values of tan β. For tan β ¼ 10, Λ values below 113 TeV are excluded. These results improve upon the previous limit in Ref. [119] by 20 TeV (15 TeV) in the low (high) tan β region.

XI. CONCLUSION
This article summarizes and extends the search for the production of electroweak SUSY particles using 20 fb −1 of ffiffi ffi s p ¼ 8 TeV pp collision data collected with the ATLAS detector at the LHC. New analyses targeting scenarios with compressed mass spectra, VBF production of charginos and neutralinos, and the direct production of stau pairs provide sensitivity to EW SUSY scenarios not optimally covered in previous publications. The new and previous results are combined to set exclusion limits in a wide range of simplified and phenomenological SUSY models. Forχ þ 1χ − 1 production withl L -mediated decays,χ AE 1 with masses up to ∼500 GeV are excluded. In theχ AE 1χ 0 2 andχ 0 2χ 0 3 scenarios withl L -mediated decays,χ AE 1 andχ 0 2 masses are excluded up to 700 GeV and 670 GeV, respectively. For all threel L -mediated decay scenarios, the value of the slepton mass is not seen to have a significant effect on the sensitivity. Exclusions are also set in pMSSM, NUHM2, and GMSB models, improving upon previous limits.

APPENDIX: CROSS-SECTION CALCULATION FOR THE SAME-SIGN CHARGINO-PAIR PRODUCTION VIA VECTOR-BOSON FUSION
The cross sections for same-sign chargino-pair production via vector-boson fusion (including radiative processes) are calculated to LO in the strong coupling constant using MADGRAPH 5-1.3.33 [57]. The default The 95% CL exclusion limit in the GMSB scenario, using the 4l analysis from Ref. [21]. The green contour corresponds to the limit from the 2SS=3l þ jets analysis from Ref. [119]. value of 99 is used for the maximum number of QCD and QED couplings. Same-sign chargino pairs are generated in association with two additional partons with jηj < 5 and no p T requirement. No jet-parton matching is performed. All SUSY particles, except for theχ AE 1 ,χ 0 1 ,χ 0 2 ,l, andν, are decoupled by setting their physical masses to ∼100 TeV. Theχ AE 1 andχ 0 2 are assumed to be mass degenerate. The sleptons are assumed to be mass degenerate with sneutrinos and have masses set halfway betweenχ AE 1 andχ 0 1 masses. Cross sections are also calculated using MADGRAPH 5-2.2.3 and are in agreement with those calculated using MADGRAPH 5-1.3.33. Details from the "proc_card.dat" are provided below. import model mssm define p¼g u c d s u ∼ c ∼ d ∼ s ∼ define j¼g u c d s u ∼ c ∼ d ∼ s ∼ define lþ¼e þ mu þ define l− ¼ e − mu − define vl¼ve vm vt define vl ∼ ¼ ve ∼ vm ∼ vt ∼ generate p p>x1þ x1þ j j @1 add process p p>x1− x1− j j @2 output −f