Search for supersymmetry in pp collisions 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} TeV in events with a single lepton, jets, and missing transverse momentum

Results are reported from a search for new physics processes in events containing a single isolated high-transverse-momentum lepton (electron or muon), energetic jets, and large missing transverse momentum. The analysis is based on a 4.98 fb−1 sample of proton–proton collisions at a center-of-mass energy of 7 TeV, obtained with the CMS detector at the LHC. Three separate background estimation methods, each relying primarily on control samples in the data, are applied to a range of signal regions, providing complementary approaches for estimating the background yields. The observed yields are consistent with the predicted standard model backgrounds. The results are interpreted in terms of limits on the parameter space for the constrained minimal supersymmetric extension of the standard model, as well as on cross sections for simplified models, which provide a generic description of the production and decay of new particles in specific, topology based final states. Electronic Supplementary Material The online version of this article (doi:10.1140/epjc/s10052-013-2404-z) contains supplementary material, which is available to authorized users.


Introduction
This paper reports results from an updated and improved search for new physics processes in proton-proton collisions at a center-of-mass energy of 7 TeV, focusing on the signature with a single isolated lepton (electron or muon), multiple energetic jets, and large missing momentum transverse to the beam direction (/ E T ). The data sample was collected by the Compact Muon Solenoid (CMS) experiment during 2011 at the Large Hadron Collider (LHC) and corresponds to an integrated luminosity of 4.98 fb −1 , roughly one hundred times larger than the sample used for our previous search [1].

Electronic supplementary material
The online version of this article (doi: 10.1140/epjc/s10052-013-2404-z) contains supplementary material, which is available to authorized users. * e-mail: cms-publication-committee-chair@cern.ch The lepton + jets + / E T signature is prominent in models based on supersymmetry (SUSY) [2][3][4][5][6][7]. In R-parityconserving models [8], SUSY particles are produced in pairs, and their decay chains end with the lightest supersymmetric particle (LSP). In some scenarios, the LSP is a neutralino ( χ 0 ), a heavy, electrically neutral, weakly interacting particle with the properties of a dark-matter candidate [9]. The presence of two such LSPs in each SUSY event typically leads to a large missing transverse momentum, depending on the details of the SUSY mass spectrum. The isolated lepton indicates a weak decay of a heavy particle, such as a W boson or a chargino ( χ ± ). Multiple jets can be produced in complex decay chains of SUSY particles. This signature arises in many SUSY models, including the constrained minimal supersymmetric extension of the standard model (CMSSM) [10,11], and in simplified models [12][13][14][15], which are based on simplified mass spectra and decays of new particles. Both of these frameworks are used to interpret the results. Searches in this or similar channels have been reported by CMS [1,16] and ATLAS [17][18][19].
Searches for SUSY particles are complicated by the presence of standard model (SM) backgrounds that can share many of the features of signal events. In the single-lepton final state, backgrounds arise primarily from the production of tt and W + jets events, with smaller contributions from Z + jets, single-top quark production, and QCD multijet events. In the event topology studied here, a large observed value of / E T in a standard model event is usually genuine, resulting from the production of one or more high-momentum neutrinos. A smaller contribution to events in the high-/ E T tail in this search can arise from the mismeasurement of jets in high cross section processes such as QCD multijet events. To determine the contributions from these backgrounds, we use methods that are primarily based on control samples in data, sometimes in conjunction with specific information from simulated event samples or from additional measurements that provide constraints on the background processes.
Three complementary methods are used to analyze the data, providing valuable cross-checks and probing different signal regions. The Lepton Spectrum (LS) method was used in the CMS single-lepton [1] and opposite-sign dilepton [20] SUSY searches performed using the 2010 data sample. It uses the observed lepton transverse momentum (p T ) spectrum and other control samples to predict the / E T distribution associated with the dominant SM backgrounds. This method is sensitive to SUSY models in which the / E T distribution is decoupled from the lepton p T spectrum, as is the case when two undetected LSPs produce a large missing transverse momentum. The Lepton-Projection Variable (L P ) method uses the L P variable, which was developed for the CMS measurement of the W polarization in W + jets events [21]. This variable, described in Sect. 6, is correlated with the helicity angle of the lepton in the W-boson rest frame. Both the L P and the LS methods take advantage of well-understood properties of the W polarization in tt and W+jets events for the background determination. The methods are complementary in that they rely on significantly different approaches to determining the backgrounds, based on different kinematic variables and different signal regions. The ANN method uses an artificial neural network discriminant built from several kinematic quantities. The ANN discriminant is then used in conjunction with / E T to define signal and sideband regions, from which the background yield is determined. A key variable in the ANN is M T , an approximate invariant mass of the system comprising the lep- This paper is organized as follows. Sections 2 and 3 describe the CMS detector and the event samples. The event preselection requirements that are common to all methods are discussed in Sect. 4. Sections 5, 6, and 7 describe the LS, L P , and ANN methods, respectively, for obtaining SM background estimates from control samples in data. The observed yields in data are compared with the background estimate obtained for each method. Systematic uncertainties are described in Sect. 8. Finally, the results, interpretation, and conclusions of the analysis are presented in Sects. 9 and 10.

The CMS detector
The CMS detector, described in detail in Ref. [22], is a multipurpose apparatus designed to study high-p T physics processes in proton-proton collisions, as well as a broad range of phenomena in heavy-ion collisions. The central element of CMS is a 3.8 T superconducting solenoid, 13 m in length and 6 m in diameter. Within the magnet are (in order of in-creasing distance from the beam pipe) high-precision silicon pixel and silicon strip detectors for charged particle tracking; a lead-tungstate crystal electromagnetic calorimeter for measurements of photons, electrons, and the electromagnetic component of jets; and a hadron calorimeter, constructed from scintillating tiles and brass absorbers, for jet energy measurements. Beyond the magnet is the muon system, comprising drift tube, cathode strip, and resistive-plate detectors interleaved with steel absorbers. Most of the detector systems are divided into subsystems that cover the central (barrel) and forward (endcap) regions. The first level of the CMS trigger consists of custom hardware processors that use information from the calorimeter and the muon system to select up to 100 kHz of the most interesting events. These events are then analyzed in the High Level Trigger (HLT) processor farm, which uses information from all CMS detector systems to reduce the event rate to about 300 Hz.
In describing the angular distribution of particles and the acceptance of the detector, we frequently make use of the pseudorapidity, η = − ln[tan(θ/2)], where the polar angle θ of the particle's momentum vector is measured with respect to the z axis of the CMS coordinate system. The z axis points along the direction of the counterclockwise-moving proton beam; the azimuthal angle φ is measured in a plane perpendicular to this axis. The separation between two momentum vectors in η-φ space is characterized by the quantity R = ( η) 2 + ( φ) 2 , which is approximately invariant under Lorentz boosts along the z axis.

Data and simulated event samples
The data samples used in the analysis were selected using triggers based on / E T , lepton p T , and the transverse momenta (p j T ) of the observed jets j . The overall level of jet activity was measured with the quantity H trigger T = j p j T , the scalar sum of jet transverse momenta satisfying p formed using the GEANT4 package [25]. A variety of Monte Carlo (MC) event generators are used to model the backgrounds. The QCD multijet samples are generated with the PYTHIA 6.4.22 [26] MC generator with tune Z2 [27]. The dominant background, tt, is studied with a sample generated using MADGRAPH 5.1.1.0 [28]. The W + jets and Z + jets processes are also simulated with MADGRAPH. Single-top (s-channel, t-channel, and tW) production is simulated with POWHEG [29]. To model the effect of multiple pp interactions per beam crossing (pileup), simulated events are generated with a nominal distribution of multiple vertices, then reweighted to match the distribution of the number of collision vertices per bunch crossing as measured in data.
The ANN method uses the LM0 model [30] to train the neural network. Because of its large cross section (54.9 pb at NLO), LM0 has already been excluded [1], but its kinematic distributions still provide a reasonably generic description of SUSY behavior with respect to the variables used in the neural network. The parameters for LM0 are m 0 = 200 GeV, m 1/2 = 160 GeV, A 0 = −400 GeV, tan β = 10, and μ > 0.
The results are interpreted in two ways: (i) as constraints on CMSSM parameter space and (ii) as constraints on cross sections for event topologies described in the framework of simplified models. In both cases, a large number of simulated event samples are required to scan over the relevant space of model parameters. For this reason, the scans are performed with the CMS fast simulation package [31], which reduces the time associated with the detector simulation.
Both the LS and L P background determination methods rely on knowledge of the W-boson polarization in W + jets and in tt events. The polarization effects are well modeled in simulated event samples, which are used in conjunction with control samples in data. The angular distribution of the (positively) charged lepton in the W + rest frame can be written as: where f +1 , f −1 , and f 0 denote the polarization fractions associated with the W-boson helicities +1, −1, and 0, respectively. The angle θ * is the polar angle of the charged lepton in the W + rest frame, measured with respect to a z axis that is aligned with the momentum direction of the W + in the top-quark rest frame. The polarization fractions thus determine the angular distribution of the lepton in the W rest frame and, together with the Lorentz boosts, control the p T distributions of the lepton and the neutrino in the laboratory frame.
The W polarization in W + jets events exhibits a more complex behavior than that in tt production. Both CMS [21] and ATLAS [34] have reported measurements of these effects, which are consistent with ALPGEN [35] and MAD-GRAPH [28] simulations predicting that the W + and W − bosons are both predominantly left-handed in W + jets events at high p T . An NLO QCD calculation [36] has demonstrated that the predicted polarization fractions are stable with respect to QCD corrections. As discussed in later sections, this detailed knowledge of the W-boson polarization provides key information for measuring the SM backgrounds using control samples in data. Table 1 summarizes the main variables and requirements used in the event preselection, which is designed to be simple and robust. Except where noted, a common set of preselection requirements is used by each of the three analysis methods. Events are required to have at least one good reconstructed primary vertex, at least three jets (L P method and ANN method) or four jets (LS method), and exactly one isolated muon or exactly one isolated electron. These basic requirements select an event sample that is dominated by genuine, single-lepton events from SM processes. The primary vertex must satisfy a set of quality requirements, including |z PV | < 24 cm and ρ PV < 2 cm, where z PV and ρ PV are the longitudinal and transverse distances of the primary vertex with respect to the nominal interaction point in the CMS detector.

Event preselection
Jets are reconstructed offline using the anti-k T clustering algorithm [37] with a distance parameter of 0.5. The particle four-vectors reconstructed by the CMS particle-flow algorithm [23,24], are used as inputs to the jet clustering algorithm. The particle-flow algorithm combines information from all CMS sub-detectors to provide a complete list of long-lived particles in the event. Corrections based on simulation are applied to the jet energies to establish a uniform response across the detector and a first approximation to the absolute energy scale [38]. Additional jet energy corrections are applied to the data using measurements of energy balance in dijet and photon + jet control samples in data. These additional corrections take into account residual differences between the jet energy scale in data and simulation. The effect of pileup was significant during much of the datataking period. Extra energy clustered into jets due to pileup is taken into account with an event-by-event correction to the jet momentum four-vectors. Jet candidates are required to satisfy quality criteria that suppress noise and spurious energy deposits in the calorimeters. The performance of jet reconstruction and the corrections are described in Ref. [38]. In this analysis, reconstructed jets are required to satisfy p T > 40 GeV and |η| < 2.4. The / E T vector is defined as the negative of the vector sum of the transverse momenta of all the particles reconstructed and identified by the particleflow algorithm.
In the muon channel, the preselection requires a single muon candidate [39] satisfying p T (μ) > 20 GeV and |η| < 2.1. Several requirements are imposed on the elements that form the muon candidate. The reconstructed track must satisfy quality criteria related to the number of hits in the pixel, strip, and muon detectors, and it must have an impact parameter d 0 in the transverse plane with respect to the beam spot satisfying |d 0 | < 0.02 cm and an impact parameter d z with respect to the primary vertex along the z direction satisfying |d z | < 1.0 cm.
To suppress background in which the muon originates from a semileptonic decay of a hadron containing a bottom or charm quark, we require that the muon candidate be spatially isolated from other energy in the event. A cone of size R = 0.3 is constructed around the initial muon momentum direction in η-φ space. The muon combined isolation variable, , is defined as the sum of the transverse energy E T (as measured in the electromagnetic and hadron calorimeters) and the transverse momentum p T (as measured in the silicon tracker) of all reconstructed objects within this cone, excluding the muon. This quantity is used to compute the combined isolation relative to the muon transverse momentum, I comb rel = I comb /p T (μ), which is required to satisfy I comb rel < 0.1. Electron candidates [40] are reconstructed by matching energy clusters in the ECAL with tracks in the silicon tracking system. Candidates must satisfy p T > 20 GeV and |η| < 2.4, excluding the barrel-endcap transition region (1.442 < |η| < 1.56). Quality and photon-conversion rejection requirements are also imposed. The relative isolation variable, defined in a manner similar to that in the muon channel, must satisfy I comb rel < 0.07 in the barrel region and I comb rel < 0.06 in the endcaps. The requirements on d 0 and d z are the same as those used in the muon channel.
The preselection requirements have a large effect on the sample composition. The lepton isolation requirement is critical for the rejection of QCD multijet processes, which have very large cross sections but are reduced to a low level by the isolation and the other preselection requirements. While many lepton candidates are produced in the semileptonic decays of hadrons containing b or c quarks, from π and K decays in flight, and from misidentification of hadrons, the vast majority of these candidates are either within or near hadronic jets. The background from W+jets events (primarily from W → eν or W → μν, but also W → τ ν) is initially also very large. This contribution is heavily suppressed by the three-or four-jet requirement. Depending on the particular signal region, either tt or W + jets production emerges as the largest contribution to the background in the sample of events with moderate to large / E T . Events with a second isolated-lepton candidate satisfying the criteria listed in Table 1 are vetoed. This requirement not only suppresses SM background, but also minimizes the statistical overlap between the event sample used in this search and those used in multilepton searches. However, tt events with dileptons can still be present, and this contribution must be determined, particularly because the presence of two neutrinos in the decay chains can result in large values of / E T . The background involving W → τ ν decays, both from tt events and from direct W production, must also be determined. To help suppress the dilepton background, the requirements on the veto leptons are somewhat looser than those on the signal lepton. For both muons and electrons, the p T threshold is p T > 15 GeV, the isolation requirement is I comb rel < 0.15, and the impact parameter requirement is |d 0 | < 0.1 cm (the d z requirement is kept the same as for the signal lepton). In addition, some of the quality requirements for both the muon and electron are loosened.
Further event selection requirements are used in the individual background estimation methods described in Sects. 5, 6, and 7. The methods use the quantity H T , which is defined as the scalar sum of the transverse momenta of particle-flow jets j with p j T > 40 GeV and |η j | < 2.4, The three background determination methods presented in the following three sections use different approaches to estimating the SM backgrounds using control samples in data. In Sect. 9, we compare the results of the different methods and make some observations about their features.

Overview of the Lepton Spectrum method
This section describes the Lepton Spectrum (LS) method, which is named for the technique used to determine the dominant background source: genuine, single-lepton processes. Such processes account for about 75 % of the total SM background in the signal regions and arise primarily from tt, single-top, and W + jets events. Their contribution to the / E T distribution is estimated by exploiting the fact that, when the lepton is produced in W-boson decay, the / E T distribution is fundamentally related to the lepton p T spectrum, unlike the / E T for many SUSY models. A more detailed description of the Lepton Spectrum method is given in the references [1,41].
Non-single-lepton backgrounds are also determined using control samples in the data. Such events arise mainly from (i) tt dilepton events, in which zero, one, or both of the leptons is a τ and (ii) tt and W + jets events with a single τ → (μ, e) decay. Background from QCD multijet events is expected from simulation to be very small. However, the uncertainties in such simulations are difficult to quantify, because the QCD multijet background in the phase space relevant to this analysis arises from extreme tails of processes with very large cross sections. We therefore use control samples in data to measure the QCD multijet background. Simulated event samples are used for the determination of the Z + jets contribution, which is estimated with sufficient precision to be below one event for most of the signal regions.

Estimation of single-lepton backgrounds
The physical foundation of the Lepton Spectrum method is that, when the lepton and neutrino are produced together in two-body W decay (either in tt or in W + jets events), the lepton p T spectrum is directly related to the / E T spectrum. The lepton and the neutrino share a common Lorentz boost from the W rest frame to the laboratory frame. As a consequence, the lepton spectrum reflects the p T distribution of the W, regardless of whether the lepton was produced in a top-quark decay or in a W + jets event. With suitable corrections, discussed below, the lepton p T spectrum can therefore be used to predict the / E T spectrum for SM single-lepton backgrounds.
The / E T distribution in many SUSY models is dominated by the presence of two LSPs. In contrast to the SM backgrounds, the / E T and lepton p T distributions in SUSY processes are therefore nearly decoupled. The / E T distribution for such models extends to far higher values than the lepton spectrum. Figure 1 shows the relationship between the lepton-p T and / E T distributions in the laboratory frame for two simulated event samples: (i) the predicted SM mixture of tt and W + jets events and (ii) the SUSY LM6 benchmark model. When taken from data, the upper-left region in Fig. 1 (top) provides the key control sample of high-p T leptons from SM processes. This region typically has very little contamination from SUSY events, which populate the high-/ E T region but have relatively low lepton p T values.
The lepton p T spectrum is measured with a muon control sample defined by the preselection criteria and the H T requirements. Unlike the signal region, no / E T requirement is applied, because even a modest one (/ E T > 25 GeV) would bias the high end of the lepton p T spectrum, which is critical for making the background prediction. Only muon events are used as a control sample, because the QCD multijet background is significant in the low-/ E T region of the electron sample. The number of events that are common to both the control sample and the signal region is small. For example, the overlap as measured in simulated tt events is 3.6 % for H T ≥ 750 GeV, / E T ≥ 250 GeV, and p T ≥ 250 GeV. Because no / E T requirement is placed on the muon control sample, a small amount of QCD background remains and must be measured and subtracted. The scaling from the muon to the electron samples is obtained by fitting their ratio in the data over the range 60 ≤ / E T ≤ 250 GeV, with systematic uncertainties evaluated by varying the fit range. The resulting correction factor is N(e)/N(μ) = 0.88 ± 0.03 ± 0.03, where the uncertainties are statistical and systematic, respectively.
To use the lepton spectrum to predict the / E T spectrum in single-lepton SM background processes, three main issues must be understood: (i) the effect of the W-boson polarization in both tt and W + jets events, (ii) the effect of the applied lepton p T threshold, and (iii) the difference between the experimental resolutions on the measurements of lepton p T and / E T . The status of theoretical and experimental knowledge of W-boson polarization in tt and in W + jets events is discussed in Sect. 3. The helicity zero polarization state results in a forward-backward symmetric angular distribution of the lepton and the neutrino in the W rest frame (with respect to the W momentum direction), leading to identical lepton and neutrino spectra in the laboratory frame. In contrast, the helicity ±1 states result in angular asymmetries that lead to somewhat different lepton and neutrino p T spectra in the laboratory frame. These effects are taken into account by applying correction factors obtained from simulation to the measured lepton spectrum, with uncertainties as described in Sect. 8.
The second key issue in the Lepton Spectrum method is the effect of the threshold (p T > 20 GeV) applied to the leptons in both the signal and control samples. Because of the anticorrelation between the lepton p T and the / E T arising from non-zero W-boson helicity states, the threshold requirement removes SM background events in the high-/ E T signal region but not the events in the control sample with high-p T muons that are used to predict the high tail of the / E T spectrum. For the tt background, this effect partially compensates for the bias from the W polarization. For W + jets events, in contrast, the polarization effects for W + and W − approximately cancel, but the lepton p T threshold shifts the predicted yield upward. Correction factors from simulation are used to account for these effects (as well as for polarization effects), which are well defined and understood.
Finally, the resolution on the reconstructed / E T is poorer than that for the lepton p T , so the / E T spectrum is somewhat broadened with respect to the prediction from the lepton spectrum. We measure / E T resolution functions in the data using QCD multijet events obtained with a set of single-jet triggers spanning the range from E T ≥ 30 GeV to E T ≥ 370 GeV. These resolution functions, or templates, quantify the / E T resolution as a function of the number of jets and the H T of the event. These templates are used to smear the measured lepton momenta. Because the templates are taken from data, they include not only the intrinsic detector resolutions, but also acceptance effects. The overall effect of the smearing is modest, changing the background prediction by 5-15 %, depending on the / E T threshold applied. The raw background predictions for the single-lepton background are corrected to account for the effects described above, as well as for the small contamination of the single-lepton control sample arising from dilepton and single-τ events with high-p T leptons. These backgrounds are measured separately, as described below. The overall correction factor is defined such that the single-lepton prediction in a given signal region in simulation matches the yield from single-lepton processes.
The predicted single-lepton background yield varies from about 150 events for the signal region with 250 ≤ / E T < 350 GeV and H T ≥ 500 GeV to about 3 events for the region with / E T ≥ 550 GeV and H T ≥ 1000 GeV. These predictions, as well as the expectations from simulation, are presented in Tables 2, 3, and 4 and discussed in more detail in Sect. 5.4. Table 2 Event yields for the Lepton Spectrum method for H T ≥ 500 GeV. The upper part of the table gives the background predictions that are based on simulated (MC) event samples and the yield for the SUSY signal points LM3 and LM6. The lower part gives the backgrounds predicted using control samples in the data (data-driven prediction). The actual yield observed in data is given at the bot-tom, with the separate muon and electron yields given in parentheses (N μ , N e ) after the total yield. The uncertainties on the background predictions are statistical and systematic. The MC yields are not used in setting limits and are included only for reference. The uncertainties on the MC yields are statistical only

Estimation of non-single-lepton backgrounds
The non-single-lepton backgrounds include dilepton events in several categories, events with W → τ ν followed by τ → decays (in either tt or W + jets events), and QCD multijet processes. These subdominant backgrounds are estimated using control samples in data, in conjunction with information from simulation. The contribution from Drell- Yan and Z + jets is very small and is estimated directly from simulation. Dilepton background events (including the τ as one of the leptons) contain at least two neutrinos, so these events can be important in the tails of the / E T distributions. These backgrounds are divided into the following categories: (i) 2 events with one lost or ignored lepton ( = e, μ), (ii) + τ events with τ → hadrons, and (iii) + τ events with τ → lepton. A lost lepton is one that is either not reconstructed or is out of the detector acceptance. An ignored lepton is one that is reconstructed but fails either the lepton-identification requirements (including isolation) or the p T threshold requirement.
To estimate the background from dilepton events with lost or ignored leptons, we compute the ratio of the combined yield of dilepton events in the ee, eμ, and μμ channels in data to the corresponding combined yield in simulated event samples. This ratio, which is 0.91 ± 0.07 for H T ≥ 500 GeV, 0.93 ± 0.15 for H T ≥ 750 GeV, and 0.87 ± 0.37 for H T ≥ 1000 GeV, is used to rescale the / E T distribution of dilepton events that appear in the signal region in simulation. (Events within 20 GeV of the nominal Z mass are excluded in the e + e − and μ + μ − channels.) This approach is used because the dilepton control sample in data is small, and using it to obtain the shapes of / E T distributions would result in large statistical uncertainties. For all / E T bins above 250 GeV, the predicted yield from this background contribution is less than 6 events, and for all / E T bins above 350 GeV, the yield is at or below 1 event. The / E T distribution associated with the reconstructed dilepton events in data is well described by the simulation.
Dilepton events can also involve τ decays, either τ → hadrons or τ → . The / E T distributions in the dilepton events in data, when suitably modified to reflect the presence of a leptonic or hadronic τ decay, provide an accurate description of the shape of the / E T distribution of these backgrounds. Thus, to estimate the shape from the τ → hadrons background, we effectively replace a lepton in a reconstructed dilepton event with a hadronic τ jet. Both hadronic and leptonic τ response functions are used, providing a probability distribution for a τ to produce a jet or a lepton with a given fraction p T (jet)/p T (τ ) or p T ( )/p T (τ ). These response functions, obtained from simulation, are computed in bins of p T (τ ). This procedure can change the total number of jets above threshold in the event, as well as other properties such as H T and / E T , which are recalculated. Simulated event samples are used to determine, for each of these processes i, the ratio r i = N i feed /N control of the number of events observed in the single-lepton channel to the number of events in the control sample, as a function of / E T . This procedure effectively normalizes all such contributions to the control samples in data. For all / E T bins above 250 GeV, the number of dilepton events with a τ → hadrons decay is predicted to be about 7 events or less and is much smaller in the higher / E T bins. The number of dilepton events with a τ → decay is predicted to be less than 3 events for all / E T bins above 250 GeV and is much smaller in the higher / E T bins.
Estimates for the τ → single-lepton backgrounds from tt and W + jets processes are based on a procedure similar to that used for the dilepton backgrounds, but in this case the single-lepton sample itself is used as the control sample. The / E T distribution obtained by applying the τ → response function to the data is rescaled by a ratio from simulation that gives the yield of τ → background events divided by the yield of events in the single-lepton control sample, as a  Tables 2, 3, and 4, which are discussed in more detail in Sect. 5.4.
Background from QCD multijet events is suppressed to a level well below 1 event in nearly all signal regions, as shown in Tables 2, 3, and 4. The QCD multijet background is determined by first defining a control sample with small missing transverse momentum (/ E T < 50 GeV) and with a lepton impact parameter relative to the beam spot |d 0 | > 0.02 cm. These requirements select a sample with little contamination from other SM processes such as tt and W + jets processes. Using this control sample, we measure the shape of the distribution in the combined relative isolation variable, I comb rel (see Sect. 4). The shape of this distribution has very little correlation to / E T or to the lepton impact parameter (d 0 ), and so can be applied in the high-/ E T signal regions. For each signal region in the data, we determine the background at low values of I comb rel by first scaling the measured QCD multijet background shape in the relative isolation variable to the high-I comb rel sideband of the signal region. The shape is then used to extrapolate the yield to the low-I comb rel signal region. In the high-/ E T signal regions, some non-QCD SM background can be present at high I comb rel , where the QCD background shape is normalized. We therefore subtract the estimated background from tt, W + jets, and Z + jets from this region. These yields are taken from simulation, with systematic uncertainties determined from a comparison with a control region in the data. Tables 2, 3, and 4 compare the background yields predicted from the control samples in data with the yields obtained directly from simulation for H T ≥ 500 GeV, H T ≥ 750 GeV, and H T ≥ 1000 GeV, respectively. We observe that the single-lepton background is the dominant contribution in all regions. The various sources of uncertainties associated with these background determinations are discussed in Sect. 9.   sum over three sources: single-lepton backgrounds (from tt, single-top, and W + jets events), dilepton background from tt, and single-τ events (from both tt and W + jets processes). The vertical span of the green bar corresponds to the statistical uncertainty on the background prediction. (The systematic uncertainties are computed in wider bins used for setting  the limits and are given in Tables 2, 3, and 4.) In each signal region, the blue histogram shows the contribution from the dilepton and single-τ backgrounds only. It is evident that the single-lepton background is dominant in all cases. The / E T distributions for the SUSY benchmark models LM3 and LM6 are overlaid (not summed) for comparison. Systematic uncertainties and the interpretation are presented in Sect. 9.

Overview of the Lepton Projection method
The Lepton-Projection (LP) method uses the difference between SM and SUSY processes in the correlation of the lepton transverse momentum and the missing transverse momentum. As previously discussed, in the SM processes the / E T corresponds to the neutrino in the decay of the W boson, either in W+jets or in tt events. The kinematics of W decays are dictated by the V-A nature of the W coupling to leptons and the helicity of the W boson, as discussed in Sect. 3. Since W bosons that are produced with high transverse momentum in W + jets events exhibit a sizable left-handed polarization, there is a significant asymmetry in the p T spectra of the neutrino and charged lepton. A smaller asymmetry is expected in W bosons from t quark (t antiquark) decays, which yield W bosons which are predominantly longitudinally polarized with smaller left-handed (right-handed) components for W + (W − ).
We have measured the fraction of the helicity states of the W boson using an angular analysis of leptonic W decays [21]. Since the total momentum of the W boson in these decays, and therefore its center-of-mass frame, cannot be accurately determined because the momentum of the neutrino along the beam axis cannot be measured, an observable that depends only on transverse quantities is used. A variable that is highly correlated with the cosine of the polar angle (with respect to the W boson flight direction) in the center-of-mass frame of the W boson is the "lepton projection variable": where p T ( ) is the transverse momentum of the charged lepton and p T (W) is the transverse momentum of the W boson. The latter quantity is obtained from the vector sum of the charged lepton transverse momentum and the missing transverse momentum in the event.
Since SUSY decay chains result in large values of / E T , and often result in relatively low values of the lepton momentum as well, the L P distribution for SUSY events tends to peak near zero, whereas W + jets and tt yield a broad range of L P values. This behavior is illustrated in Fig. 4, which compares the L P distribution from both SM processes and from two representative SUSY benchmark points (LM3 and LM6).
In the L P method, two regions in L P are defined: the region with L P < 0.15 is used as the signal region; the region with L P > 0.3 is used as the control region, i.e., a sample that is depleted in the signal expected and is instead dominated by SM processes. These regions are selected using simulated event samples of W + jets, Z + jets, and tt, that are collectively referred to as electroweak (EWK) processes in what follows, as well as with simulated SUSY events with SUSY particle masses near the region currently under exploration.

Background estimation in the L P method
The key ingredient of the analysis is the estimate of the number of events in the signal region from the SM processes. We define a translation factor, which is the ratio of the number of events in the signal and control regions for the EWK processes. The translation factor is obtained from MC simulation of the EWK processes, and the uncertainties on this factor are included in the systematic uncertainty of the background estimate. In the case of muons, where the background from QCD multijets is negligible, the total number of events predicted from SM processes in the signal region, N pred SM (L P < 0.15), can be determined directly from the number of events observed in the data in the control region, N data (L P > 0.3): In the case of the electrons, the presence of events from QCD multijet processes necessitates an independent evaluation of this background prior to the application of the translation factor for EWK processes. The number of events estimated with this method is then compared to the number of events observed in the data in the signal region, N data (L P < 0.15), for indications of an excess of events over the SM expectation. The analysis is performed in different regions of the event mass scale. To characterize the latter without affecting the correlation of the charged lepton and the neutrino in SM events, the scalar sum of the lepton transverse momentum and the missing transverse momentum, S lep T , is used: In order to make the search optimization less dependent on the unknown energy scale of a new physics signal, the analysis is performed in disjoint ranges of S lep T and the results in these ranges are combined. In addition, the selection is also binned in a second dimension, the H T variable, defined in Eq. (2).
As indicated in Table 1, the event selection used in this analysis is slightly different from the corresponding one in the LS analysis. To increase the sensitivity to SUSY decays, this analysis requires three or more jets. While this results in a significant increase in W + jets events, the additional SM background is mostly concentrated in the control region in L P .
The event yields in the muon and electron channels, as predicted from simulation, are shown in Table 5. As discussed previously, the dominant backgrounds to the lepton plus jets and / E T signature arise from the production and decay of W + jets and tt. The production of single W bosons in association with jets, and with large transverse momenta, is in general the larger of the two, especially at lower jet multiplicities. The majority of the tt background arises from semileptonic tt decays, with fully leptonic tt decays in which a lepton is either ignored or not reconstructed contributing about 20 % of the total tt background.
A source of background, which is not listed in Table 5, stems from QCD multijet events in which a jet is misreconstructed as a lepton. The simulation indicates that the magnitude of this background is small in the control region and negligible in the signal region. Nevertheless, since the uncertainties in simulating these backgrounds can be significant, we use control data samples to estimate the background in the muon and electron channels.
To estimate the background from QCD multijets in the muon final state, we use the relative combined isolation, We employ an additional control data sample, which is specially selected to be enriched in QCD multijets, to determine the ratio of multijet events at low values of the relative isolation. Using this ratio and the number of multijet events expected in the control region of the sample passing the preselection requirements, we estimate the background from multijet events in the signal region to be always smaller than 1 % of the EWK backgrounds. This level of background is negligible and is thus ignored in what follows.
The main sources of electrons in QCD multijet events are misidentified jets and photon conversions. This background is expected to be more substantial than the corresponding one in the muon sample, and its estimate exhibits a large dependence on the details of the simulation. For this reason, we estimate this background from control samples in data. The method relies on the inversion of one or more of the electron identification requirements in order to obtain a sample of anti-selected events, which is dominated by jets misidentified as electrons. We find that the inversion of the requirements on the spatial matching of the calorimeter cluster and the charged-particle track in pseudorapidity and azimuth leaves the relative fraction of the different background sources in QCD multijets unchanged. Moreover, to increase the number of events in this control sample, the requirements on d 0 and d z are removed, while the isolation requirement is loosened. These changes to the event selection have a negligible effect on the L P distribution in the data. In the simulated event samples, it is found that the L P distribution from the control sample events provides a good description of the corresponding distribution from QCD background passing all selection requirements.
The L P distribution obtained with this control sample is used as a template to fit, along with the L P distribution from EWK processes, the L P distribution in the data. In this fit, the EWK template is taken from simulation. This approach, which provides a template obtained from data for the QCD contamination, was applied in the measurement of the polarization of high-p T W bosons [21]. The fit is performed in the control region (L P > 0.3), where the possible presence of signal is highly suppressed. The numbers of QCD and EWK events obtained by the fit are used to estimate the total SM contamination in the signal region (L P < 0.15). The method for estimating the number of SM events expected in the signal region is applied in each range of S lep T and H T . The method for estimating the SM expectation in the signal region is checked using two different control samples, where both the fit and signal regions have a negligible expected SUSY yield. The first sample is defined as all events satisfying the preselection requirements but confined to low values of S lep T : 150 < S lep T < 250 GeV. The method described above is employed to predict the number of events expected in the signal region for both muons and electrons. This prediction is found to be fully consistent with the number of events observed in the data in the signal region. The results of the fits and the yields of QCD and EWK events in the region of low S lep T (<250 GeV) are displayed in Fig. 5 for the electron and muon samples. As can be seen in Fig. 5, the QCD contamination in the signal region, L P < 0.15, is negligible, as expected, since low values of L P favor events with low-p T leptons and high / E T . The second sample, used only for events with muons, is collected with a separate trigger without any requirements on H T or / E T . The muon transverse momentum threshold is raised to p T (μ) > 35 GeV, while the H T threshold is lowered from 500 GeV to 200 GeV and the jet multiplicity requirement is reversed, to be fewer than three jets. Given these requirements on H T and on the jet multiplicity, this control sample is dominated by SM processes. It is found that the estimated background agrees well with the number of events seen in the signal region L P < 0.15.

Results of the L P method
The L P distributions in three ranges of S lep T , are displayed in Fig. 6 for muons (top) and electrons (bottom). Tables 6 and 7 list the numbers of events observed and the number of events expected from all SM processes as presented above, in the signal region, for the muon and electron channels, respectively. The predictions, along with the numbers of events observed in each range of S lep T and H T , are also displayed graphically in Fig. 8 for muons and in Fig. 7 for electrons. The uncertainties quoted in Table 7 correspond to the statistical uncertainty of the fit, while the predictions displayed in Fig. 7 include the total statistical and systematic uncertainty.
All estimates of the total contribution expected from SM processes in the various bins in (S lep T , H T ) are consistent with the numbers of events observed in the data, with no visible excess from a potential SUSY signal. The result is interpreted as a limit in SUSY parameter space in the context of the CMSSM in Sect. 9.

Overview of the method
The Artificial Neural Network (ANN) method uses a multivariate analysis to combine several event characteristics, other than / E T , into a single variable that distinguishes signal from background. Signal events then preferentially populate a signal region in the two-dimensional plane of the ANN output (z ANN ) and / E T , and the sidebands in this plane provide an estimate of the residual background.
Four input variables drive the ANN. The first two are n jets , the number of jets with p T > 40 GeV, and H T , the scalar sum of the p T of each jet with p T > 40 GeV. The SUSY signal typically has heavy particles decaying via complex cascades, and as such, is likely to produce more jets and larger H T than SM backgrounds. The third variable is φ(j 1 , j 2 ), the angle between the two leading p T jets in the transverse plane, which makes use of the greater likelihood that the two highest p T jets are produced back-toback in SM than in SUSY events. The final variable is M T , the transverse mass of the lepton and / E T system. In tt and W + jets events, the lepton and / E T generally arise from the decay of a W boson, and as a result, M T peaks near the W boson mass, with larger values arising only when there are  additional neutrinos from τ or semileptonic decays. By contrast, in SUSY events, M T tends to be greater than the W mass because of / E T due to undetected LSPs. Figure 9 shows the distributions of these variables for simulated SM and SUSY events. The most powerful input variable is M T ; n jets and H T also have considerable discriminating power. The φ(j 1 , j 2 ) variable is weaker, but it still improves the sensitivity of the search. Lepton p T also discriminates between the SM and SUSY, but it is not included in the ANN because its strong correlation with / E T in the SM would spoil the background estimate. Additional variables either do little to improve sensitivity or introduce a correlation between z ANN and / E T . The input variables have similar distributions in the muon and electron channels, so we choose to train the ANN on the two channels combined, and use the same ANN for both. In general, the SM simulation describes the data adequately apart from a possible small structure near 130 GeV in the M T distribution. Reweighting the simulation to match the M T distribution in data does not affect the results of the analysis.
The ANN infrastructure uses standard ROOT utilities [42]. During training, weights are determined that minimize the root-mean-square deviation of background events from zero and signal events from unity. For the SUSY parameter space under study, our sensitivity depends only mildly on the details of the signal sample that trains the ANN. Specifically, for LM points 0 through 13 [30], the sensitivity is comparable (less than 30 % variation) whether the ANN is trained on LM0, LM6 or LM9, even though these three training samples have rather different characteristics. We select LM0 for training because it gives the best overall performance. The SM simulation provides the background sample. Figure 10 compares the distributions of z ANN for data and SM simulation for all events surviving the preselection. The two distributions are consistent within the uncertainties. The SM contribution is concentrated at small values of z ANN , while the LM3 and LM6 SUSY distributions, which are also shown, extend to high values of z ANN where the SM is suppressed.
We define two signal regions in the two-dimensional / E T and z ANN plane. One region, referred to as the "low-/ E T " signal region, has z ANN > 0.4 and 350 < / E T < 500 GeV, while The muon and electron channels are combined the other, the "high-/ E T " signal region, has the same z ANN range, but / E T > 500 GeV. The high-/ E T signal region minimizes the probability that the expected background fluctuates up to a LM6 signal when signal contamination is taken into account. We observe 10 events in the low-/ E T signal region and 1 event in the high-/ E T signal region.

Background estimation using the ANN sidebands
The sidebands in the two dimensional plane of / E T and z ANN provide a strategy for estimating the background. The signal and sideband regions are shown in Fig. 11 and are denoted A, B, C, and D for the low-/ E T signal region and A, B , C, and D for the high-/ E T signal region. The choice of boundaries for the sideband regions balances the competing needs of statistics and insensitivity to signal contamination against preserving similar event compositions in the signal and sideband regions.
The predicted yield in region D is given by where N i is the yield in region i, and the predicted yield in region D is defined similarly. This procedure is equivalent  D and D ), normalized using the yields in regions A and C. We test this estimation procedure using SM simulation: Fig. 12 (top) shows that the / E T distributions for low and high z ANN are similar.
If a signal is present, it enters primarily in the signal regions D and D , but there are also significant contributions relative to the SM in regions B and B , somewhat increasing the predicted backgrounds in D and D . This effect is accounted for in the final results. Table 8 summarizes the event yields in the sideband subtraction regions for the various components of the SM background. The W + jets and tt dominate in all the regions, though their relative proportion varies. The W + jets events are most important at low z ANN since M T , which largely drives z ANN , tends to peak near the W-boson mass. Because the W bosons (and hence their daughters) can be highly boosted, these events extend to very high values of / E T . As seen in Fig. 10, tt events are more likely to have high values of z ANN than are W + jets events; this is because of the presence of dilepton tt events, in which both W bosons (from the top quark pair) decay leptonically, but only one lepton is identified (dilepton ( )), giving large M T . There is also a small contribution from events in which the lepton comes from the decay of a τ produced from a top quark decay, with the other top quark decaying either leptonically (dilepton (τ → )) or hadronically (single τ ). The remaining small backgrounds come from single-top-quark, QCD multijet and Z + jets events.
There are too few events in the simulated QCD multijet and Z + jets samples to populate the high / E T regions (B, B , D and D ). For the results quoted in Table 8 for QCD multijet and Z + jets events, we employ an extrapolation technique based on loosening the z ANN and / E T requirements. The extrapolated numbers for all the regions are consistent with   . Because the SM simulation appears to describe the data well, with, for example, consistent exponential decay constants describing the / E T distributions in the ANN sidebands, we choose to scale the background prediction of the data by κ. The uncertainty in the background from the relative cross sections of SM processes and other effects is quantified in Sect. 8. In the low-/ E T signal region, we expect 9.5 ± 2.2 (stat.) events, and in the high-/ E T signal region 0.7 ± 0.5 (stat.) events. The observed yields are 10 and 1 events, respectively, consistent with the background prediction.

Systematic uncertainties
Systematic uncertainties affect both the background estimates and the signal efficiencies. The sources of systematic uncertainty in the background predictions vary among the three methods, both because the final event selections differ and because the background estimation methods themselves differ. The systematic uncertainties stem from lack of perfect knowledge of the detector response and from uncertainties in the properties of the SM backgrounds. Common uncertainties for all methods are described in Sect. 8.1, while details that are specific to each method are given in Sects. 8.2, 8.3, and 8.4 for the LS, L P , and ANN methods, respectively. Tables 9, 10, and 11 list the main uncertainties associated with each method. The systematic uncertainties affecting the signal efficiency and luminosity, which are largely common to all methods, are described in Sect. 8.5.

Common uncertainties in the background predictions
The jet energy scale (JES) and its effect on / E T in the event can affect the H T and / E T distributions and can also lead to differences between the lepton p T spectrum and / E T spectrum. To understand the effects of energy-scale variations, we vary the jet energy scale as a function of p T and η by amounts determined in independent studies of jet energy scale uncertainties [38], and corresponding to 2 GeV or less for jets with p T > 40 GeV, and then recompute H T and / E T . We also vary the energy scale of "unclustered" calorimeter deposits by 10 % to take into account energy not clustered into jets (this effect is very small).
The uncertainty in the lepton efficiency accounts for differences between data and simulation and uncertainties in the trigger efficiencies. The lepton efficiencies are studied using a sample of lepton pairs with invariant mass close to the Z peak, in which one lepton satisfies tight selection criteria, and the second, reconstructed with relaxed criteria, serves as a probe of the tighter reconstruction and isolation requirements ("tag-and-probe" method [43]). Discrepancies between the data and simulation for electrons are maximal at low p T (10 % effect at around 20 GeV), and we reweight events as a function of lepton p T to quantify the effect. The total lepton efficiency in data is described by simulation with an accuracy of 3 %. Studies of the trigger that separately , and lepton requirements show that the lepton inefficiencies dominate, and amount to 2 % to 3 % for leptons that are reconstructed successfully offline. Muon p T scale uncertainties are obtained from the study of the q/p T (transverse curvature with sign given by the electric charge q) distribution of muons in Z events in data. By comparing the q/p T distribution of positive and negative muons it is possible to quantify the amount of bias in the measurement of q/p T .
The relative amount of tt and W + jets background affects each analysis method through corrections obtained from simulation. The contributions from tt and W + jet have not Table 10 Sources of systematic uncertainties for the L P method and their effects on the background prediction in bins of S lep T for the muon and electron channels. The full list of systematic uncertainties are given for the range 500 < H T < 750 GeV, and the total uncertainties are shown for the two ranges 750 < H T < 1000 GeV and H T > 1000 GeV. The total uncertainty is the individual uncertainties summed in quadrature  been specifically measured in the narrow region of phase space studied in this analysis and their relative contribution must be evaluated. The tt cross section is validated using an algorithm based on the reconstructed top-quark masses for both the hadronic and the leptonic top-quark decays. The uncertainty in the tt cross section is determined by comparing yields in data and simulation after a selection based on top mass variables. The W + jets cross section is validated by comparing event yields between data and simulation in Z + jets events in a dedicated dilepton event selection with similar kinematics. We assign an uncertainty to the W + jets cross section based on the agreement of the data and simulation in the Z + jets sample. Using the uncertainties obtained for the tt and W + jets cross sections, we probe different rel-ative contributions of tt and W + jets events in our sample and the effect on our background predictions. Uncertainties in the polarization fraction for the W boson, either in tt or W + jets events, must be taken into account. For the W polarization in tt events, the theoretical uncertainties are very small (see Sect. 3) and have negligible effect on the background predictions. The W polarization in W + jets events, which is described in more detail in Sect. 3, is more complicated than in tt production. In this case, we consider the effect of conservative variations of the helicity fractions in bins of W-boson p T and η with respect to the theoretical NLO calculations [36].
For the dilepton tt background, tt ( ), the uncertainties are evaluated somewhat differently for the different methods. In the L P and ANN methods this background is evaluated together with the same control sample as for the main single-lepton background prediction. Uncertainties in the prediction can arise from finite detector acceptance, inefficient lepton identification, and cross section uncertainties. In the LS method the dilepton tt background is not predicted using the single-lepton background prediction and separate control samples must be used. Thus the uncertainties for the dilepton tt background are estimated separately and described in the next section.
The small residual QCD multijet background is probed by inverting the requirement on I comb rel or the electron selection criteria to obtain QCD dominated control samples. Contamination from tt and W + jet events in these control samples must be considered and the uncertainties on their cross sections are the dominant uncertainty for these methods.
The Z + jets contribution to the signal regions is very small and uncertainties on this background prediction come from lepton efficiency and cross section uncertainties. In addition, for the LS method there is a small Z+jets contamination to the single-lepton control sample, which must be subtracted, and lepton efficiency and cross-section uncertainties are considered for this as well.

Lepton Spectrum method background prediction uncertainty
For the LS method, the systematic uncertainties for each of the different background predictions from control samples in data (1 , dilepton, 1 τ , QCD, and Z + jets) are included in Tables 2, 3, and 4. To determine the systematic uncertainties for the largest source of background, 1-events (arising from tt, W + jets, and single-top processes), we evaluate deviations for the / E T -dependent correction factor, which is determined from simulation and applied to the 1-background prediction (see Sect. 5.2). Table 9 gives a breakdown of the contributions of the systematic uncertainties for the 1-prediction in bins of / E T and for H T > 750 GeV. The uncertainties in the 1-prediction for the H T > 500 GeV and H T > 1 TeV signal regions are similar to those listed in Table 9. The largest source of uncertainty arises from the potential difference in the muon p T and the / E T scales, because the muon p T spectrum is used to predict the / E T spectrum. The statistical uncertainties in the correction factors (denoted as K-factors in Table 9) for the 1-method are slightly smaller than the combined systematic uncertainty of the correction factor. Table 9 does not include an uncertainty from jet resolution effects because this is taken into account by the smearing of the lepton p T spectra by QCD multijet / E T templates (described in Sect. 5.2). For the purposes of setting limits, the total systematic uncertainty in the 1-background prediction is treated as correlated across all bins in / E T . Tables 2, 3, and 4 also list the non-single-lepton backgrounds, which account for about 25 % of the total, with a relative uncertainty of 5-10 % in the lowest-/ E T bin and about 30 % in the highest-/ E T bin. For the dilepton prediction of lost and ignored leptons (described in Sect. 5.3) the main sources of systematic uncertainty arise from the lepton reconstruction and identification efficiencies and the top-quark p T spectrum. The uncertainties on the lepton efficiencies are described in Sect. 8.1, and the uncertainty associated with the top-quark p T spectrum is determined from varying the fraction of events in the tail of this distribution in simulation in a manner consistent with the uncertainty in this tail as observed in data. This uncertainty is then propagated through the background determination procedure.

Lepton Projection method background prediction uncertainty
For the L P method, the estimate of the total number of events expected from SM processes in the signal region, N pred SM (L P < 0.15), relies on the knowledge of the translation factor, R CS , as well as the number of events observed in the control region, subtracted for the QCD background, N data (L P > 0.3). There are, therefore, two sources of uncertainty in this estimate: uncertainties in the number of events from EWK processes in the control region and uncertainty in R CS . The relative change on the predicted background from each source of systematic uncertainty is listed in Table 10 for both muons and electrons. The largest uncertainty for high S lep T bins is the statistical uncertainty in the data in the control region. The second largest uncertainty comes from the JES uncertainty. The effect from the JES uncertainty is larger in the electron channel, since the JES affects also the shape of the L P distribution used in the fit of the control region. The uncertainty in the resolution of the measurement of the hadronic energy recoiling against the lepton and / E T is evaluated conservatively by smearing the total recoil energy in simulation by an additional 7.5 % along the direction of the recoil and by 3.75 % in the direction orthogonal to the recoil. This decreases the resolution more than 10 % for the high recoils (above 250 GeV) of the signal region and thus covers the difference between data and simulation.

ANN method background prediction uncertainty
For the ANN method, the systematic uncertainty in the background prediction is dominated by the statistics of the simulation, which probes for bias in the background estimation. Another important uncertainty comes from the p T spectrum of the W boson in W + jets events, since it affects the / E T distribution of these events, which preferentially populate the z ANN sideband. To assess the impact, we reweight the p T spectrum of W boson events, using the differences in the p T spectra of Z bosons in data and simulation as a guide. This uncertainty is driven by the statistics of the Z + jets sample. The relative proportions of W+jets and tt events differ in the z ANN signal and sideband regions so the background prediction depends on their relative cross sections. Those tt events with two leptons in the final state, only one of which is observed, have large / E T and are the source of most SM events in the signal region. In addition to the tt cross section, this background depends on lepton acceptance and identification inefficiencies. Additional sources of systematic uncertainty are the hadronic and leptonic energy scales. Table 11 summarizes these uncertainties.

Signal efficiency and other multiplicative uncertainties
The systematic uncertainty in the signal yield arises from the uncertainty in the signal efficiency. In general, this un-certainty is correlated across / E T or S lep T bins. The JES component of the signal efficiency uncertainty is computed separately for each model point in CMSSM and simplified model parameter space and is correlated with the JES uncertainty in the single-lepton background prediction. The systematic uncertainties in the signal efficiency associated with lepton reconstruction and the trigger amount to 3 %. The uncertainty in the integrated luminosity is 2.2 % [44]. The systematic uncertainty in the signal efficiency, not including the JES component, is 6 % for each of the analyses.

Results and interpretation
The LS, L P , and ANN methods each yield SM background predictions that are consistent with the number of events observed in data. We therefore proceed to set exclusion limits on SUSY model parameters. All limits are computed using the modified-frequentist CLs method [45] with a one-sided profile likelihood test statistic. To interpret the absence of an observed signal, three complementary approaches are used.

Constraints on CMSSM parameter space
First, we scan over models in the CMSSM and determine whether the number of events predicted at each model point in parameter space can be excluded by the measurements. This procedure relies on the fact that the CMSSM parameter space can be described with just five parameters, and we fix three of them to commonly used values (A 0 = 0 GeV, μ > 0, tan β = 10). Each model point has a complete SUSY particle spectrum and a well defined cross section, which typically involves several production subprocesses. The CMSSM simulated samples are initially generated using leading-order cross sections. At each point in CMSSM parameter space, the predicted yields for each production subprocess (e.g., gg → g g) are corrected using the NLO cross sections discussed in Ref. [46]. Using the observed yield in data and the predicted background, we determine whether the CMSSM yield for the particular model point can be excluded at the 95 % confidence level (CL). This procedure is complicated by the fact that the control regions in data could potentially be contaminated by signal events. This effect is taken into account for each model by removing the expected contribution to the predicted background arising from signal contamination of the control regions. Figures 14, 15, and 16 show the CMSSM exclusion region [47] for the three background estimation methods, evaluated in the m 1/2 vs. m 0 plane, with the values of the remaining CMSSM parameters fixed at tan β = 10, A 0 = 0 GeV, and μ > 0. Figure 17 displays all of the results together. The excluded regions are below the plotted curves, corresponding to SUSY particle masses below certain values. For reference, the plots display curves of constant gluino and squark masses. The lines of constant gluino mass are approximately horizontal with m( g) ≈ 2.5 m 1/2 . Lines of constant squark mass are strongly curved in the m 1/2 vs. m 0 plane. At low m 0 , the analyses exclude gluinos with masses up to about 1.3 TeV, but the sensitivity falls with increasing m 0 . To determine the one standard deviation (σ ) theoretical uncertainty on the observed limit, the signal yields are recomputed after changing each of the process-dependent SUSY production cross sections at each model point by ±1σ of their uncertainty arising from the parton distribution functions and renormalization and factorization scales [46].

Constraints on simplified model parameter space
The second approach to interpretation is based on the use of simplified models [14,15], which provide a more generic description of new physics signatures. Such models do not include a full SUSY particle spectrum, but instead include only the states needed to describe a particular set of decay chains of interest. Rather than excluding a model, the procedure is to calculate cross section upper limits on a given topological signature. (Such cross section limits can, however, be converted into limits on particle masses within the assumptions of the particular model.) Because simplified models do not describe a full SUSY spectrum, the number of free parameters is small. Furthermore, the parameters are simply the masses of the SUSY particles, in contrast to the grand-unified-theory-scale parameters used in the CMSSM. An advantage of simplified models is that, as a consequence, certain relationships between particle masses that arise with the CMSSM no longer hold, and the spectra can be much more generic.
We consider the "Topology 3 weakino" (T3w) simplified model, which involves the production of two gluinos and their decay via the mechanism shown in Fig. 18. One gluino is forced to decay into two quark jets plus the LSP ( χ 0 ) via the three-body decay g → qq χ 0 , while the other gluino decays via g → qq χ ± , followed by χ ± → W ± χ 0 . The W ± boson can then decay leptonically. The T3w model is specified by masses of the gluino, the LSP ( χ 0 ), and an intermediate chargino ( χ ± ). We calculate cross section limits as a function of M( g), assuming a fixed value for the LSP mass M( χ 0 ) = 50 GeV and setting the chargino mass according to M( χ ± ) = 0.5(M( χ 0 ) + M( g)). The nominal production cross section for the gluino pair production mechanism is given in Ref. [46]. Figure 19 shows the cross sections excluded by each method for this model. The limits fluctuate significantly at low M( g) because of the low signal efficiency in this region.   The data can be interpreted using a third approach, which is applicable to models that do not fall within the scope of either the CMSSM or the simplified model discussed in this section. A model builder can investigate the sensitivity of any one of the three methods presented in this paper to a given signal hypothesis by applying the event selection requirements listed in Table 1, together with the final requirements that define the signal regions. We provide a simple efficiency model for the most important observables used in the event selections. The efficiency model can then be applied to a basic (PYTHIA) simulation of the signal process.
The efficiency model is based on parametrizations of the efficiencies for the event selection requirements with respect to the main reconstruction objects and quantities, such as H T , / E T , and lepton p T . The efficiency of the analysis for a given model can be estimated by applying these individual reconstruction efficiencies, which are given as a function of the most important parameter (such as lepton p T ), to the corresponding kinematic distributions in the model. This procedure would then yield an estimate for the number of signal events from the model. Finally, the sensitivity of the analysis to the model can be obtained by comparing the yield of signal events obtained in this manner with the background yields given in this paper. Kinematic correlations (which can be model dependent) are not taken into account, but this approach nonetheless provides a first approximation to the sensitivity.
The efficiencies for each analysis object are described using "turn-on" curves, which are simply error functions, where x represents the variable most relevant for the reconstruction of the particular object. The error function is parametrized in terms of the plateau efficiency, plateau ; the turn-on threshold, x thresh ; and the characteristic width of the turn-on region, σ . These parameters are obtained by fitting simulated event samples as a function of the true (generated) value.
The selection efficiency associated with the lepton reconstruction, identification, and isolation requirements is estimated as a function of lepton p T by considering muons and electrons (including those from τ decay) generated in the PYTHIA-simulated hard-scattering process. The lepton isolation requirement has a large effect on the efficiency, which consequently depends on the number of jets in the event. To reduce the model dependence arising from this effect, two categories of leptons are considered. First, we assign zero efficiency to leptons that are within R < 0.4 of a quark or gluon with p T > 40 GeV in the hard-scattering process. The efficiency for the remaining leptons is described by a turn-on curve whose parameters are listed in Table 12. The efficiencies are specified for both the lepton selection and for the lepton veto.
The number of jets and the resulting H T value for each event are computed using information available at the generator level. The same clustering algorithm used to reconstruct jets in the data is applied to the generator-level particles. The resulting generator-level jets are required to satisfy R > 0.3 with respect to the leptons described above. The / E T variable is estimated at the generator level from the transverse momenta of neutrinos and any new weakly interacting particles, such as the χ 0 . The parametrizations of the efficiency turn-on curves for the H T and / E T requirements are listed in Tables 13 and 14, respectively. For the requirements used with the LS method, the information given in these tables generally reproduces the efficiency from full simulation to within about 15 %.
In the L P method, the variables L P and S  Table 12 Efficiency-model parameters for lepton efficiencies as a function of x ≡ p T . The leptons are required to lie within the fiducial region and must satisfy the p T thresholds specified in Table 1 Lepton plateau   j 2 ), and M T . The efficiency turn-on curve for z ANN > 0.4 is approximated by the parameter values plateau = 0.98, x thresh = 0.41, and σ = 0.1.
With these additional procedures, the emulation of the efficiencies for the L P and ANN methods is found to be accurate to within ∼15 %, as for the LS method.

Summary
Using a sample of proton-proton collisions at √ s = 7 TeV corresponding to an integrated luminosity of 4.98 fb −1 , we have performed a search for an excess of events with a single, isolated high-p T lepton, at least three jets, and large missing transverse momentum. To provide a robust and redundant determination of the SM backgrounds, three methods are used, each of which relies primarily on control samples in the data.
The Lepton Spectrum (LS) method exploits the relationship between two key observables, the lepton p T distribution and the / E T distribution. In the dominant SM background processes, which have a single, isolated lepton, this connection arises from the fact that the lepton and neutrino are produced together in the two-body decay of the W boson, regardless of whether the W is produced in tt or W + jets events. In many SUSY models, however, the / E T is associated with the production of two neutralinos, which decouples / E T from the lepton p T spectrum. Smaller backgrounds arising from tt dilepton events, from τ → decays in tt or W + jets events, and from QCD multijet processes are also estimated using control samples in the data. In the sample investigated with this method, at least four jets are required, which helps to suppress the background from W + jets events. Nine signal regions are considered, specified by three thresholds on H T and three bins of / E T . The observed yields in each region are consistent with the background estimates based on control samples in the data.
The Lepton Projection (L P ) method exploits information on the W-boson polarization in tt and W + jets events. The dimensionless L P variable itself is sensitive to the helicity angle of the lepton from W decay, but it also provides discrimination between signal and background through the ra- SUSY-like events. The S lep T variable maps out a diagonal line in the plane of lepton p T vs. / E T and reflects the W transverse momentum for the boosted W boson. The L P distributions are studied in bins of S lep T , and H T , and at least three jets are required. In each signal region, the data are in agreement with expectations from the SM.
The artificial neural network (ANN) method provides a means to obtain the / E T distribution of background events in data by constructing a neural network variable z ANN , which has a very small correlation with / E T . This variable also provides strong discrimination between signal and background events, so that the background regions do not suffer from large signal contamination in the models considered. A key element of the z ANN variable is the transverse mass of the lepton-/ E T system, but additional variables, such as the number of observed jets, play a role as well. In the ANN analysis, no excess of events is observed in the signal regions with respect to the SM background prediction.
Because these methods probe extreme kinematic regions of the background phase space, the use of redundant approaches provides confidence in the results. Although the LS and L P methods both make use of information on the Wboson polarization in the background, they are based on different kinematic variables and have different signal regions. The LS method breaks the background into several pieces (single lepton, τ → , dilepton, and QCD) and provides a direct background prediction for the / E T distribution. In contrast, the L P method defines a powerful kinematic variable that is used to obtain a global background prediction by extrapolating an overall background shape from a control region into the signal region. The ANN method similarly uses a global approach to estimating the background. The neuralnet variable incorporates information used in neither of the other two methods.
The results from each method are interpreted in the context of both the CMSSM and a so-called simplified model, T3w, which has a minimal SUSY particle spectrum. The CMSSM limits exclude gluino masses up to approximately 1.3 TeV in the part of the parameter space in which m 0 < 800 GeV, but the bound gradually weakens for larger values of m 0 . For the T3w simplified model, we obtain cross section upper limits as a function of gluino mass. Finally, we provide an approximate model of our signal efficiency that can be used in conjunction with a simple PYTHIA simulation to determine whether other models can be probed by these data.