Search for supersymmetry in events with one lepton and multiple jets in proton-proton collisions at s p = 13 TeV

A search for supersymmetry is performed in events with a single electron or muon in proton-proton collisions at a center-of-mass energy of 13 TeV. The data were recorded by the CMS experiment at the LHC and correspond to an integrated luminosity of 2 . 3 fb − 1 . Several exclusive search regions are defined based on the number of jets and b -tagged jets, the scalar sum of the jet transverse momenta, and the scalar sum of the missing transverse momentum and the transverse momentum of the lepton. The observed event yields in data are consistent with the expected backgrounds from standard model processes. The results are interpreted using two simplified models of supersymmetric particle spectra, both of which describe gluino pair production. In the first model, each gluino decays via a three-body process to top quarks and a neutralino, which is associated with the observed missing transverse momentum in the event. Gluinos with masses up to 1.6 TeVare excluded for neutralino masses below 600 GeV. In the second model, each gluino decays via a three-body process to two light quarks and a chargino, which subsequently decays to a W boson and a neutralino. The mass of the chargino is taken to be midway between the gluino and neutralino masses. In this model, gluinos with masses below 1.4 TeV are excluded for neutralino masses below 700 GeV.


I. INTRODUCTION
Supersymmetry (SUSY) [1-8] is a well-motivated theoretical framework that postulates new physics beyond the standard model (SM). Models based on SUSY can address several open questions in particle physics, e.g. the cancellation of quadratically divergent loop corrections when calculating the squared mass of the Higgs boson. In Rparity [9] conserving SUSY models, the lightest SUSY particle (LSP) is stable and can be a viable dark matter candidate. An inclusive search for SUSY in the singlelepton channel was performed with 13 TeV data recorded in 2015 by the CMS experiment at the CERN LHC, corresponding to an integrated luminosity of 2.3 fb −1 . Similar searches were performed in 7 TeV [10][11][12] and in 8 TeV [13][14][15] data by the CMS and ATLAS experiments. First results in the single-lepton final state at 13 TeV are also available from both collaborations [16][17][18]. In this paper, we present a search for gluino pair production designed to be sensitive to a variety of SUSY models.
In this analysis, the main backgrounds arise from W þ jets events and top quark-antiquark (tt þ jets) events, which also lead to W-boson production. In W þ jets events, or in tt þ jets events with a single leptonic W-boson decay, the missing transverse momentump miss T , defined as the negative vector sum of the transverse momenta of all reconstructed particles in the event, provides a measurement of the neutrino transverse momentum. The quantitỹ p l T þp miss T , wherep l T is the lepton transverse momentum vector, corresponds to the transverse momentum of the W boson in background events of this type. We also define the magnitude of the missing transverse momentum, E miss T ¼ jp miss T j, and the sum L T ¼ p l T þ E miss T , where p l T is the magnitude ofp l T . A key analysis variable is the azimuthal angle ΔΦ, measured in the plane perpendicular to the beams, betweeñ p l T andp l T þp miss T . In background events with a single W-boson decay, ΔΦ corresponds to the azimuthal angle between the transverse momentum vectors of the charged lepton and the W boson. In such events, the distribution of ΔΦ falls rapidly and has a maximum value determined by the mass and transverse momentum of the W boson. The higher the boost of the W boson, the smaller the maximum value of ΔΦ. In SUSY events corresponding to our signal models, however, E miss T typically receives a large contribution from the missing momentum of the two neutralino LSPs. As a consequence, the ΔΦ distribution in signal events is roughly uniform. The main backgrounds can therefore be suppressed by rejecting events with a small value of ΔΦ. The primary remaining background arises from tt þ jets production, where both W bosons decay into a charged lepton and a neutrino, with one lepton being not well identified or falling outside the detector acceptance. This background populates the high region of ΔΦ.
Since many models of gluino pair production lead to final states with a large number of jets, the signal-tobackground ratio is very small in regions with low jet multiplicity. We therefore restrict the search to regions of large jet multiplicity and use low jet multiplicity regions, dominantly populated by events from SM processes, to estimate the background. Exclusive search regions are characterized by the number of jets (n jet ), the number of b-tagged jets (n b ), the scalar sum of the transverse momenta p T of the jets (H T ), and L T .
The results are interpreted in terms of simplified models [19][20][21][22] of gluino pair production. In the first model, designated T1tttt and shown in Fig. 1 (left), gluinos are pair produced and subsequently undergo three-body decays to tt þχ 0 1 , whereχ 0 1 is the lightest neutralino. In the second model, termed T5qqqqWW and shown in Fig. 1 (right), the gluinos undergo three-body decays to a quark-antiquark pair (qq) from the first or second generation and a chargino (χ AE 1 ). The chargino mass is taken to be m~χAE 1 ¼ 0.5ðmg þ m~χ0 1 Þ. The chargino then decays to a W boson and theχ 0 1 , where the W boson can be virtual, depending on the mass difference between the chargino and the lightest neutralino.
The organization of this paper is as follows. Section II describes the CMS detector. The event reconstruction and selection are discussed in Secs. III and IV, respectively. The background estimations are given in Sec. V. An overview of the main systematic uncertainties is presented in Sec. VI. The results are discussed and interpreted in Sec. VII, and a summary is given in Sec. VIII.

II. THE CMS DETECTOR
The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. A silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter (HCAL), each composed of a barrel and two end-cap sections, reside within the solenoid volume. Forward calorimeters extend the pseudorapidity (η) [23] coverage provided by the barrel and end-cap detectors. Muons are measured in the range jηj < 2.4, with detection planes made using three technologies: drift tubes, cathode strip chambers, and resistive plate chambers.
The silicon tracker measures charged particles within the range jηj < 2.5. Isolated particles with transverse momenta p T ¼ 100 GeV, emitted at jηj < 1.4, have track resolutions of 2.8% in p T , and 10 ð30Þμm in the transverse (longitudinal) impact parameter [24]. The ECAL and HCAL measure energy depositions in the range jηj < 3, with quartz fiber and steel forward calorimeters extending the coverage to jηj < 5. When information from the various detector systems is combined, the resulting jet energy resolution is typically 15% at 10 GeV, 8% at 100 GeV, and 4% at 1 TeV [25]. The momentum resolution for electrons with p T ≈ 45 GeV from Z → ee decays ranges from 1.7% for electrons that do not shower in the barrel region to 4.5% for electrons that shower in the end caps [26]. Matching muons to tracks measured in the silicon tracker yields relative transverse momentum resolutions for muons with 20 < p T < 100 GeV of 1.3%-2.0% in the barrel, and less than 6% in the end caps. The p T resolution in the barrel is below 10% for muons with p T up to 1 TeV [27]. The CMS trigger system consists of two levels, where the first level (L1), composed of custom hardware processors, uses information from the calorimeters and muon detectors to select the most interesting events in a fixed time interval of less than 4μs. The high-level trigger (HLT) processor farm further decreases the event rate from around 100 kHz to less than 1 kHz, before data storage.
A more detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [23].

III. EVENT RECONSTRUCTION AND SIMULATION
All objects in the event are reconstructed using the particle-flow event reconstruction algorithm [28,29], which reconstructs and identifies each individual particle through an optimized combination of information from the various elements of the CMS detector. The energy of electrons is determined from a combination of the electron momentum at the primary interaction vertex as determined by the tracker, the energy of the corresponding ECAL cluster, and the energy sum of all bremsstrahlung photons spatially compatible with originating from the electron track [26]. Electron candidates are required to satisfy identification criteria designed to suppress contributions from misidentified jets, photon conversions, and electrons from heavyflavor quark decays. Muons are reconstructed using a stand-alone muon track in the muon system serving as a seed to find a corresponding track in the silicon detector [27]. Additional criteria include requirements on the track and hit parameters. Events are vetoed if additional electrons or muons with looser identification requirements are found. Depending on the mass difference between the chargino (χ AE 1 ) and the neutralino (χ 0 1 ), the W boson can be virtual.
The degree of isolation of a lepton from other particles provides a strong indication of whether it was produced in a hadronic jet, such as a jet resulting from the fragmentation of a b quark, or in the leptonic decay of a W boson or other heavy particle. Lepton isolation is quantified by performing a scalar sum of the transverse momenta of all particles that lie within a cone of specified size around the lepton momentum vector, excluding the contribution of the lepton itself. To maintain high efficiency for signal events, which typically contain a large number of jets from the SUSY decay chains, we use a p T -dependent cone radius R ¼ ð0.2; 10 GeV=p T ½GeV; 0.05Þ for ðp T < 50 GeV;50 GeV < p T < 200 GeV;p T > 200 GeVÞ, respectively. The isolation variable is defined as a relative quantity, I rel , by dividing this scalar sum by the p T of the lepton. For selected muons or electrons, we require I rel < 0.2 and I rel < 0.1, respectively, while for additional leptons used in the event veto, we require I rel < 0.4. When computing the isolation variable, an areabased correction is applied to remove the contribution of particles from additional proton-proton collisions within the same or neighboring bunch crossings (pileup).
The energy of charged hadrons is determined from a combination of their momenta measured in the tracker and the matching ECAL and HCAL energy depositions, corrected for zero-suppression effects in the readout electronics, and for the response function of the calorimeters to hadronic showers. Finally, the energy of neutral hadrons is obtained from the corresponding corrected ECAL and HCAL energies.
Jets are clustered with the anti-k T algorithm [30] with a distance parameter of 0.4 [25], as implemented in the FASTJET package [31]. Jet momentum is determined as the vectorial sum of all particle momenta in the jet. An offset is subtracted from the jet energies to take into account the contribution from pileup [32]. Jet energy corrections are obtained from simulation and are confirmed with in situ measurements of the energy balance in dijet and photon þ jet events [25]. Additional selection criteria are applied to each event to remove spurious jetlike features originating from isolated noise patterns in certain HCAL regions.
To identify jets originating from b quarks, we use an inclusive combined secondary vertex tagger (CSVv2) [33,34], which employs both secondary vertex and trackbased information. The working point is chosen to have about 70% b-tagging efficiency and a 1.5% light-flavor misidentification rate [35]. Double counting of objects is avoided by not considering jets that lie within a cone of radius 0.4 around a selected lepton.
While the main backgrounds are determined from data, as described in Sec. V, simulated events are used to validate the techniques and to estimate extrapolation factors as needed. In addition, some smaller backgrounds are estimated entirely from simulation. The leading-order (LO) MADGRAPH5 [36] event generator, using the NNPDF3.0LO [37] parton distribution functions (PDFs), is used to simulate tt þ jets, W þ jets, Z þ jets, and multijet events.
Single-top quark events in the t-channel and the tW process are generated using the next-to-leading order (NLO) POWHEGv1.0 [38][39][40][41][42] program, and in the s-channel process, as well as for ttW and ttZ production, using NLO MADGRAPH5_AMC@NLO [43]. All signal events are generated with MADGRAPH5, with up to two partons in addition to the gluino pair. Both programs use the NNPDF3.0NLO [37] PDF. The gluino decays are based on a pure phase-space matrix element [44], with signal production cross sections [45][46][47][48][49] computed at NLO plus next-to-leading-logarithm (NLL) accuracy.
Showering and hadronization of all partons is performed using the PYTHIA 8.2 ackage [44]. Pileup is generated for some nominal distribution of the number of proton-proton interactions per bunch crossing, which is weighted to match the corresponding distribution in data. The detector response for all backgrounds is modeled using the GEANT4 [50] package, while for the signal, the CMS fast simulation program [51] is used to reduce computation time. The fast simulation has been validated against the detailed GEANT4-based simulation for the variables relevant for this search, and efficiency corrections based on measurements in data are applied.

IV. TRIGGER AND EVENT SELECTION
The events are selected with an L1 trigger requiring H T > 150 GeV, followed by HLT requirements of H T > 350 GeV (online reconstruction) and at least one isolated lepton (an electron or muon) satisfying p T > 15 GeV. A trigger efficiency of 94 AE 1% is observed in the kinematic regime of the analysis, defined by lepton p T > 25 GeV and H T > 500 GeV, where the trigger efficiency reaches its maximum.
The electron or muon candidate is required to have a minimum p T of 25 GeV. Events with additional electrons or muons with p T > 10 GeV, satisfying the criteria for vetoed leptons, are rejected. Jets are selected with p T > 30 GeV and jηj < 2.4. In all search regions, we require at least five jets, where the two highest-p T jets must satisfy p T > 80 GeV.
To separate possible new-physics signals from background, we use the L T variable, which is defined as the scalar sum of the lepton p T and the missing transverse energy E miss T , and reflects the leptonic energy scale of the event. A minimum L T of 250 GeV is required, such that the analysis is not only sensitive to events with high E miss T , but also to signal events with very small E miss T , but higher lepton p T . An additional kinematic quantity important for the search is given by the hadronic energy scale of the event H T . A cutflow for the benchmark signal models is given in Table I.
After imposing the minimum requirements on L T and H T , several search regions are defined in bins of n jet , n b , L T , and H T , where n jet and n b are the numbers of jets and b-tagged jets, respectively. Defining search bins in b-jet multiplicity enables the analysis to target specific event topologies and to separate them from SM backgrounds. The phase space is divided into exclusive [0, 1, 2, ≥ 3] b-tagged jet categories when defining search bins, with a minimum b-jet p T of 30 GeV.
All search bins with at least one b-tagged jet, called "multi-b" bins in the following, are sensitive to the T1tttt model, while the search bins requiring zero b-tagged jets, called "zero-b" bins, are sensitive to the T5qqqqWW model. The baseline selection and the background estimation method differ for these two b-tag categories. For T1tttt, we expect a large number of jets and find in simulation that the n jet distribution peaks at eight jets for most mass points. We require at least six jets for the multi-b analysis and define two independent categories with 6-8 and ≥9 jets. For the zero-b analysis, where the investigated simplified T5qqqqWW model has fewer jets, we require, in the search region, 5, 6-7, or ≥8 jets. Depending on the specific SUSY particle masses, the hadronic event activity varies. To accommodate this, we define search bins in H T . Figure 2 shows the H T distributions for the multi-b and the zero-b selection. To exploit the strong separation power associated with the L T variable, we divide the search region into four bins in L T , such that sufficient statistical accuracy is given in each control bin to predict the background in the corresponding search bin.
After these selections, the main backgrounds are leptonically decaying W þ jets and semileptonic tt events. These backgrounds, both of which contain one lepton and one neutrino (from the W boson decay) in the final state, are mostly located at small ΔΦ values due to the correlation between the lepton and the neutrino. Therefore, the region with large ΔΦ is defined as the search region, while the events with small ΔΦ are used as the control sample. I. Expected event yields for SUSY signal benchmark models, normalized to 2.3 fb −1 . The baseline selection corresponds to all requirements up to and including the requirement on L T . The last two lines are exclusive for the zero-b and the multi-b selection, respectively. The events are corrected with scale factors to account for differences in the lepton identification and isolation efficiencies, trigger efficiency, and the b-tagging efficiency between simulation and data.  and several signal points are overlaid for illustration, but without stacking. Overflows are included in the last bin. The label DY refers to qq → Z=γ Ã → l þ l − events, and QCD refers to multijet events. The event yields for the benchmark models have been scaled up by a factor of 10. The ratio of data to simulation is given below each of the panels. All uncertainties are statistical only. Figure 3 shows the ΔΦ distributions for the zero-b and multi-b search regions. The ratio of the background event yield in the search region to that in the control region is determined in the corresponding signal-depleted sideband regions, which have smaller values of n jet , as discussed in Sec. V. Since the angle between the W boson and the lepton depends on the W momentum, being smaller for W bosons with higher boost, the ΔΦ requirement for the signal region is chosen depending on L T , which is a measure of the W boson p T . For the zero-b analysis, ΔΦ is required to be larger than 1.0 for most regions except for those with large L T , where the requirement is relaxed to 0.75, while the multi-b analysis has a relaxed ΔΦ requirement of 0.75 and 0.5 for medium-and high-L T regions, respectively.
In total, we define 30 search bins in the multi-b analysis and 13 search bins in the zero-b analysis, as described in detail in Table II.

V. BACKGROUND ESTIMATION
The dominant backgrounds in this search are from tt þ jets and W þ jets events, whose contributions vary with the multiplicity of b-tagged jets and the kinematic region in H T and L T . To determine these backgrounds, we define two regions for each bin in L T , H T , and n b : the search region (SR) with large values of ΔΦ, and the control region (CR) with low values of ΔΦ, with the separation requirement depending on the L T value, as shown in Table II. We further divide each of these bins into low-n jet sideband (SB) and high-n jet main band (MB) regions.
About 10%-15% of the SM background events in the CR are expected to be multijet events (denoted in the following as QCD) and are predicted as described in Sec. V C. Since the multijet background is negligible in the SR, it is subtracted from the number of background events in the CR when calculating the transfer factor R data CS to extrapolate from CR (low-ΔΦ) to SR (high-ΔΦ). This transfer factor R data CS is determined from data in the low-n jet SB regions, separately for each L T , H T , and n b search region: where N SB data ðSRÞ is the number of events in the low-n jet SB high-ΔΦ signal region, N SB data ðCRÞ the number of events in the low-n jet SB low-ΔΦ control region, and N SB QCD pred ðCRÞ the predicted number of QCD multijet events in the SB CR.
In the regions with one b tag and four or five jets, about 80% tt þ jets events and 15%-20% W þ jets and single top quark events are expected, while in all other multi-b regions, tt background is completely dominant. Because only a single SM background dominates in the multi-b analysis, just one R CS factor is needed for each L T , H T , and n b range. In the zero-b bins, the contributions from W þ jets and tt þ jets are roughly equal. Here, an extension of the multi-b strategy is employed, which takes into account differences in the R CS values for these two backgrounds.
An overview of the ðn jet ; n b Þ regions used in this analysis, as discussed in detail in Secs. VA-V C, is given in Table III  and several signal points are overlaid for illustration, but without stacking. The wider bins are normalized to a bin width of 0.1. The label DY refers to qq → Z=γ Ã → l þ l − events, and QCD refers to multijet events. The event yields for the benchmark models have been scaled up by a factor of 10. The ratio of data to simulation is given below each of the panels.
A. Estimate of the leading backgrounds for n b ≥ 1 For the multi-b analysis, the SB region, where R CS is determined, is required to have four or five jets, while the MB region must satisfy n jet ∈ ½6 − 8 or n jet ≥ 9. To account for possible differences in this extrapolation from SB to MB as a function of jet multiplicity, we apply multiplicative correction factors κ EW , determined from simulation. The predicted number N MB pred ðSRÞ of background events in each MB SR is then given by with Here R data CS is determined from Eq. (1), N MB data ðCRÞ is the number of data events in the CR of the MB region, and N MB QCD pred ðCRÞ is the predicted number of multijet events in the MB. The label EW refers to all backgrounds other than multijets. The residual difference of the values of R CS between the SB and MB regions is evaluated in simulation as the correction factor κ EW given by Eq. (3), where R MC CS ðMB; EWÞ is the R CS in a search MB region from simulation and R MC CS ðSB; EWÞ is the R CS in the corresponding SB region in simulation for the EW background.
The κ EW factor is determined separately for each search bin, except that an overall κ EW factor is applied for the n b ≥ 2 search bins with the same H T and L T , since the κ EW factors are found to be nearly independent of n b . Similarly, R CS at very high H T is determined jointly across all three n b bins to increase the number of events, as the overall uncertainty of the background prediction for several of the search bins is dominated by the statistical uncertainty of the yield in the SR of the side band.
The value of R CS for the total background is equal to the sum of the R CS values of each background component, weighted with the relative contributions of the components. For semileptonic tt and W þ jets events, which contain both one neutrino from the hard interaction, R CS typically has values of 0.01 to 0.04, depending on the search bin. In events with more than one neutrino, e.g. in tt events in which both W bosons decay leptonically, R CS is higher, with values of around 0.5. This is visible in Fig. 3, where at high ΔΦ a large fraction of events is due to dileptonic tt þ jets background, while the low-ΔΦ region is dominated by events with only one neutrino. A larger R CS is also expected for events with three neutrinos, such as ttZ, when the tt system decays semileptonically and the Z boson decays to

Analysis
Multi-b analysis Zero-b analysis MB MB two neutrinos. The influence of these latter processes is small since their relative contribution to the background is minor. Most of the SRs with six or more jets are dominated by semileptonic tt events, and therefore this background dominates the total R CS value of ≈0.05. As the R CS for dileptonic tt events is an order of magnitude larger than for semileptonic tt events, a slight change in composition in the CR from low-to high-n jet multiplicity translates into κ EW slightly different from unity. This change in the dileptonic tt contribution is accounted for by assigning an uncertainty on the n jet extrapolation based on a dileptonic control sample in data, as discussed in Sec. VI.
B. Estimate of the leading backgrounds for n b = 0 For search bins in which b-tagged jets are vetoed, the background contributions from W þ jets and tt þ jets events are estimated by applying the R CS method separately to each of the two components. This strategy implies the use of two sidebands enriched in W þ jets and tt þ jets events, respectively. We write the total background in each search region n SR jet (with a ΔΦ requirement as shown in Table II) as where the predicted yields of W þ jets and tt þ jets background events are denoted by N SR W and N SR tt , respectively. Additional backgrounds from rare sources are estimated from simulation and denoted by N SRðMCÞ other . The expected number of events for each of the background components can be described by where N CR data is the total number of events in the CR of the MB region and f i is the relative yield of component i. The relative contributions of the two components are determined by a fit of templates obtained from simulation to the n b multiplicity distribution in the CR of the MB region. The contribution of the QCD multijet background in the CR is fixed to the yield estimated from data as described in Sec. V C. The contribution of other rare background components is obtained from simulation as well, as is done in the SR. Uncertainties in these two components are propagated as systematic uncertainties to the final prediction. Examples of these fits are shown in Fig. 4.
The two R CS values, for W þ jets and tt þ jets, are measured in two different low-n jet SB regions. For the tt þ jets estimate, a sideband with the requirements 4 ≤ n jet ≤ 5 and n b ¼ 1 is used. The value of R tt CS is then given by The correction factors κ b and κ tt are determined from simulation. The factor κ b corrects for a potential difference of R tt CS between samples with zero or one b jet and for the small contributions of backgrounds other than tt þ jets or QCD multijet events. The factor κ tt corrects for a residual dependence of R tt CS on n jet , in analogy to the κ EW factor defined in Sec. VA. Both values, κ b and κ tt , are close to unity, and statistical uncertainties from the simulation are propagated to the predicted yields.
Similarly, the W þ jets contribution is estimated using R CS values from a sideband with 3 ≤ n jet ≤ 4 and n b ¼ 0.
With respect to the SB used for the estimate of R tt CS , a lower jet multiplicity is chosen in order to limit the contamination from tt þ jets events. Only the muon channel is used since it has a negligible contamination from QCD multijet events, contrary to the electron channel. A systematic uncertainty is derived from simulation to cover potential differences between the μ and the combined e and μ samples. The 3 ≤ n jet ≤ 4 (250 ≤ L T < 350 GeV, H T ≥ 500 GeV, ΔΦ < 1) and (right) 6 ≤ n jet ≤ 7 (250 ≤ L T < 350 GeV, H T ≥ 750 GeV, ΔΦ < 1) in data (muon channel). The solid lines represent the templates scaled according to the fit result (blue for tt þ jets, green for W þ jets, turquoise for QCD, and red for the remaining backgrounds), the dashed line shows the sum after the fit, and the points with error bars represent data. value of R W CS is given by Again, the factor κ W corrects for a residual dependence of R W CS on the jet multiplicity. The raw value of R data CS measured in the SB has to be corrected for the contamination of tt þ jets events. The tt þ jets yields are subtracted in the numerator and denominator according to The event yields N CR data and N SR data are measured in the SB CRs and SRs. The fraction of tt þ jets events f tt is again obtained by a fit to the n b multiplicity in the SB CR. The R CS value for tt þ jets in this SB is obtained from simulation.
Systematic uncertainties are assigned to κ tt and κ W according to the difference between the R CS values in the sideband and the result of a linear fit over the full range of n jet . The uncertainties vary from 3% to 43% for κ tt and from 1% to 49% for κ W . The two sources are treated as being independent.

C. Estimate of the multijet background
Multijet events enter this analysis mostly when reconstructed electrons originate from misidentified jets or from photon conversion in the inner detector. This background is estimated from the yield of "antiselected" electron candidates in each region, which pass looser identification and isolation requirements, and fail the tighter criteria for selected electrons. These events are scaled by the ratio of jets and photons that pass the tight electron identification requirements to the number of antiselected electron candidates in a multijet-enriched control sample with no b-tagged jets and three or four other jets. The assumption is that this sample is devoid of genuine prompt electrons. The estimation method was introduced previously [10,52], and it relies on the L P variable: For the dominant SM backgrounds, tt þ jets and W þ jets, the distribution of L P is a well-understood consequence of the W boson polarization and falls from 0 to 1. In contrast, the distribution of L P for multijet events peaks near L P ¼ 1.
The ratio of selected to antiselected electron candidates is obtained from a fit to the L P distribution in bins of L T . The shape of the QCD multijet contribution used in the fit is taken from the antiselected sample, while the shape of all other contributions is taken from simulation, as the behavior due to W polarization is well understood. The ratios are found to be in the range 0.1-0.2.
In principle, the background estimation with the R CS method requires knowledge of the multijet contribution in the SR and CR separately. Since the multijet background estimation is performed inclusively with respect to ΔΦ, an R CS factor for multijet events is determined as well. In practice, since the resulting R CS values are all found to be below 2%, the multijet contamination is negligible for the SR. Therefore, the previously described R CS method takes into account only the QCD multijet contribution in the CR, as written in Eq. (1). For the muon channel, the contribution from QCD multijet background is typically of the order of 1% of the total background. To estimate this contribution, a procedure similar to the one outlined above is applied and assigned a 100% uncertainty.

VI. SYSTEMATIC UNCERTAINTIES
Systematic uncertainties either influence κ, and thereby the predictions for the background, or modify the expected signal yield.
The main systematic uncertainty on the background arises from the extrapolation of R CS from the low n jet region, where it is measured, to the MB regions of higher jet multiplicities, where it is applied. Therefore, a systematic uncertainty on R CS is determined in a dedicated control region with dileptonic events. The ratio of the semileptonic to dileptonic tt þ jets final states for different numbers of reconstructed jets is of major importance since the total R CS is based on the fraction of the two channels and their corresponding R CS values, which differ significantly in tt þ jets events. To ensure that the data are described well by simulation, a high-purity dilepton tt þ jets control sample is selected from the data by requiring two leptons of opposite charge. For same-flavor leptons, it is also required that the invariant mass of the lepton pair be more than 10 GeV away from the Z boson mass peak. To study the behavior of the dileptonic events in the single-lepton selection, one of the two leptons is removed from the event. Since these "lost leptons" are principally from τ → hadrons þ ν decays, we replace the removed lepton by a jet with 2=3 of the original lepton's p T to accommodate for the missing energy due to the neutrino from the τ decay, and we recalculate the L T , ΔΦ, and H T values of the now "singlelepton" event. In order to maximize the number of events, no ΔΦ requirement is applied, and all events are used twice, with each reconstructed lepton being considered as the lost lepton. We refer to the samples produced using this procedure as the dilepton CRs.
A key test is performed by comparing the jet multiplicity distribution in the sample resulting from single-lepton baseline selection (excluding the SRs) with the corresponding simulated event sample, and by comparing the dilepton CRs with the corresponding simulated event sample. Both comparisons show the same trend, a slight overprediction by simulation of the rate of high jet multiplicity events. The ratio of event yields in data-to-simulation is computed for each comparison, and the two ratios are then divided to see whether the behavior in data relative to simulation is the same in both pairs of samples. This double ratio is consistent with unity within statistical uncertainty. The systematic uncertainty in the description of the n jet distribution in simulation is determined from this double ratio and is mainly due to the statistical uncertainty of the data samples, which is within 8%-40%, and therefore larger than the observed slope of the double ratio vs n jet .
The remaining uncertainties are smaller than the one from the dileptonic tt þ jets fraction. In particular, the applied jet energy scale (JES) factors are varied up and down according to their uncertainty [25] as a function of jet p T and η, and these changes are propagated to E miss T . The scale factors applied to the efficiencies for the identification of b-quark jets and for the misidentification of c-quark, light-quark, or gluon jets are also varied up and down according to their uncertainties [34]. Uncertainties for the efficiency of lepton reconstruction and identification are handled in the same way. For pileup, a 5% uncertainty in the inelastic cross section [53] is used to obtain its impact on the uncertainty in the pileup. In a few bins with a low number of simulated events, the reweighting leads to a large uncertainty. All these uncertainties apply to both the background prediction and the signal yield. The luminosity is measured with the pixel cluster counting method, and the absolute luminosity scale calibration is derived from an analysis of van der Meer scans performed in August 2015, resulting in an uncertainty of 2.7% [54].
The W þ jets and tt þ jets cross sections are changed by 30% [55] to cover possible biases in the estimation of the background composition in terms of W þ jets vs tt þ jets events, which would lead to a slight change in the κ value. These changes have only a small impact on the zero-b analysis, where the relative fraction of the two processes is determined from a fit. Also, the following changes in the simulation are performed, with differences between the values of κ in the reweighted and original samples defining the uncertainties. Motivated by measurements at ffiffi ffi s p ¼ 8 TeV, simulated tt þ jets events are reweighted by a factor ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi T ÞÞ, to improve the modeling of the top quark p T spectrum [56]. The reweighting preserves the normalization of the sample, and the difference relative to the results obtained with the unweighted sample is assigned as a systematic uncertainty. The polarization of W bosons is varied by reweighting events by the factor wðcos θ Ã Þ ¼ 1 þ αð1 − cos θ Ã Þ 2 , where θ Ã is the angle between the charged lepton and W boson in the W boson rest frame. In W þ jets events, we take α to be 0.1, guided by the theoretical uncertainty and measurements found in Refs. [52,[57][58][59]. For tt þ jets events, we take α ¼ 0.05. For W þ jets events, where the initial state can have different polarizations for W þ vs W − bosons, we take, as the uncertainty, the larger change in κ resulting from reweighting only the W þ bosons in the sample, and from reweighting all W bosons. The ttV cross section is varied by 100%. The systematic uncertainty in the multijet estimation depends on n jet and n b , and ranges from 25% to 100%.
For the zero-b analysis, an additional systematic uncertainty is applied, based on linear fits of R CS as a function of n jet as described in Sec. V B, and a 50% cross-section uncertainty is used for all backgrounds other than W þ jets, tt þ jets, ttV, and multijets.
For the signal, an uncertainty in initial-state radiation (ISR) is applied, based on the p T of the gluino-gluino system, which corresponds to a 15% uncertainty at p T between 400 and 600 GeV, and 30% at larger p T . This uncertainty is based on measurements of ISR in Z þ jets and tt þ jets events [16,60]. The factorization and renormalization scales are each changed by a factor of 0.5 and 2. Uncertainties in the signal cross section are also taken into account. The impact of the systematic uncertainties in the total background prediction for the multi-b and zero-b analyses is summarized in Table IV. While the systematic uncertainty is determined for each signal point, the uncertainties typical for most signals are summarized, for illustration, in Table V.

VII. RESULTS AND INTERPRETATION
The backgrounds for all SRs are determined, as described previously, in different SB regions with lower jet or b-jet multiplicities. The results of the background prediction and the observed data are shown in Table VI and Fig. 5 for the multi-b events. In this figure, the outline of the filled histogram represents the total number of background events from the prediction. For illustration, the relative amount of tt þ jets, W þ jets, and of other backgrounds is shown as well, based on the fractions estimated in simulation. Table VII and Fig. 6 show the results for the zero-b events. Here, the filled histogram represents the predictions from data for tt þ jets and W þ jets events, and for the remaining backgrounds, where the latter include the multijet prediction determined from data and rare backgrounds taken from simulation. The data agree with SM expectations, and no excess is observed.
To set limits, separate likelihood functions, one for the multi-b analysis and one for the zero-b analysis, are constructed from the Poisson probability functions for all four data regions (the CRs and SRs in the SB as well in the MB) to determine the background in the MB SR. In addition, the κ values from simulation are included to correct any residual differences between the SB and MB regions, with uncertainties incorporated through lognormal constraints. The estimated contribution from multijet events in the two CRs is also included. A possible signal contamination is taken into account by including signal terms in the fit for both the sideband and the control regions. For the zero-b analysis, the relative contributions of W þ jets and tt þ jets events as determined in the fits to the n b distribution in the CRs are treated as external measurements. The correlation between the W þ jets and tt þ jets yields introduced by these fits is taken into account. A profile likelihood ratio in the asymptotic approximation [61] is used as the test statistic. Limits are then calculated at the 95% confidence level (CL) using the asymptotic CL s criterion [62,63].
The cross-section limits obtained for the T1tttt model using the multi-b analysis, and for the T5qqqqWW model using the zero-b analysis, are shown in Fig. 7 as a function of mðgÞ and mðχ 0 1 Þ, assuming branching fractions of 100% as shown in Fig. 1. Using thegg pair production cross section calculated at next-to-leading order within next-to-leading-logarithmic accuracy, exclusion limits are set as a function of the ðmg; m~χ0 1 Þ mass hypothesis.  The data are shown by black points with error bars, while the total SM background predictions are shown by a grey line, with a hatched region representing its uncertainty. For illustration, the relative fraction of the different SM background contributions, as determined from simulation, is shown by the stacked, colored histograms, whose total normalization is set by the total background yields obtained from the control samples in the data. The expected event yields for two T1tttt SUSY benchmark models are shown by open histograms, each of which is shown stacked on the total background prediction. The vertical dashed and dotted lines separate different n jet and L T bins, respectively, as indicated by the x-axis labels. Lower panel: The ratio of the yield observed in data to the predicted background yield is shown for each bin. The error bars on the data points indicate the combined statistical and systematic uncertainty in the ratio. The grey hatched region indicates the uncertainty on the ratio that arises from the uncertainty on the background prediction.

VIII. SUMMARY
A search for supersymmetry has been performed with 2.3 fb −1 of proton-proton collision data recorded by the CMS experiment at ffiffi ffi s p ¼ 13 TeV in 2015. The data are analyzed in several exclusive categories, differing in the number of jets and b-tagged jets, the scalar sum of all jet transverse momenta, and the scalar sum of the missing transverse momentum and the transverse momentum of the lepton. The main background is significantly reduced by requiring a large azimuthal angle between the directions of the momenta of the lepton and of the reconstructed W boson. No significant excess is observed, and the results are interpreted in terms of two simplified models that describe gluino pair production. For the simplified model T1tttt, in which each gluino decays through an off-shell top squark to a tt pair and the lightest neutralino, gluino masses up to 1.6 TeV are excluded for neutralino masses below 600 GeV. Neutralino masses below 850 GeV can be excluded for a gluino mass up to 1.4 TeV. Similar to Ref. [16], these results extend the limits obtained from the 8 TeV searches [13][14][15] by about 250 GeV. The second simplified model T5qqqqWW also contains gluino pair production, with the gluinos decaying to first or second generation squarks and a chargino, which then

95% CL upper limit on cross section [pb]
FIG. 7. Cross-section limits at a 95% C.L. for the (left) T1tttt and (right) T5qqqqWW models, as a function of the gluino and LSP masses. In T5qqqqWW, the pair-produced gluinos each decay to a quark-antiquark pair of the first or second generation (qq), and a chargino (χ AE 1 ) with its mass taken to be m~χAE 1 ¼ 0.5ðm~g þ m~χ0 1 Þ. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the AE1σ uncertainty bands related to the theoretical (experimental) uncertainties. decays to a W boson and the lightest neutralino. The chargino mass in this decay chain is taken to be m~χAE 1 ¼ 0.5ðmg þ m~χ0 1 Þ. In this model, gluino masses below 1.4 TeV are excluded for neutralino masses below 700 GeV. For a gluino mass of 1.3 TeV, neutralinos with masses up to 850 GeV can be excluded. These results improve existing limits [17] on the neutralino mass in this channel for gluino masses between 900 GeV and 1.4 TeV.

ACKNOWLEDGMENTS
We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centers and personnel of the Worldwide LHC Computing Grid for delivering, so effectively, the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: BMWFW and FWF  [56] CMS Collaboration, Measurement of the differential cross section for top quark pair production in pp collisions at ffiffi ffi s p ¼ 8 TeV, Eur. Phys. J. C 75, 542 (2015).