Search for pair-produced third-generation squarks decaying via charm quarks or in compressed supersymmetric scenarios in pp collisions at √s = 8 TeV with the ATLAS detector

Aad, G.; et al., [Unknown]; Aben, R.; Angelozzi, I.; Beemster, L.J.; Bentvelsen, S.; Berge, D.; Bobbink, G.J.; Bos, K.; Boterenbrood, H.; Butti, P.; Castelli, A.; Colijn, A.P.; de Jong, P.; de Nooij, L.; Deigaard, I.; Deluca, C.; Deviveiros, P.O.; Dhaliwal, S.; Ferrari, P.; Gadatsch, S.; Geerts, D.A.A.; Hartjes, F.; Hessey, N.P.; Hod, N.; Igonkina, O.; Kluit, P.; Koffeman, E.; Lee, H.C.; Linde, F.; Mahlstedt, J.; Mechnich, J.; Oussoren, K.P.; Pani, P.; Salek, D.; Valencic, N.; van den Wollenberg, W.; van der Deijl, P.C.; van der Geer, R.; van der Graaf, H.; van der Leeuw, R.; van Vulpen, I.; Verkerke, W.; Vermeulen, J.C.; Vreeswijk, M.; Weits, H. DOI 10.1103/PhysRevD.90.052008 Publication date 2014 Document Version Final published version Published in Physical Review D. Particles, Fields, Gravitation, and Cosmology

Search for pair-produced third-generation squarks decaying via charm quarks or in compressed supersymmetric scenarios in pp collisions at ffiffi s p ¼ 8 TeV with the ATLAS detector G. Aad et al. * (ATLAS Collaboration) (Received 3 July 2014;published 24 September 2014) Results of a search for supersymmetry via direct production of third-generation squarks are reported, using 20.3 fb −1 of proton-proton collision data at ffiffi ffi s p ¼ 8 TeV recorded by the ATLAS experiment at the LHC in 2012. Two different analysis strategies based on monojetlike and c-tagged event selections are carried out to optimize the sensitivity for direct top squark-pair production in the decay channel to a charm quark and the lightest neutralino (t 1 → c þχ 0 1 ) across the top squark-neutralino mass parameter space. No excess above the Standard Model background expectation is observed. The results are interpreted in the context of direct pair production of top squarks and presented in terms of exclusion limits in the (m~t 1 , m~χ0 1 ) parameter space. A top squark of mass up to about 240 GeV is excluded at 95% confidence level for arbitrary neutralino masses, within the kinematic boundaries. Top squark masses up to 270 GeV are excluded for a neutralino mass of 200 GeV. In a scenario where the top squark and the lightest neutralino are nearly degenerate in mass, top squark masses up to 260 GeV are excluded. The results from the monojetlike analysis are also interpreted in terms of compressed scenarios for top squark-pair production in the decay channelt 1 → b þ ff 0 þχ 0 1 and sbottom pair production withb 1 → b þχ 0 1 , leading to a similar exclusion for nearly mass-degenerate third-generation squarks and the lightest neutralino. The results in this paper significantly extend previous results at colliders.

I. INTRODUCTION
Supersymmetry (SUSY) [1][2][3][4][5][6][7][8][9] is a theoretically favored candidate for physics beyond the Standard Model (SM). It naturally solves the hierarchy problem and provides a possible candidate for dark matter in the Universe. SUSY enlarges the SM spectrum of particles by introducing a new supersymmetric partner (sparticle) for each particle in the SM. In particular, a new scalar field is associated with each left-and right-handed quark state, and two squark mass eigenstatesq 1 andq 2 result from the mixing of the scalar fields. In some SUSY scenarios, a significant mass difference between the two eigenstates in the bottom squark and top squark sectors can occur, leading to rather light sbottom b 1 and stopt 1 mass states, where the sbottom and stop are the SUSY partners of the SM bottom and top quarks, respectively. In addition, naturalness arguments suggest that the third-generation squarks should be light with masses below 1 TeV [10,11]. In a generic supersymmetric extension of the SM that assumes R-parity conservation [12][13][14][15][16], sparticles are produced in pairs and the lightest supersymmetric particle (LSP) is stable. In this paper the LSP is assumed to be the lightest neutralino [17] (χ 0 1 ).
For a mass difference Δm ≡ m~t 1 − m~χ0 1 > m t and depending on the SUSY parameters and sparticle mass hierarchy, the dominant decay channels are expected to bẽ t 1 → t þχ 0 1 ort 1 → b þχ AE 1 , where the latter decay mode involves charginos (χ AE 1 ) that subsequently can decay into the lightest neutralino via W ðÃÞ emission, leading to a fourbody decayt 1 → b þ ff 0 þχ 0 1 , where ff 0 denotes a pair of fermions (see Fig. 1). If the chargino is heavier than the stop and m W þ m b < Δm < m t , the dominant decay mode is expected to be the three-body Wbχ 0 1 decay. Several searches on 7 TeV data have been carried out in these decay channels in zero-, one-, and two-lepton final states [18][19][20][21] and have been extended using 8 TeV data [22][23][24][25].
In the scenario for which Δm < m W þ m b , the four-body decay mode above competes with the stop decay to a charm quark and the LSP (t 1 → c þχ 0 1 ), which proceeds via a loop decay (see Fig. 1). The corresponding final state is characterized by the presence of two jets from the hadronization of the charm quarks and missing transverse momentum (p miss T denoting its magnitude by E miss T ) from the two undetected LSPs. However, given the relatively small mass difference (Δm), both the transverse momenta of the two charm jets and the E miss T are low, making it very difficult to extract the signal from the large multijet background. In this study, the event selection makes use of the presence of initial-state radiation (ISR) jets to identify signal events. In this case, the squark-pair system is boosted leading to larger E miss T . As an example, for a stop with a mass of 200 GeV and Δm of 5 GeV, about 18% of the events have E miss T > 150 GeV and a jet with p T > 150 GeV. Two different approaches are used to maximize the sensitivity of the analysis across the different Δm regions. A "monojetlike" analysis is carried out, where events with low jet multiplicity and large E miss T are selected, that is optimized for small Δm (Δm ≤ 20 GeV). For Δm ≥ 20 GeV, the charm jets receive a large enough boost to be detected. In addition to the requirements on the presence of ISR jets, the identification of jets containing the decay products of charm hadrons (c tagging) is used, leading to a "c-tagged" analysis that further enhances the sensitivity to the SUSY signal in the region m~t 1 > 200 GeV and Δm ≥ 20 GeV. Results for searches in this channel have been previously reported by collider experiments [26][27][28]. In addition to the decay channelt 1 → c þχ 0 1 , the monojetlike results are reinterpreted in terms of the search for stop pair production witht 1 → b þ ff 0 þχ 0 1 and small Δm. In such a scenario, the decay products of the top squark are too soft to be identified in the final state, and the signal selection relies on the presence of an ISR jet.
In the case of sbottom pair production, assuming a SUSY particle mass hierarchy such that the sbottom decays exclusively asb 1 → b þχ 0 1 (see Fig. 1), the expected signal for direct sbottom pair production is characterized by the presence of two energetic jets from the hadronization of the bottom quarks and large missing transverse momentum from the two LSPs in the final state. Results on searches in this channel at colliders have been reported [21,23,[29][30][31]. In this study, the monojetlike results are also reinterpreted in terms of the search for sbottom pair production withb 1 → b þχ 0 1 in a compressed scenario (small sbottom-neutralino mass difference) with two soft b jets and an energetic ISR jet in the final state.
The paper is organized as follows. The ATLAS detector is described in the next section. Section III provides details of the simulations used in the analysis for background and signal processes. Section IV discusses the reconstruction of jets, leptons, and the E miss T , while Sec. V describes the event selection. The estimation of background contributions and the study of systematic uncertainties are discussed in Secs. VI and VII. The results are presented in Sec. VIII, and are interpreted in terms of the search for stop and sbottom pair production. Finally, Sec. IX is devoted to the conclusions.

II. EXPERIMENTAL SETUP
The ATLAS detector [32] covers almost the whole solid angle around the collision point with layers of tracking detectors, calorimeters, and muon chambers. The ATLAS inner detector has full coverage [33] in ϕ and covers the pseudorapidity range jηj < 2.5. It consists of a silicon pixel detector, a silicon microstrip detector, and a straw tube tracker that also measures transition radiation for particle identification, all immersed in a 2 T axial magnetic field produced by a solenoid.
High-granularity liquid-argon (LAr) electromagnetic sampling calorimeters, with excellent energy and position resolution, cover the pseudorapidity range jηj < 3.2. The hadronic calorimetry in the range jηj < 1.7 is provided by a scintillator-tile calorimeter consisting of a large barrel and two smaller extended barrel cylinders, one on either side of the central barrel. In the end caps (jηj > 1.5), LAr hadronic calorimeters match the outer jηj limits of the end cap electromagnetic calorimeters. The LAr forward calorimeters provide both the electromagnetic and hadronic energy measurements, and extend the coverage to jηj < 4.9.
The muon spectrometer measures the deflection of muon tracks in the large superconducting air-core toroid magnets in the pseudorapidity range jηj < 2.7, using separate trigger and high-precision tracking chambers. Over most of the η range, a precise measurement of the track coordinates in the principal bending direction of the magnetic field is provided by monitored drift tubes. At large pseudorapidities, cathode strip chambers with higher granularity are used in the innermost plane over 2.0 < jηj < 2.7. The muon trigger system covers the pseudorapidity range jηj < 2.4.

III. MONTE CARLO SIMULATION
Monte Carlo (MC) simulated event samples are used to assist in computing detector acceptance and reconstruction efficiencies, determine signal and background contributions, and estimate systematic uncertainties on the final results.
Samples of simulated W þ jets and Z þ jets events are generated using   [34], including leading-order (LO) matrix elements for up to five partons in the final state and using massive b=c quarks, with CT10 [35] parton distribution functions (PDFs) and its own model for FIG. 1 (color online). Diagrams for the pair production of top squarks with the decay modest 1 → c þχ 0 1 ort 1 → b þ ff 0 þχ 0 1 , and the pair production of sbottom squarks with the decay modẽ b 1 → b þχ 0 1 . In one case, the presence of a jet from initial-state radiation is also indicated for illustration purposes.
hadronization. Similar samples are generated using the   [36] generator and are employed to assess the corresponding modeling uncertainties. The MC predictions are initially normalized to next-to-next-to-leadingorder (NNLO) predictions according to DYNNLO [37,38] using MSTW2008 NNLO PDF sets [39].
The production of top-quark pairs (tt) is simulated using the POWHEG-R2129 [40] MC generator. ALPGEN and   [41] MC simulated samples are used to assess tt modeling uncertainties. Single top production samples are generated with POWHEG for the s and Wt channels and MC@NLO is used to determine systematic uncertainties, while ACERMC-V3.8 [42] is used for single top production in the t channel. Finally, samples of tt production associated with additional vector bosons (tt þ W and tt þ Z processes) are generated with   [43]. In the case of POWHEG and MADGRAPH, parton showers are implemented using   [44], while   [45] interfaced to JIMMY [46] is used for the ALPGEN and MC@NLO generators. A top-quark mass of 172.5 GeV and the CTEQ6L1 PDFs are used. The Perugia 2011C [47] and AUET2B [48] tunes for the underlying event are used for the tt, single top, and tt þ W=Z processes, respectively. The cross section prediction at NNLO þ NNLL (next-tonext-to-leading-logarithm) accuracy, as determined by Top þ þ2.0 [49], is used in the normalization of the tt [50] sample. An approximate NLO þ NNLL prediction is used for the Wt [51] process and NLO cross sections are considered for tt þ W and tt þ Z processes.
Diboson samples (WW, WZ, and ZZ production) are generated using SHERPA using massive b=c quarks, with CT10 PDFs, and are normalized to NLO predictions [52]. Additional samples are generated with HERWIG to assess uncertainties. Finally, Higgs boson production including ZH, WH, and ttH processes is generated using  with CTEQ6L1 PDFs.
Stop pair production witht 1 → c þχ 0 1 is modeled with MADGRAPH with one additional jet from the matrix element. The showering is done with PYTHIA-6 and using the AUET2B tune for the underlying event, which involves CTEQ6L1 PDFs. Samples are produced with stop masses between 100 and 400 GeV andχ 0 1 masses between 70 and 390 GeV. The Δm step size increases with Δm from 2 to 30 GeV and the maximum Δm considered is 82 GeV. The region Δm < 2 GeV is not considered since in this regime the stop can become long-lived leading to the signature studied in Ref. [54]. Similarly, MC simulated samples are produced separately fort 1 → b þ ff 0 þχ 0 1 andb 1 → b þ χ 0 1 processes across the stop-neutralino and sbottomneutralino mass planes. In the case of thet 1 → b þ ff 0 þ χ 0 1 process, samples are produced with stop masses in the range between 100 and 300 GeV and Δm that varies between 10 and 80 GeV. For sbottom pair production with b 1 → b þχ 0 1 , samples are produced with sbottom masses in the range between 100 and 350 GeV andχ 0 1 masses in the range between 1 and 340 GeV, with an sbottom-neutralino mass difference that varies between 10 and 50 GeV. Signal cross sections are calculated to NLO in the strong coupling constant, adding the resummation of soft gluon emission at next-to-leading-logarithmic (NLO þ NLL) accuracy [55][56][57]. The nominal cross section and the uncertainty are taken from an envelope of cross-section predictions using different PDF sets and factorization and renormalization scales, as described in Ref. [58].
Differing pileup (multiple proton-proton interactions in the same or neighboring bunch crossings) conditions as a function of the instantaneous luminosity are taken into account by overlaying simulated minimum-bias events generated with PYTHIA-8 onto the hard-scattering process and reweighting them according to the distribution of the mean number of interactions observed. The MC generated samples are processed either with a full ATLAS detector simulation [59] based on GEANT4 [60] or a fast simulation based on the parametrization of the response of the electromagnetic and hadronic showers in the ATLAS calorimeters [61] and a simulation of the trigger system. The results based on fast simulation are validated against fully simulated samples. The simulated events are reconstructed and analyzed with the same analysis chain as for the data, using the same trigger and event selection criteria discussed in Sec. V.

IV. RECONSTRUCTION OF PHYSICS OBJECTS
Jets are reconstructed from energy deposits in the calorimeters using the anti-k t jet algorithm [62] with the distance parameter (in η-ϕ space) ΔR ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi ðΔηÞ 2 þ ðΔϕÞ 2 p set to 0.4. The measured jet transverse momentum (p T ) is corrected for detector effects, including the noncompensating character of the calorimeter, by weighting energy deposits arising from electromagnetic and hadronic showers differently. In addition, jets are corrected for contributions from pileup, as described in Ref. [63]. Jets with corrected p T > 20 GeV and jηj < 2.8 are considered in the analysis. In order to remove jets originating from pileup collisions, central jets (jηj < 2.4) with p T < 50 GeV and with charged-particle tracks associated to them must have a jet vertex fraction (JVF) above 0.5, where the JVF is defined as the ratio of the sum of transverse momentum of matched tracks that originate from the primary vertex to the sum of transverse momentum of all tracks associated with the jet. The presence of leptons (muons or electrons) in the final state is used in the analysis to define control samples and to reject background contributions in the signal regions (see Secs. V and VI). Muon candidates are formed by combining information from the muon spectrometer and inner tracking detectors as described in Ref. [64] and are required to have p T > 10 GeV, jηj < 2.4, and ΔR > 0.4 with respect to any jet with p T > 20 GeV. The latter requirement is increased to 30 GeV in the case of the monojetlike analysis. This increases the efficiency for the selection of real muons from W boson decays. It also avoids biases in the muon selection due to the presence of low-p T jets with large pileup contributions affecting the Wð→ μνÞ þ jets events, as determined by simulations. This is particularly relevant for the monojetlike analysis since, as described in Sec. VI, the Wð→ μνÞ þ jets control samples in data are used to constrain the irreducible Zð→ ννÞ þ jets background contribution in the signal regions. In addition, muons are required to be isolated: the sum of the transverse momenta of the tracks not associated with the muon in a cone of radius ΔR ¼ 0.2 around the muon direction is required to be less than 1.8 GeV.
Electron candidates are initially required to have p T > 10 GeV and jηj < 2.47, and to pass the medium electron shower shape and track selection criteria described in Ref. [65] and reoptimized for 2012 data. Overlaps between identified electrons and jets in the final state are resolved. Jets are discarded if their separation ΔR from an identified electron is less than 0.2. The electrons separated by ΔR between 0.2 and 0.4 from any remaining jet are removed. In the monojetlike analysis, electrons are selected with p T > 20 GeV in both the control and signal regions. The use of the same p T threshold in the control and signal regions minimizes the impact from lepton reconstruction and identification uncertainties on the final results. The 20 GeV p T requirement together with the monojetlike selection also applied to define the control regions brings the background from jets misidentified as electrons to negligible levels without the need for electron isolation requirements. As detailed in Secs. V and VI, slightly different requirements on the lepton p T are applied in the c-tagged analysis to define signal regions and background control samples. In this case, the electrons are required to have p T > 10 GeV and p T > 20 GeV for signal and control samples, respectively, and to be isolated: the total track momentum not associated with the electron in a cone of radius 0.2 around the electron candidate is required to be less than 10% of the electron p T . In the c-tagged analysis, the use of a tighter electron veto in the signal regions, compared to that in the monojetlike analysis, contributes to the reduction of the sizable background from top-quark-related processes.
E miss T is reconstructed using all energy deposits in the calorimeter up to a pseudorapidity jηj < 4.9 and without including information from identified muons in the final state. Clusters associated with either electrons or photons with p T > 10 GeV and those associated with jets with p T > 20 GeV make use of the corresponding calibrations for these objects. Softer jets and clusters not associated with these objects are calibrated using both calorimeter and tracking information [66].
Jets are tagged as containing the decay products of charm hadrons (c tagging) via a dedicated algorithm using multivariate techniques. It combines information from the impact parameters of displaced tracks and topological properties of secondary and tertiary decay vertices reconstructed within the jet. The algorithm provides three probabilities: one targeted for light-flavor quarks and gluon jets (P u ), one for charm jets (P c ), and one for b-quark jets (P b ). From these probabilities, anti-b and anti-u discriminators are calculated: and used for the selected jets in the final state. Figure 2 shows the distributions of the anti-b and anti-u discriminators for the first-and the third-leading jets (sorted in decreasing jet p T ), respectively. The data are compared to MC simulations for the different SM processes, separated by jet flavor [67], and the data-driven multijet background prediction (see Sec. VI C), and include the signal preselection defined in Sec. V without applying the tagging requirements. Good agreement is observed between data and simulations. Two operating points specific to c tagging are used. The medium operating point [log ðP c =P b Þ > −0.9, log ðP c =P u Þ > 0.95] has a c-tagging efficiency of ≈20%, and a rejection factor of ≈8 for b jets, ≈200 for light-flavor jets, and ≈10 for τ jets. The loose operating point ½log ðP c =P b Þ > −0.9] has a c-tagging efficiency of ≈95%, with a factor of 2.5 rejection of b jets but without any significant rejection for light-flavor or τ jets. The efficiencies and rejections are quoted for jets with 30 GeV< p T < 200 GeV and jηj < 2.5 in simulated tt events, and reach a plateau at high jet p T . The c-tagging efficiency is calibrated using data with the method described in Ref. [68] for 7 TeV collisions. This method makes use of a jet sample enriched in charmquark-initiated jets containing a D Ãþ meson identified in the D 0 ð→ K − π þ Þπ þ decay mode [69]. The same calibration method applied to the 8 TeV data leads to reduced uncertainties. The standard calibration techniques are used for the b-jet [70,71] and light-jet [72] rejections: a data-to-simulation multiplicative scale factor of about 0.9, with a very moderate jet p T dependence, is applied to the simulated heavy-flavor tagging efficiencies in the MC samples. The total uncertainty for the c-tagging efficiency varies between 20% at low p T and 9% at high p T and includes uncertainties on the heavy-flavor content of the charm-quark jet enriched sample and on the b-tagging scale factors; uncertainties on the D Ãþ mass fit; uncertainties on the jet energy scale and resolution; and uncertainties on the extrapolation of the results to inclusive charm-quark jets. Similarly, data-to-simulation multiplicative scale factors of order 1.5 are applied to the simulated efficiency for tagging light jets (mistags). They are determined with a precision in the range between 20% and 40% depending on jet p T and η.

V. EVENT SELECTION
The data sample considered in this paper was collected with tracking detectors, calorimeters, muon chambers, and magnets fully operational, and corresponds to a total integrated luminosity of 20.3 fb −1 . The uncertainty on the integrated luminosity is 2.8%, and it is estimated, following the same methodology detailed in Ref. [73], from a preliminary calibration of the luminosity scale derived from beam-separation scans performed in November 2012. The data were selected online using a trigger logic that selects events with E miss T above 80 GeV, as computed at the final stage of the three-level trigger system of ATLAS [74]. With respect to the final analysis requirements, the trigger selection is fully efficient for E miss T > 150 GeV, as determined using a data sample with muons in the final state. Table I summarizes the different event selection criteria applied in the signal regions. The following preselection criteria are applied.
(i) Events are required to have a reconstructed primary vertex consistent with the beamspot envelope and having at least five associated tracks; when more than one such vertex is found, the vertex with the largest summed p 2 T of the associated tracks is chosen.
(ii) Events are required to have E miss T > 150 GeV and at least one jet with p T > 150 GeV and jηj < 2.8 (jηj < 2.5) in the final state for the monojetlike (c-tagged) selection.
(iii) Events are rejected if they contain any jet with p T > 20 GeV and jηj < 4.5 that presents a charged fraction [75], electromagnetic fraction in the calorimeter, or sampling fraction inconsistent with the requirement that they originate from a proton-proton collision [76]. Additional requirements based on the timing and the pulse shape of the cells in the calorimeter are applied to suppress coherent noise and electronic noise bursts in the calorimeter producing anomalous energy deposits [77], which have a negligible effect on the signal efficiency. (iv) Events with isolated muons with p T > 10 GeV are vetoed. Similarly, events with electrons with p T > 20 GeV (p T > 10 GeV) are vetoed in the monojetlike (c-tagged) selection.

A. Monojetlike selection
The monojetlike analysis targets the region in which the stop and the lightest neutralino are nearly degenerate in mass so that the jets from the charm-quark fragmentation (c jets) are too soft to be identified. Stop pair production events are then characterized by large E miss T and a small number of jets, and can be identified via the presence of an energetic jet from initial-state radiation. A maximum of three jets with p T > 30 GeV and jηj < 2.8 in the event are allowed. An additional requirement on the azimuthal separation of Δϕðjet; p miss T Þ > 0.4 between the missing transverse momentum direction and that of each of the selected jets is imposed. This requirement reduces the

B. c-tagged selection
The kinematics of the charm jets from the stop decays depend mainly on Δm. As Δm decreases, the p T of the charm jets become softer and it is more likely that other jets from initial-state radiation have a higher transverse momentum than the charm jets. As a consequence, the stop signal is expected to have relatively large jet multiplicities and a c-tagged jet can be found among any of the subleading jets. An optimization of the c-tagged selection criteria is performed across thet andχ 0 1 mass plane to maximize the sensitivity to a SUSY signal. In the c-tagged analysis, the events are required to have at least four jets with p T > 30 GeV, jηj < 2.5, and Δϕðjet; p miss T Þ > 0.4. A veto against b jets is applied to the selected jets in the event by using a loose c-tag requirement. In addition, at least one of the three subleading jets is required to be c tagged using the medium criteria. The leading jet is required to have p T > 290 GeV and two separate signal regions, here denoted by C1 and C2, are defined with E miss T > 250 GeV and E miss T > 350 GeV, respectively. The tighter requirement on E miss T for the C2 signal region targets models with larger stop and neutralino masses.
The W=Z þ jets backgrounds are estimated using MC event samples normalized using data in control regions. The simulated W=Z þ jets events are reweighted to data as a function of the generated p T of the vector boson, following a procedure similar to that in Ref.
[78] based on the comparison of data and simulation in an event sample enriched in Z þ jets events, which is found to improve the agreement between data and simulation. The weights applied to the simulation result from the comparison of the reconstructed boson p T distribution in data and SHERPA MC simulation in W þ jets and Z þ jets control samples where the jet and E miss T preselection requirements (see Table I) have been applied. The weights are defined in several bins in boson p T . Due to the limited number of data events at large boson p T , an inclusive last bin with boson p T > 400 GeV is used. The uncertainties of the reweighting procedure are taken into account in the final results. The top-quark background contribution to the monojetlike analysis is very small and is determined using MC simulated samples. In the case of the c-tagged analysis, the top-quark background is sizable, as it is enhanced by the jet multiplicity and c-tag requirements, and is estimated using MC simulated samples normalized in a topquark-enriched control region. The simulated tt events are reweighted based on the measurement in the data [79], indicating that the differential cross section as a function of the p T of the tt system is softer than that predicted by the MC simulation.
The normalization factors for W=Z þ jets and tt background contributions are extracted simultaneously using a global fit to all control regions and include systematic uncertainties, to properly take into account correlations. The remaining SM backgrounds from tt þ W=Z, single top, diboson, and Higgs processes are determined using Monte Carlo simulated samples, while the multijet background contribution is extracted from data. Finally, the potential contributions from beam-related background and cosmic rays are estimated in data using jet timing information and are found to be negligible.
In the following subsections, details on the definition of W=Z þ jets and tt control regions and on the data-driven determination of the multijet background are given. This is followed by a description of the background fits and the validation of the resulting background estimations.

A. W=Z þ jets background
In the monojetlike analysis, control samples in data, orthogonal to the signal regions, with identified electrons or muons in the final state and with the same requirements on the jet p T , subleading jet vetoes, and E miss T are used to determine the W=Z þ jets electroweak background contributions from data. A Wð→ μνÞ þ jets control sample is defined using events with a muon with p T > 10 GeV and W transverse mass [80] in the range 30 GeV < m T < 100 GeV. Similarly, a Z=γ Ã ð→ μ þ μ − Þ þ jets control sample is selected, requiring the presence of two muons with invariant mass in the range 66 GeV < m μμ < 116 GeV. The E miss T -based online trigger used in the analysis does not include muon information in the E miss T calculation. This allows the Wð→ μνÞ þ jets and Z=γ Ã ð→ μ þ μ − Þ þ jets control samples to be collected with the same trigger as for the signal regions. Finally, a Wð→ eνÞ þ jets-dominated control sample is defined with an electron candidate with p T > 20 GeV. The E miss T calculation includes the contribution of the energy cluster from the identified electron in the calorimeter, since Wð→ eνÞ þ jets processes contribute to the background in the signal regions when the electron is not identified. In the Wð→ μνÞ þ jets and Z=γ Ã ð→ μ þ μ − Þþ jets control regions, the E miss T does not include muon momentum contributions, motivated by the fact that these control regions are used to estimate the irreducible Zð→ ννÞ þ jets background in the signal regions.
The definition of the control regions in the c-tagged analysis follows closely that of the monojetlike approach with differences motivated by the background composition and the contribution from heavy-flavor jets. A tighter cut of 81 GeV < m μμ < 101 GeV is used to define the Z=γ Ã ð→ μ þ μ − Þ þ jets control sample, as required to further reject tt contamination. This is complemented with a corresponding Z=γ Ã ð→ e þ e − Þ þ jets control sample, with the same mass requirements, for which the energy clusters associated with the identified electrons are then removed from the calorimeter. The Z=γ Ã ð→ e þ e − Þ þ jets control sample is collected using a trigger that selects events with an electron in the final state. As in the monojetlike case, in the Wð→ eνÞ þ jets control region the E miss T calculation includes the contribution from the identified electron. The electron also contributes to the number of jets in the final state, since the presence of a misidentified electron in the signal region can potentially affect the c-tagging results. The c-tagging and the heavy-flavor composition are two of the major uncertainties (of the order of 10%-30%) in the c-tagged selection and the same tagging criteria as used in the signal selection are therefore applied to the Wð→ μνÞ þ jets, Wð→ eνÞ þ jets, Z=γ Ã ð→ μ þ μ − Þ þ jets, and Z=γ Ã ð→ e þ e − Þ þ jets control regions. Since this reduces significantly the selection efficiency related to these control regions, the kinematic selections on the leading jet p T and E miss T are both reduced to 150 GeV, where the trigger selection still remains fully efficient. This introduces the need for a MC-based extrapolation of the normalization factors, as determined using data at relatively low-leading jet p T and E miss T , to the signal regions. This extrapolation is tested in dedicated validation regions as described in Sec. VI E.
Monte Carlo-based transfer factors determined from the SHERPA simulation and including the boson p T reweighting explained above are defined for each of the signal selections to estimate the different electroweak background contributions in the signal regions. As an example, in the case of the dominant Zð→ ννÞ þ jets background process in the monojetlike selection, its contribution to a given signal region N Zð→ννÞ signal is determined using the Wð→ μνÞ þ jets control sample in data according to where N MCðZð→ννÞÞ signal denotes the background predicted by the MC simulation in the signal region, and N data Wð→μνÞ;control , N MC Wð→μνÞ;control , and N non-W Wð→μνÞ;control denote, in the control region, the number of Wð→ μνÞ þ jets candidates in data and MC simulation, and the non-Wð→ μνÞ background contribution, respectively. The N non-W Wð→μνÞ;control term refers mainly to top-quark and diboson processes, but also includes contributions from other W=Z þ jets processes. The transfer factors for each process [e.g., the last term in Eq. (2)] are defined as the ratio of simulated events for the process in the signal region over the total number of simulated events in the control region.
In the monojetlike analysis, the Wð→ μνÞ þ jets control sample is used to define transfer factors for Wð→ μνÞ þ jets and Zð→ ννÞ þ jets processes. As discussed in Secs. VI D and VII, the use of the Wð→ μνÞ þ jets control sample to constrain the normalization of the Zð→ ννÞ þ jets process translates into a reduced uncertainty on the estimation of the main irreducible background contribution, due to a partial cancellation of systematic uncertainties and the statistical power of the Wð→ μνÞ þ jets control sample in data, about 7 times larger than the Z=γ Ã ð→ μ þ μ − Þ þ jets control sample. The Wð→ eνÞ þ jets control sample is used to constrain Wð→ eνÞ þ jets, Wð→ τνÞ þ jets, The c-tagged analysis follows a similar approach to determine the normalization factors for each of the W=Z þ jets background contributions. However, in this case the Zð→ ννÞ þ jets, Z=γ Ã ð→ e þ e − Þ þ jets, and Z=γ Ã ð→ μ þ μ − Þ þ jets normalization factors are extracted from the combined Z=γ Ã ð→ l þ l − Þ þ jets (l ¼ e; μ) control sample, motivated by the fact that these processes involve identical heavy-flavor production mechanisms. Simulation studies indicate a very similar heavy-flavor composition in the control and signal regions. Figure 3 shows, for the M1 monojetlike kinematic selection and in the different control regions, the distributions of the E miss T and the leading-jet p T in data and MC simulations. The MC predictions include data-driven normalization factors as a result of the use of transfer factors from the control to signal regions discussed above. Similarly, the distributions for events in the W=Z þ jets control regions of the c-tagged selection are shown in Fig. 4. Altogether, the MC simulation provides a good description of the shape of the measured distributions for both the monojetlike and c-tagged selections in the different control regions.

B. Top-quark background
The background contribution from top-quark-related production processes to the monojetlike selection is small and is entirely determined from MC simulations. In the case of the c-tagged analysis, single top and tt þ W=Z processes are directly taken from MC simulations and the tt MC predictions are normalized to the data in a separate control region. The tt background contribution is dominated by events with hadronic τ-lepton decays and ISR jets in the final state. A tt control sample is selected with two opposite-charge leptons (ee, μμ, or eμ configurations) in the final state, the same selection criteria for jet multiplicity and c tagging as in the signal region, and relaxed E miss T > 150 GeV and leading jet p T > 150 GeV requirements. In order to reduce the potential Z=γ Ã ð→ e þ e − Þ þ jets and Z=γ Ã ð→ μ þ μ − Þ þ jets contamination in the tt control sample, ee and μμ events with a dilepton invariant mass within 15 GeV of the nominal Z boson mass are rejected. Figure 5 compares the distributions for data and simulation in the tt control region. The MC simulation provides a good description of the shape of the measured distributions.

C. Multijets background
The multijet background with large E miss T mainly originates from the misreconstruction of the energy of a jet in the calorimeter and to a lesser extent is due to the presence of neutrinos in the final state from heavy-flavor decays. In this analysis, the multijet background is determined from data, using a jet smearing method as described in Ref. [81], which relies on the assumption that the E miss T of multijet events is dominated by fluctuations in the jet response in the detector that can be measured in the data. Different response functions are used for untagged and heavy-flavor tagged jets. For the M1 monojetlike and C1 c-tagged analyses, the multijet background constitutes about 1% of the total background, and is negligible for the other signal regions.

D. Background fits
The use of control regions to constrain the normalization of the dominant background contributions from Zð→ ννÞ þ jets, W þ jets (and tt in the case of the c-tagged analysis) reduces significantly the relatively large theoretical and experimental systematic uncertainties, of the order of 20%-30%, associated with purely MCbased background predictions in the signal regions. A complete study of systematic uncertainties is carried out in the monojetlike and c-tagged analyses, as detailed in Sec. VII. To determine the final uncertainty on the total background, all systematic uncertainties are treated as nuisance parameters with Gaussian shapes in a fit based on the profile likelihood method [82], that takes into account correlations among systematic variations. The fit takes also into account cross contamination between different background sources in the control regions.
A simultaneous likelihood fit to the Wð→ μνÞ þ jets, Wð→ eνÞ þ jets, Z=γ Ã ð→ l þ l − Þ þ jets, and tt control  Tables II-IV for the  monojetlike selections, and in Table V for the c-tagged analysis. As the tables indicate, the W=Z þ jets background predictions receive multiplicative normalization factors that vary in the range between 1.1 and 0.9 for the monojetlike analysis, depending on the process and the kinematic selection, and between 0.8 and 0.9 for the c-tagged analyses. In the c-tagged analysis, the tt background predictions are normalized with a scale factor 1.1 for both the C1 and C2 selections.

E. Validation of the background determination
In the monojetlike analysis, the control regions are defined using the same requirements for E miss T , leading jet p T , event topologies, and jet vetoes as in the signal regions, such that no extrapolation in E miss T and jet p T is needed from the control to signal regions. The agreement between data and background predictions is confirmed in a low-p T validation region defined using the same monojetlike selection criteria with E miss T and leading jet p T limited to the range 150-220 GeV.
In the case of the c-tagged analysis, for which the control regions are defined with lower thresholds on the leading jet p T and E miss T compared to those of the signal regions, the Wð→ μνÞ þ jets, Wð→ eνÞ þ jets, Z=γ Ã ð→ l þ l − Þ þ jet, and tt yields fitted in the control regions are then validated in dedicated validation regions (here denoted by V1-V5). The definition of the validation regions is presented in Table VI and is such that there is no overlap of events with the control and signal regions. The validation regions V1-V4 differ from TABLE V. Data and background predictions in the W=Z þ jets and tt control regions before and after the fit is performed for the c-tagged selection. The background predictions include both the statistical and systematic uncertainties. The individual uncertainties are correlated, and do not necessarily add in quadrature to the total background uncertainty.   One medium c tag among jets 2-4 (2-3) for V1-V4 (V5) Three (two) loose c tags acting as b veto, for other 3 (2) jets for V1-V4 (V5) the signal regions only on the thresholds imposed on the E miss T and leading jet p T . In the case of V5, the same requirements as one of the signal regions on E miss T and leading jet p T are imposed but the number of jets is limited to be exactly three. Similar to the transfer factors from the control to signal regions, transfer factors from the control to the validation regions are also defined based on MC simulation. The same experimental systematic uncertainties are evaluated and taken into account in the extrapolation. These transfer factors are subject to the modeling uncertainties of the simulation, which are also applied in the validation regions. Hence, the extrapolation to the validation regions is identical to that of the signal regions. Table VII presents the comparison between data and the scaled MC predictions in the validation regions and Fig. 6 presents the E miss T and leading jet p T distributions for the V3 to V5 regions. Good agreement, within uncertainties, is observed between data and predictions demonstrating a good understanding of the background yields.

VII. SYSTEMATIC UNCERTAINTIES AND BACKGROUND FITS
In this section the impact of each source of systematic uncertainty on the total background prediction in the signal regions, as determined via the global fits explained in Sec. VI D, is discussed separately for monojetlike and c-tagged selections. Finally, the experimental and theoretical uncertainties on the SUSY signal yields are discussed.

A. Monojetlike analysis
Uncertainties on the absolute jet and E miss T energy scale and resolution [63] translate into an uncertainty on the total background that varies between 1.1% for M1 and 1.3% for M3. Uncertainties related to jet quality requirements and pileup description and corrections to the jet p T and E miss T introduce a 0.2% to 0.3% uncertainty on the background predictions. Uncertainties on the simulated lepton identification and reconstruction efficiencies, energy/momentum scale, and resolution translate into a 1.2% and 0.9% uncertainty in the total background for M1 and M3 selections, respectively.
Variations of the renormalization/factorization and parton-shower matching scales and PDFs in the SHERPA W=Z þ jets background samples translate into a 1% to 0.4% uncertainty in the total background. Variations within uncertainties in the reweighting procedure for the simulated W and Z p T distributions introduce less than a 0.2% uncertainty on the total background estimates.
Model uncertainties, related to potential differences between W þ jets and Z þ jets final states, affecting the normalization of the dominant Zð→ ννÞ þ jets and the small Z=γ Ã ð→ τ þ τ − Þ þ jets and Z=γ Ã ð→ e þ e − Þ þ jets background contributions, as determined in the Wð→ μνÞþ jets and Wð→ eνÞ þ jets control regions, are studied in detail. This includes uncertainties related to PDFs and renormalization/factorization scale settings, the parton-shower parameters, and the hadronization model used in the MC simulations, and the dependence on the lepton reconstruction and acceptance. As a result, an additional 3% uncertainty on the Zð→ ννÞ þ jets, Z=γ Ã ð→ τ þ τ − Þ þ jets, and Z=γ Ã ð→ e þ e − Þ þ jets contributions is included for all the selections. Separate studies using parton-level predictions for W=Z þ jet production, as implemented in MCFM-6.8 [83], indicate that NLO strong corrections affect the Wð→ μνÞ þ jetsto-Zð→ ννÞ þ jets ratio by less than 1% in the E miss T and leading jet p T kinematic range considered. In addition, the effect from NLO electroweak corrections on the  W þ jets-to-Z þ jets ratio is taken into account [84][85][86]. Dedicated parton-level calculations are performed with the same E miss T and leading jet p T requirements as in the M1 to M3 signal regions. The studies suggest an effect on the W þ jets-to-Z þ jets ratio that varies between about 2% for M1 and 3% for M2 and M3, although the calculations suffer from large uncertainties, mainly due to the limited knowledge of the photon PDFs inside the proton. In this analysis, these results are conservatively adopted as an additional uncertainty on the Zð→ ννÞ þ jets, Z=γ Ã ð→ τ þ τ − Þ þ jets, and Z=γ Ã ð→ e þ e − Þ þ jets contributions. Altogether, this translates into an uncertainty on the total background that varies from 1.9% and 2.1% for the M1 and M2 selections, respectively, to about 2.6% for the M3 selection.
Theoretical uncertainties on the predicted background yields for top-quark-related processes include uncertainties on the absolute tt, single top, and tt þ Z=W cross sections; uncertainties on the MC generators and the modeling of parton showers employed (see Sec. III); variations in the set of parameters that govern the parton showers and the amount of initial-and final-state soft gluon radiation; and uncertainties due to the choice of renormalization and factorization scales and PDFs. This introduces an uncertainty on the total background prediction that varies between 1.6% and 1.0% for the M1 and M3 selections, respectively. Uncertainties on the diboson contribution are estimated in a similar way and translate into an uncertainty on the total background in the range between 0.7% and 1.3%. A conservative 100% uncertainty on the multijet background estimation is adopted, leading to a 1% uncertainty on the total background for the M1 selection. Finally, statistical uncertainties related to the data control regions and simulation samples lead to an additional uncertainty on the final background estimates in the signal regions that vary between 1.2% for M1 and 1.4% for M3 selections. Other uncertainties related to the trigger efficiency and the determination of the total integrated luminosity [73] are also included, which cancel out in the case of the dominant background contributions that are determined using datadriven methods, leading to a less than 0.3% uncertainty on the total background.

B. c-tagged analysis
In the c-tagged analysis, the jet energy scale uncertainty translates into a 0.3% to 2.2% uncertainty in the final background estimate. Uncertainties related to the loose and medium c tag introduce a 2.8% and 2.5% uncertainty on the background yield for the C1 and C2 selections, respectively. Uncertainties related to the jet energy resolution, soft contributions to E miss T , modeling of multiple pp interactions, trigger and lepton reconstruction, and identification (momentum and energy scales, resolutions, and efficiencies) translate into about a 1.2% (1.4%) uncertainty for the C1 (C2) selection. Variations of the renormalization/ factorization and parton-shower matching scales and PDFs in the SHERPA W=Z þ jets background samples translate into a 3.0% and 3.3% uncertainty in the total background for the C1 and C2 selections, respectively. Uncertainties in the reweighting of the simulated W and Z p T distributions, affecting the extrapolation of the MC normalization factors from the control to the signal regions, introduce a less than 0.6% uncertainty in the final background estimates. In the c-tagged analysis, the Z þ jets and W þ jets background is enriched in heavy-flavor jets produced in association with the vector boson and the same heavy-flavor processes are present in the signal region and the V þ jets control regions. Theoretical uncertainties on the background predictions for top-related processes and diboson contributions are computed following the same prescriptions as in the monojetlike analysis and constitute the dominant sources of systematic uncertainty. In the case of top-related processes, this translates into an uncertainty on the total background prediction of 5.2% and 5.0% for the C1 and C2 selections, respectively. Similarly, the uncertainties on the diboson contributions lead to an uncertainty on the total background of 5.5% (11.5%) for the C1 (C2) selection. The limited number of SM MC events and data events in the control regions lead to an additional uncertainty of 3.0% (4.4%) for the C1 (C2) signal region. Finally, a conservative 100% uncertainty on the multijet background contribution in the control and signal regions is also adopted, which translates into a 0.4% and 0.9% uncertainty on the total background for the C1 and C2 selections, respectively.

C. Signal systematic uncertainties
Different sources of systematic uncertainty on the predicted SUSY signals are considered. Experimental uncertainties related to the jet and E miss T reconstruction, energy scales, and resolutions introduce uncertainties in the signal yields in the range 3% to 7% and 10% to 27% for the monojetlike and c-tagged analyses, respectively, depending on the stop and neutralino masses considered. In the c-tagged analysis, uncertainties on the simulated c-tagging efficiencies for loose and medium tags introduce 9% to 16% uncertainties in the signal yields. In addition, a 2.8% uncertainty on the integrated luminosity is included. Uncertainties affecting the signal acceptance times efficiency (A × ε) related to the generation of the SUSY samples are determined using additional samples with modified parameters. This includes uncertainties on the modeling of the initial-and final-state gluon radiation, the choice of renormalization/factorization scales, and the parton-shower matching scale settings. Altogether this translates into an uncertainty on the signal yields that tends to increase with decreasing Δm and varies between 8% and 12% in the monojetlike analyses, and between 17% and 38% in the c-tagged selections, depending on the stop and neutralino masses. Finally, uncertainties on the predicted SUSY signal cross sections include PDF uncertainties, variations on the α s ðM Z Þ value employed, as well as variations of the renormalization and factorization scales by factors of 2 and 0.5. Altogether, this results in a total theoretical uncertainty on the cross section that varies between 14% and 16% for stop masses in the range between 100 and 400 GeV.

VIII. RESULTS AND INTERPRETATION
The data and the expected background predictions for the monojetlike and c-tagged analyses are summarized in Table VIII. Good agreement is observed between the data and the SM predictions in each case. The SM predictions for the monojetlike selections are determined with a total uncertainty of 2.9%, 3.2%, and 4.6% for the M1, M2, and M3 signal regions, respectively, which include correlations between uncertainties on the individual background contributions. Similarly, the SM predictions for the c-tagged analyses are determined with a total uncertainty of 10% for C1 and 14% for C2 selections. Figure 7 shows the measured leading jet p T and E miss T distributions for the monojetlike selections compared to the background predictions. Similarly, Fig. 8 presents the leading jet p T , E miss T , and jet multiplicity distributions for the c-tagged selections. For illustration purposes, the distributions of two different SUSY scenarios for stop pair production in thet 1 → c þχ 0 1 decay channel with stop masses of 200 GeV and neutralino masses of 125 and 195 GeV are included.
The agreement between the data and the SM predictions for the total number of events in the different signal regions is translated into 95% confidence level (C.L.) upper limits on the visible cross section σ × A × ε using the CL s modified frequentist approach [87], considering the systematic uncertainties on the SM backgrounds, and assuming there is no signal contamination in the control regions. The upper limits are derived from pseudoexperiments and from an asymptotic approximation [82], which gives similar results. For the monojetlike analysis, values of σ × A × ε in the range between 96 and 9.6 fb are excluded at 95% C.L. In the case of the c-tagged analysis, visible cross sections above 1.76 and 0.95 fb, for the C1 and the C2 selections, respectively, are excluded at 95% C.L., as shown in Table IX. A. Stop pair production witht 1 → c þχ 0   1 The results are then translated into exclusion limits on the pair production of top squarks witht 1 → c þχ 0 1 (BR ¼ 100%) as a function of the stop mass for different neutralino masses. Expected and observed 95% C.L. exclusion limits are set using the CL s approach, for which a simultaneous fit to the signal and control regions is performed including statistical and systematic uncertainties. Uncertainties on the signal acceptance times efficiency, the background predictions, and the luminosity are considered, and correlations between systematic uncertainties on signal and background predictions are taken into account. The fit accounts for any potential contamination of signal events in the control regions which a priori has been estimated to be very small. In addition, observed limits are computed using AE1σ variations on the theoretical predictions for the SUSY cross sections. For each SUSY point considered, observed and expected limits are computed separately for the different monojetlike and c-tagged analyses, and the one with the best expected limit is adopted as the nominal result. Finally, the 95% C.L. observed limits corresponding to the −1σ variations on the SUSY theoretical cross sections are then quoted. Figure 9 shows the results separately for the monojetlike and c-tagged analyses, illustrating their complementary regions of sensitivity. As anticipated, the monojetlike selections drive the exclusion limits at very low Δm for which the M2 and M3 signal regions enhance the sensitivity to large stop and neutralino masses. The c-tagged results determine the exclusion limits in the rest of the plane. Figure 10 presents the combined results. Masses for the stop up to 240 GeV are excluded at 95% C.L. for arbitrary neutralino masses, within the kinematic boundaries. For neutralino masses of about 200 GeV, stop masses below 270 GeV are excluded at 95% C.L. In the compressed scenario with the stop and neutralino nearly degenerate in mass, the exclusion extends up to stop masses of 260 GeV. The region with Δm < 2 GeV is not considered in the exclusion since in this regime the stop could become long-lived. These results significantly extend previous exclusion limits [27,28] on the stop and neutralino masses in this channel.

B. Stop and sbottom pair production
The monojetlike results are also interpreted in terms of exclusion limits on the stop pair production in the four-body decay modet 1 → b þ ff 0 þχ 0 1 (BR ¼ 100%) and the sbottom pair production withb 1 → b þχ 0 1 (BR ¼ 100%), using the same CL s approach as explained above. As already mentioned, this is particularly relevant in a mass-degenerate scenario in which the decay products of the squarks are too soft to be identified in the final state, and  and leading jet p T distributions for the C1 selection before the cut in the variable shown (as indicated by the vertical arrows) is applied. In the case of the E miss T distribution, the cuts corresponding to the C1 and C2 selections are both indicated. (Bottom) Measured leading jet p T and jet multiplicity for the C2 selection. The data are compared to the SM predictions. For illustration purposes, the distributions of two different SUSY scenarios are included. The error bands in the ratios include both the statistical and systematic uncertainties on the background predictions.  ) shows the 95% C.L. upper limit on the number of signal events, given the expected number (and AE1σ on the expectation) of background events. The CL B value, i.e. the confidence level observed for the background-only hypothesis, and the p 0 values, which represent the probability of the background alone to fluctuate to the observed numbers of events or higher, are also reported. The p 0 values are truncated at 0.5 if the number of observed events is below the number of expected events. The limits derived using an asymptotic approximation instead of pseudoexperiments are given in parentheses. the signal selection relies on the presence of an ISR jet. Figure 11 shows the expected and observed 95% C.L. exclusion limits as a function of the stop and neutralino masses for thet 1 → b þ ff 0 þχ 0 1 decay channel. For Δm ∼ m b , stop masses up to 255 GeV are excluded at 95% C.L. Top squarks with mass of about 150 and 200 GeV are excluded for m b < Δm < 50 GeV and m b < Δm < 35 GeV, respectively.

Signal region
Finally, Fig. 12 presents the expected and observed 95% C.L. exclusion limits as a function of the sbottom and neutralino masses for theb 1 → b þχ 0 1 decay channel, compared to previous results. In the scenario with (BR ¼ 100%) as determined separately for the monojetlike (left) and the c-tagged (right) selections. The observed (red line) and expected (blue line) upper limits from this analysis are compared to previous results from Tevatron experiments [27,28], and from LEP [26] experiments at CERN with squark mixing angle θ ¼ 0°. The dotted lines around the observed limit indicate the range of observed limits corresponding to AE1σ variations on the NLO SUSY cross-section predictions. The shaded area around the expected limit indicates the expected AE1σ ranges of limits in the absence of a signal. A band for Δm < 2 GeV indicates the region in the phase space for which the stop can become long-lived. . Exclusion plane at 95% C.L. as a function of stop and neutralino masses for the decay channel t 1 → c þχ 0 1 (BR ¼ 100%). The observed (red line) and expected (blue line) upper limits from this analysis are compared to previous results from Tevatron experiments [27,28], and from LEP [26] experiments at CERN with squark mixing angle θ ¼ 0°. The dotted lines around the observed limit indicate the range of observed limits corresponding to AE1σ variations on the NLO SUSY cross-section predictions. The shaded area around the expected limit indicates the expected AE1σ ranges of limits in the absence of a signal. A band for Δm < 2 GeV indicates the region in the phase space for which the stop can become long-lived.  11 (color online). Exclusion plane at 95% C.L. as a function of stop and neutralino masses for the decay channel t 1 → b þ ff 0 þχ 0 1 (BR ¼ 100%). The dotted lines around the observed limit indicate the range of observed limits corresponding to AE1σ variations on the NLO SUSY cross-section predictions. The shaded area around the expected limit indicates the expected AE1σ ranges of limits in the absence of a signal. A band for m~t 1 − m~χ0 1 − m b < 2 GeV indicates the region in the phase space for which the stop can become long-lived. mb 1 − m~χ0 1 ∼ m b , this analysis extends the 95% C.L. exclusion limits up to an sbottom mass of 255 GeV.

IX. CONCLUSIONS
In summary, this paper presents results of a search for stop pair production in the decay channelt 1 → c þχ 0 1 using 20.3 fb −1 of proton-proton collision data at ffiffi ffi s p ¼ 8 TeV recorded with the ATLAS experiment at the LHC. Two different analysis strategies based on monojetlike and c-tagged event selections are carried out that optimize the sensitivity across the stop-neutralino mass plane. Good agreement is observed between the data and the SM predictions. The results are translated into 95% C.L. exclusion limits on the stop and neutralino masses. A stop mass of about 240 GeV is excluded at 95% confidence level for m~t 1 -m~χ0 1 < 85 GeV, as the maximum mass difference in which the decay modet 1 → c þχ 0 1 dominates. Stop masses up to 270 GeV are excluded for a neutralino mass of 200 GeV. In a scenario with the stop and the lightest neutralino nearly degenerate in mass, stop masses up to 260 GeV are excluded. The results from the monojetlike analysis are also reinterpreted in terms of stop pair production in the four-body decay channelt 1 → b þ ff 0 þ χ 0 1 and sbottom pair production withb 1 → b þχ 0 1 , leading to a similar exclusion for the mass-degenerate scenario. The results in this paper significantly extend previous results [23,[26][27][28][29][30] at colliders. FIG. 12 (color online). Exclusion plane at 95% C.L. as a function of sbottom and neutralino masses for the decay channel b 1 → b þχ 0 1 (BR ¼ 100%). The observed (red line) and expected (blue line) upper limits from this analysis are compared to previous results from CDF [29], D0 [30], and ATLAS [23]. For the latter, the area below the dashed-dotted line is excluded. The dotted lines around the observed limit indicate the range of observed limits corresponding to AE1σ variations on the NLO SUSY cross-section predictions. The shaded area around the expected limit indicates the expected AE1σ ranges of limits in the absence of a signal. A band for m~b 1 − m~χ0 1 − m b < 2 GeV indicates the region in the phase space for which the sbottom can become long-lived.
[75] The charged fraction is defined as f ch ¼ P p track;jet T =p jet T , where P p track;jet T is the scalar sum of the transverse momenta of tracks associated with the primary vertex within a cone of radius ΔR ¼ 0.4 around the jet axis, and p jet T is the jet transverse momentum as determined from calorimetric measurements.
[80] The transverse mass m T is defined by the lepton (l) and neutrino (ν) p T and direction as m T ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2p l T p ν T ð1 − cosðϕ l − ϕ ν ÞÞ p , where the ðx; yÞ components of the neutrino momentum are taken to be the same as the corresponding p miss