ATLAS Run 1 searches for direct pair production of third-generation squarks at the Large Hadron Collider

This paper reviews and extends searches for the direct pair production of the scalar supersymmetric partners of the top and bottom quarks in proton–proton collisions collected by the ATLAS collaboration during the LHC Run 1. Most of the analyses use 20 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{fb}^{-1}}$$\end{document}fb-1 of collisions at a centre-of-mass energy of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s} = 8$$\end{document}s=8 TeV, although in some case an additional \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4.7\ {\mathrm{fb}^{-1}}$$\end{document}4.7fb-1 of collision data at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}= 7$$\end{document}s=7 TeV are used. New analyses are introduced to improve the sensitivity to specific regions of the model parameter space. Since no evidence of third-generation squarks is found, exclusion limits are derived by combining several analyses and are presented in both a simplified model framework, assuming simple decay chains, as well as within the context of more elaborate phenomenological supersymmetric models.

The scalar partners of the right-handed and left-handed chiral components of the top-quark state (t R andt L respectively) are among these particles. In many supersymmetric models, the large Yukawa coupling of the top quark to the Higgs sector makes the Higgs boson mass sensitive to the masses of the scalar top (referred to as stop in the following) states, such that, to avoid fine tuning, their masses are often required to be light. Thet R andt L components mix to form the mass eigenstatest 1 andt 2 ,t 1 being defined as the lighter of the two. The scalar superpartner of the lefthanded chiral component of the bottom quark (b L ) belongs to the same weak isospin doublet as thet L , hence they usually share the same supersymmetry-breaking mass parameter: a light stop can therefore imply the existence of a light scalar bottom. The lightest sbottom mass eigenstate is referred to asb 1 .
The ATLAS and CMS collaborations have searched for direct production of stops and sbottoms [16][17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32][33][34][35] using about 4.7 fb −1 of data from the proton-proton collisions produced by the LHC at √ s = 7 TeV and 20 fb −1 at √ s = 8 TeV. These searches have found no evidence of third-generation squark signals, leading to exclusion limits in many SUSY models. The aim of this paper is to sum-marise the sensitivity of the ATLAS experiment to R-parityconserving 1 [38][39][40][41][42] models including the direct pair production of stops and sbottoms using the full √ s = 8 TeV proton-proton collision dataset collected during Run 1 of the LHC. 2 The third-generation squarks are assumed to decay to the stable lightest supersymmetric particle (LSP) directly or through one or more intermediate stages. The analyses considered are those previously published by the ATLAS collaboration on the topic, together with new ones designed to increase the sensitivity to scenarios not optimally covered so far. A wide range of SUSY scenarios are studied by combining different analyses to improve the global sensitivity.
The paper is organised as follows: Sect. 2 briefly reviews the expected phenomenology of third-generation squark production and decay; Sect. 3 reviews the general analysis approach followed by the ATLAS collaboration for SUSY searches; Sects. 4 and 5 present the exclusion limits obtained in specific models by combining the results of several analyses. Two different types of models have been considered: simplified models, where the third-generation squarks are assumed to decay into typically one or two different final states, and more complex phenomenological supersymmetric models, where the stop and sbottom have many allowed decay channels. Conclusions are drawn in Sect. 6.
For the sake of brevity, the body of the paper provides no details of the ATLAS detector and object reconstruction, of the analyses used in the limit derivation, or of how the signal Monte Carlo simulation samples were generated. However, a comprehensive set of appendices is provided to supply additional information to the interested reader. Appendix A briefly summarises the layout of the ATLAS detector and the general principles used in the reconstruction of electrons, muons, jets, jets containing b-hadrons (b-jets), and the missing transverse momentum vector p miss T (whose magnitude is referred to as E miss T ). Appendix B discusses the analyses used to derive the exclusion limits presented in Sects. 4 and 5. The analyses that have already been published are only briefly reviewed, while those presented for the first time in this paper are discussed in detail. Appendix C provides further details of a combination of analyses which is performed for the first time in this paper. Finally, Appendix D provides details about the generation and simulation of the signal Monte Carlo samples used to derive the limits presented. 1 It is also assumed that the decay of the third-generation squarks is prompt: long-lived and metastable stops/sbottoms are discussed elsewhere [36,37]. 2 The analysis exploiting the measurement of the tt cross section discussed in this paper also uses 4.7 fb −1 of proton-proton collisions at √ s = 7 TeV.

Third-generation squark phenomenology
The cross section for direct stop pair production in protonproton collisions at √ s = 8 TeV as a function of the stop mass as calculated with PROSPINO [43,44] is shown in Fig. 1a. It is calculated to next-to-leading order accuracy in the strong coupling constant, adding the resummation of soft gluon emission at next-to-leading-logarithmic accuracy (NLO+NLL) [45][46][47]. In this paper, the nominal cross section and its uncertainty are taken from an envelope of cross-section predictions using different parton distribution function (PDF) sets and factorisation and renormalisation scales described in Ref. [44]. The difference in cross section between the sbottom and stop pair production is known to be small [46], hence the values of Fig. 1a are used for both.
Searches for direct production of stops and sbottoms by the ATLAS collaboration have covered several possible finalstate topologies. The experimental signatures used to identify these processes depend on the masses of the stop or sbottom, on the masses of the other supersymmetric particles they can decay into, and on other parameters of the model, such as the stop and sbottom left-right mixing and the mixing between the gaugino and higgsino states in the chargino-neutralino sector.
Assuming that the lightest supersymmetric particle is a stable neutralino (χ 0 1 ), and that no other supersymmetric particle plays a significant role in the sbottom decay, the decay chain of the sbottom is simplyb 1 → bχ 0 1 (Fig. 2a). A significantly more complex phenomenology has to be considered for the stop, depending on its mass and on thẽ χ 0 1 mass. Figure 1b shows the three main regions in the mt 1 -mχ0 1 plane that are taken into account. They are identified by different values of m(t 1 ,χ 0 1 ) = mt 1 − mχ0 1 . In the region where m(t 1 ,χ 0 1 ) > m t , the favoured decay ist 1 → tχ 0 1 (Fig. 2b). The region where m W +m b < m(t 1 ,χ 0 1 ) < m t is characterised by the three-body decay 3 (t 1 → W bχ 0 1 through an off-shell top quark, Fig. 2c). The region where the value of m(t 1 ,χ 0 1 ) drops below m W + m b , sees the four-body decaỹ t 1 → b f f χ 0 1 , (where f and f indicate generic fermions coming from the decay of an off-shell W boson, Fig. 2d) competing with the flavour-changing decay 4t 1 → cχ 0 1 of Fig. 2e; the dominant decay depends on the details of the supersymmetric model chosen [50].
If the third-generation squark decay involves more SUSY particles (other than theχ 0 1 ), then additional dependencies on 3 In scenarios that depart from the minimal flavour violation assumption, flavour-changing decays liket 1 → cχ 0 1 ort 1 → uχ 0 1 could have a significant branching ratio up to m(t 1 ,χ 0 1 ) ∼ 100 GeV [48]. 4 The decayt 1 → uχ 0 1 , in the assumption of minimal flavour violation [49], is further suppressed with respect tot 1 → cχ 0 1 by corresponding factors of the CKM matrix. gb 1b1 → tχ ± 1 tχ ± 1 ; hb 1b1 → bχ 0 2 bχ 0 2 . The diagrams do not show "mixed" decays, in which the two pair-produced third-generation squarks decay to different final states SUSY parameters arise. For example, if the lightest chargino (χ ± 1 ) is the next-to-lightest supersymmetric particle (NLSP), then the stop tends to have a significant branching ratio fort 1 → bχ ± 1 (Fig. 2f), or, for the sbottom,b 1 → tχ ± 1 if kinematically allowed (Fig. 2g). The presence of additional particles in the decay chain makes the phenomenology depend on their masses. Several possible scenarios have been considered, the most common ones being the gauge-universality inspired mχ± 1 = 2mχ0 1 , favoured, for example, in mSUGRA/CMSSM models [51][52][53][54][55][56]; other interpretations include the case of a chargino almost degenerate with the neutralino, a chargino almost degenerate with the squark, or a chargino of fixed mass. Another possible decay channel considered for the sbottom isb 1 → bχ 0 2 → bhχ 0 1 (Fig. 2h), which occurs in scenarios with a large higgsino component of the two lightest neutralinos.  3 Diagrams oft 2 decays considered as simplified models: ã t 2t2 →t 1 Zt 1 Z ; bt 2t2 →t 1 ht 1 h; ct 2t2 → tχ 0 1 tχ 0 1 . The diagrams do not show "mixed" decays, in which the two pair-produced third-generation squarks decay to different final states. The decayt 2 → γt 1 is not an allowed process Despite the lower production cross section and similar final states tot 1 , the heavier stop state (t 2 ) pair production has also been studied: the search for it becomes interesting in scenarios where the detection oft 1 pair production becomes difficult (for example if m(t 1 ,χ 0 1 ) ∼ m t ). The diagrams of the investigated processes are shown in Fig. 3.
Two types of SUSY models are used to interpret the results in terms of exclusion limits. The simplified model approach assumes that either a stop or a sbottom pair is produced and that they decay into well-defined final states, involving one or two decay channels. Simplified models are used to optimise the analyses for a specific final-state topology, rather than the complex (and model-dependent) mixture of different topologies that would arise from a SUSY model involving many possible allowed production and decay channels. The sensitivity to simplified models is discussed in Sect. 4.
More complete phenomenological minimal supersymmetric extensions of the Standard Model (pMSSM in the following [57]) are also considered, to assess the performance of the analyses in scenarios where the stop and sbottom typically have many allowed decay channels with competing branching ratios. Three different sets of pMSSM models are considered, which take into account experimental constraints from LHC direct searches, satisfying the Higgs boson mass and dark-matter relic density constraints, or additional constraints arising from considerations of naturalness. The sensitivity to these models is discussed in Sect. 5.

General discussion of the analysis strategy
The rich phenomenology of third-generation supersymmetric particles requires several event selections to target the wide range of possible topologies. A common analysis strategy and common statistical techniques, which are extensively described in Ref. [58], are employed.
Signal regions (SR) are defined, which target one specific model and SUSY particle mass range. The event selection is optimised by relying on the Monte Carlo simulation of both the Standard Model (SM) background production processes and the signal itself. The optimisation process aims to maximise the expected significance for discovery or exclusion for each of the models considered.
For each SR, multiple control regions (CR) are defined: they are used to constrain the normalisation of the most relevant SM production processes and to validate the MC predictions of the shapes of distributions of the kinematic variables used in the analysis. The event selection of the CRs is mutually exclusive with that of the SRs. It is, however, chosen to be as close as possible to that of the signal region while keeping the signal contamination small, and such that the event yield is dominated by one specific background process.
A likelihood function is built as the product of Poisson probability functions, describing the observed and expected number of events in the control and signal regions. The observed numbers of events in the various CRs and SRs are used in a combined profile likelihood fit [59] to determine the expected SM background yields for each of the SRs. Systematic uncertainties are treated as nuisance parameters in the fit and are constrained with Gaussian functions with standard deviation equal to their value. The fit procedure takes into account correlations in the yield predictions between different regions due to common background normalisation parameters and systematic uncertainties, as well as contamination from SUSY signal events, when a particular model is considered for exclusion.
The full procedure is validated by comparing the background predictions and the shapes of the distributions of the key analysis variables from the fit results to those observed in dedicated validation regions (VRs), which are defined to be orthogonal to, and kinematically similar, to the signal regions, with low potential contamination from signal.
After successful validation, the observed yields in the signal regions are compared to the prediction. The profile likelihood ratio statistic is used first to verify the SM backgroundonly hypothesis, and, if no significant excess is observed, to exclude the signal-plus-background hypothesis in specific signal models. A signal model is said to be excluded at 95 % confidence level (CL) if the CL s [60,61] of the profile likeli-hood ratio statistics of the signal-plus-background hypothesis is below 0.05.
Several publications, targeting specific stop and sbottom final-state topologies, were published by the ATLAS collaboration at the end of the proton-proton collision run at √ s = 8 TeV, using a total integrated luminosity of about 20 fb −1 . Each of these papers defined one or more sets of signal regions optimised for different simplified models with different mass hierarchies and decay modes for the stop and/or sbottom. A few additional signal regions, focusing on regions of the parameter space not well covered by existing analyses have been defined since then. All signal regions that are used in this paper are discussed in detail in Appendix B, while Table 1 introduces their names and the targeted models. Each analysis is identified by a short acronym defined in the second column of Table 1. The signal region names of previously published analyses are retained, but, to avoid confusion and to ease the bookkeeping, the analysis acronym is prepended to their names. For example, SRA1 from the t0L analysis of Ref. [16], which is a search for stop pair production in channels with no leptons in the final state, is referred to as t0L-SRA1.

Interpretations in simplified models
The use of simplified models for analysis optimisation and result interpretation has become more and more common in the last years. The attractive feature of this approach is that it focuses on a specific final-state topology, rather than on a complex (and often heavily model-dependent) mixture of several different topologies: only a few SUSY particles are assumed to be produced in the proton-proton collision -often just one type -and only a few decay channels are assumed to be allowed. In the remainder of this section, several exclusion limits derived in different supersymmetric simplified models are presented. Details about how the MC signal samples used for the limit derivations were produced are available in Appendix D.

Stop decays with no charginos in the decay chain
A first series of simplified models is considered. It includes direct stop pair production as the only SUSY production process, and assumes that no supersymmetric particle other than thet 1 itself and the LSP, taken to be the lightest neutralinõ χ 0 1 , is involved in the decay. Under this assumption, there is little model dependence left in the stop phenomenology, as discussed in Sect. 2. The stop decay modes are defined mainly by the mass separation m(t 1 ,χ 0 1 ) between the stop and the neutralino, as shown in Fig. 1b. The corresponding diagrams are shown in Fig. 2.  Table 1 and discussed in Appendix B, or by their combination. A detailed discussion of which analysis is relevant in each range of m(t 1 ,χ 0 1 ) follows. m(t 1 ,χ 0 1 ) < m W + m b This kinematic region is characterised by the presence of two competing decays: the flavour-violating decayt 1 → cχ 0 1 (Fig. 2e) and the fourbody decayt 1 → b f f χ 0 1 (Fig. 2d). Which one of the two becomes dominant depends on the model details, in particular on the mass separation between the stop and the neutralino, and on the amount of flavour violation allowed in the model [50]. Several analyses have sensitivity in this region of the mt 1 −mχ0 1 plane. The monojet-like signal regions (tc-M1-3) dominate the sensitivity in the region with m(t 1 ,χ 0 1 ) m b , regardless of the decay of the stop pair, which goes undetected: their selection is based on the presence of an initial-state radiation (ISR) jet recoiling against the stop-pair system, which is assumed to be invisible. At larger values of m(t 1 ,χ 0 1 ), signal regions requiring the presence of a c-tagged jet (tc-C1-2) complement the monojet-like signal regions by targeting thet 1 → cχ 0 1 decay. Limits on four-body decays can be set using signal regions which include low transverse momentum electrons and muons (t1L-bCa_low and WW).
The limits reported in Fig. 4 for these values of m all assume that the branching ratio of the stop decay into either t 1 → cχ 0 1 ort 1 → b f f χ 0 1 is 100 %. However, this assumption can be relaxed, and exclusion limits derived as a function of the branching ratio of thet 1 → cχ 0 1 decay, BR(t 1 → cχ 0 1 ), assuming that BR(t 1 → cχ 0 1 ) + BR(t 1 → b f f χ 0 1 ) = 1. Two different scenarios, with m(t 1 ,χ 0 1 ) = 10, 80 GeV, are considered. The first compressed scenario is characterised by lowp T stop decay products, and the set of signal regions which have sensitivity is the tc-M, independently of the decay of the stop. In the second scenario, the phase space available for thet 1 decay is larger, and the full set of tc-M, tc-C, t1L-bCa_low, t1L-bCa_med and WW-SR selections have different sensitivity, depending on BR(t 1 → cχ 0 1 ). The cross-section limit is derived by combining the analyses discussed above. The SR giving the lowest expected exclusion CL s for each signal model and for each value of BR(t → cχ 0 1 ) is chosen. Figure 5 shows the result of these combinations. For m(t 1 ,χ 0 1 ) = 10 GeV, the sensitivity is completely dominated by the tc-M signal regions, hence no significant dependence on BR(t → cχ 0 1 ) is observed. In this case, stop masses up to about 250 GeV are excluded. For m(t 1 ,χ 0 1 ) = 80 GeV, the sensitivity is dominated by the tc-C signal regions at high values of BR(t → cχ 0 1 ). For lower values of BR(t → cχ 0 1 ), the "soft-lepton" and WW signal regions both become competitive, the latter yielding a higher sensitivity at smaller g3b SR-0 -4j-A, SR-0 -4j-B, SR-0 -4j-C, Gluino-mediatedt 1 andb 1 production, Strongly produced final states with two same-sign or three leptons [63] SS3L SR3b, SR0b, SR1b, Generic gluino and squark production,b 1 → tχ ± 1 SR3Llow, S3Lhigh Spin correlation in tt production events [64] SC -t 1 → tχ 0 1 with mt 1 ∼ m t + mχ0 1 tt production cross section [65] x s e c -t 1 → tχ 0 1 , three-body decay   Fig. 4 Summary of the ATLAS Run 1 searches for direct stop pair production in models where no supersymmetric particle other than thẽ t 1 and theχ 0 1 is involved in thet 1 decay. The 95 % CL exclusion limits are shown in the mt 1 -mχ0 1 mass plane. The dashed and solid lines show the expected and observed limits, respectively, including all uncertainties except the theoretical signal cross-section uncertainty (PDF and scale). Four decay modes are considered separately with a branching ratio of 100 %: m W + m b < m(t 1 ,χ 0 1 ) < m t In this case, the three-body decay of Fig. 2c is dominant. The signal regions that are sensitive to this decay are the dedicated signal region defined in the analysis selecting one-lepton final states (the t1L-3body) and the combination of several signal regions from the analysis selecting two-lepton final states, the t2L. The exclusion limits shown in Fig. 4 assume BR(t 1 → bWχ 0 1 ) = 1. The WW signal regions are found to be sensitive to the kinematic region separating the three-body from the four-body stop decay region. m(t 1 ,χ 0 1 ) ∼ m t In this case, the neutralinos are produced with low p T , and the kinematic properties of the signal are similar to those of SM tt production. Exclusion limits in this region were obtained by two analyses performing precision SM measurements. The first one is the measurement of the tt inclusive production cross section σ tt . Limits ont 1 pair production were already set in Ref. [65], which measured σ tt in the different-flavour, opposite-sign channel eμ. They were derived assuming at 1 decay into an on-shell top quark,t 1 → tχ 0 1 . An extension of the limits into the threebody stop decay is discussed in Appendix B.1. For a massless neutralino, the analysis excludes stop masses from about 150 GeV to about m t . The limit deteriorates for higher neutralino masses, mainly because of the softer b-jet spectrum and the consequent loss in acceptance. The second analysis considered is that of the top quark spin correlation (SC) which considers SM tt production with decays to final states containing two leptons (electrons or muons). The shape and normalisation of the distribution of the azimuthal angle between the two leptons is sensitive to the spin of the produced particles,   The pink lines, whose darkness indicate the value of BR(t → cχ 0 1 ) according to the legend, indicate the observed limit on the production cross section hence it allows the analysis to differentiate between stop pair and tt production. The limit obtained is shown in the bottom middle (dark orange) of the inset of Fig. 4. A small region of m(t 1 ,χ 0 1 ) ≈ 180 GeV is excluded with this measurement assuming a small neutralino mass. m(t 1 ,χ 0 1 ) > m t In this kinematic region, the decayt 1 → tχ 0 1 (see Fig. 2b) is dominant. The best results in this region are obtained by a statistical combination of the results of the multijet (t0L) and one-lepton (t1L) analyses. They both have dedicated signal regions targeting this scenario and the expected sensitivity is comparable for the two analyses. The number of required leptons makes the two signal regions mutually exclusive.
To maximise the sensitivity to thet 1 → tχ 0 1 decays a statistical combination of the t0L and t1L signal regions is performed. The details of the combination are given in Appendix C and the final limit is shown in Fig. 4 by the largest shaded region (yellow). The expected limit on the stop mass is about 50 GeV higher at low mχ0 1 than in the individual analyses. The observed limit is increased by roughly the same amount and stop masses between 200 and 700 GeV are excluded for small neutralino masses. 5 A similar combination is performed to target a scenario where the stop can decay ast 1 → tχ 0 1 with branching ratio x and ast 1 → bχ ± 1 with branching ratio 1 − x. Assuming gauge universality, the mass of the chargino is set to be twice that of the neutralino. Neutralino masses below 50 GeV are not considered, to take into account limits on the lightest chargino mass obtained at LEP [66][67][68][69][70]. The exclusion limits are derived for x = 75, 50, 25 and 0 %. 6 Regardless of the branching ratio considered, it is always assumed that mt 1 > m t + mχ0 1 and mt 1 > m b + mχ± 1 , such that the two decayst → tχ 0 1 andt → bχ ± 1 are both kinematically allowed. A statistical combination, identical to the one described above, is used for x = 75 %. For smaller values of x, no combined fit is performed, as the sensitivity is dominated by the t1L analysis almost everywhere: rather either the t0L or the t1L analysis is used, depending which one gives the smaller expected CL s value. Figure 6 shows the result of the combination in the  5 This result holds if the top quark produced in thet 1 decay has a righthanded chirality. The dependence of the individual limits on the top quark chirality is discussed in Refs. [16,17]. 6 A value of x = 0 % is in fact not achievable in a real supersymmetric model. Nevertheless, this value has been considered as the limiting case of a simplified model.

Stop decays with a chargino in the decay chain
In the pMSSM, unless the higgsino-gaugino mass parameters are related by M 1 μ, M 2 , the mass difference between the lightest neutralino and the lightest chargino cannot be too large. The mass hierarchy mχ0 is, hence, well motivated, leading to the decay chain shown in Fig. 2f.
If additional particles beside the stop and the lightest neutralino take part in the stop decay, the stop phenomenology quickly becomes complex. Even if the chargino is the only other relevant SUSY particle, the stop phenomenology depends on the chargino mass, on the stop left-right mixing, and on the composition of the neutralino and chargino in terms of bino, wino and higgsino states. Figure 7 shows the exclusion limits obtained by the analyses listed in Table 1 and discussed in Appendix B if a branching ratio of 100 % fort → bχ ± 1 is assumed. The exclusion limits are presented in a number of mt 1 -mχ0 1 planes, each characterised by a different hypothesis on the chargino mass. For all scenarios considered, the chargino is assumed to decay asχ ± 1 → W ( * )χ 0 1 , where the ( * ) indicates a possibly virtual W boson. m(χ ± 1 ,χ 0 1 ) = 5, 20 GeV This scenario assumes that the difference in mass between the lightest chargino and the neutralino is small (Fig. 7a), which is a rather common feature of models where, for example, the LSP has a large wino or higgsino component. Two hypotheses have been considered,

Observed limits
Expected limits All limits at 95% CL ) = 10 GeV t1L, t2L, WW Fig. 7 Summary of the ATLAS Run 1 searches for direct stop pair production in models where the decay modet 1 → bχ ± 1 withχ ± 1 → W * χ 0 1 is assumed with a branching ratio of 100 %. Various hypotheses on thẽ t 1 ,χ ± 1 , andχ 0 1 mass hierarchy are used. Exclusion limits at 95 % CL are shown in thet 1 −χ 0 1 mass plane. The dashed and solid lines show the expected and observed limits, respectively, including all uncertainties except the theoretical signal cross-section uncertainty (PDF and scale). Wherever not superseded by any √ s = 8 TeV analysis, results obtained by analyses using 4.7 fb −1 of proton-proton collision data taken at √ s = 7 TeV are also shown, with the corresponding reference. The four plots correspond to interpretations of a the b0L and t1L softlepton analyses in two scenarios ( m(χ ± 1 ,χ 0 1 ) = 5 GeV in light green and m(χ ± 1 ,χ 0 1 ) = 20 GeV in dark green), for a total of four limits; b the b0L, t1L and t2L analyses in scenarios with a fixed chargino mass mχ± For m(χ ± 1 ,χ 0 1 ) = 20 GeV, the signal efficiencies of the b0L signal regions decrease because of the lepton and jet veto applied. The t1L signal regions with soft leptons, instead, gain in sensitivity, profiting from the higher transverse momentum of the fermions from the off-shell W decay produced in the chargino decay. mχ± 1 = 106, 150 GeV This scenario (Fig. 7b) assumes a fixed chargino mass. The SR yielding the lowest expected exclusion CL s for this scenario depends on the value of m(χ ± 1 ,χ 0 1 ). For m(χ ± 1 ,χ 0 1 ) < 20 GeV, the b0L signal regions provide the best sensitivity; for larger values of m(χ ± 1 ,χ 0 1 ), the t1L and t2L signal regions provide better sensitivity because of the same mechanism as in the m(χ ± 1 ,χ 0 1 ) = 5, 20 GeV scenario above. The exclusion extends up to about 600 GeV for small values of m(χ ± 1 ,χ 0 1 ). A region of the parameter space with mt 1 up to about 260 GeV and mχ0 1 between 100 GeV and mχ± 1 is not yet excluded. m(t 1 ,χ ± 1 ) = 10 GeV The fourth scenario ( Fig. 7d) assumes a rather compressedt 1 −χ ± 1 spectrum. The region at low mt 1 and large mχ0 1 is characterised by low mass separations between all particles involved, and it is best covered by the t1L-bCc_diag, the t1L soft lepton, and the WW signal regions. At larger values of the stop mass, the leptons emitted in theχ ± 1 decay have larger p T , and the t2L signal regions provide the best sensitivity. m˜t 1 = 300 GeV The final scenario considered is one where the stop mass is fixed at 300 GeV, and the exclusion limits are expressed in the mχ± 1 -mχ0 1 plane. In the case of the compressed scenario, corresponding to a small mass difference m(χ ± 1 ,χ 0 1 ), the fermions from the W ( * ) decay can escape detection and only the two b-jets and E miss T would be identified in the final state. Thus, the b0L signal regions are expected to have a large sensitivity in this case, while for larger values of m(χ ± 1 ,χ 0 1 ), the lepton can be observed, yielding a final-state signature investigated by the t1L softlepton signal region. A combination of the b0L and t1L signal regions is performed by choosing, for each point of the plane, the SR giving the lowest CL s for expected exclusion. The result, reported in Fig. 8, shows that a large portion of the plane is excluded, with the exception of a region where the mass separations between thet 1 , theχ ± 1 and theχ 0 1 are small.
[GeV] Summarising, in the simplified models witht 1 → bχ ± 1 → bW ( * )χ 0 1 , stop masses up to 450-600 GeV are generally excluded. Scenarios where m(t 1 ,χ 0 1 ) is small are particularly difficult to exclude and in these compressed scenarios, stop masses as low as 200 GeV are still allowed (Fig. 7b). A small unexcluded area is also left for a small region around (mt Fig. 7c), where the sensitivity of the analyses is poor because the signal kinematics are similar to SM tt production.

Limits on pair production oft 2
Although the pair production oft 1 has a cross section larger than that oft 2 , and although the decay patterns of the two particles can be similar, it can be convenient to search for the latter in regions where the sensitivity to the former is limited. This is the case, for example, in the region where Fig. 4, where the separation oft 1 pair production from SM top quark pair production is difficult. The t2t1Z and t2t1h analyses are designed to detectt 2 pair production in this region of the mt 1 − mχ0 1 plane, followed by the decayst 2 →t 1 Z andt 2 →t 1 h. The Higgs boson h is assumed to have a mass of 125 GeV and SM branching ratios.
The exclusion limits were first derived in a scenario in which the pair-producedt 2 decays either throught 2 → Zt 1 with a branching ratio of 100 % (Fig. 3a), or throught 2 → ht 1 (again with a branching ratio of 100 %; Fig. 3b). In both cases, thet 1 is assumed to decay throught 1 → tχ 0 1 , and its mass is set to be 180 GeV above that of the neutralino (assumed to be the LSP), which is the region not excluded in Fig. 4. The final state contains two top quarks, two neutralinos, and either two Z or two h bosons. Exclusion limits at 95 % CL in the scenario wheret 2 pair production is assumed, followed by the decayt 2 → Zt 1 (blue) ort 2 →t 1 h (red) and then byt 1 → tχ 0 1 with a branching ratio of 100 %, as a function of thet 2 andχ 0 1 mass. Thet 1 mass is determined by the relation mt The dashed lines indicate the expected limit and the solid lines indicate the observed limit Figure 9 shows the exclusion limits for the t2t1h and the t2t1Z analyses. In both cases, a limit on mt 2 is set at about 600 GeV for a massless neutralino. In the case of at 2 decay through a Higgs boson, the limit covers neutralino masses lower than in the case of the decay through a Z boson.
The assumption on the branching ratio of thet 2 has also been relaxed, and limits have been derived assuming that the three decayst 2 → Zt 1 ,t 2 → ht 1 andt 2 → tχ 0 1 (

Sbottom decays
Under the assumption that no supersymmetric particle takes part in the sbottom decay apart from the lightest neutralino, the sbottom decays asb 1 → bχ 0 1 with a branching ratio of 100 % (Fig. 2a). The final state arising from sbottom pair production hence contains two b-jets and E miss T . The b0L signal 7 For the combination of the t0L and t1L analyses, the limits extracted for thet 1 → tχ 0 1 decay with branching ratio of 100 % have simply been rescaled by appropriate factors depending on the branching ratio oft 2 → tχ 0 1 considered here.
regions were explicitly optimised to be sensitive to this scenario. In case of a mass degeneracy between the sbottom and the neutralino, the general consideration that the monojet-like tc-M selection is almost insensitive to the details of the decay of the produced particles still holds: the tc-M signal regions offer the best sensitivity for scenarios where mb 1 ∼ mχ0 1 . Figure 11 shows the limits of the tc and b0L analyses on the mb If other supersymmetric particles enter into the decay chain, then multiple decay channels would be allowed. Similarly to the stop, the case in which other neutralinos or charginos have a mass below the sbottom is well motivated. The branching ratios of the sbottom to the different decay channels depend on the supersymmetric particle mass hierarchy, on the mixing of the left-right components of the sbottom, and on the composition of the charginos and neutralinos in terms of bino, wino, and higgsino states.
An exclusion limit is derived under the assumption that the sbottom decays with a branching ratio of 100 % intõ b 1 → tχ ± 1 (Fig. 2g). The chargino is assumed to decay throughχ ± 1 → W ( * )χ 0 1 with a branching ratio of 100 %. The final state is a complex one, and offers many handles for background rejection: it potentially contains up to ten jets, two b-jets, and up to four leptons. The limits of Fig. 12a, shown in the mb 1 − mχ0 1 plane, were obtained by using the three-lepton signal regions SS3L, either fixing the mass of the neutralino to mχ0 1 = 60 GeV or by making the assumption that mχ± 1 = 2mχ0 1 . In the two scenarios considered, sbottom masses up to about 440 GeV are excluded, with a mild dependency on the neutralino mass.
The last case considered is one where the pair-produced sbottoms decay throughb 1 → bχ 0 2 , followed by the decay of χ 0 2 into aχ 0 1 and a SM-like Higgs boson h (Fig. 2h). The final state contains up to six b-jets, four of which are produced by the two Higgs bosons decays. Since multiple b-jets are present in the final state, the three-b-jets signal regions (g3b) are used to place limits in this model.
The limit, derived as a function of mb 1 and mχ0 2 assuming a fixed neutralino mass ofχ 0 1 = 60 GeV, is shown in Fig. 12b. Sbottom masses between about 300 and 650 GeV are excluded forχ 0 2 masses above 250 GeV.

Interpretations in pMSSM models
The interpretation of the results in simplified models is useful to assess the sensitivity of each signal region to a specific topology. However, this approach fails to test signal regions on the complexity of the stop and sbottom phenomenology Expected t2t1Z Observed t2t1h Expected t2t1h Observed t0/t1L comb. Expected t0/t1L comb.   Fig. 11 Observed (solid lines) and expected (dashed lines) 95 % CL limits on sbottom pair production where the sbottom is assumed to decay asb 1 → bχ 0 1 with a branching ratio of 100 %. The purple lines refer to the limit of the tc analysis, while the blue lines refer to the b0L analysis that appears in a realistic SUSY model. To this extent, the signal regions are used to derive exclusion limits in the context of specific pMSSM models.

ATLAS
The pMSSM [57] is obtained from the more general MSSM by making assumptions based on experimental results:

-No new source of CP violation beyond the Standard
Model. New sources of CP violation are constrained by experimental limits on the electron and neutron electric dipole moments. -No flavour-changing neutral currents. This is implemented by requiring that the matrices for the sfermion masses and trilinear couplings are diagonal. -First-and second-generation universality. The soft-SUSYbreaking mass parameters and the trilinear couplings for the first and second generation are assumed to be the same based on experimental data from, e.g., the neutral kaon system [71].
With the above assumptions, and with the choice of a neutralino as the LSP, the pMSSM adds 19 free parameters on top of those of the SM. The complete set of pMSSM parameters is shown in Table 2.
A full assessment of the ATLAS sensitivity to a scan of the 19-parameters space has been performed in Ref. [72]. Here, a set of additional hypotheses are made, to focus on the sensitivity to a specific, well-motivated set of models with enhanced third generation squark production: -The common masses of the first-and second-generation squarks have been set to a multi-TeV scale, making these quarks irrelevant for the processes studied at the energies investigated in this paper. This choice is motivated by the    [85]. The stop sector is therefore completely fixed, given the Higgs boson mass, the value of X t and one of the two stop mass parameters. 8 -The trilinear couplings A b in the sbottom sector are found to have limited impact on the phenomenology, and are therefore set to zero. -The gluino mass parameter M 3 is set such to evade LHC constraints on gluino-pair production.
These assumptions reduce the number of additional free parameters of the model to the mass parameters of the electroweak sector (μ, M 1 , M 2 ) and two of the three third-generation squark mass parameters (mq L3 , mt R , mb R ). All the assumptions made either have a solid experimental basis, or are intended to simplify the interpretation in terms of direct production of stops and sbottoms (as, for example, the assumption on the slepton mass parameters).
Three types of models have been chosen, that, by implementing in different ways constraints arising from naturalness arguments and the dark-matter relic density measurement, further reduce the number of parameters to be scanned over. They are described below, and summarised in Table 3 8 In particular, a minimum value of M S ∼ 800 GeV is allowed if the maximal mixing condition X t /M S = √ 6 is realised.
together with additional information on the most relevant production and decay channels.

Naturalness-inspired pMSSM
The model is inspired by naturalness criteria, which require a value of μ in the range of a few hundred GeV, favour stop masses below one TeV, place weak constraints on the gluino mass and give no constraints on the mass of other SUSY particles [86]. The exclusion limits are determined as a function of the higgsino mass parameter μ and the left-handed squark mass parameter mq L3 . The parameter mq L3 is scanned in the range 350 GeV < mq L3 < 900 GeV. The parameter μ is scanned in the range 100 GeV < μ < mq L3 − 150 GeV, where the lower bound is determined by limits on the chargino mass arising from LEP [66 -70]. The right-handed stop mass parameter mt R and the stop mixing parameter X t are determined by choosing the maximal mixing scenario X t /M S = √ 6 and by the requirement of having a Higgs boson mass of about 125 GeV. The other squark and slepton masses, as well as the bino mass parameter M 1 , are set to 3 TeV. The wino mass parameter M 2 is set such that M 2 = 3μ. The gluino mass parameter M 3 is set to 1.7 TeV.
With this choice of the model parameters, the spectrum is characterised by two light neutralinos χ 0 1 ,χ 0 2 and one chargino χ ± 1 , all with masses of the order of μ, a lightb 1 with a mass of the order of mq L3 , and a lightt 1 with mass of the order of mq L3 up to mq L3 ∼ 700 GeV (the constraint on M S does not allow the mass oft 1 to increase beyond about 650 GeV). The production processes considered are direct pair production ofb 1 andt 1 with similar masses. Because of the abundance of light higgsino states, many different decays can occur.
Well-tempered neutralino pMSSM The models are designed to loosely satisfy dark-matter thermal-relic density constraints (0.09 < c h 2 < 0.15, where h is the Hubble constant), while keeping fine tuning (defined as in Ref. [87]) to less than 1 %. The exclusion limits are determined as a function of M 1 and mq L3 , or M 1 and mt R , with μ ∼ −M 1 in both cases to satisfy the dark-matter constraints through the presence of well-tempered neutralinos [88]. The constraints on the Higgs boson mass are satisfied in a way similar to the naturalness-inspired pMSSM model above. All other parameters are the same as in the naturalness-inspired pMSSM model. These models tend to have three neutralinos and two charginos with masses lower thant 1 orb 1 , giving rise to a diverse phenomenology.
h/Z-enriched pMSSM These models are defined such that Higgs and Z bosons are produced abundantly in the SUSY particles' decay chains. The assumption of M 1 = 100 GeV ensures the presence of a bino-like neutralino LSP, while M 3 = 2.5 TeV ensures that direct gluino production is highly Exclusion limits for these pMSSM models are determined by combining many of the SRs defined for the searches discussed in this paper (t0L, t1L, tb, 9 t2t1Z, g3b, tc). For each 9 The tb signal region, discussed in detail in Appendix B.2.3, implement a one-lepton selection, designed to be sensitive to final states containing a top quark, a b-quark and E miss T . It complements the selections of the t0L and t1L signal regions targeting tt E miss T final states.  Figure 13 shows the exclusion limit for the naturalnessinspired set of pMSSM models based on the t0L, t1L and tb analyses. The t0L and t1L analyses have a similar expected sensitivity. These SRs were optimised assuming a 100 % BR fort 1 → tχ 0 1 ort 1 → bχ ± 1 , while for these pMSSM models, the stop decays tot 1 → tχ 0 1 ,t 1 → bχ ± 1 and t 1 → bχ 0 2 with similar branching ratios (and the sbottom to bothb 1 → bχ 0 1 andb 1 → tχ ± 1 ). The tb signal regions, discussed in detail in Appendix B.2.3, are designed to be sensitive to final states containing a top quark, a b-quark and missing transverse momentum and address such mixeddecay scenarios by requiring a lower jet multiplicity.
The signal regions that dominate the sensitivity are the tb, t0L-SRC1 and t1L-bCd_bulk at low values of mq L3 , and tb, t0L-SRA1, t0L-SRA2 and t1L-tNbC_mix at intermediate and high values of mq L3 . The excluded region for models with mq L3 ∼ 900 GeV and μ ∼ 150 GeV is due to the saturation of mt 1 at high mq L3 values: to satisfy the Higgs boson mass constraint requires M S ∼ 800 GeV, hence mt 1 at mq L3 ∼ 900 GeV is smaller than that at mq L3 ∼ 800 GeV. The large fluctuations of the observed limit with respect to the expected one are due to transitions between different signal regions providing the best expected exclusion in different regions of the plane. Figure 14a, b show the exclusion limit obtained for the set of pMSSM models with well-tempered neutralinos as a function of mq L3 and mt R , respectively. In both cases, the exclusion is largely dominated by the t0L analysis. For  Fig. 14a, the signal region dominating the sensitivity at low mq L3 is t0L-SRC1, while at higher mq L3 values t0L-SRA1 and t0L-SRA2 dominate the sensitivity. The drop in sensitivity at mq L3 = 410 GeV, M 1 = 260 GeV is due to the opening of thet 1 → tχ 0 2 andt 1 → tχ 0 3 transition, kinematically suppressed for smaller values of the difference mq L3 − M 1 . Such decays introduce more intermediate states in the decay, effectively reducing the transverse momenta of the final state objects. The large fluctuations of the observed limit are again due to transitions between different signal regions. For Fig. 14b, the sensitivity is entirely dominated by the various t0L-SRC. The difference in sensitivity between these two scenarios is due to the presence of both a stop and Finally, Fig. 15a, b show the exclusion limit obtained for the set of h/Z -enriched pMSSM models. These models yield large b-jet multiplicities to the final state through direct sbottom decays, top-quark decays andχ 0 2 → h/Zχ 0 1 . The exclusion is dominated by the t0L and g3b analyses for Fig. 15a and by and the t0L analysis for Fig. 15b.
More informations about the limits obtained, including the SLHA files for the points mentioned in Table 3, can be found in Refs.

Conclusions
The search programme of the ATLAS collaboration for the direct pair production of stops and sbottoms is summarised and extended by new analyses targeting scenarios not optimally covered by previously published searches. The paper is based on 20 fb −1 of proton-proton collisions collected at the LHC by ATLAS in 2012 at a centre-of-mass energy √ s = 8 TeV. Exclusion limits in the context of simplified models are presented. In general, stop and sbottom masses up to several hundred GeV are excluded, although the exclusion limits significantly weaken in the presence of compressed SUSY mass spectra or multiple allowed decay chains. Three classes of pMSSM models, based on general arguments of Higgs boson mass naturalness and compatibility with the observed dark-matter relic density have also been studied and exclusion limits have been set. Large regions of the considered parameter space are excluded.
We thank CERN for the very successful operation of the LHC, as well as the support staff from our institutions without whom ATLAS could not be operated efficiently.

A The ATLAS detector and object reconstruction
The ATLAS detector [91] consists of inner tracking devices surrounded by a superconducting solenoid, electromagnetic and hadronic calorimeters and a muon spectrometer immersed in a toroidal magnetic field. The inner detector (ID), in combination with a superconducting solenoid magnet with a central field of 2 T, provides precision tracking and momentum measurements of charged particles in a pseudorapidity 10 range |η| < 2.5. The ID consists of a silicon pixel detector, a silicon microstrip detector and a straw tube tracker (|η| < 2.0) that also provides transition radiation measurements for electron identification. A high-granularity electromagnetic calorimeter system, with acceptance covering |η| < 3.2, uses liquid argon (LAr) as the active medium. A scintillator-tile calorimeter provides hadronic coverage for |η| < 1.7. The end-cap and forward regions, spanning 1.5 < |η| < 4.9, are instrumented with LAr electromagnetic and hadronic calorimeters. The muon spectrometer has separate trigger and high-precision tracking chambers which provide trigger coverage for |η| < 2.4 and muon identification and momentum measurements for |η| < 2.7.
The data sample used in this analysis was taken during the period from March to December 2012 with the LHC operating at a pp centre-of-mass energy of √ s = 8 TeV. 11 Following requirements based on beam, detector conditions and data quality, the complete dataset corresponds to an inte- 10 ATLAS uses a right-handed system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upward. Cylindrical coordinates (r, φ) are used in the transverse plane, φ being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2). The distance R in the η-φ space is defined as R = ( η) 2 + ( φ) 2 . 11 The limits derived using a measurement of the tt production cross section discussed in Sect. 4.1 and Appendix B.1 also uses 4.7 fb −1 of pp collisions data collected at a centre-of-mass energy of √ s = 7 TeV. grated luminosity of 20.3 fb −1 , with an associated uncertainty of 2.8 %. The uncertainty is derived following the same methodology as that detailed in Ref. [92]. Events used in the analyses presented in this paper were selected using the ATLAS three-level trigger following different chains based on the signatures being considered. A common set of cleaning cuts, aimed at rejecting events heavily contaminated by non-collision backgrounds, or events containing badly measured or fake jets is applied to all analyses. The experimental signature of third-generation supersymmetric particles includes the production of b-jets in association with missing transverse momentum and possibly additional jets and charged leptons. Different signatures are investigated in this paper to gain sensitivity to a variety of possible topologies arising from the production and decay of stops and sbottoms. Different event selections share common definitions of the final reconstructed objects, which are detailed in the remainder of this Appendix. Analysis-specific departures from those definitions are detailed for each case in Appendix B or in the specific analysis paper.
The reconstructed primary vertex [93] is required to be consistent with the luminous region and to have at least five associated tracks with p T > 400 MeV; when more than one such vertex is found, the vertex with the largest summed p 2 T of the associated tracks is chosen.
Jets are constructed from three-dimensional clusters of noise-suppressed calorimeter cells [94] using the anti-k t algorithm [95][96][97] with a distance parameter R = 0.4 and calibrated with a local cluster weighting algorithm [98]. An area-dependent correction is applied for energy from additional proton-proton collisions based on an estimate of the pileup activity in a given event using the method proposed in Ref. [99]. Jets are calibrated as discussed in Ref.
[100] and required to have p T > 20 GeV and |η| < 4.5. Events containing jets arising from detector noise, cosmic-ray muons, or other non-collision sources are removed from consideration [100].
Jets arising from a b-quark fragmentation and within the acceptance of the inner detector (|η| < 2.5) are identified with an algorithm that exploits both the track impact parameters and secondary vertex information [101]; this algorithm is based on a neural network using the output weights of the IP3D, JetFitter+IP3D, and SV1 algorithms (defined in Refs. [102,103]). A lower cut on the output of the neural network defines the b-tagged jets. Three different working points are used, with a nominal efficiency of 60, 70 and 80 % as evaluated on simulated top quark pair production events. The corresponding rejection factors against jets originating from light (c) quarks are 25 (3), 135 (5) and 600 (8).
Electrons are reconstructed from energy clusters in the electromagnetic calorimeter matched to a track in the inner detector [104] and are required to have |η| < 2.47. Several criteria, including calorimeter shower shape, quality of the match between the track and the cluster, and the amount of transition radiation emitted in the TRT detector, are used to define three selections with decreasing efficiency and increasing purity, named respectively 'loose', 'medium' and 'tight' [104]. These three electron selections are used throughout this paper in the definitions of various signal and control regions. Muons, which are identified either as a combined track in the muon spectrometer and inner detector systems, or as an inner detector track matched with a muon spectrometer track segment [105,106], are required to have |η| < 2.4.
Electrons and muons (generically referred to by the symbol ) are usually required to have transverse momentum p T > 10 GeV. For specific scenarios with compressed mass spectra, lowp T leptons are expected and the p T threshold is lowered to 6 GeV for muons and to 7 GeV for electrons.
The missing transverse momentum p miss T (with magnitude E miss T ) is the negative vector sum of the p T measured in the clusters of calorimeter cells, which are calibrated according to their associated reconstructed object (e.g. jets and electrons), and the p T of the muons. Calorimeter cells not associated with any reconstructed object are also used in the calculation of p miss T . The missing transverse momentum from the tracking system (denoted by p miss,track T , with magnitude E miss,track T ) is computed from the vector sum of the reconstructed inner detector tracks with p T > 500 MeV and |η| < 2.5, associated with the primary vertex in the event.

B Analyses used in the paper
Several signal regions are used in this paper, either standalone or in combination with others, to derive exclusion limits in the many models considered. This Appendix provides a review of the already published analyses and a more extended documentation of the signal regions not previously published.

B.1 Review of already published signal regions
The discussion of analyses that have already been published is reduced to a summary for the sake of brevity. Table 1 provides a reference to the papers where full details of the signal, control and validation region selections, together with the strategies adopted for the estimation of the background processes are found.

Multijet final states (t0L)
The analysis is designed to be sensitive to final states arising from all-hadronic decays of directly pair-produced stops [16]. Two sets of signal regions were optimised to maximise the sensitivity to topologies arising fromt 1 → tχ 0 1 decays, assumed to happen with a branching ratio of one. The first set of signal regions, named t0L-SRA, assumes that both top quark hadronic decays can be fully resolved by indentifying the six final-state jets. The SM background [dominated by tt and Z + heavy flavour (HF) jets production] is rejected based on the presence of two hadronic systems consistent with top quarks and large E miss T . The second set of signal regions, named t0L-SRB targets a similar scenario, but aims at topologies where the top quarks have a large boost, and some of the decay products are merged into a single jet. The event selection is designed to select final states with a maximum of five R = 0.4 anti-k t jets, to be mutually exclusive with t0L-SRA, and relies on the presence of R = 0.8 and R = 1.2 anti-k t jets containing the hadronic decay products of the two top quarks. The jet masses, the transverse mass of the E miss T and the nearest b-jet, and other variables are used to discriminate against the dominant SM tt, Z + HF jets and W + HF jets production background processes.
Finally, a third set of signal regions, named t0L-SRC, is designed to increase the analysis sensitivity to the decayt 1 → bχ ± 1 . The presence of the intermediate chargino state tends to decrease the jet multiplicity: these signal regions require five anti-k t jets with R = 0.4, and base the signal selection on a set of transverse mass variables aimed at rejecting the dominant SM tt production process.

One-lepton final states (t1L)
The large number of signal regions defined in this analysis stems from the variety and complexity of the possible stop final states considered [17]. All signal regions are characterised by the presence of one lepton, a second-lepton veto, a minimum of two jets and large E miss T . A first set of four signal regions (t1L-tN) were optimised assuming a branching ratio of 100 % for the decaỹ t → tχ 0 1 . These signal regions aim at having sensitivity to different m(t,χ 0 1 ), in particular t1L-tN_diag targets scenarios with small m(t 1 ,χ 0 1 ) and makes use of the shape information of the E miss T and m T distributions. 12 The t1L-tN_boost SR targets models with the largest m(t,χ 0 1 ), where the top quark produced by the stop decay has a large boost and large-R jets are used to reconstruct the top quark decays.
The decayt →χ ± 1 b introduces additional degrees of freedom in the decay. The final-state kinematics is largely driven by the mass separation between the stop and the chargino m(t,χ ± 1 ), and by that between the chargino and the neutralino m(χ ± 1 ,χ 0 1 ). Several signal regions, identified by the 12 The transverse mass m T of the lepton with transverse momentum p T and the missing transverse momentum vector p miss T with magnitude E miss T is defined as and it is extensively used in one-lepton final states to reject SM background processes containing a W boson decaying leptonically. prefix t1L-bC were designed and optimised depending on the mass hierarchy and, consequently, on the different kinematics of the lepton and b-jets.
The four signal regions t1L-bCa_low, t1L-bCa_med, t1L-bCb_med1 and t1L-bCb_high target small values of m(χ ± 1 ,χ 0 1 ) and have the common feature of making use of a dedicated soft-lepton selection: muons and electrons are identified down to a p T threshold of 6 and 7 GeV, respectively, requiring a special treatment for the estimate of possible background processes arising from lepton misidentification. They are collectively referred to as "soft-lepton" signal regions. Both t1L-bCa signal regions require a hard ISR jet to boost the stop pair system and produce a sizeable E miss T . The t1L-bCb targets large values of m(t 1 ,χ ± 1 ) and exploits the presence of two relatively hard b-jets in the event.
The signal region t1L-bCc_diag targets a mass hierarchy complementary to that of the t1L-bCb. The small value of m(t,χ ± 1 ) gives rise to soft b-jets that go undetected, hence b-tagged jets are vetoed for this region.
Topologies arising from scenarios where both m(t,χ ± 1 ) and m(χ ± 1 ,χ 0 1 ) are sizeable are targeted by the three t1L-bCd regions: they all require four jets in the final state, are characterised by different b-jet multiplicities, and apply different selections on the E miss T , m T and am T2 13 variables. A veto on additional isolated tracks and τ lepton candidates identified with loose criteria helps to suppress the dominant SM background from dileptonic tt decays.
The last two signal regions listed in Table 1, t1L-3body and t1L-tNbC_mix, were optimised for two additional possible scenarios. If m(t,χ 0 1 ) < m t and the mass hierarchy or the model parameters suppress the decay through a chargino, then the dominant stop decay ist → bWχ 0 1 , through an off-shell top quark (three-body decay). The dedicated signal region relies on the shape information from the m T and am T2 variable distributions. Finally, t1L-tNbC_mix is designed to recover sensitivity in scenarios where the stop is assumed to decay with similar probabilities to tχ 0 1 and bχ ± 1 : the selection aims to reject the dominant dileptonic tt background by making use of the topness [109] variable.

Two-lepton final states (t2L)
If the SUSY mass hierarchy forbids the presence of sleptons in the stop decay chain, final states containing two leptons (e or μ) and a large amount of E miss T would arise from stop pair production. The main background is given by SM processes containing two W bosons in the final state (mainly tt and W W ) [18]. To discrimi-nate the stop signal from the SM background, the stransverse mass variable m T2 [107,108] is used. The stransverse mass, computed using the two leptons as visible particles and the missing transverse momentum vector, exhibits a kinematical end-point at m W for most SM processes. Because of the presence of additional E miss T due to the LSP, the end-point for a SUSY signal can be at larger values, depending on the mass separation between the particles involved in the decay. The analysis is optimised assumingt 1 →χ ± 1 b with BR = 100 % and m(χ ± 1 ,χ 0 1 ) > m W , but it is also sensitive to the threebody decay mode of the stop. To derive exclusion limits, five signal regions (t2L) have been defined, requiring different jet multiplicities and different m T2 thresholds. A selection requiring two b-jets and based on m T2 computed using them as visible particles is sensitive to the chargino decay mode with m(t 1 ,χ ± 1 ) > m t . Finally, a multivariate discriminant is built which targets thet 1 → tχ 0 1 decay mode.

Final states from compressed stop decays (tc)
If the difference in mass between the stop and the neutralino is smaller than the W boson mass, then the only possible decay channels aret →χ 0 1 c ort → W * b, where the decay products of the off-shell W * would, in general, be soft. This analysis [19] has defined two sets of signal regions, both optimised for thẽ t →χ 0 1 c decay. A common preselection requires the presence of a highp T jet, large E miss T and applies a lepton veto. The first set of signal regions named tc-M, targets scenarios with the stop mass almost degenerate with the neutralino mass, and applies a selection that exploits a monojet-like signature arising from the presence of an ISR jet. Three different signal regions have been designed, characterised by increasing thresholds on the leading jet p T and E miss T . The second set of signal regions, named tc-C, targets less compressed scenarios, and exploits the presence of jets originating from the fragmentation of c-quarks in the final state. A dedicated c-tagging algorithm was used to reject the dominant SM background processes arising mostly from tt and Z → νν (produced in association with heavy-flavour jets) production. As in the case of the tc-M signal regions, different thresholds on the leading jet p T and on E miss T are used to identify a looser and a tighter tc-C region.

Final states with a Z boson (t2t1Z)
A Z boson can be emitted in the decay oft 2 →t 1 Z , producing final states with large lepton multiplicities. It can be useful to look fort 2 (rather thant 1 ) production if, for example, the mass oft 1 is very close to the sum of the top quark and neutralino masses, which would lead tot 1 pair production final states difficult to distinguish from SM tt production. Models are investigated with m(t 1 ,χ 0 1 ) = 180 GeV with the decayt 1 → tχ 0 1 . The final state would contain, beyond the Z boson, several jets arising from thet 1 decay. Similar final states can be obtained in GMSB models where the Z boson is emitted in theχ 0 1 →G Z decay if the gravitinoG is the LSP and the neutralino the NLSP.
This analysis [20] defines five different signal regions divided into two sets. The first set, named t2t1Z-SR2, requires two same-flavour leptons whose invariant mass is consistent with that of a Z boson, m Z , and at least one b-tagged jet. The three signal regions are characterised by the different selection thresholds applied to the E miss T , to the transverse momentum of the dilepton system p T ( ) and to the jet multiplicity. The second set of signal regions, named t2t1Z-SR3, requires three leptons, two of which must form an oppositesign same-flavour pair whose invariant mass is consistent with m Z . Both signal regions require at least five jets, among which at least one has to be b-tagged. The two signal regions are characterised by the different selection thresholds applied to p T ( ) and to the leading lepton p T .
The event selection of b0L-SRA requires large E miss T , exactly two b-jets and vetoes the presence of additional jets; the rejection of the SM ttproduction background is carried out by making use of the contransverse mass [110] of the two b-jets. Its distributions shows a kinematical end-point at about 135 GeV for tt production, while extending to higher values for the signal.
A selection relying on the presence of an ISR jet is instead needed if the third-generation squark mass is almost degenerate with that of the neutralino. This is the purpose of b0L-SRB, which selects a hard, non-b-tagged leading jet recoiling against the squark pair system. The selection includes the requirement of two b-tagged jets, a veto on additional hadronic activity, and the presence of large E miss T . [62] is designed to search for gluino-mediated sbottom and stop production in events with no leptons or one lepton (electron or muon) in the final state. However, it was found to have sensitivity for directb 1 production followed bỹ b 1 →χ 0 2 b →χ 0 1 hb, where h is the SM Higgs boson with mass m h = 125 GeV, and also sensitivity to some of the pMSSM models considered in this paper. Such final states are characterised by a large multiplicity of b-jets both iñ gg →t 1t1 tt andgg →b 1b1 bb where there are up to four b-jets in the final state.

Final states with three b-jets (g3b) This analysis
Three sets of signal regions have been designed to target different mass hierarchies of the gluino-mediated sbottom and stop production models. All signal regions have at least four jets with p T > 30 GeV, three identified b-jets, large E miss T and a large m eff , defined as the scalar sum of the p T of the jets and E miss T .
Strongly produced final states with two same sign or three leptons (SS3L) Final states containing many leptons or same-sign (SS) leptons can arise from the pair production of gluinos and squarks, when the produced particles decay to the LSP through multiple intermediate stages, or when several top quarks appear as part of the decay chain. The analysis was developed for the gluino-mediated stop production processgg →t 1t1 tt followed byt 1 → tχ 0 1 , which can yield final states containing up to four leptons, including SS pairs. Similar final states arise from the sbottom decayb 1 → tχ ± 1 , which are studied in this paper.
This analysis [63] concentrates on final states containing either three leptons or a SS lepton pair produced in association with many jets. Five signal regions (identified by the prefix SS3L) are defined, which are characterised by different light-and heavy-flavour jet multiplicities, high selection thresholds on E miss T and m eff , and different thresholds on the transverse mass of the lepton with the highest transverse momentum and the E miss T .

Spin correlation in tt production events (SC) If the mass of thet 1 is such that mt
ics are similar to that of Standard Model tt production. One possible approach is to derive exclusion limits on the stop mass by performing SM precision measurements. This analysis has measured the azimuthal angle difference between the two leptons arising from the dileptonic tt decay [64].
The events are required to contain, beside the two leptons, at least two additional jets, one of which is required to be btagged. In events containing two leptons of the same flavour, the Z production background is suppressed by applying a selection on the dilepton invariant mass. The distribution of the azimuthal angle between the two leptons is sensitive to the spin correlations of the tt system: it is hence used to extract limits on possible contaminations from direct scalar top production events.
tt production cross section (xsec) The measurement of the tt production cross section using events containing two different-flavour leptons eμ and b-tagged jets is used in Ref. [65] to extract limits on the direct pair production of t 1 with mass close to the top quark. The assumed decay is t 1 → tχ 0 1 . The tt production cross section σ tt is obtained by using the equations where N 1 and N 2 are the number of events with two different flavour leptons having exactly one or two b-tagged jets, respectively, L is the integrated luminosity, eμ the efficiency for a tt event to pass the lepton selection, b is the probability of having a b-jet within acceptance and for it to be tagged, C b is a correlation coefficient which is close to unity, and N bkg 1 and N bkg 2 are the number of events with one or two b-tagged jets from SM events different from tt production. The values of σ tt and b are extracted from the data by solving the two simultaneous Eqs. (2) and (3), avoiding the need to estimate b from simulation.
Stop-pair production events with mt 1 > m t + mχ0 1 have similar eμ and b-jet kinematics to SM tt production events, so the fitted value of b in a combined sample is compatible with that from tt production events alone, and the fitted cross section corresponds closely to the sum of tt and stop-pair production cross sections. Limits on stop pair production are extracted by calculating 95 % CL limits on the stop pair production signal strength μ (defined as the ratio of the obtained stop cross section to the theoretical prediction) based on the comparison of the measured cross section with that predicted for SM tt production events alone. A 95 % CL signal strength smaller than unity for a given signal point implies its exclusion.
This interpretation, which made use of collision data with both √ s = 7 and 8 TeV, is extended here to the three-body decayt 1 → W bχ 0 1 . The main difference with respect to the scenario considered in Ref. [65] is that the three-body decay tends to yield b-jets with lower p T , leading to a fitted b for the combined sample which is different from that expected for tt events alone. The limits obtained are summarised in Fig. 16 for a neutralino mass of 1 GeV. A 95 % CL limit that excludes stop masses below 175 GeV is obtained. The figure also shows the effect on the limit of a "sneaky top squark" scenario [111]: the presence of at 1 with mass similar to that of the top quark could bias the measurement of the top-quark mass itself. The bias in the top-mass measurement introduced by the existence of at 1 with mass mt 1 = 170 GeV depends on the analysis technique and channel, and was evaluated to be at most 1 GeV for the two-and three-dimensional template techniques used in the ATLAS top mass measurement in the lepton+jets channel [112]. The effect of a potential bias of 1 and 2.5 GeV on the top-mass measurement was studied by recalculating the observed 95 % CL limit on μ when reducing the predicted SM tt production cross section from the baseline value of m t = 172.5 ± 1.0 GeV to those obtained for top mass central values of 173.5 and 175 GeV. The corresponding limit on the stop mass is reduced by about 5 and 15 GeV, respectively.
The dependence of the exclusion limits on the neutralino mass was studied and found to be important: the effect of an increasing neutralino mass is to decrease the p T of the b-jets, The sensitivity of the tt cross-section measurement tot 1 pair production assuming a branching ratio of 100 % intõ t 1 → bχ ± 1 , followed byχ ± 1 → W ( * )χ 0 1 with mχ± 1 = 2mχ0 1 was also investigated. The presence of the intermediate chargino state tends to lower the p T of the leptons and of the b-jets significantly, hence decreasing both eμ and b . No exclusion limit can be derived for this scenario. Summarising, the limits on stop pair production obtained in Ref. [65] have been extended by considering the stop threebody decay. Stop masses between 150 GeV and m t can be excluded for a neutralino mass of 1 GeV. The exclusion holds provided that any bias in the top-quark mass measurement by a nearby stop is not significant. Studies indicate that this potential bias would affect the limit on the stop mass by less than 5 GeV.

B.2 Description of the new signal regions
New analyses were developed to target topologies and regions of the SUSY parameter space not well covered by previously published signal regions. They are identified throughout this paper and in Table 1 with the acronym WW, tb and t2t1h. Their contribution to the exclusion limits derived both in simplified and pMSSM models is outlined in Sects. 4 and 5 respectively. In this Appendix, further details about these analyses are provided for the interested reader. Additional informations about selection efficiencies, sensitivities of the different signal regions and individual limit plots, please refer to Refs. [89,90].

B.2.1 Final states with two leptons at intermediate values of m T2 (WW)
The measurement of the production cross section of nonresonant W W pairs in the two-lepton channel at the LHC [113][114][115] has given rise to theoretical speculations [116][117][118] which interpret the possible excess as due to the production of a light stop. The mass hierarchy favoured by these speculations includes at 1 with mass around 200 GeV, aχ ± 1 degenerate with it, and mχ± 1 − mχ0 1 of a few tens of GeV: possible hadronic decay products of thet 1 → bχ ± 1 transition would have low p T and would allow the events to survive the tight jet-veto selections applied in the SM cross-section measurement. Dedicated signal regions, defined by requiring two different-flavour opposite-sign leptons in the final states, are designed to have maximum sensitivity to such scenarios. The approach is also sensitive to scenarios where the stop decays predominantly through the three-bodyt 1 → bWχ 0 1 or four-bodyt 1 → b νχ 0 1 decay. MC simulated events are used to model the signal and to describe all backgrounds that produce two prompt leptons from W , Z or h decay. For processes whose predicted yield in the signal regions is small, or whose topology resembles very closely that of the signal, making it hard to define a proper control region, the background estimate is fully based on MC simulation. For tt, Z + jets and W W production processes, which are the dominant backgrounds, the acceptance of the signal regions selection is estimated with MC simulation, while the normalisation is estimated in dedicated control regions. The MC samples used are the same as in Ref. [18].
The identification criteria for electrons, muons and jets follow the strategy defined in Appendix A: baseline electrons, which are used in the estimation of the fake-lepton background, are selected by applying the "medium" identification criteria. Signal electrons are identified using the "tight" criteria, and they are further required to be isolated. Signal muons correspond to baseline muons with an additional calorimeterand track-based isolation requirement applied. Jets that have |η| < 2.5 and p T > 20 GeV are used for the event selection, although all jets up to |η| < 4.5 are retained for the computation of the missing transverse momentum.
Candidate stop production events, preselected by the same trigger and data quality requirements used in Ref. [18], are further required to contain one electron and one muon of opposite charge, with an invariant mass m > 20 GeV. The leading (in p T ) and next-to-leading leptons are required to have p T > 25 GeV and p T > 20 GeV, respectively.
At this stage of the selection, the background is dominated by production of top-quark pairs and Z → τ τ , followed by W W and W t production.
A requirement of m T2 > 20 GeV, where m T2 is the stransverse mass of the two leptons, strongly reduces the Z → τ τ background, which is expected to have a kinematical end-point at m T2 = m τ . The ratio R 1 of the E miss T and the effective mass, defined as the scalar p T sum of the E miss T , the leptons and the jets, is useful in suppressing the tt background, which is typically characterised by a larger hadronic activity than in signal events. The selection chosen is R 1 > 0.3 + m eff (with m eff in TeV).
After the above selections, the SM background is dominated by W W production. Two differences between this process and the stop pair production signal are further exploited: firstly, the W W production is dominated by quark-antiquark scattering, while stop pair production is mostly initiated by gluon-gluon processes, and secondly the stop pair production signal has four invisible (two neutralinos and two neutrinos) and two undetected (the two b-jets) objects, while the W W process has only two. The first difference implies a higher longitudinal boost of the system emerging from the hard scattering in signal events than in background events. The variable was defined in Ref. [119], and it is an estimator of the boost. The second difference implies a higher E miss T for signal events. This is exploited by making use of Finally, the variable cos θ b , the cosine of the angle between the direction of motion of one of the two leptons and the beam axis in the centre-of-mass frame of the two visible leptons [119], is sensitive to the spin of the produced particles, hence it provides additional rejection power against the W W production process.
A set of seven signal regions were optimised for the discovery of stop pair production, with the stop decaying either ast 1 →χ ± 1 b with a branching ratio of 100 % (assuming mt 1 − mχ± 1 < 10 GeV), or ast 1 → bW ( * )χ 0 1 . The definitions of the signal regions are shown in Table 4.
The background from non-prompt leptons originating from heavy-quark decays or from photon conversions in the signal regions, or from hadrons misidentified as leptons (collectively referred to as fake leptons in the following), is estimated as in Ref. [18].
Specific control regions, whose event yield is expected to be dominated by each of these production processes, are defined and included in the fit to constrain the normalisation parameters. The control region CRT for tt production is defined by changing the following selections with respect to the signal regions: m T2 > 35 GeV, R 1 < 0.3. Its purity is 92 %. The CR for W W production (CRW) is defined by m T2 > 35 GeV, X > 0.04, and has a purity of 72 %. Finally, the CR for Z + jets (CRZ) is defined by m T2 < 20 GeV, 30 GeV < m < 80 GeV, with a purity of 86 %. The normalisation factors of the W W , tt, Z + jets production processes (μ W W , μ tt and μ Z respectively) are determined by a combined profile likelihood fit. When testing the signalplus-background hypothesis for rejection, the fit takes automatically into account the signal contamination in the control regions. For signal scenarios considering light (mt 1 < 150 GeV) stops decaying throught 1 → bW ( * )χ 0 1 , the signal contamination becomes so large that μ W W becomes unrealistically low. For such cases the fit is performed excluding CRW and taking the normalisation of the W W background from MC simulation.
Systematic uncertainties, affecting both the modelling of the detector response (detector-related systematic uncertainties) and the theoretical prediction of the cross sections and acceptances of the background processes (theory-related systematic uncertainties) affect the predicted rates in the signal regions. Their classification and estimation follows closely those defined in Ref. [18]. A few differences, discussed in the following, exist on the estimation of the theory-related uncertainties. The total uncertainty on the yield of the W W production process is composed of three terms: the uncertainty on the NLO hard-scattering calculation is taken to be the difference between the prediction of POWHEG and aMC@NLO both using PYTHIA for the parton shower; the uncertainty addressing the choice of the parton-shower model is estimated as the difference of the aMC@NLO predictions showered either with HERWIG or PYTHIA; the uncertainty due to the choice of the renormalisation and factorisation scale is evaluated by changing the scales independently by a factor of two or one-half and taking the maximum difference. The estimated relative uncertainties on the signal region yields are about 6 % in SR1-SR4 and SR6; 11 % in SR7 and 29 % in SR5. Similar comparisons performed on the W Z and Z Z process yield uncertainties ranging from 30 to 45 % depending on the signal region considered. Additional systematic uncertainties are assigned to the small expected yields from Z + jets production (80 %), W t (50-100 % depending on the SR considered), and non-prompt lepton background. The values of the normalisation factors obtained when performing the fit to the control regions only are shown in Table 5.
The overall predictions of the fit are compared to the data in dedicated validation region that are kinematically close to the signal region. They are defined by applying the preselection requirements of Table 4 with the additional selections shown in Table 6. The m T2 distribution in WW-VR2 and WW-VRT is shown in Fig. 17.
For all signal regions, the expected background yield is dominated by production of W W (35 % in SR1 to 66 % in SR4). Other important background processes are Z + jets in SR1 (20 %), non-prompt leptons in SR2 (12 %), tt in all other SR, with contributions of about 10 %. The distributions of E miss T and m T2 in the signal region WW-SR3 are shown in Fig. 18. Table 6 Summary of the validation regions used in the WW analysis. The preselection requirements of Table 4 are also applied in all three validation regions Table 7 compares the predicted and observed numbers of events in each of the signal regions. No excess above the SM prediction is observed, hence the results are first used to derive model-independent 95 % CL exclusion limits on the minimum number of events beyond the Standard Model in the signal region assuming no signal contamination in the control regions, and then to extract limits on σ vis = σ × ×A, where σ is the cross section for non-SM processes, is the selection efficiency and A is the selection acceptance. These limits are also reported in Table 7. Finally, 95 % CL exclusion limits are derived in specific supersymmetric models of direct pair production of stops. The first exclusion limit (Fig. 19a) is derived in a model where the stop is assumed to decay ast 1 → bχ ± 1 with a branching ratio of 100 %, followed by the decay of the chargino into the neutralino, assumed to be the stable LSP, throughχ ± 1 → W ( * )χ 0 1 . The chargino mass is assumed to satisfy the relation mχ± 1 = mt 1 − 10 GeV, and the limit is derived in the mt 1 -mχ0 1 plane. Stop masses up to about 250 GeV are excluded, almost independently of the neutralino mass. The second limit is derived in a model where thet 1 decays through its three-body or four-body decay (depending on its mass and on that of the neutralino) intõ t 1 → b νχ 0 1 with a branching ratio of 100 %, under the assumption that the decay happens through an off-shell top quark and an on-or off-shell W boson. The limit is shown in Fig. 19b and fills a gap between the exclusions of the t2L and t1L analyses.

B.2.2 Final states containing two top quarks and a Higgs boson (t2t1h)
If the lightest stop has a mass such that m(t 1 ,χ 0 1 ) ∼ m t , the sensitivity of the searches for the production of at 1 pair is greatly reduced. One of the approaches followed is to search for direct pair production oft 2 instead. This is the strategy used, for example, by the t2t1Z analysis, whose signal regions were optimised to detect the decay of a pair-producedt 2 followed by the decayt 2 → Zt 1 .
Inspired by the search for a SM Higgs boson produced in association with a top quark pair, a search was developed and optimised for the decayt 2 → ht 1 , where the Higgs boson is assumed to have SM properties, and thet 1 is assumed to decay ast 1 → tχ 0 1 with a BR of 100 %. The final state is hence characterised by a large jet multiplicity, by the presence of many b-jets from the top quark and Higgs boson decays and by E miss T associated with the presence of neutrinos from semileptonic decays of the top quark and of neutralinos.  The selection of electrons, muons, jets and b-jets follows the principles outlined in Appendix A. The specific choices made for the p T and pseudorapidity thresholds and working points of the final-state objects, as well as the trigger selection, are the same as those in Ref. [120]. The selection requires the presence of exactly one electron or muon with p T > 25 GeV, E miss T > 50 GeV, at least six jets with p T > 25 GeV and |η| < 2.5, of which at least two are required to be b-tagged. The working point chosen for the b-tagging is such that the efficiency to tag b-jets (evaluated on a MC sample of tt production) is about 70 %.
The modelling of the production of tt pairs in association with heavy flavour (tt+HF) is of key relevance in this analysis. A detailed categorisation of tt+HF is made for the purpose of comparisons with different generators and of the propagation of systematic uncertainties on the different heavy-flavour components. The categorisation is also used to reweight the different flavour components of the tt+jets background to obtain a better modelling. These categorisation and reweighting procedures are discussed in detail in Ref. [120]. In particular, the tt + bb component, which is simulated with POWHEG, is reweighted to a full NLO calculation [121] performed in SHERPA 1.4.1+OpenLoops [122,123]. The reweighting is done at generator level using a number of kinematic variables such as the top quark p T , tt system p T , R and p T of the dijet system not coming from the top-quark decay. A different reweighting is applied to the tt + cc and tt+ lightjets components, which is based on the ratio of the differential cross sections at √ s = 7 TeV obtained in data and simulation as a function of the top quark p T and tt system p T [124].
The selected events are categorised into different channels, depending on the number of b-tagged jets (two,  Fig. 20, where the background estimate both before and after the fit is compared to the data in the high-m T region. The dominant post-fit uncertainties are those on the absolute normalisation of the tt + bb and tt + cc processes.
The full list of detector systematic uncertainties considered, discussed in detail in Ref. [120], includes, beside a total uncertainty of 2.8 % on the integrated luminosity, systematic uncertainties on the identification efficiency and energy scale uncertainty of the leptons, reconstruction efficiency and energy scale and resolution uncertainties for jets, b-tagging efficiency and mis-tag rate uncertainties. Further modelling uncertainties are considered, which include, beside production cross-section uncertainties for W/Z +jets, single top and tt, dedicated uncertainties on the NLO calculation of the tt +bb process and on the modelling of the tt +cc component.
No significant excess above the expected background is observed, hence 95 % CL limits are derived in a model wheret 2 production is assumed, followed by the decaỹ t 2 →t 1 h (with a branching ratio of 100 %) andt 1 → tχ 0 1 (again with a branching fraction of 100 %). 14 The limit is derived as a function of thet 2 andχ 0 1 masses, under the assumption that m(t 1 ,χ 0 1 ) = 180 GeV, and it is presented in Sect. 4.3.

B.2.3 Final states containing two b-jets, a charged lepton, and missing transverse momentum (tb)
Several phenomenological models, where bothχ ± 1 andχ 0 1 are lighter than the stop (or the sbottom), allow for thet 1 → tχ 0 1 , t 1 → bχ ± 1 andb 1 → bχ 0 1 ,b 1 → tχ ± 1 decay channels to be open with competing branching ratios. Naturalness arguments require the higgsino mass parameter μ to be smaller than a few hundred GeV, while they impose virtually no constraint on the bino and wino mass parameters M 1 and M 2 . If μ M 1 , M 2 , then the lightest chargino and neutralino masses are both of the order of μ and hence m(χ ± 1 ,χ 0 1 ) is small. Therefore, pair production of stops can lead to Both the leptonic and hadronic decays of the top quark have been studied, and the leptonic channel was found to give a better sensitivity to the signal models of interest. The dominant SM background processes in the signal regions are semileptonic tt and single top production. The SM background is evaluated using a combination of Monte Carlo and partially data-driven techniques. Events are selected online by a trigger requiring the presence of one electron or muon. The online selection thresholds are such that the plateau efficiency is reached for lepton transverse momenta of 25 GeV.
The identification criteria for electrons, muons, jets and E miss T follow the principles outlined in Appendix A. In particular, electrons and muons are required to be isolated: the scalar p T sum of tracks in a cone R = 0.2 around the electron (muon) is required to be smaller than 10 % of the electron transverse momentum (1.2 GeV). The electron or muon track is excluded from the sum. The b-tagging algorithm is used at an operating point with 70 % efficiency in simulated top-quark pair production events. Signal regions are defined as detailed in Table 8, requiring one and only one electron or muon, two b-tagged jets and a large E miss T . Three of the SRs, labelled tb-SRin have no additional jet veto applied, while one of them (tb-SREx1) has a veto requirement on the number of jets (N xjets ) with p T >50 GeV in addition to the two leading b-tagged jets. The final SR optimisation is performed by using selections on the momenta of the objects, the m T and the m eff variables. In addition, the following kinematic variables are used in the event selection: φ b min : the minimum azimuthal distance between the closest b-tagged jet and the E miss T . This variable is used to remove multijet backgrounds with a cut of φ b min > 0.4. m eff : the scalar sum of the p T of the two b-jets (with p T > 25 GeV and |η| < 2.8 ) plus at most one light jet (with p T > 25 GeV and |η| < 2.5) and the E miss T . The number of light jets, n, included in this sum depends on the signal region under study, although n = 1 was mostly used. is also defined to indicate the minimum between m b (1) and m b (2). am T2 : the asymmetric stransverse mass [108] is a kinematic variable which can be used to separate processes in which two decays giving missing transverse momentum occur. It is defined as follows: where p T (v i ) are reconstructed transverse momentum vectors, q T represent the missing transverse momenta from the two decays, with a total missing transverse momentum, p T , and χ is a free parameter representing the unknown neutralino mass, which is assumed to be zero in the calculation. The am T2 variable is calculated with different choices for p T (v 1 ) and p T (v 2 ), depending on the value of m b (n) (n = 1, 2), the invariant mass of the n th b-tagged jet b n and the lepton: The optimisation is carried out using both a pMSSM signal model and simplified models where m(χ ± 1 ,χ 0 1 ) = 5 or 10 GeV. In the case of the pMSSM model, additional non-btagged jets are expected in the final state via the production of other SUSY particles, hence the optimisation points to SRs with no requirement on the N xjets variable (tb-SRIn). In the case of the simplified models, additional jets come only from initial-or final-state radiation, and as a consequence a strict selection on N xjets is applied as in the selection tb-SREx1.
The main SM backgrounds are top-pair production, W production in association with heavy-flavour jets and singletop production. The MC cross section is used to normalise the single-top background and all the other minor SM backgrounds, such as Z +jets, diboson production, tt+W and tt+Z . The normalisation factors of the tt and W + jets backgrounds are determined by a combined profile-likelihood fit. Specific control regions, whose event yield is expected to be dominated by each of these production processes, are defined and included in the fit to constrain the normalisation parameters. The tt control regions (CRT) are defined by inverting the selection on am T2 , requiring am T2 <160 (180) GeV for the inclusive (exclusive) SRs. The purity of the tt process in the CRTs is in excess of 95 %. The W + jets control regions (CRW) are defined by requiring m T <120 GeV. For the control regions corresponding to the tb-SRIn, events with one b-tagged jet are included in the CRW. Top quark pair production dominates the CRWs, with a W + jets purity of 30 % or better. The normalisation factors μ W and μ tt are presented in Table 9. The background model is then validated using validation regions, where little signal contamination is expected.
The distributions of the variable am T2 in the four SRs are shown in Fig. 21 together with the expected distribution from some of the signal models used to optimise the analysis. Table 10 compares the predicted and observed numbers of events in each of the signal regions. No excess above the where σ is the cross section for non-SM processes, is the selection efficiency and A is the selection acceptance. All these limits are also reported in Table 10.
Since the number of events observed agrees with the SM predictions, 95 % CL exclusion limits are derived in specific supersymmetric models of direct pair production of stops. Simplified models were simulated with the two decays t 1 → tχ 0 1 ,t 1 → bχ ± 1 each having a 50 % BR for values of m(χ ± 1 ,χ 0 1 ) = 5, 20 GeV. Furthermore, by using a weighted combination of these simplified models with models corresponding to a 100 % BR in eithert 1 → tχ 0 1 ort 1 → bχ ± 1 , limits can be obtained for any value of the stop BR. Figure 22 shows the exclusion limits for BR(t 1 → tχ 0 1 ) = 25, 50 and 75 % for the two values of m(χ ± 1 ,χ 0 1 ) considered. Finally, 95 % CL exclusion limits are also derived for a natural pMSSM model and are presented in Fig. 23

C Further details of the statistical combination of the t0L and t1L signal regions
This section provides additional details on the combination of the t0L and t1L signal regions targeting scenarios in which the stop decays into eithert 1 → tχ 0 1 or the mixed case wherẽ t 1 → tχ 0 1 andt 1 → bχ ± 1 are both allowed, as discussed in Sect. 4.1 The statistical combination of the two analyses is performed by running the combined fit simultaneously on the control and signal regions of the two analyses. The detector systematic uncertainties are treated as correlated by using, for each of the uncertainties considered, a single nuisance parameter. The supersymmetric signal parameter strength used is the same for the two analyses, while the normalisation parameters for the background processes are kept independent in each analysis. 15 The nuisance parameters associated with modelling uncertainties of the various processes are also kept independent.
The control regions of the two analyses are not mutually exclusive: the events that belong to both a CR of t0L and one of t1L are, at most about 2 % of the total number of events of the t0L CR. The strategy adopted is to remove them from the corresponding t0L CR for the combination. It has been verified that such removal does not affect the individual results of the t0L analysis.

1.58
For each combination performed, the fit setup is validated by checking that the background normalisation parameters obtained are compatible with those obtained separately by the two analyses, by verifying that no additional constraint on the nuisance parameters is introduced with respect to the individual fits, and by checking that no artificial correlation is introduced between any of the fit parameters.
The 95 % CL limit derived from the combination is shown in Fig. 24, where the combined limit is compared to the individual limits obtained by the t0L and t1L analyses independently.

D Signal generation details
Several SUSY models are considered throughout this paper. This section provides the details of how these signal models are generated. For all SUSY models discussed below, the detector response is simulated by passing the generated events through a detector simulation [125] based on GEANT4 [126] or through a fast simulation using a para-metric response to the showers in the electromagnetic and hadronic calorimeters [127] and GEANT4-based simulation elsewhere. All samples are produced with a varying number of simulated minimum-bias interactions overlaid on the hard-scattering event to account for multiple pp interactions in the same or nearby bunch crossings (pileup). The simulation is reweighted to match the number of minimum bias interactions in data, which varies between approximately 10 and 30 interactions in each bunch crossing. Corrections are applied to the simulated samples to account for differences between data and simulation for the trigger and reconstruction efficiencies, momentum scale and resolution of the final-state objects, including the efficiency of identifying jets originating from the fragmentation of b-quarks, together with the probability for mis-tagging light-flavour and charm quarks.

Simplified models
The signal samples for the scenario where both stops decay to a top quark and a neutralino are generated using Herwig++ 2.5.2 [128] interfaced to PYTHIA 6.426 [129]. The neutralino is fixed to be a pure bino, enhancing the decay of thet R component of t 1 to a right-handed top quark. Signal samples where the two stops decay ast 1 → bχ ± 1 are generated with Mad-Graph 5.1.4.8 [130]. For models where the W boson is onshell, thet 1 decay is treated by MadGraph, while if the W is off-shell, PYTHIA is used to decay thet 1 . In these samples, thet 1 is assumed to be mostly at L , and the chargino is assumed to decay throughχ ± 1 → W ( * ) b with a branching ratio of 100 %. Several assumptions about the chargino masses are considered as described in the body of the paper.
Models in which the stop is assumed to decay either ast 1 → tχ 0 1 ort 1 → bχ ± 1 with different branching ratios are obtained by appropriately weighting three samples: one where both stops decay throught 1 → tχ 0 1 , a second one where both stops decay throught 1 → bχ ± 1 , and a third one, where one of the two pair-produced stops decays ast 1 → tχ 0 1 and the other one decays ast 1 → bχ ± 1 . This last sample is generated with MadGraph, thet 1 is assumed to be a maximal mixing oft L andt R . The mass of the chargino in this sample satisfies the gauge-universality relation mχ± The three-body stop decay samples are generated with Herwig++, which performs the matrix element calculation of the three-body decay. The four-body stop decay samples are generated with MadGraph.
For all samples considered, the mass of the bottom quark is fixed to 4.8 GeV and its width is assumed to be zero.
The samples where both stops decay ast 1 → cχ 0 1 are generated with MadGraph, with one additional parton from the matrix element. Similarly to the case of the limit derived as a function of the stop branching ratio intot 1 → tχ 0 1 , the samples with both stops decaying ast 1 → cχ 0 1 and those where both stops decay through the four-body decay are appropriately weighted and combined with a third sample where one stop decays ast 1 → cχ 0 1 and the other decays through the four-body decay to produce a sample of arbitrary branching ratio intot 1 → cχ 0 1 (assuming that t 1 → cχ 0 1 and the four-body decay are the only possible stop decays). Such mixed samples are also generated with MadGraph.
Sbottom pair production samples are also all produced with MadGraph interfaced to PYTHIA, and no more than one additional parton is added to the matrix element. The PDF set used for all signal samples is CTEQ6L1 [131].