1 Introduction

Weak-scale supersymmetry (SUSY) [19] is an extension to the Standard Model (SM) that provides a solution to the instability of the scalar SM sector with respect to new high-scale physics. For each known boson or fermion, SUSY introduces a particle with identical quantum numbers except for a difference of half a unit of spin. In the framework of a generic R-parity conserving minimal supersymmetric extension of the SM (MSSM) [1014], SUSY particles are produced in pairs and the lightest supersymmetric particle (LSP) is stable. In a large variety of models, the LSP is the lightest neutralino, \(\tilde{\chi}^{0}_{1}\), which is only weakly interacting. The scalar partners of right-handed and left-handed quarks, \(\tilde{q}_{R}\) and \(\tilde{q}_{L}\), mix to form two mass eigenstates, \(\tilde{q}_{1}\) and \(\tilde{q}_{2}\), with \(\tilde{q}_{1}\) defined to be the lighter one. In the case of the supersymmetric partner of the top quark (\(\tilde {t}\), stop), large mixing effects can lead to one stop mass eigenstate, \(\tilde {t}_{1}\), that is significantly lighter than the other squarks. Depending on the SUSY particle mass spectrum, stop pair production and decay can result in final states topologically similar to \(t\bar{t}\) events.

In this Letter, a search for direct stop pair production is presented in \(\sqrt{s}=7~\mathrm{TeV}\) proton–proton collisions recorded with the ATLAS detector at the Large Hadron Collider, considering a SUSY particle mass hierarchy such that \(m_{t} > m_{\tilde{t}_{1}} > (m_{\tilde{\chi}^{\pm}_{1}} +m_{b})\) and the \(\tilde{t}_{1}\) decays exclusively via \(b+\tilde{\chi}^{\pm}_{1}\). The mass of all other supersymmetric particles are set to be above 2 TeV, and large stop gauge mixing results in \(m_{\tilde{t}_{2}} \gg m_{\tilde {t}_{1}}\) so that only \(\tilde{t}_{1}\) pair production is considered. The stop is predominantly right-handed, but this has little effect on the acceptance and efficiency for the final interpretation. The chargino (\(\tilde{\chi}^{\pm}_{1}\)) mass is set to 106 GeV (above the present exclusion limit of 103.5 GeV [15]) and it is assumed to decay through a virtual W boson (\(\tilde{\chi}^{\pm}_{1} \rightarrow W^{*} \tilde{\chi}^{0}_{1} \)). The choice of chargino mass is identical to that used in a previous study reported by the CDF experiment [16], thus allowing easy comparison of the CDF and ATLAS results. Stops within a mass range between 110 GeV and 160 GeV would be produced with relatively large cross sections—between 245 pb and 41 pb. In this search, dilepton final states (=e, μ) are considered. Although these events could contribute to an anomaly in the measured \(t\bar{t}\) cross section, the relative contribution would be small due to the low transverse momenta of the visible decay products. Events are required to contain at least one energetic jet, large missing transverse momentum (\(E^{\mathrm{miss}}_{\mathrm{T}}\)) and low transverse momenta (p T) leptons, to target the light stop final state.

By targeting very light top squarks, this analysis is complementary to other direct stop searches recently presented by the ATLAS experiment [17, 18].

2 The ATLAS detector

The ATLAS detector [19] is a multi-purpose particle physics detector with a forward-backward symmetric cylindrical geometry and nearly 4π coverage in solid angle.Footnote 1 It contains four superconducting magnet systems, which comprise a solenoid surrounding the inner tracking detector (ID), and the barrel and two end-cap toroids equipping a muon spectrometer. The ID consists of a silicon pixel detector, a silicon microstrip detector (SCT), and a transition radiation tracker (TRT). In the pseudorapidity region \(\left|\eta\right| < 3.2\), high-granularity liquid-argon (LAr) electromagnetic (EM) sampling calorimeters are used. An iron/scintillator tile calorimeter provides coverage for hadron detection over \(\left|\eta\right| < 1.7\). The end-cap and forward regions, spanning \(1.5~<~\left|\eta\right|~<~4.9\), are instrumented with LAr calorimeters for both EM and hadronic measurements. The muon spectrometer surrounds the calorimeters and consists of a system of precision tracking chambers (|η|<2.7), and detectors for triggering (|η|<2.4).

3 Simulated event samples

Monte Carlo (MC) simulated event samples are used to develop and validate the analysis procedure and to evaluate the SM backgrounds in the signal region. Production of top quark pairs is simulated with MC@NLO 4.01 [2022], using a top-quark mass of 172.5 GeV. Samples of W(→ℓν) and Z/γ (→ℓℓ), produced with accompanying jets (of both light and heavy flavour), are obtained with ALPGEN 2.14 [23]. Diboson (WW, WZ, ZZ) production is simulated with HERWIG 6.520 [24] and single top production with MC@NLO 4.01. Fragmentation and hadronisation for the ALPGEN 2.14 and MC@NLO 4.01 samples are performed with HERWIG 6.520, using JIMMY 4.31 [25] for the underlying event. Expected diboson yields are normalised using next-to-leading-order (NLO) QCD predictions obtained with MCFM [26, 27]. The top-quark contribution is normalised to approximate next-to-next-to-leading-order (NNLO) calculations [28]. The inclusive W and Z/γ production cross sections are normalised to the NNLO cross sections obtained using FEWZ [29]. ALPGEN 2.14 and POWHEG  [30] samples are used to assess the systematic uncertainties associated with the choice of generator for \(t\bar{t}\) production, and AcerMC [31] samples are used to assess the uncertainties associated with initial- and final-state radiation (ISR/FSR) [32]. The choice of the parton distribution functions (PDFs) depends on the generator. CT10 [33] sets are used for all MC@NLO samples. MRST LO** [34] sets are used with HERWIG and PYTHIA, and CTEQ6L1 [35] with ALPGEN 2.14. The stop production models are simulated using PYTHIA 6.425 [36]. Signal cross sections are calculated to next-to-leading order in the strong coupling constant, including the resummation of soft gluon emission at next-to-leading-logarithmic accuracy (NLO + NLL) [3739]. An envelope of cross-section predictions is defined using the 68 % C.L. ranges of the CTEQ6.6 (including the α S uncertainty) and MSTW 2008 NLO [40] PDF sets, together with independent variations of the factorisation and renormalisation scales by factors of two or one half. The nominal cross-section value is taken to be the midpoint of the envelope and the uncertainty assigned is half the full width of the envelope, following the PDF4LHC recommendations [41]. All MC samples are produced using a GEANT4-based [42] detector simulation [43]. The effect of multiple proton–proton collisions from the same or different bunch crossings is incorporated into the simulation by overlaying additional PYTHIA minimum bias events onto hard-scattering events. Simulated events are weighted to match the distribution of the mean number of interactions per bunch crossing observed in data.

4 Data and event selection

The analysis uses the full 2011 proton-proton collision data sample. After applying the beam, detector and data-quality requirements, the data sample corresponds to a total integrated luminosity of 4.7 fb−1. Events were triggered using a combination of single and double lepton triggers. The single electron triggers vary with the data-taking period, and the tightest of these has an efficiency of ∼97 % for electrons with p T>25 GeV. The single muon trigger used for all data-taking periods reaches an efficiency plateau of ∼75 % (∼90 %) in the barrel (end-caps) for muons with p T>20 GeV. All efficiencies are quoted with respect to reconstructed leptons, passing the baseline lepton definitions. The double lepton triggers reach similar plateau efficiencies, but at lower p T thresholds (greater than 17 GeV for electrons passing the dielectron trigger, and greater than 12 GeV for muons passing the dimuon trigger; for the electron-muon trigger the thresholds are 15 and 10 GeV for electrons and muons, respectively). If a lepton has an offline p T above the single lepton trigger plateau threshold in a given event, the relevant single lepton trigger is used. Double lepton triggers are used for events with no such lepton. An exception to this rule is applied in the μμ channel. In this case when one lepton has p T>20 GeV and the second p T>12 GeV, a logical OR of both triggers is used to recover efficiency.

Jet candidates are reconstructed using the anti-k t jet clustering algorithm [44] with a radius parameter of 0.4. The input to this algorithm is the three-dimensional energy clusters seeded by calorimeter cells with energy significantly above the noise resulting from the electronics and additional proton–proton interactions (calorimeter clusters). The jet candidate energies are corrected for the effects of calorimeter non-compensation, inhomogeneities and energy loss in material in front of the calorimeter, by using p T- and η-dependent calibration factors based on MC simulations and validated with extensive test-beam and collision-data studies [45]. Furthermore, the reconstructed jet is modified such that the jet direction points to the primary vertex, defined as the vertex with the highest summed track \(p_{\mathrm{T}}^{2}\). Only jet candidates with corrected transverse momenta p T>20 GeV and |η|<4.5 are subsequently retained. Jets likely to have arisen from detector noise or cosmic rays are rejected [45]. Electron candidates are required to have p T>10 GeV, |η|<2.47, and pass the “medium” shower shape and track selection criteria of Ref. [46]. Muon candidates are reconstructed using either a full muon spectrometer track matched to an ID track, or a muon spectrometer segment matched to an extrapolated ID track [47]. They must be reconstructed with sufficient hits in the pixel, SCT and TRT detectors. They are required to have p T>10 GeV and |η|<2.4.

Following object reconstruction, overlaps between candidate jets and leptons are resolved. Any jet candidate lying within a distance \(\Delta R = \sqrt{(\Delta \eta)^{2} + (\Delta \phi )^{2}}=0.2\) of an electron is discarded. Subsequently, any electron or muon candidate remaining within a distance ΔR=0.4 of any surviving jet candidate is discarded.

The measurement of the missing transverse momentum \(\mathbf{p}_{\mathrm{T}}^{\mathrm{miss}}\), and its magnitude \(E_{\mathrm{T}}^{\mathrm{miss}}\), is based on the transverse momenta of all electrons, muons and jets as described above, and of all calorimeter clusters with |η|<4.5 not associated to such objects.

Following overlap removal, electrons are further required to have p T>17 GeV and to pass the “tight” [46] quality criteria, which places additional requirements on the ratio of calorimetric energy to track momentum, and the fraction of high-threshold hits in the TRT. Electrons are also required to be isolated: the p T sum of tracks above 1 GeV within a cone of size ΔR=0.2 around each electron candidate (excluding the electron candidates themselves) is required to be less than 10 % of the electron p T. Muons must have p T>12 GeV and must be isolated: the p T sum of tracks within a cone of size ΔR=0.2 around the muon candidate is required to be less than 1.8 GeV. Jets are subject to the further requirements p T>25 GeV, |η|<2.5 and a “jet vertex fraction”Footnote 2 higher than 0.75.

The top background measurement described below uses a b-tagging algorithm [48], which exploits the topological structure of weak b- and c-hadron decays inside a candidate jet to identify jets containing a b-hadron decay. The nominal b-tagging efficiency, computed from \(t\bar{t}\) MC events, is on average 60 %, with a misidentification (mis-tag) rate for light-quark/gluon jets of less than 1 %. To correct small differences in the b-tagging efficiency observed in the simulation with respect to the data, a scale factor is applied to all simulated samples.

During part of the data-taking period, a localised electronics failure in the electromagnetic calorimeter created a dead region (Δη×Δϕ≈1.4×0.2). For jets in this region, a correction to their energy is made using the energy depositions in the neighbouring cells, and is propagated to \(E^{\mathrm{miss}}_{\mathrm{T}}\). If the energy correction exceeds 10 GeV or 10 % of the \(E^{\mathrm{miss}}_{\mathrm{T}}\), the event is discarded. Events with reconstructed electrons in the calorimeter dead region are also rejected.

Events are subject to the following requirements. The primary vertex in the event must have at least five associated tracks and each event must contain exactly two selected leptons (electrons or muons) of opposite sign. Both of these leptons must additionally satisfy the full list of signal lepton requirements, and the dilepton invariant mass, m ll, must be greater than 20 GeV across all flavour combinations. In addition, events in the signal region must have at least one jet with p T>25 GeV, \(E_{\mathrm{T}}^{\mathrm{miss}}>20~\mbox{GeV}\), missing transverse momentum significanceFootnote 3 \(E^{\mathrm{miss},\mathrm{sig}}_{\mathrm{T}}>7.5~\mbox{GeV}^{1/2}\) to reject multijet events, and leading lepton p T<30 GeV (to provide further rejection of the dominant dileptonic \(t\bar{t}\) background). Events in the ee and μμ channels are subject to a further requirement on the dilepton invariant mass to reject events arising from Z production and decay. This selection, summarised in Table 1, has a low signal efficiency, but strong background rejection. The main factor in the efficiency loss is the lowest lepton p T requirement needed to reach the efficiency plateau of the dilepton triggers. The kinematic acceptance varies between 0.06 % and 0.3 % for a neutralino mass of 55 GeV as the stop mass varies between 112 GeV and 180 GeV, and between 0.1 % and 0.003 % for a stop mass of 140 GeV as the neutralino mass varies between 1 GeV and 95 GeV (the detector efficiency for these points, defined as the efficiency for reconstructing events that already enter the kinematic acceptance, is ∼40 %).

Table 1 Signal region, top control region and Z control region requirements in each flavour channel. The Z veto rejects events with m ll>81 GeV and m ll<101 GeV

5 Background estimation

The dominant SM background, after the signal selection requirements, arises from \(t\bar{t}\) events where both top quarks decay leptonically, with the next most significant background being Z/γ +jets. Single top, W+jets, diboson and multijet events give much smaller expected contributions.

The fully leptonic \(t\bar{t}\) background in the signal region is obtained by extrapolating the number of \(t\bar{t}\) events measured in a suitable control region (CR), after correcting for contamination from non-\(t\bar{t}\) events, into the signal region (SR). This extrapolation, detailed in Eq. (1), uses the ratio of the number of simulated \(t\bar{t}\) events in the signal region to those in the control region:

$$ (N_{t\bar{t}})_{\mathrm{SR}}= \bigl[ (N_{\mathrm{data}})_{\mathrm{CR}}-(N_{\mathrm{non}\mbox{-}t\bar{t},\mathrm{MC}})_{\mathrm{CR}} \bigr] \frac{(N_{t\bar{t},\mathrm{MC}})_{\mathrm{SR}}}{(N_{t\bar{t},\mathrm{MC}})_{\mathrm{CR}}} $$
(1)

The CR is designed to give an event sample dominated by top events, whilst minimising signal contamination. It is further chosen to be kinematically similar to the signal region to minimise systematic uncertainties due to extrapolation. Selection requirements for the top control region are summarised in Table 1. In this analysis, models with small stop–chargino mass difference are considered, and hence soft b-jets are expected in the signal events which are not efficiently tagged. By requiring a b-jet in the top control region a high-purity sample of top events is obtained. The signal contamination in the considered models is typically of the order of a few per cent, rising to 30 % for models with \(m_{\tilde{\chi}^{0}_{1}}=1~\mbox{GeV}\) and high \(m_{\tilde{t}_{1}}\). The percentage of SM, non-\(t\bar{t}\) events in the CR is less than 5 % across all channels. The resulting \(t\bar{t}\) background contributions are consistent with the expected MC yields in all channels within the uncertainties. Signal contamination is taken into account when setting the exclusion limit in the next section by including, for each signal model, the expected signal yield in the top control region in the \((N_{\mathrm{non}\mbox{-}t\bar{t},\mathrm{MC}})_{\mathrm{CR}}\) term in Eq. (1).

The contribution from Z/γ +jets events to the signal region (from ee and μμ events) is evaluated in a similar way. Data are used to obtain the normalisation of the Z/γ background in a suitable CR and MC is used to extrapolate from CR to SR using an equation analogous to Eq. (1). This method is used separately for each of the ee and μμ channels (with selection requirements for the Z CR as summarised in Table 1), whereas the contribution to (including those from Z/γ ττ) is taken directly from the MC simulation due to the limited number of events in the CR. The contamination from non-Z/γ +jets SM events in the CR is less than 5 %, and the signal contamination less than 4 %. The resulting Z background contributions are consistent with the expected MC yields in the ee and μμ channels within the uncertainties. The effect of signal contamination of the Z control region on the final exclusion limit can be neglected to a very good approximation.

Single top, W+jets (including heavy-flavour contributions) and diboson backgrounds are evaluated in the signal region directly from the MC simulation. The estimated contribution from W+jets has been cross-checked using a data-driven technique (an extension of the “template fit”, described below), and found to be in good agreement.

The tight requirement on \(E^{\mathrm{miss},\mathrm{sig}}_{\mathrm{T}}\) heavily suppresses the multijet background. A data-driven template fit technique is used to verify that this background is small, and to assign an uncertainty on the yield in the signal region. The isolation requirements on the electrons and muons are reversed to enhance the multijet content of selected events. The requirements are inverted in the signal region, prior to application of the \(E^{\mathrm{miss,sig}}_{\mathrm{T}}\) requirement. The shape of the \(E^{\mathrm{miss,sig}}_{\mathrm{T}}\) distribution in data for this inverted selection (after subtracting the dominantly electroweak background using the MC simulation) is then compared to the equivalent distribution in data for the “normal” isolation requirements in order to validate that inverting the lepton isolation does not distort the shape of the distribution. The “normal” and “inverted” shapes were found to agree very closely for the full range of distributions considered in the analysis. The inverted \(E^{\mathrm{miss},\mathrm{sig}}_{\mathrm{T}}\) distribution is then renormalised to match the distribution after nominal isolation requirements. Passing this correctly normalised template through the remaining requirements gives the multijet yield in the signal region. It is found to be small in all channels, making up less than 2 % of the total background.

6 Systematic uncertainties

The total systematic uncertainty on the expected background in the combined flavour channel (the sum of ee, and μμ events) is 9.8 %, and is dominated by the uncertainties on the two largest backgrounds (dileptonic \(t\bar{t}\) and Z+jet events). The largest source of systematic uncertainty on the \(t\bar{t}\) background evaluation is the uncertainty on the jet energy scale (JES), with smaller contributions coming from the jet energy resolution (JER) uncertainty [45], the theory and MC modelling uncertainties (using the prescriptions described in Ref. [49]), the systematic uncertainties on the b-tagging efficiency [48], and the uncertainty arising from the limited numbers of MC and data events. Uncertainties [46, 50, 51] in lepton reconstruction and identification (momentum and energy scales, resolutions and efficiencies) give smaller contributions.

The primary source of uncertainty on the Z/γ +jets background estimate in the combined flavour channel is the jet energy resolution uncertainty, with smaller contributions coming from the statistical and jet energy scale uncertainties. Theoretical uncertainties on the Z/γ +jets background are investigated by varying the PDF and renormalisation scales. An uncertainty on the luminosity of 3.9 % [52, 53] is included in the systematic uncertainty calculation for backgrounds taken directly from the MC simulation. The dominant uncertainties on these backgrounds are the jet energy scale and statistical uncertainties. The systematic uncertainty on the multijet yield is obtained by varying the range in which the template fit is performed, and using the maximum deviation of the final yield to assign the uncertainty.

In the considered \(m_{\tilde{\chi}^{0}_{1}}\)\(m_{\tilde{t}_{1}}\) mass plane the theoretical uncertainty on each of the signal cross sections is approximately 16 %. These arise from considering the cross-section envelope defined using the 68 % C.L. ranges of the CTEQ6.6 and MSTW 2008 NLO PDF sets, and independent variations of the factorisation and renormalisation scales (see Sect. 3). Further uncertainties on the numbers of predicted signal events arise from the JES uncertainty (7–15 %), the JER uncertainty (1–7 %), the luminosity uncertainty (3.9 %), the uncertainties on calorimeter energy clusters used to calculate \(E^{\mathrm{miss}}_{\mathrm{T}}\) (2–6 %), the statistical uncertainty from finite MC event samples (4–20 %) and smaller contributions from uncertainties on lepton reconstruction and identification, where the quoted ranges display the maximum variation observed using all signal models considered in this analysis.

7 Results and interpretation

Table 2 shows the data observations in the signal regions in each flavour channel, and in the combined flavour channel, along with the evaluated background contributions. Good agreement is observed across all channels, and the absence of evidence for light scalar top production allows a limit to be set on the visible cross section for non-SM physics, σ vis=σ×ϵ×A, for which this analysis has an efficiency ϵ and acceptance A. The limits are calculated using the modified frequentist CLs prescription [54] by comparing the number of observed events in data with the SM and SM-plus-signal expectations.

Table 2 The expected and observed numbers of events in the signal region for each flavour channel. In the combined flavour column (“all”), the statistical uncertainty (first uncertainty quoted, includes the MC and data statistical errors) on the various background estimates have each been added in quadrature whilst the systematic uncertainties (second uncertainty quoted) have been combined taking into account the correlations between background sources. Observed and expected upper limits at 95 % confidence level on the visible cross section σ vis=σ×A×ϵ are also shown. The expected signal yields and statistical uncertainties on the yields are quoted for the two mass points illustrated in the figures

All systematic uncertainties and their correlations are taken into account via nuisance parameters using a profile likelihood technique [55]. In Fig. 1, the leading lepton p T distributions in the ee and μμ channels are illustrated along with the \(E_{\mathrm{T}}^{\mathrm{miss}}\) and \(E^{\mathrm{miss},\mathrm{sig}}_{\mathrm{T}}\) distributions of the data and simulated events in the signal region (with the background normalisations set to their nominal values).

Fig. 1
figure 1

The leading electron and muon p T distributions in the same flavour (aee and (bμμ channels, before the requirement on the leading lepton p T (which is marked by the dashed vertical line), and (c) the \(E^{\mathrm{miss}}_{\mathrm{T}}\) distribution and (d\(E^{\mathrm{miss},\mathrm{sig}}_{\mathrm{T}}\) distribution after all signal region requirements. The data and evaluated background components are shown. The hashed band indicates the total experimental uncertainty on the expectation. The dashed lines give the expectations for signal models with stop masses of 112 GeV and 160 GeV, and a neutralino mass of 55 GeV. The last histogram bins in (a) and (b) include the integrals of all events with p T>300 GeV. The final bin in (d) includes all events with an \(E^{\mathrm{miss,sig}}_{\mathrm{T}}\) of at least 15 GeV1/2. The bottom panels show the ratio of the data to the expected background (points) and the systematic uncertainty on the background (hashed area)

The observed data yield is in good agreement with the SM prediction in the combined flavour channel given in Table 2.

The results in the combined channel are used to place exclusions at 95 % confidence level in the \(m_{\tilde{t}_{1}}\)\(m_{\tilde{\chi}_{1}^{0}}\) mass plane, using the CLs method. The resulting 95 % confidence level expected (dashed) and observed (solid) limits are shown in Fig. 2. Neutralino masses down to 1 GeV are considered, since there is no LEP limit on the neutralino mass in the MSSM for the case that the lightest neutralino is predominantly bino in nature. A bino-dominated lightest neutralino is favoured by the recent LHC Higgs search results.

Fig. 2
figure 2

95 % exclusion limit in the \(m_{\tilde{t}_{1}}-m_{\tilde{\chi }_{1}^{0}}\) mass plane, with \(m_{\tilde{\chi}_{1}^{\pm}}=106~\mathrm{GeV}\). The dashed and solid lines show the 95 % C.L. expected and observed limits, respectively, including all uncertainties except for the theoretical signal cross-section uncertainty (PDF and scale). The band around the expected limit shows the ±1σ result. The dotted ±1σ lines around the observed limit represent the results obtained when moving the nominal signal cross section up or down by the theoretical uncertainty. Illustrated also is the region excluded at the 95 % C.L. by CDF [16], where the lowest neutralino mass considered was 44 GeV, indicated by the horizontal dotted line

The observed limits represent a significant extension of the CDF limit [16] for a chargino mass of 106 GeV to smaller chargino minus neutralino mass difference (the ATLAS limit extends up to a neutralino mass of 70 GeV for a stop mass of 130 GeV, whilst the CDF limit extends up to a neutralino mass of 46 GeV).

The limit on the stop mass for neutralino masses of 45 GeV (135 GeV) is comparable to the equivalent CDF limit. Increasing the chargino mass by 15 GeV leads to a modest shift of the exclusion limit to higher values of the neutralino mass, with the reach in stop mass being enhanced to a lesser degree due to the falling stop production cross section. For example, for a model with \(m_{\tilde{t}_{1}}=130~\mbox{GeV}\), \(m_{\tilde{\chi }^{\pm}_{1}}=120~\mbox{GeV}\) and \(m_{\tilde{\chi}^{0}_{1}}=60~\mbox{GeV}\), the value of A×ϵ (0.032 %) is slightly higher than for the equivalent model with \(m_{\tilde{\chi}^{\pm}_{1}}=106~\mbox{GeV}\) (0.026 %).

8 Conclusions

A search for light top squarks has been performed in the dilepton final state. SM backgrounds have been evaluated using a combination of data-driven techniques and MC simulation. Good agreement is observed between data and the SM prediction in all three flavour channels. The results are interpreted in the \(m_{\tilde{t}_{1}}\)\(m_{\tilde{\chi }_{1}^{0}}\) plane with the chargino mass set to 106 GeV, and with the assumption that the decay \(\tilde{t}_{1}\rightarrow b\tilde{\chi}^{\pm}_{1}\) occurs 100 % of the time, followed by decay via a virtual W (\(\tilde {\chi}^{\pm}_{1}\rightarrow W^{*}\tilde{\chi}^{0}_{1}\)) with an 11 % branching ratio (per flavour channel) to decay leptonically. A lower limit at 95 % confidence level is set on the stop mass in this plane using the combined flavour channel. This excludes stop masses up to 130 GeV (for neutralino masses between 1 GeV and 70 GeV). This limit exceeds that set by the CDF Collaboration for the same scenario [16].