"Search for supersymmetry with razor variables in pp collisions at sqrt(s) = 7 TeV"

The razor approach to search for R-parity conserving supersymmetric particles is described in detail. Two analyses are considered: an inclusive search for new heavy particle pairs decaying to final states with at least two jets and missing transverse energy, and a dedicated search for final states with at least one jet originating from a bottom quark. For both the inclusive study and the study requiring a bottom-quark jet, the data are examined in exclusive final states corresponding to all-hadronic, single-lepton, and dilepton events. The study is based on the data set of proton-proton collisions at sqrt(s) = 7 TeV collected with the CMS detector at the LHC in 2011, corresponding to an integrated luminosity of 4.7 inverse femtobarns. The study consists of a shape analysis performed in the plane of two kinematic variables, denoted MR and R2, that correspond to the mass and transverse energy flow, respectively, of pair-produced, heavy, new-physics particles. The data are found to be compatible with the background model, defined by studying event simulations and data control samples. Exclusion limits for squark and gluino production are derived in the context of the constrained supersymmetric model (CMSSM) for simplified-model spectra (SMS). Within the CMSSM parameter considered, squark and gluino masses up to 1350 GeV are at 95% confidence level, depending on the model parameters. For SMS scenarios, the direct production of pairs of top or bottom squarks is excluded for masses as high as 400 GeV. The razor approach to search for R – parity conserving supersymmetric particles is described in detail. Two analyses are considered: an inclusive search for new heavy particle pairs decaying to final states with at least two jets and missing transverse energy, and a dedicated search for final states with at least one jet originating from a bottom quark. For both the inclusive study and the study requiring a bottom-quark jet, the data are examined in exclusive final states corresponding to all-hadronic, single-lepton, and dilepton events. The study is based on the data set of proton-proton collisions at ﬃﬃﬃ s p ¼ 7 TeV collected with the CMS detector at the LHC in 2011, corresponding to an integrated luminosity of 4 . 7 fb − 1 . The study consists of a shape analysis performed in the plane of two kinematic variables, denoted M R and R 2 , that correspond to the mass and transverse energy flow, respectively, of pair-produced, heavy, new-physics particles. The data are found to be compatible with the background model, defined by studying event simulations and data control samples. Exclusion limits for squark and gluino production are derived in the context of the constrained minimal supersymmetric standard model (CMSSM) and also for simplified-model spectra (SMS). Within the CMSSM parameter space considered, squark and gluino masses up to 1350 GeVare excluded at 95% confidence level, depending on the model parameters. For SMS scenarios, the direct production of pairs of top or bottom squarks is excluded for masses as high as 400 GeV.

In this paper, we present the detailed methodology of an inclusive search for SUSY based on the razor kinematic variables [26,27]. A summary of the results of this search, based on 4.7 fb −1 of pp collision data at ffiffi ffi s p ¼ 7 TeV collected with the CMS detector at the LHC, can be found in Ref. [28]. The search is sensitive to the production of pairs of heavy particles, provided that the decays of these particles produce significant E miss T . The jets in each event are cast into two disjoint sets, referred to as "megajets." The razor variables M R and R 2 , defined in Sec. II, are calculated from the four-momenta of these megajets event by event, and the search is performed by determining the expected distributions of SM processes in the twodimensional (M R , R 2 ) razor plane. A critical feature of the razor variables is that they are computed in the approximate center-of-mass frame of the produced superpartner candidates.
The megajets represent the visible part of the decay chain of pair-produced superpartners, each of which decays to one or more visible SM particles and one stable, weakly interacting lightest SUSY particle (LSP), here taken to be the lightest neutralinoχ 0 1 . In this framework the reconstructed products of the decay chain of each originally produced superpartner are collected into one megajet. Every topology can then be described kinematically by the simplest example of squark-antisquark production with the direct two-body squark decayq → qχ 0 1 , denoted a "dijet plus E miss T " final state, to which the razor variables strictly apply.
The strategy and execution of the search is summarized as follows: (1) Events with two reconstructed jets at the hardwarebased first level trigger (L1) are processed by a dedicated set of algorithms in the high-level trigger (HLT). From the jets and leptons reconstructed at the HLT level, the razor variables M R and R 2 are calculated and their values are used to determine whether to retain the event for further off-line processing. A looser kinematic requirement is applied for events with electrons or muons, due to the smaller rate of SM background for these processes. The correspondence between the HLT and off-line reconstruction procedures allows events of interest to be selected more efficiently than is possible with an inclusive multipurpose trigger. (2) In the off-line environment, leptons and jets are reconstructed, and a tagging algorithm is applied to identify those jets likely to have originated from a bottom-quark jet (b jet). (3) The reconstructed objects in each event are combined into two megajets, which are used to calculate the variables M R and R 2 . Several baseline kinematic requirements are applied to reduce the number of misreconstructed events and to ensure that only regions of the razor plane where the trigger is efficient are selected. (4) Events are assigned to final-state "boxes" based on the presence or absence of a reconstructed lepton. This box partitioning scheme allows us to isolate individual SM background processes based on the final-state particle content and kinematic phase space; we are able to measure the yield and the distribution of events in the (M R , R 2 ) razor plane for different SM backgrounds. Events with at least one tagged b jet are considered in a parallel analysis focusing on a search for the superpartners of third-generation quarks. In total, we consider 12 mutually exclusive final-state boxes: dielectron events (ELE-ELE), electron-muon events (ELE-MU), dimuon events (MU-MU), single-electron events (ELE), single-muon events (MU), and events with no identified electron or muon (HAD), each inclusive or with a b-tagged jet. (5) For each box we use the low (M R , R 2 ) region of the razor plane, where negligible signal contributions are expected, to determine the shape and normalization of the various background components. An analytic model constructed from these results is used to predict the SM background over the entire razor plane. (6) The data are compared with the prediction for the background in the sensitive regions of the razor plane and the results are used to constrain the parameter space of SUSY models. This paper is structured as follows. The definition of the razor variables is given in Sec. II. The trigger and off-line event selection are discussed in Sec. III. The features of the signal and background kinematic distributions are described in Sec. IV. In Sec. V we describe the sources of SM background, and in Sec. VI the analytic model used to characterize this background in the signal regions. Systematic uncertainties are discussed in Sec. VII. The interpretation of the results is presented in Sec. VIII in terms of exclusion limits on squark and gluino production in the context both of the constrained minimal SUSY model (CMSSM) [29][30][31] and for some simplified model spectra (SMS) [32][33][34][35][36]. Section IX contains a summary. For the CMSSM, exclusion limits are provided as a function of the universal scalar and fermion mass values at the unification scale, respectively denoted m 0 and m 1=2 . For the SMS, limits are provided in terms of the masses of the produced SUSY partner and the LSP.

II. THE RAZOR APPROACH
The razor kinematic variables are designed to be sensitive to processes involving the pair-production of two heavy particles, each decaying to an unseen particle plus jets. Such processes include SUSY particle production with various decay chains, the simplest example of which is the pair production of squarks, where each squark decays to a quark and the LSP, with the LSP assumed to be stable and weakly interacting. In processes with two or more undetected energetic final-state particles, it is not possible to fully reconstruct the event kinematics. Event by event, one cannot make precise assignments of the reconstructed finalstate particles (leptons, jets, and undetected neutrinos and LSPs) to each of the original superpartners produced. For a given event, there is not enough information to determine the mass of the parent particles, the subprocess center-ofmass energy ffiffi f s p , the center-of-mass frame of the colliding protons, or the rest frame of the decay of either parent particle. As a result, it is challenging to distinguish between SUSY signal events and SM background events with energetic neutrinos, even though the latter involve different topologies and mass scales. It is also challenging to identify events with instrumental sources of E miss T that can mimic the signal topology.
The razor approach [26,27] addresses these challenges through a novel treatment of the event kinematics. The key points of this approach are listed below.
(i) The visible particles (leptons and jets) are used to define two megajets, each representing the visible part of a parent particle decay. The megajet reconstruction ignores details of the decay chains in favor of obtaining the best correspondence between a signal event candidate and the presumption of a pair-produced heavy particle that undergoes two-body decay. (ii) Lorentz-boosted reference frames are defined in terms of the megajets. These frames approximate, event by event, the center-of-mass frame of the signal subprocess and the rest frames of the decays of the parent particles. The kinematic quantities in these frames can be used to extract the relevant SUSY mass scales. (iii) The razor variables, M R , M R T , and R ≡ M R T =M R , are computed from the megajet four-momenta and the E miss T in the event. The M R variable is an estimate of an overall mass scale, which in the limit of massless decay products equals the mass of the heavy parent particle. It contains both longitudinal and transverse information, and its distribution peaks at the true value of the new-physics mass scale. The razor variable M R T is defined entirely from transverse information: the transverse momenta (p T ) of the megajets and the E miss T . This variable has a kinematic end point at the same underlying mass scale as the M R mean value. The ratio R quantifies the flow of energy in the plane perpendicular to the beam and the partitioning of momentum between visible and invisible particles. (iv) The shapes of the distributions in the (M R , R 2 ) plane are described for the SM processes. Razor variable distributions exhibit peaks for most SM backgrounds, as a result of turn-on effects from trigger and selection thresholds as well as of the relevant heavy mass scales for SM processes, namely the topquark mass and the W and Z boson masses. However, compared with signals involving heavier particles and new-physics sources of E miss T , the SM distributions peak at smaller values of the razor variables. For values of the razor variables above the peaks, the SM background distributions (and also the signal distributions) exhibit exponentially falling behavior in the (M R , R 2 ) plane. Hence, the asymptotic behavior of the razor variables is determined by a combination of the parton luminosities and the intrinsic sources of E miss T . The multijet background from processes described by quantum chromodynamics (QCD), which contains the smallest level of intrinsic E miss T amongst the major sources of SM background, has the steepest exponential falloff. Backgrounds with energetic neutrinos from W=Z boson and top-quark production exhibit a slower falloff and resemble each other closely in the asymptotic regime. Thus, razor signals are characterized by peaks in the (M R , R 2 ) plane on top of exponentially falling SM background distributions. Any SUSY search based on razor variables is then more similar to a "bump hunt," e.g., a search for heavy resonances decaying to two jets [37], than to a traditional SUSY search. This justifies the use of a shape analysis, based on an analytic fit of the background, as described in Sec. VI.

A. Razor megajet reconstruction
The razor megajets are defined by dividing the reconstructed jets of each event into two partitions. Each partition contains at least one jet. The megajet four-momenta are defined as the sum of the four-momenta of the assigned jets. Of all the possible combinations, the one that minimizes the sum of the squared-invariant-mass values of the two megajets is selected. In simulated event samples, this megajet algorithm is found to be stable against variations in the jet definition and it provides an unbiased description of the visible part of the two decay chains in SUSY signal events. The inclusive nature of the megajets allows an estimate of the SM background in the razor plane.
Reconstructed leptons in the final state can be included as visible objects in the reconstruction of the megajets, or they can be treated as invisible, i.e., as though they are escaping weakly interacting particles [26]. For SM background processes such as WðlνÞ þ jets, the former choice yields more transversely balanced megajets and lower values of R. If the leptons are treated as invisible in these processes, the E miss T corresponds to the entire W boson p T value, similar to the case of ZðννÞ þ jets events.

B. Razor variables
To the extent that the reconstructed pair of megajets accurately reflects the visible portion of the underlying parent particle decays, the kinematics of the event are equivalent to that of the pair production of heavy squarks q 1 ,q 2 , withq i → q iχ 0 1 , whereχ 0 1 denotes the LSP and q i denotes the visible products of the decays as represented by the megajets.
The razor analysis approximates the unknown center-ofmass and parent particle rest frames with a razor frame defined unambiguously from measured quantities in the laboratory frame. Two observables M R and M R T estimate the heavy mass scale M Δ . Consider the two visible fourmomenta written in the rest frame of the respective parent particles: whereû q i (i ¼ 1; 2) is a unit three-vector and m q i represents the mass corresponding to the megajet, e.g., the top-quark mass fort → tχ 0 1 . Here we have parametrized the magnitude of the three-momenta by the mass scale M Δ , where In the limit of massless megajets we then have Þ=mq and the four-momenta reduce to The razor variable M R is defined in terms of the momenta of the two megajets by wherep q i is the momentum of megajet q i (i ¼ 1; 2) and p q i z is its component along the beam direction. For massless megajets, M R is invariant under a longitudinal boost. It is always possible to perform a longitudinal boost to a razor frame where p q 1 z þ p q 2 z vanishes, and M R becomes just the scalar sum of the megajet threemomenta added in quadrature. For heavy particle production near threshold, the three-momenta in this razor frame do not differ significantly from the three-momenta in the actual parent particle rest frames. Thus, for SUSY signal events, M R is an estimator of M Δ , and for simulated samples we find that the distribution of M R indeed peaks around the true value of M Δ . This definition of M R is improved with respect to the one used in Ref. [26], to avoid configurations where the razor frame is unphysical.
The razor observable M R T is defined as wherep T q i is the transverse momentum of megajet q i (i ¼ 1; 2) and p q i T is the corresponding magnitude; similarly,Ẽ miss T is the missing transverse momentum in the event and E miss T its magnitude. Given a global estimator M R and a transverse estimator M R T , the razor dimensionless ratio is defined as For signal events, M R T has a maximum value (a kinematic end point) at M Δ , so R has a maximum value of approximately 1. Thus, together with the shape of M R peaking at M Δ , this behavior is in stark contrast with, for example, QCD multijet background events, whose distributions in both M R and R 2 fall exponentially. These properties allow us to identify a region of the two-dimensional (2D) razor space where the contributions of the SM background are reduced while those of signal events are enhanced.

C. SUSY and SM in the razor plane
The expected distributions of the main SM backgrounds in the razor plane, based on simulation, are shown in Fig. 1,  along with the results from the CMSSM low-mass benchmark model LM6 [38], for which M Δ ¼ 831 GeV. The peaking behavior of the signal events at M R ≈ M Δ , and the exponential falloff of the SM distributions with increasing M R and R 2 , are to be noted. For both signal and background processes, events with small values of M R are suppressed because of a requirement that there be at least two jets above a certain threshold in p T (Sec. III E).
In the context of SMS, we refer to the pair production of squark pairsq,q Ã , followed byq → qχ 0 1 , as "T2" scenarios [39], where theq Ã state is the charge conjugate of theq state. Figure 2(a) shows a diagram for heavy-squark pair production. The distributions of M R and R 2 for different LSP masses are shown in Figs. 2(b) and 2(c). Figure 2(d) shows the distribution of signal events in the razor plane. The colored bands running from top left to bottom right show the approximate SM background constant-yield contours. The associated numbers indicate the SM yield suppression relative to the reference line marked "1." Based on these kinematic properties, a 2D analytical description of the SM processes in the (M R , R 2 ) plane is developed.

A. The CMS apparatus
A hallmark of the CMS detector [40] is its superconducting solenoid magnet, of 6 m internal diameter, providing a field of 3.8 T. The silicon pixel and strip tracker, the crystal electromagnetic calorimeter (ECAL), and the brass/scintillator hadron calorimeter (HCAL) are contained within the solenoid. Muons are detected in gasionization detectors embedded in the steel flux-return yoke, based on three different technologies: drift tubes, resistive plate chambers, and cathode strip chambers (CSCs). The ECAL has an energy resolution better than 0.5% above 100 GeV. The combination of the HCAL and ECAL provides jet energy measurements with a resolution ΔE=E ≈ 100%= ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi E=GeV p ⊕ 5%. The CMS experiment uses a coordinate system with the origin located at the nominal collision point, the x axis pointing towards the center of the LHC ring, the y axis pointing up (perpendicular to the plane containing the LHC ring), and the z axis along the counterclockwise beam direction. The azimuthal angle, ϕ, is measured with respect to the x axis in the ðx; yÞ plane, and the polar angle, θ, is defined with respect to the z axis. The pseudorapidity is η ¼ − ln½tanðθ=2Þ.
For the data used in this analysis, the peak luminosity of the LHC increased from 1 × 10 33 cm −2 s −1 to over 4 × 10 33 cm −2 s −1 . For the data collected between ð1-2Þ × 10 33 cm −2 s −1 , the increase was achieved by increasing the number of bunches colliding in the machine, keeping the average number of interactions per crossing at about seven. For the rest of the data, the increase in the instantaneous luminosity was achieved by increasing the number and density of the protons in each bunch, leading to an increase in the average number of interactions per crossing from around 7 to around 17. The presence of multiple interactions per crossing was taken into account in the CMS Monte Carlo (MC) simulation by adding a random number of minimum bias events to the hard interactions, with the multiplicity distribution matching that in data.

B. Trigger selection
The CMS experiment uses a two-stage trigger system, with events flowing from the L1 trigger at a rate up to 100 kHz. These events are then processed by the HLT computer farm. The HLT software selects events for storage and off-line analysis at a rate of a few hundred Hz. The HLT algorithms consist of sequences of off-linestyle reconstruction and filtering modules.
The 2010 CMS razor-based inclusive search for SUSY [26] used triggers based on the scalar sum of jet p T , H T , for hadronic final states and single-lepton triggers for leptonic final states. Because of the higher peak luminosity of the LHC in 2011, the corresponding triggers for 2011 had higher thresholds. To preserve the high sensitivity of the razor analysis, CMS designed a suite of dedicated razor triggers, implemented in the spring of 2011. The total integrated luminosity collected with these triggers was 4.7 fb −1 at ffiffi ffi s p ¼ 7 TeV.
The razor triggers apply thresholds to the values of M R and R driven by the allocated bandwidth. The algorithms used for the calculation of M R and R are based on calorimetric objects. The reconstruction of these objects is fast enough to satisfy the stringent timing constraints imposed by the HLT.
Three trigger categories are used: hadronic triggers, defined by applying moderate requirements on M R and R for events with at least two jets with p T > 56 GeV; electron triggers, similar to the hadronic triggers, but with looser requirements for M R and R and requiring at least one electron with p T > 10 GeV satisfying loose isolation criteria; and muon triggers, with similar M R and R requirements and at least one muon with jηj < 2.1 and p T > 10 GeV. All these triggers have an efficiency of ð98 AE 2Þ% in the kinematic regions used for the off-line selection.
In addition, control samples are defined using several nonrazor triggers. These include prescaled inclusive hadronic triggers, hadronic multijet triggers, hadronic triggers based on H T , and inclusive electron and muon triggers.

C. Physics object reconstruction
Events are required to have at least one reconstructed interaction vertex [41]. When multiple vertices are found, the one with the highest scalar sum of charged track p 2 T is taken to be the event interaction vertex. Jets are reconstructed off-line from calorimeter energy deposits using the infrared-safe anti-k T [42] algorithm with a distance parameter R ¼ 0.5. Jets are corrected for the nonuniformity of the calorimeter response in energy and η using corrections derived from data and simulations and are required to have p T > 40 GeV and jηj < 3.0 [43]. To match the trigger requirements, the p T of the two leading jets is required to be greater than 60 GeV. The jet energy scale uncertainty for these corrected jets is 5% [43]. The E miss T is defined as the negative of the vector sum of the transverse energies (E T ) of all the particles found by the particle-flow algorithm [44].
Electrons are reconstructed using a combination of shower shape information and matching between tracks and electromagnetic clusters [45]. Muons are reconstructed using information from the muon detectors and the silicon tracker and are required to be consistent with the reconstructed primary vertex [46].
The selection criteria for electrons and muons are considered to be tight if the electron or muon candidate is isolated, satisfies the selection requirements of Ref. [47], and lies within jηj < 2.5 and jηj < 2.1, respectively. Loose electron and muon candidates satisfy relaxed isolation requirements.

D. Selection of good quality data
The 4.7 fb −1 integrated luminosity used in this analysis is certified as having a fully functional detector. Events with various sources of noise in the ECAL or HCAL detectors are rejected using either topological information, such as unphysical charge sharing between neighboring channels, or timing and pulse shape information. The last requirement exploits the difference between the shapes of the pulses that develop from particle energy deposits in the calorimeters and from noise events [48]. Muons produced from proton collisions upstream of the detector (beam halo) can mimic proton-proton collisions with large E miss T and are identified using information obtained from the CSCs. The geometry of the CSCs allows efficient identification of beam halo muons, since halo muons that traverse the calorimetry will mostly also traverse one or both CSC end caps. Events are rejected if a significant amount of energy is lost in the masked crystals that constitute approximately 1% of the ECAL, using information either from the separate readout of the L1 hardware trigger or by measuring the energy deposited around the masked crystals. We select events with a well-reconstructed primary vertex and with the scalar P p T of tracks associated to it greater than 10% of the scalar P p T of all jet transverse momenta. These requirements reject 0.003% of an otherwise good inclusive sample of proton-proton interactions (minimum bias events).

E. Event selection and classification
Electrons enter the megajet definition as ordinary jets. Reconstructed muons are not included in the megajet grouping because, unlike electrons, they are distinguished from jets in the HLT. This choice also allows the use of WðμνÞ þ jets events to constrain and study the shape of ZðννÞ þ jets events in fully hadronic final states.
The megajets are constructed as the sum of the fourmomenta of their constituent objects. After considering all possible partitions into two megajets, the combination is selected that has the smallest sum of megajet squaredinvariant-mass values.
The variables M R and R 2 are calculated from the megajet four-momenta. The events are assigned to one of the six final state boxes according to whether the event has zero, one, or two isolated leptons, and according to the lepton flavor (electrons and muons), as shown in Table I. The lepton p T , M R , and R 2 thresholds for each of the boxes are chosen so that the trigger efficiencies are independent of M R and R 2 .
The requirements given in Table I determine the full analysis regions of the (M R , R 2 ) plane for each box. These regions are large enough to allow an accurate characterization of the background, while maintaining efficient triggers. To prevent ambiguities when an event satisfies the selection requirements for more than one box, the boxes are arranged in a predefined hierarchy. Each event is uniquely assigned to the first box whose criteria the event satisfies. Table I shows the box-filling order followed in the analysis.
Six additional boxes are formed with the requirement that at least one of the jets with p T > 40 GeV and jηj < 3.0 be tagged as a b jet, using an algorithm that orders the tracks in a jet by their impact parameter significance and discriminates using the track with the second-highest significance [49]. This algorithm has a tagging efficiency of about 60%, evaluated using b jets containing muons from semileptonic decays of b hadrons in data, and a misidentification rate of about 1% for jets originating from u, d, and s quarks or from gluons, and of about 10% for jets coming from c quarks [49]. The combination of these six boxes defines an inclusive event sample with an enhanced heavy-flavor content.

IV. SIGNAL AND STANDARD MODEL BACKGROUND MODELING
The razor analysis is guided by studies of MC event samples generated with the PYTHIA v6.426 [50] (with Z2 tune) and MADGRAPH v4.22 [51] programs, using the CTEQ6 parton distribution functions (PDFs) [52]. Events generated with MADGRAPH are processed with PYTHIA [50] to provide parton showering, hadronization, and the underlying event description. The matrix element/parton shower matching is performed using the approach described in Ref. [53]. Generated events are processed with the GEANT4 [54] based simulation of the CMS detector, and then reconstructed with the same software used for data.
The simulation of the tt, W þ jets, Z þ jets, single-top (s, t, and t-W channels), and diboson samples is performed using MADGRAPH. The events containing top-quark pairs are generated accompanied by up to three extra partons in the matrix-element calculation [55]. Multijet samples from QCD processes are produced using PYTHIA.
To generate SUSY signal MC events in the context of the CMSSM, the mass spectrum is first calculated with the SOFTSUSY program [56] and the decays with the SUS-HIT [57] package. The PYTHIA generator is used with the SUSY Les Houches Accord (SLHA) interface [58] to generate the events. The generator-level cross sections and the K factors for the next-to-leading-order (NLO) cross sections are computed using PROSPINO [59].
We also use SMS MC simulations in the interpretation of the results. In an SMS simulation, a limited set of hypothetical particles is introduced to produce a given topological signature. The amplitude describing the production and decay of these particles is parametrized in terms of the particle masses. Compared with the constrained SUSY models, SMS provide benchmarks that focus on one final-state topology at a time, with a broader variation in the masses determining the final-state kinematics. The SMS are thus useful for comparing search strategies as well as for identifying challenging areas of parameter space where search methods may lack sensitivity. Furthermore, by providing a tabulation of both the signal acceptance and the 95% confidence level (C.L.) exclusion limit on the signal cross section as a function of the SMS mass parameters, SMS results can be used to place limits on a wide variety of theoretical models beyond SUSY. TABLE I. Definition of the full analysis regions for the mutually exclusive boxes, based on the M R and R 2 values, and, for the categories with leptons, on their p T value, listed according to the hierarchy followed in the analysis, the ELE-MU (HAD) being the first (last).

Lepton boxes
The considered SMS scenarios produce multijet final states with or without leptons and b-tagged jets [39]. While the SUSY terminology is employed, interpretations of SMS scenarios are not restricted to SUSY scenarios.
In the SMS scenarios considered here, each produced particle decays directly to the LSP and SM particles through a two-body or three-body decay. Simplified models that are relevant to inclusive hadronic jets þ E miss T analyses are gluino pair production with the direct three-body decaỹ g → qqχ 0 1 (T1), and squark-antisquark production with the direct two-body decayq → qχ 0 1 (T2). For b-quark enriched final states, we have considered two additional gluino SMS scenarios, where each gluino is forced into the three-body decayg → bbχ 0 1 with 100% branching fraction (T1bbbb), or where each gluino decays throughg → ttχ 0 1 (T1tttt). For b-quark enriched final states we also consider SMS that describe the direct pair production of bottom or top squarks, with the two-body decaysb → bχ 0 1 (T2bb) and t → tχ 0 1 (T2tt). Note that first-generationqq production (unlikeqq Ã production) is not part of the simplified models used for the interpretation of the razor results, even though it is often the dominant process in the CMSSM for low values of the scalar-mass parameter m 0 . This is because of the additional complication that the production rate depends on the gluino mass. However, the acceptance forqq production is expected to be somewhat higher than forqq Ã , so the limits from T2 can be conservatively applied toqq production with analogous decays.
For each SMS, simulated samples are generated for a range of masses of the particles involved, providing a wider spectrum of mass spectra than allowed by the CMSSM. A minimum requirement of Oð100 GeVÞ on the mass difference between the mother particle and the LSP is applied, to remove phase space where the jets from superpartner decays become soft and the signal is detected only when it is given a boost by associated jet production. By restricting attention to SMS scenarios with large mass differences, we avoid the region of phase space where accurate modeling of initial-and final-state radiation from quarks and gluons is required, and where the description of the signal shape has large uncertainties.
The production of the primary particles in each SMS is modeled with SUSY processes in the appropriate decoupling limit of the other superpartners. In particular, forqq Ã production, the gluino mass is set to a very large value so that it has a minimal effect on the kinematics of the squarks. The mass spectrum and decay modes of the particles in a specific SMS point are fixed using the SLHA input files, which are processed with PYTHIA v6.426 with Tune Z2 [60,61] to produce signal events as an input to a parametrized fast simulation of the CMS detector [62], resulting in simulated samples of reconstructed events for each choice of masses for each SMS. These samples are used for the direct calculation of the signal efficiency, and together with the background model are used to determine the 95% C.L. upper bound on the allowed production cross section.

V. STANDARD MODEL BACKGROUNDS IN THE (M R , R 2 ) RAZOR PLANE
The distributions of SM background events in both the MC simulations and the data are found to be described by the sum of exponential functions of M R and R 2 over a large part of the (M R , R 2 ) plane. Spurious instrumental effects and QCD multijet production are challenging backgrounds due to difficulties in modeling the high p T and E miss T tails. Nevertheless, these event classes populate predictable regions of the (M R , R 2 ) plane, which allows us to study them and reduce their contribution to negligible levels. The remaining backgrounds in the lepton, dilepton, and hadronic boxes are processes with genuine E miss T due to energetic neutrinos and charged leptons from vector boson decay, including W bosons from top-quark and diboson production. The analysis uses simulated events to characterize the shapes of the SM background distributions, determine the number of independent parameters needed to describe them, and to extract initial estimates of the values of these parameters. Furthermore, for each of the main SM backgrounds a control data sample is defined using ≈250 pb −1 of data collected at the beginning of the run. These events cannot be used in the search, as the dedicated razor triggers were not available. Instead, events in this run period were collected using inclusive nonrazor hadronic and leptonic triggers, thus defining kinematically unbiased data control samples. We use these control samples to derive a data-driven description of the shapes of the background components and to build a background representation using statistically independent data samples; this is used as an input to a global fit of data selected using the razor triggers in a signal-free region of the (M R , R 2 ) razor plane.
The two-dimensional probability density function P j ðM R ; R 2 Þ describing the R 2 versus M R distribution of each considered SM process j is found to be well approximated by the same family of functions F j ðM R ; R 2 Þ: where k j , M 0 R;j , and R 2 0;j are free parameters of the background model. After applying a baseline selection in the razor kinematic plane, M R > M min R and R 2 > R 2 min , this function exhibits an exponential behavior in R 2 (M R ), when integrated over M R (R 2 ): The fact that the function in Eq. (5) depends on R 2 and not simply on R motivates the choice of R 2 as the kinematic variable quantifying the transverse imbalance. The values of M 0 R;j , R 0;j , k j , and the normalization constant are floated when fitting the function to the data or simulation samples.
The function of Eq. (5) describes the QCD multijet, the lepton þ jets (dominated by W þ jets and tt events), and the dilepton þ jets (dominated by tt and Z þ jets events) backgrounds in the simulation and data control samples. The initial filtering of the SM backgrounds is performed at the trigger level and the analysis proceeds with the analytical description of the SM backgrounds.

A. QCD multijet background
The QCD multijet control sample for the hadronic box is obtained using events recorded with prescaled jet triggers. The trigger used in this study requires at least two jets with average uncorrected p T thresholds of 60 GeV. The QCD multijet background samples provide ≳95% of the events with low M R , allowing the study of the M R shapes with different thresholds on R 2 , which we denote R 2 min . The study was repeated using data sets collected with many jet trigger thresholds and prescale factors during the course of the 2011 LHC data taking, with consistent results.
The M R distributions for events satisfying the HAD box selection in this multijet control data sample are shown for different values of the R min threshold in Fig. 3(a). The M R distribution is exponentially falling, except for a turn-on at low M R resulting from the p T threshold requirement on the jets entering the megajet calculation. The exponential region of these distributions is fitted for each value of R 2 min to extract the absolute value of the coefficient in the exponent, denoted S. The value of S that maximizes the likelihood in the exponential fit is found to be a linear function of R 2 min , as shown in Fig 3(b). Fitting S to the form S ¼ −a − bR 2 min determines the values of a and b. The R 2 min distributions are shown for different values of the M R threshold in Fig. 4(a). The R 2 distribution is exponentially falling, except for a turn-on at low R 2 . The exponential region of these distributions is fitted for each value of M min R to extract the absolute value of the coefficient in the exponent, denoted by S 0 . The value of S 0 that maximizes the likelihood in the exponential fit is found to be a linear function of M min R as shown in Fig. 4 R determines the values of c and d. The slope d is found to be equal to the slope b to within a few percent, as seen from the values of these parameters listed in Figs. 3(b) and 4(b), respectively. The equality d ¼ b is essential for building the 2D probability density function that analytically describes the R 2 versus M R distribution, as it reduces the number of possible 2D functions to the function given in Eq. (5). Note that in Eq. (5) the k j parameters are the b j , d j parameters used in the description of the SM backgrounds.

B. Lepton þ jets backgrounds
The major SM backgrounds with leptons and jets in the final state are ðW=ZÞ þ jets, tt, and single-top-quark production. These events can also contain genuine E miss T . In both the simulated and the data events in the MU and ELE razor boxes, the M R distribution is well described by the sum of two exponential components. One component, which we denote the "first" component, has a steeper slope than the other, "second" component, i.e., jS 1 j > jS 2 j, and thus the second component is dominant in the high-M R region. The relative normalization of the two components is considered as an additional degree of freedom. Both the S 1 and S 2 values, along with their relative and absolute normalizations, are determined in the fit. The M R distributions are shown as a function of R 2 min in Fig. 5 for the zero b-jet MU data, which is dominated by W þ jets events. The dependence of S 1 and S 2 on R 2 min is shown in Fig. 6.
The corresponding results from simulation are shown in Figs. 7 and 8. It is seen that the values of the slope parameters b 1 and b 2 from simulation, given in Fig. 8, agree within the uncertainties with the results from data, given in Fig. 6.
The R 2 distributions as a function of M min R for the data are shown in Fig. 9 for the MU box with the requirement of zero b-tagged jets. The S 0 1 and S 0 2 parameters characterizing the exponential behavior of the first and second WðμνÞ þ jets components are shown in Fig. 10. The corresponding    results from simulation are shown in Figs. 11 and 12. The results for the slopes d 1 and d 2 from simulation, listed in Fig. 12, are seen to be in agreement with the measured results, listed in Fig. 10. Furthermore, the extracted values of d 1 and d 2 are in agreement with the extracted values of b 1 and b 2 , respectively. This is the essential ingredient to build a 2D template for the (M R ,R 2 ) distributions, starting with the function of Eq. (5).
The corresponding distributions for the tt MC simulation with ≥1 b-tagged jet are presented in Appendix A, for events selected in the HAD box.

C. Dilepton backgrounds
The SM contributions to the ELE-ELE and MU-MU boxes are expected to be dominated by Z þ jets events, and the SM contribution to the ELE-MU box by tt events, all at the level of ≳95%. We find that the M R distributions as a function of R 2 min , and the R 2 distribution as a function of M min R , are independent of the lepton-flavor combination for both the ELE-ELE and MU-MU boxes, as determined using simulated ttð2l2ν þ jetsÞ events. In addition, the asymptotic second component is found to be process independent.

VI. Background model and fits
As described earlier, the full 2D SM background representation is built using statistically independent data control samples. The parameters of this model provide the input to the final fit performed in the fit region (FR) of the data samples, defining an extended, unbinned maximum likelihood (ML) fit with the ROOFIT fitting package [63]. The fit region is defined for each of the razor boxes as the region of low M R and small R 2 , where signal contamination is expected to have negligible impact on the shape fit. The 2D model is extrapolated to the rest of the (M R , R 2 ) plane, which is sensitive to new-physics signals and where the search is performed.
For each box, the fit is conducted in the signal-free FR of the (M R , R 2 ) plane; their definition can be found in Figs. 15, 17, 19, 21, 23, and 25. These regions are used to provide a full description of the SM background in the entire (M R , R 2 ) plane in each box. The likelihood function for a given box is written as [64] where N is the total number of events in the FR region of the box, the sum runs over all the SM processes relevant for that box, and the N j are normalization parameters for each SM process involved in the considered box. We find that each SM process in a given final-state box is well described in the (M R , R 2 ) plane by the function P j defined as where the first (F 1st j ) and second (F 2nd j ) components are defined as in Eq. (5), and f j 2 is the normalization fraction of the second component with respect to the total. When fitting this function to the data, the shape parameters of each F j ðM R ; R 2 Þ function, the absolute normalization, and the relative fraction f j 2 are floated in the fit. Studies of simulated events and fits to data control samples with either a b-jet requirement or a b-jet veto indicate that the parameters corresponding to the first components of these backgrounds (with steeper slopes at low M R and R 2 ) are box dependent. The parameters describing the second components are box independent, and at the current precision of the background model, they are identical between the dominant backgrounds considered in these final states.
We validate the choice of the background shape by use of a sample of tt MC simulated events corresponding to an integrated luminosity of 10 fb −1 . Besides being the dominant background in the ≥1 b-tag search, tt events are the dominant background for the inclusive search for large values of M R and R 2 . The result for the HAD box in the inclusive razor path is shown in Fig. 13 expressed as the projection of the 2D fit on M R and R 2 . As the same level of agreement is found in all boxes both in the inclusive and in the ≥1 b-tagged razor path, we proceed to fit all the SM processes with this shape.

A. Fit results and validation
The shape parameters in Eq. (5) are determined for each box via the 2D fit. The likelihood of Eq. (8) is multiplied by Gaussian penalty terms [65] to account for the uncertainties of the shape parameters k j , M 0 R;j , and R 2 0;j . The central values of the Gaussians are derived from analogous 2D fits in the low-statistics data control sample. The penalty terms pull the fit to the local minimum closer to the shape derived from the data control samples. Using pseudoexperiments, we verified that this procedure does not bias the determination of the background shape. As an example, the k j parameter uncertainties are typically ∼30%. Additional background shape uncertainties due to the choice of the functional form were considered and found to be negligible, as discussed in Appendix B.
The result of the ML fit projected on M R and R 2 is shown in Fig. 14 for the inclusive HAD box. No significant discrepancy is observed between the data and the fit model for any of the six boxes. In order to establish the compatibility of the background model with the observed data set, we define a set of signal regions (SRi) in the tail of the SM background distribution. Using the 2D background model determined using the ML fit, we derive the distribution of the expected yield in each SRi using pseudoexperiments, accounting for correlations and uncertainties in the parameters describing the background model. In order to correctly account for the uncertainties in the parameters describing the background model and their correlations, the shape parameters used to generate each pseudoexperiment data set are sampled from the covariance matrix returned by the ML fit performed on the actual data set. The actual number of events in each data set is drawn from a Poisson distribution centered on the yield returned by the covariance matrix sampling. For each pseudoexperiment data set, the number of events in the SRi is found. For each of the SRi, the distribution of the number of events derived by the pseudoexperiments is used to calculate a two-sided p-value (as shown for the HAD box in Fig. 15), corresponding to the probability of observing an equal or less probable outcome for a counting experiment in each signal region.   Fig. 15) and projected into the full analysis region. The full uncertainty in the total background prediction is drawn in these projections, including the one due to the variation of the background shape parameters and normalization.    Fig. 19) and projected into the full analysis region. The histograms are described in the text.  corresponding results for the ≥1 b-tagged samples are presented in Appendix C.

VII. SIGNAL SYSTEMATIC UNCERTAINTIES
We evaluate the impact of systematic uncertainties on the shape of the signal distributions, for each point of each SUSY model, using the simulated signal event samples. The following systematic uncertainties are considered, with the approximate size of the uncertainty given in parentheses: (i) PDFs (up to 30%, evaluated point by point) [66]; (ii) jet-energy scale (up to 1%, evaluated point by point) [43]; (iii) lepton identification, using the "tag-and-probe" technique based on Z → ll events [67] (l ¼ e; μ, 1% per lepton). In addition, the    Fig. 25) and projected into the full analysis region. The histogram is described in the text. following uncertainties, which affect the signal yield, are considered: (i) luminosity uncertainty [68] (2.2%); (ii) theoretical cross section [69] (up to 15%, evaluated point by point); (iii) razor trigger efficiency (2%); (iv) lepton trigger efficiency (3%). An additional systematic uncertainty is considered for the b-tagging efficiency [49] (between 6% and 20% in p T bins). We consider variations of the function modeling, the signal uncertainty (log normal versus Gaussian), and the binning, and find negligible deviations in the results. The systematic uncertainties are included using the best-fit shape to compute the likelihood values for each pseudoexperiment, while sampling the same pseudoexperiment from a different function, derived from the covariance matrix of the fit to the data. This procedure is repeated for both the background and signal probability density functions.

VIII. INTERPRETATION OF THE RESULTS
In order to evaluate exclusion limits for a given SUSY model, its parameters are varied and an excluded cross section at the 95% C.L. is associated with each configuration of the model parameters, using the hybrid version of the CL s method [70][71][72], described below. For each box, we consider the test statistic given by the logarithm of the likelihood ratio ln Q ¼ ln½Lðs þ bjH i Þ=LðbjH i Þ, where H i ði ¼ 0; 2Þ is the hypothesis under test: H 1 (signal plus background) or H 0 (background only). The likelihood function for the background-only hypothesis is given by Eq. (8). The likelihood corresponding to the signal-plus-background hypothesis is written as where σ is the signal cross section, i.e., the parameter of interest; L is the integrated luminosity; ϵ is the signal acceptance times efficiency; and P S ðM R;i ; R 2 i Þ is the two-dimensional probability density function for the signal, computed numerically from the distribution of simulated signal events. The signal and background shape parameters, and the normalization factors L and ϵ, are the nuisance parameters.
For each analysis (inclusive razor or inclusive b-jet razor) we sum the test statistics of the six corresponding boxes to compute the combined test statistic.
The distribution of ln Q is derived numerically with an MC technique. The values of the nuisance parameters in the likelihood are randomized for each iteration of the MC generation, to reflect the corresponding uncertainty. Once the likelihood is defined, a sample of events is generated according to the signal and background probability density functions. The value of ln Q for each generated sample is then evaluated, fixing each signal and background parameter to its expected value. This procedure corresponds to a numerical marginalization of the nuisance parameters.
Given the distribution of ln Q for the background-only and the signal-plus-background pseudoexperiments, and FIG. 27. The diagrams corresponding to the SMS models considered in this analysis. the value of ln Q observed in the data, we calculate CL sþb and 1 − CL b [70]. From these values, CL s ¼ CL sþb =CL b is computed for that model point. The procedure is independently applied to each of the two analyses (inclusive razor and inclusive b-jet razor).
The CMSSM model is studied in the (m 0 , m 1=2 ) plane, fixing tan β ¼ 10, A 0 ¼ 0, and sgnðμÞ ¼ þ1. A point in the plane is excluded at the 95% C.L. if CL s < 0.05. The result obtained for the inclusive razor analysis is shown in Fig. 26(a). The shape of the observed exclusion curves reflects the changing relevant SUSY strong-production processes across the parameter space, with squark-antisquark and gluino-gluino production dominating at low and high m 0 , respectively. The observed limit is less constraining than the median expected limit at lower m 0 due to a local excess of events at large R 2 in the hadronic box.
For large values of m 0 , boxes with leptons in the final state have a sensitivity comparable to that of the hadronic boxes, as cascade decays of gluinos yield leptons production.  show the CMSSM exclusion limits based on the HAD box only and on the leptonic boxes only.
The results are also interpreted as cross section limits on a number of simplified models [32][33][34][35][36] where a limited set of hypothetical particles and decay chains are introduced to produce a given topological signature. For each model studied, we derive the maximum allowed cross section at the 95% C.L. as a function of the mass of the produced particles (gluinos or squarks, depending on the model) and the LSP mass, as well as the exclusion limit corresponding to the SUSY cross section. We study several SMS benchmark scenarios [39]: The solid gray region indicates model points where the selection efficiency is found to have dependence on ISR modeling in the simulation of signal events above a predefined tolerance; no interpretation is presented for these model points.
(i) gluino-gluino production with four light-flavor jetsþE miss T in the superpartner decays, T1 in Fig. 27. (ii) squark-antisquark production with two light-flavor jets þ E miss T in the superpartner decays, T2 in Fig. 27. (iii) gluino-gluino production with four b jets þ E miss T in the superpartner decays, T1bbbb in Fig. 27. (iv) squark-antisquark production with two b jets þ E miss T in the superpartner decays, T2bb in Fig. 27. (v) gluino-gluino production with four top quarks þ E miss T in the superpartner decays, T1tttt in Fig. 27. (vi) squark-antisquark production with two top quarks þ E miss T in the superpartner decays, T2tt in Fig. 27. In all cases, additional jets in the final state can arise from initial-and final-state radiation (ISR and FSR), simulated by PYTHIA. We show in Figs. 28 and 29 the excluded cross section at 95% C.L. as a function of the mass of the produced particle (gluinos or squarks, depending on the model) and the LSP mass, as well as the exclusion curve corresponding to the NLO þ NLL SUSY cross section [74][75][76][77][78], where NLL indicates the next-to-leading logarithmic. A result is not quoted for the region of the SMS plane in which the signal efficiency strongly depends on the ISR and FSR modeling (gray area), as a consequence of the small mass difference between the produced superpartner and the LSP and the consequent small p T for the jets produced in the cascade.
In Fig. 30, we present a summary of the 95% C.L. excluded largest parent mass for various LSP masses in each of the simplified models studied, showing separately the results from the inclusive razor analysis and the inclusive b-jet razor analysis. A comparison of the razor results with those obtained from other approaches is given in Ref. [79].  The solid gray region indicates model points where the selection efficiency is found to have dependence on ISR modeling in the simulation of signal events above a predefined tolerance; no interpretation is presented for these model points.

IX. SUMMARY
Using a data sample of ffiffi ffi s p ¼ 7 TeV proton-proton collisions collected by the CMS experiment at the LHC in 2011, corresponding to an integrated luminosity of 4.7 fb −1 , we have performed a search for pair-produced supersymmetric particles such as squarks and gluinos in the razor-variable plane. A 2D shape description of the relevant standard model processes determined from data control samples and validated with simulated events has been used, and no significant excess over the background expectations has been observed. The results are presented as a 95% C.L. limit in the (m 0 , m 1=2 ) CMSSM parameter space. We exclude squark and gluino masses up to 1350 GeV for mðqÞ ∼ mðgÞ, while for mðqÞ > mðgÞ we exclude gluino masses up to 800 GeV. For simplified models, we exclude gluino masses up to 1000 GeV, and first-and secondgeneration squark masses up to 800 GeV. The direct production of top or bottom squarks is excluded for squark masses up to 400 GeV.

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: BMWF and FWF  Figure 31 shows the M R distribution as a function of R 2 min for tt MC events with ≥1 b-tagged jets in the HAD box. The S 1 and S 2 parameters characterizing the exponential behavior of the first and second WðμνÞ þ jets FIG. 30 (color online). Summary of the 95% C.L. excluded largest parent mass for each of the simplified models studied, for various LSP masses. The results from the b-jet razor analysis are shown immediately below those from the inclusive razor analysis for each of the four categories of events indicated.   components are shown in Fig. 32. The corresponding distributions for R 2 , and for the S 0 1 and S 0 2 parameters, are shown in Figs. 33 and 34, respectively. The conclusions derived from the data and MC studies of Sec. V hold also for tt MC events.

APPENDIX B: ALTERNATIVE BACKGROUND SHAPE ANALYSIS
In order to quantify a systematic uncertainty associated with the choice of the fit function, we first generalize our 2D function to allow for deviations from the exponential behavior, once projected onto M R or R 2 . To do this, we (i) identify a set of functions that describe the data, (ii) use one as a default description, (iii) use the rest to quantify the systematic variation, (iv) randomly choose one of the three functions when generating the pseudoexperiments used to set limits, and (v) use the nominal function when evaluating the likelihood.
For a 1D fit of the M R distribution, an obvious choice is where β ≠ 1 accounts for deviations from the exponential function. In this analysis, we need a 2D function of M R and R 2 that allows us to measure the deviation from the nominal shape on the projections. For this purpose, we introduce a generalization of the razor 2D function: which has the two following properties: where with M min R and R 2 min respectively the thresholds applied on M R and R 2 before projecting onto R 2 and M R . Using this function to evaluate systematic uncertainties corresponds to the 2D generalization needed here. We proceed as follows: (i) we repeat the fit in the fit region of each box, using F SYS ðM R ; R 2 Þ rather than FðM R ; R 2 Þ for the second component of the background model (the one that extrapolates to the signal region), with n floated in the fit. We determine n fit AE σ n in this fit. (ii) we assign an allowed range to the difference n − 1 taking the larger of n fit − 1 and σ n , which we refer to as ½n min ; n max . (iii) we repeat the fit in the fit region fixing n to first to n min and then to n max and we take these fits as the alternative background descriptions. In particular, we find that the fit returns values of n fit that are very close to n. Following the prescription outlined above, we take the fit uncertainty as the shift in n.
The main conclusion of the study is that the systematic uncertainty in the choice of the function is already covered by the large uncertainty in the fit parameters and that the effect corresponds to an increase of about 15% in the 68% C.L. range, once this contribution is summed in quadrature with the already quoted uncertainty.
As an example, we present the results of the above procedure for the bins in the HAD box. Figure 35 shows the fit result with n floated in the full region of the HAD box, projected onto M R and R 2 . The quality of the fit is similar to that of the nominal procedure. We find n ¼ 0.96 AE 0.04. We then take n min ¼ 0.96 and n max ¼ 1.04. We show in    Table II the bin-by-bin background prediction for the nominal fit and the two alternative fits. We use a finer binning than the one used to compute the p-values in the nominal analysis. For comparison, we also show the values obtained with n floated in the fit. For all cases, we quote the predicted background as the center of the 68% probability range and the associated uncertainty corresponds to half the range. The range is defined by integrating the background distribution (derived from the pseudoexperiments) using the probability value as the ordering algorithm. Similar results are obtained for all boxes.

APPENDIX D: GUIDE ON EMULATING THE RAZOR ANALYSIS FOR ADDITIONAL STUDIES
In this appendix, we provide a guide to facilitate use of the razor analysis results for the interpretation of signal scenarios not considered here. We assume the existence of an event generator that can simulate LHC collisions for a given theoretical model. We also assume that this event generator is interfaced to a parton shower simulation, such that a list of produced particles at the generator level is available. The procedure described in this appendix represents a simplification of the analysis, giving conservative limits within the AE1 standard deviation band of the nominal result.
The following classes of stable particles are relevant to this analysis: (i) invisible particles (neutrinos and any weakly interacting stable new particles, for example the LSP in SUSY models); (ii) electrons; (iii) muons; (iv) all    other stable electrically charged SM particles; and (v) all other stable electrically neutral SM particles. It is possible to emulate the razor analysis as follows: (i) all the visible stable particles are clustered into generator-level jets using the anti-k T algorithm with a distance parameter of 0.5. (ii) the generator-level E miss T is computed as E miss T;Gen ¼ − P p p p T , where the sum runs over all the visible stable particles p.    (iii) the detector resolution is applied to electrons and muons according to a simplified Gaussian resolution function. The RMS of the Gaussian smearing depends on the η and p T values of the lepton, as well as its flavor. Similarly, the E miss T and jet momenta are smeared according to a Gaussian response model. (iv) the detector efficiency is applied to electrons and muons generating unweighted events from the reconstruction efficiency, interpreted as a probability (see Appendix D 1). The efficiency depends on the η and p T values of the lepton, its flavor, and its generator-level isolation, as computed from the stable particles in the event. (v) the analysis selection and box classification is applied. This procedure allows us to estimate the R 2 versus M R distribution for a signal model and the efficiency in each box. This is the information that is needed to associate a 95% C.L. upper limit to a given input model. The procedure matches the full simulation of CMS to within 20% and in general provides a result that is yet closer to the CMS full simulation. The result is in general conservative, since the computation of the upper limit starts from a simplified binned likelihood, which reduces the sensitivity to a signal. This procedure is not expected to correctly simulate the special case of slowly moving electrically charged particles (e.g., staus). The remainder of this appendix describes each step of the razor emulation in more detail, including the calculation of the exclusion limit.

Emulation of reconstructed electrons and muons
The emulation of reconstructed electrons and muons consists of two independent steps: the accounting for the detector resolution and for the reconstruction efficiency.
The effects of detector resolution can be incorporated through a Gaussian smearing of the genuine p T of a given lepton, while the lepton η and ϕ can be considered to be unaffected by the detector resolution. The generated lepton is then replaced by the reconstructed one, having the same flight direction with a p T value randomly extracted according to a Gaussian distribution centered at p Gen T and with σðp Gen T Þ taken from Fig. 48. Any lepton outside the two η ranges considered in Fig. 48 should be discarded from the analysis.
To account for the reconstruction efficiency of a given lepton, the generator-level isolation is computed as follows: where the sum runs over all the stable charged and neutral visible particles p within a distance ΔR ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi ðΔηÞ 2 þ ðΔϕÞ 2 p < 0.5 from the lepton. Figure 49 shows the reconstruction probability versus the generated electron p T (before accounting for the detector resolution) for three ranges of GenIso in the ECAL barrel (jηj < 1.4442) and end cap (1.5660 < jηj < 2.5000) regions. Different values are obtained for the tight and the loose electrons used to define the boxes.
Similarly, the reconstruction efficiency for the tight muons is shown in Fig. 50. The reconstruction of loose muons can be considered to be fully efficient for muons with p T > 10 GeV, since no isolation requirement is applied.
Once the lepton reconstruction probability is found, the detector efficiency effects can be imposed numerically: the lepton is rejected if a uniformly distributed random number in the range [0,1] is found to be larger than the reconstruction efficiency.

T
The reconstruction of jets and E miss T can be emulated by applying a Gaussian resolution to the generator-level quantities. We show in Fig. 51 the dependence of the Gaussian σ jet on the jet p T (for the two relevant bins of η) and the E miss T . The dependence on η or other quantities can be safely neglected. One should apply the resolution function to all the reconstructed jets and to the E miss T and then impose the acceptance selection on the reconstructed jets.

Building the 2D templates
Once detector effects have been accounted for, jets are clustered in two megajets. The razor variables can be computed from the four-momenta of the two megajets. Figure 52 (left and middle panels) shows the M R and R 2 distributions for a sample of pair-produced gluinos of mass 800 GeV, where each gluino decays to a tt pair and a LSP of mass 300 GeV, obtained with the CMS fast simulations program and with the emulation described in this appendix. The efficiencies obtained for the six boxes are compared in Fig. 52 (right panel).

Evaluating the exclusion limit
The exclusion limit can be computed from the 2D signal templates and the box efficiencies, starting with the observed yield and the expected background. We consider  a simplified likelihood obtained by defining bins in the (M R , R 2 ) plane. Each bin i requires the observed yield n i and the expected backgroundb i AE δ i computed by integrating the background model and taking into account the uncertainty in shape. The likelihood in a given box is then written as where ϵ i is the signal efficiency in that bin, L is the luminosity, and σ is the signal cross section; log N ðb i jb i ; δ i Þ is the log-normal distribution describing the uncertainty in the background. log N ðρj1; δ ρ Þ is the distribution describing the uncertainty in the signal efficiency. A value δ ρ ∼ 0.20 (including the uncertainty in the integrated luminosity) is large enough to account for the use of a simplified detector emulation and the typical systematic uncertainty quoted in the analysis. Once this uncertainty is included, the uncertainty in the luminosity can be neglected to a good level of precision. Similarly, the total likelihood can be written as In this case, the signal systematic parameter ρ is common to the six boxes. A Bayesian upper limit (UL) on the cross section can then be computed assuming a flat prior distribution in σ:  An implementation of this simplified limit calculator is provided in the supplemental material [80] together with the values of n, b, and δ for each bin in each box. Figure 53 shows the limit on the T2tt and T1tttt models, obtained by applying the simplified procedure described in this appendix. We generate a sample of SUSY events using the PYTHIA 8 [81] program, scanning the two SMS planes. We then emulate the detector effects as described in this appendix to derive the efficiency and the (M R , R 2 ) signal probability density functions. We use this information to compute the excluded cross section for each point in the SMS plane.