A search for new phenomena in pp collisions at √ s = 13 TeV in final states with missing transverse momentum and at least one jet using the α T variable CMS

A search for new phenomena is performed in final states containing one or more jets and an imbalance in transverse momentum in pp collisions at a centre-of-mass energy of 13[Formula: see text]. The analysed data sample, recorded with the CMS detector at the CERN LHC, corresponds to an integrated luminosity of 2.3[Formula: see text]. Several kinematic variables are employed to suppress the dominant background, multijet production, as well as to discriminate between other standard model and new physics processes. The search provides sensitivity to a broad range of new-physics models that yield a stable weakly interacting massive particle. The number of observed candidate events is found to agree with the expected contributions from standard model processes, and the result is interpreted in the mass parameter space of fourteen simplified supersymmetric models that assume the pair production of gluinos or squarks and a range of decay modes. For models that assume gluino pair production, masses up to 1575 and 975[Formula: see text] are excluded for gluinos and neutralinos, respectively. For models involving the pair production of top squarks and compressed mass spectra, top squark masses up to 400[Formula: see text] are excluded.


Introduction
The standard model (SM) of particle physics is successful in describing a wide range of phenomena, although it is widely believed to be only an effective approximation of a more complete theory that supersedes it at a higher energy scale. Supersymmetry (SUSY) [1][2][3][4] is a modification to the SM that extends its underlying space-time symmetry group. For each boson (fermion) in the SM, a fermionic (bosonic) superpartner, which differs in spin by one-half unit, is introduced.
Experimentally, SUSY is testable through the prediction of an extensive array of new observable states (of unknown masses) [5,6]. In the minimal supersymmetric extension to e-mail: cms-publication-committee-chair@cern.ch the SM [6], the gluinos g, light-and heavy-flavour squarks q, b, t, and sleptons are, respectively, the superpartners to gluons, quarks, and leptons. An extended Higgs sector is also predicted, as well as four neutralino χ 0 1,2,3,4 and two chargino χ ± 1,2 states that arise from mixing between the higgsino and gaugino states, which are the superpartners of the Higgs and electroweak gauge bosons. The assumption of Rparity conservation [7] has important consequences for cosmology and collider phenomenology. Supersymmetric particles are expected to be produced in pairs at the LHC, with heavy coloured states decaying, potentially via intermediate SUSY states, to the stable lightest SUSY particle (LSP). The LSP is generally assumed to be the χ 0 1 , which is weakly interacting and massive. This SUSY particle is considered to be a candidate for dark matter (DM) [8], the existence of which is supported by astrophysical data [9]. Hence, a characteristic signature of R-parity-conserving coloured SUSY production at the LHC is a final state containing an abundance of jets, possibly originating from top or bottom quarks, accompanied by a significant transverse momentum imbalance, p miss T . The proposed supersymmetric extension of the SM is also compelling from a theoretical perspective, as the addition of superpartners to SM particles can modify the running of the gauge coupling constants such that their unification can be achieved at a high energy scale [10][11][12]. A more topical perspective, given the recently discovered Higgs boson [13][14][15], is the possibility that scale-dependent radiative corrections to the Higgs boson mass from loop processes can be largely cancelled through the introduction of superpartners, thus alleviating the gauge hierarchy problem [16,17]. Alternatively, these radiative corrections can be accommodated through an extreme level of fine tuning of the bare Higgs boson mass. A "natural" solution from SUSY, with minimal fine-tuning, implies that the masses of the χ 0 1 , third-generation squarks, and the gluino are at or near the electroweak scale [18].
The lack of evidence to date for SUSY has also focused attention on regions of the natural parameter space with sparse experimental coverage, such as phenomenologically well motivated models for which both the t and the χ 0 1 are light and nearly degenerate in mass [19][20][21][22][23][24][25][26][27]. This class of models, with "compressed" mass spectra, typically yield SM particles with low transverse momenta ( p T ) from the decays of SUSY particles. Hence, searches rely on the associated production of jets, often resulting from initial-state radiation (ISR), to achieve experimental acceptance.
This paper presents an inclusive search for new-physics phenomena in hadronic final states with one or more energetic jets and an imbalance in p miss T , performed in protonproton (pp) collisions at a centre-of-mass energy √ s = 13 TeV. The analysed data sample corresponds to an integrated luminosity of 2.3 ± 0.1 fb −1 collected by the CMS experiment. Earlier searches using the same technique have been performed in pp collisions at both √ s = 7 [28][29][30] and 8 TeV [31,32] by the CMS Collaboration. The increase in the centre-of-mass energy of the LHC, from √ s = 8 to 13 TeV, provides a unique opportunity to search for the characteristic signatures of new physics at the TeV scale. For example, the increase in √ s leads to a factor 35 increase in the parton luminosity [33] for the pair production of coloured SUSY particles, each of mass 1.5 TeV, which were beyond the reach of searches performed at √ s = 8 TeV by the ATLAS [34,35] (and references therein) and CMS [36][37][38][39][40][41][42][43] Collaborations. Several searches in this final state, interpreted within the context of SUSY, have already provided results with the first data at this new energy frontier [44][45][46][47][48][49][50].
Two important features of this search are the application of selection criteria with low thresholds, in order to maximise signal acceptance, and the categorisation of candidate signal events according to multiple discriminating variables for optimal signal extraction over a broad range of models. The search is based on an examination of the number of reconstructed jets per event, the number of these jets identified as originating from bottom quarks, and the scalar and vector p T sums of these jets. These variables provide sensitivity to the different production mechanisms (squark-squark, squarkgluino, and gluino-gluino) of massive coloured SUSY particles at hadron colliders, third-generation squark signatures, and both large and small mass splittings between the parent SUSY particle and the LSP. However, the search is sufficiently generic and inclusive to provide sensitivity to a wide range of SUSY and non-SUSY models that postulate the existence of a stable, only weakly interacting, massive particle. In addition to the jets+ p miss T topology, the search considers final states containing a "monojet" topology, which is expected to improve the sensitivity to DM particle production in pp collisions [51,52].
The dominant background process for a search in all-jet final states is multijet production, a manifestation of quantum chromodynamics (QCD). An accurate estimate of this background is difficult to achieve, given the lack of precise theoretical predictions for the multijet production cross section and kinematic properties. Hence, this search adopts a strategy that employs several variables to reduce the multijet contribution to a low level with respect to other SM backgrounds, rather than estimate a significant contribution with high precision.
The search is built around two variables that are designed to provide robust discrimination against multijet events. A dimensionless kinematic variable α T [28,53] is constructed from jet-based quantities and provides discrimination between genuine sources of p miss T from stable, weakly interacting particles such as neutrinos or neutralinos that escape the detector, and instrumental sources such as the mismeasurements of jet energies. The φ * min [28] variable exploits azimuthal angular information and also provides strong rejection power against multijet events, including rare energetic events in which neutrinos carry a significant fraction of the energy of a jet due to semileptonic decays of heavy-flavour mesons. Very restrictive requirements on the α T and Δφ * min variables are employed in this search to ensure a low level of contamination from the multijet background.
The organisation of this paper is as follows. Sections 2 and 3 describe, respectively, the CMS apparatus and the simulated event samples. Sections 4 and 5 describe the event reconstruction and selection criteria used to identify candidate signal events and control region samples. Section 6 provides details on the estimation of the multijet and all other SM backgrounds. Finally, the search results and interpretations, in terms of simplified SUSY models, are described in Sects. 7 and 8, and summarised in Sect. 9.

The CMS detector
The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing an axial magnetic field of 3.8 T. The bore of the solenoid is instrumented with several particle detection systems. A silicon pixel and strip tracker measures charged particles within the pseudorapidity range |η| < 2.5. A lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter (HCAL), each composed of a barrel and two endcap sections, extend over a range |η| < 3.0. Outside the bore of the solenoid, forward calorimeters extend the coverage to |η| < 5.0, and muons are measured within |η| < 2.4 by gasionisation detectors embedded in the steel flux-return yoke outside the solenoid. A two-tier trigger system selects pp collision events of interest. The first level of the trigger system, 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 processor farm further decreases the event rate from around 100 kHz to less than 1 kHz, before data storage. The CMS detector is nearly hermetic, which allows for momentum balance measurements in the plane transverse to the beam axis. 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. [54].

Simulated event samples
The search relies on multiple event samples, in data or generated from Monte Carlo (MC) simulations, to estimate the contributions from SM backgrounds, as described in Sect. 5. The SM backgrounds for the search are QCD multijet, top quark-antiquark (tt), and single top production, and the associated production of jets and a vector boson (W, Z → νν). Residual contributions from other processes, such as WW, WZ, ZZ (diboson) production and the associated production of tt and a W or Z boson, are also considered. Other processes, such as Drell-Yan (qq → Z/γ * → + − ) and γ + jets production, are also relevant for some control regions, defined in Sect. 5.5.
The MadGraph5_amc@nlo 2.2.2 [55] event generator code is used at leading-order (LO) accuracy to produce samples of W + jets, Z + jets, γ + jets, tt, and multijet events. The same code is used at next-to-leading-order (NLO) accuracy to generate samples of single top quarks (both s-and t-channel production), WZ, ZZ, ttW, and ttZ events. The NLO powheg v2 [56,57] generator is used to describe WW events and the tW-channel production of single top quark events. The simulated samples are normalised according to production cross sections that are calculated with NLO and next-to-NLO precision [55,[57][58][59][60][61][62], or with LO precision in the case of multijet and γ +jets production. The Geant4 [63] package is used to simulate the detector response.
Event samples for signal models involving gluino or squark pair production in association with up to two additional partons are generated at LO with MadGraph5_a mc@nlo, and the decay of the SUSY particles is performed with pythia 8.205 [64]. Inclusive, process-dependent, signal production cross sections are calculated with NLO plus nextto-leading-logarithmic (NLL) accuracy [33,[65][66][67][68][69]. The theoretical systematic uncertainties are typically dominated by the parton density function (PDF) uncertainties, evaluated using the CTEQ6.6 [70] and MSTW2008 [71] PDFs. The detector response for signal models is provided by the CMS fast simulation package [72].
The NNPDF3.0 LO and NLO [73] parton distribution functions (PDFs) are used, respectively, with the LO and NLO generators described above. The pythia program with the CUETP8M1 underlying event tune [74] is used to describe parton showering and hadronisation for all simulated samples. To model the effects of multiple pp collisions within the same or neighbouring bunch crossings (pileup), all simulated events are generated with a nominal distribution of pp interactions per bunch crossing and then reweighted to match the pileup distribution as measured in data. On average, approximately fifteen different pp collisions, identifiable via their primary interaction vertex, are reconstructed per event. Finally, (near-unity) corrections to the normalisation of the simulated samples for the γ +jets, W(→ μν)+jets, tt, Z(→ μμ) + jets, and Z(→ νν) + jets processes are derived using data sidebands to the control regions.

Event reconstruction
Global event reconstruction is provided by the particle flow (PF) algorithm [75,76], designed to identify each particle using an optimised combination of information from all detector systems. In this process, the identification of the particle type (photon, electron, muon, charged hadron, neutral hadron) plays an important role in the determination of the particle direction and energy.
Among the vertices reconstructed within 24 (2) cm of the detector centre parallel (perpendicular) to the beam axis, the primary vertex (PV) is assigned to be the one with the largest sum of charged particle (track) p 2 T values. Chargedparticle tracks associated with reconstructed vertices from pileup events are not considered by the PF algorithm as part of the global event reconstruction.
Photon candidates [77] are identified as ECAL energy clusters not linked to the extrapolation of any track to the ECAL. The energy of photons is directly obtained from the ECAL measurement, corrected for contributions from pileup events. Various quality-related criteria must be satisfied in order to identify photons with high efficiency while minimising the misidentification of electrons and associated bremsstrahlung, jets, or ECAL noise as photons. The criteria include the following: the shower shape of the energy deposition in the ECAL must be consistent with that expected from a photon, the energy detected in the HCAL behind the photon shower must not exceed 5% of the photon energy, and no matched hits in the pixel tracker must be found.
Electron candidates [78] are identified as a track associated with an ECAL cluster compatible with the track trajectory, as well as additional ECAL energy clusters from potential bremsstrahlung photons emitted as the electron traverses material of the silicon tracker. The energy of electrons is determined from a combination of the track momentum at the main interaction vertex, the corresponding ECAL cluster energy, and the energy sum of all bremsstrahlung photons associated with the track. The quality criteria required for electrons are similar to those for photons, with regards to the ECAL shower shape and the relative contributions to the total energy deposited in the ECAL and HCAL. Additional requirements are also made on the associated track, which consider the track quality, energy-momentum matching, and compatibility with the PV in terms of the transverse and longitudinal impact parameters.
Muon candidates [79] are identified as a track in the silicon tracker consistent with either a track or several hits in the muon system. The track and hit parameters must satisfy various quality-related criteria, described in Ref. [79]. The energy of muons is obtained from the corresponding track momentum.
Charged hadrons are identified as tracks not classified as either electrons or muons. The energy of charged hadrons is determined from a combination of the track momentum and the corresponding ECAL and HCAL energies, corrected for contributions from pileup events and the response function of the calorimeters to hadronic showers. Neutral hadrons are identified as HCAL energy clusters not linked to any chargedhadron trajectory, or as ECAL and HCAL energy excesses with respect to the expected charged-hadron energy deposit. The energy of neutral hadrons is obtained from the corresponding corrected ECAL and HCAL energy.
Photons are required to be isolated from other activity in the event, such as charged and neutral hadrons, within a cone ΔR = √ (Δφ) 2 + (Δη) 2 = 0.3 around the photon trajectory, corrected for contributions from pileup events and the photon itself. Electrons and muons are also required to be isolated from other reconstructed particles in the event, primarily to suppress background contributions from semileptonic heavy-flavour decays in multijet events. The isolation I mini rel is defined as the scalar p T sum of all charged and neutral hadrons, and photons, within a cone around the lepton direction, divided by the lepton p T . The "mini" cone radius is dependent on the lepton p T , primarily to identify with high efficiency the collimated daughter particles of semileptonically decaying Lorentz-boosted top quarks, according to the following: R = 0.2 and 0.05 for, respectively, p T < 50 GeV and p T > 200 GeV, and R = 10 GeV/ p T for 50 < p T < 200 GeV. The variable I mini rel excludes contributions from the lepton itself and pileup events. The isolation for electrons and muons is required to satisfy, respectively, I mini rel < 0.1 and 0.2 for the signal region and nonleptonic control sample selection criteria. A tighter definition of muon isolation I μ rel is used for the definition of control regions that are required to contain at least one muon. The variable I μ rel is determined identically to I mini rel except that a cone of fixed radius R = 0.4 is assumed.
Electron and muon candidates identified by the PF algorithm that do not satisfy the quality criteria or the I mini rel isolation requirements described above, as well as charged hadrons, are collectively labelled as "single isolated tracks" if they are isolated from neighbouring tracks associated to the PV. The isolation I track rel is defined as the scalar p T sum of tracks (excluding the track under consideration) within a cone ΔR < 0.3 around the track direction, divided by the track p T . The requirement I track rel < 0.1 is imposed. Jets are clustered from the PF candidate particles with the infrared-and collinear-safe anti-k t algorithm [80], operated with a distance parameter of 0.4. The jet momentum is determined as the vectorial sum of all particle momenta in the jet, and is found in the simulation to be within 5-10% of its true momentum over the whole p T spectrum and detector acceptance. Jet energy corrections, to account for pileup [81] and to establish a uniform relative response in η and a calibrated absolute response in p T , are derived from the simulation, and are confirmed with in situ measurements using the energy balance in dijet and photon+jet events [82]. The jet energy resolution is typically 15% at 10 GeV, 8% at 100 GeV, and 4% at 1 TeV, compared to about 40, 12, and 5% obtained when the calorimeters alone are used for jet clustering. All jets are required to satisfy loose requirements on the relative composition of their particle constituents to reject noise in the calorimeter systems or failures in event reconstruction.
Jets are identified as originating from b quarks using the combined secondary vertex algorithm [83]. Control regions in data [84] are used to measure the probability of correctly identifying jets as originating from b quarks (b tagging efficiency), and the probability of misidentifying jets originating from light-flavour partons (u, d, s quarks or gluons) or a charm quark as a b-tagged jet (the light-flavour and charm mistag probabilities). A working point is employed that yields a b tagging efficiency of 65%, and charm and light-flavour mistag probabilities of approximately 12 and 1%, respectively, for jets with p T that is typical of tt events.
An estimator of p miss T is given by the projection on the plane perpendicular to the beams of the negative vector sum of the momenta of all candidate particles in an event [85], as determined by the PF algorithm. Its magnitude is referred to as E miss T .

Event selection
The kinematic selection criteria used to define the signal region, containing a sample of candidate signal events, as well as a number of control regions in data, are described below. The criteria are based on the particle candidates defined by the event reconstruction algorithms described in Sect. 4.

Common preselection criteria
A number of beam-and detector-related effects, such as beam halo, reconstruction failures, spurious detector noise, or event misreconstruction due to detector inefficiencies, can lead to events with anomalous levels of activity. These rare events, which can exhibit large values of E miss T , are rejected with high efficiency by applying a range of dedicated vetoes [85,86].
In order to suppress SM processes with genuine p miss T from neutrinos, events containing an isolated electron or muon that satisfies p T > 10 GeV and |η| < 2.5 are vetoed. Events containing an isolated photon with p T > 25 GeV and |η| < 2.5 are also vetoed, in order to select only multijet final states. Furthermore, events containing a single isolated track satisfying p T > 10 GeV and |η| < 2.5 are vetoed in order to reduce the background contribution from final states containing hadronically decaying tau leptons. Each jet j i considered by this search is required to satisfy p j i T > 40 GeV and |η j i | < 3. The number of jets within this experimental acceptance is labelled henceforth as n jet . The highest p T jet in the event is required to have p j 1 T > 100 GeV and |η j 1 | < 2.5. The second-highest p T jet in the event is used to categorise events, as described in Sect. 5.2. If the jet satisfies p j 2 T > 100 GeV, then this category of events is labelled "symmetric" and targets primarily topologies resulting from pair-produced SUSY particles. If the jet satisfies 40 < p j 2 T < 100 GeV then the event is labelled as "asymmetric," and if there exists no second jet with p j 2 T > 40 GeV, the event is labelled monojet. The asymmetric and monojet topologies target models involving the direct production of stable, weakly interacting, massive particles. The mass scale of the physics processes being probed is characterised by the scalar p T sum of the jets, defined as The magnitude of the vector p T sum of these jets, defined by H miss T |, is used to identify events with significant imbalance in p miss T . Events are vetoed if any jet satisfies p T > 40 GeV and |η| > 3 to ensure that jets reconstructed in the forward regions of the detector do not contribute significantly to H miss T . The dimensionless variable H miss T /E miss T is used to remove events that contain several jets with transverse momenta below the jet p T thresholds but an appreciable vector p T sum so as to contribute significantly to H miss T relative to E miss T . This background is typical of multijet events, which is suppressed by requiring H miss T /E miss T < 1.25. The requirement is imposed as part of the common preselection criteria used to define all control samples to minimise potential systematic biases associated with the simulation modelling for this variable. A high efficiency is maintained for SM or new-physics processes that produce unobserved particles, which are characterised by large values of p miss  Table 1, used to define a sample of all-jet events characterised by high jet activity and appreciable p miss T .

Event categorisation
Events selected by the common preselection criteria are categorised according to n jet , the number of b-tagged jets n b , and H T . Nine categories in n jet are employed: the monojet topology (n jet = 1) and four n jet bins (2, 3, 4, ≥5) for each of the asymmetric and symmetric topologies. Events are also categorised by n b (0, 1, 2, ≥3), where n b is bounded from above by n jet , resulting in 32 categories in terms of both n jet and n b . For each (n jet , n b ) category, events are binned according to H T : four 50 GeV bins at low jet activity in the range 200 < H T < 400 GeV, two 100 GeV bins in the range 400 < H T < 600 GeV, one bin covering the region 600 < H T < 800 GeV, and a final open bin for H T > 800 GeV. These categorisations are summarised in Table 1. The H T binning scheme is adapted independently per (n jet , n b ) category by removing or merging bins to satisfy a threshold on the minimum number of data events in the control regions, which are used to estimate SM backgrounds, provide checks, and validate assumptions within the methods. The lower bounds of the first and final (open) bins in H T are summarised in Table 2. In summary, the search employs a categorisation scheme for events that results in 191 bins, defined in terms of n jet , n b , and H T .

Signal region
For events satisfying the common preselection criteria described above, the multijet background dominates over all other SM backgrounds. Several variables are employed to reduce the multijet contribution to a low level with respect to other SM backgrounds.
The dimensionless kinematic variable α T [28,53], defined in Eq. (2) below, is used to provide discrimination against multijet events that do not contain significant p miss T or that contain large p miss T only because of p T mismeasurements, while retaining sensitivity to new-physics events with significant p miss T . The α T variable depends solely on the transverse component of jet four-momenta and is intrinsically robust against the presence of jet energy mismeasurements in multijet systems. For events containing only two jets, α T is defined as α T = E  T and back-to-back jets ( φ = π ), and in the limit in which the momentum of each jet is large compared with its mass, the value of α T is 0.5. For an imbalance in the E T of back-to-back jets, α T is reduced to a value <0.5, which gives the variable its intrinsic robustness. Values sig-nificantly greater than 0.5 are observed when the two jets are not back-to-back and recoil against p miss T from weakly interacting particles that escape the detector.
The definition of the α T variable can be generalised for events with more than two jets [28]. The mass scale for any process is characterised through the scalar sum of the jet transverse energies, defined as where N jet is the number of jets with E T above a predefined threshold. (The definition of E T should be contrasted with that of H T , the scalar p T sum of the jets.) For events with three or more jets, a pseudo-dijet system is formed by combining the jets in the event into two pseudo-jets. The E T for each of the two pseudo-jets is given by the scalar E T sum of its contributing jets. The combination chosen is the one that minimises E T , defined as the difference between these sums for the two pseudo-jets. This clustering criterion assumes a balancedevent hypothesis, which provides strong separation between SM multijet events and events with genuine p miss T . The α T definition can be generalised to: When jet energies are mismeasured, or there are neutrinos from heavy-flavour quark decays, the magnitudes of H miss T and E T are highly correlated. This correlation is much weaker for R-parity-conserving SUSY events, where each of the two decay chains produces an undetected LSP.
Multijet events populate the region α T < 0.5 and the α T distribution is characterised by a sharp edge at 0.5, beyond which the multijet event yield falls by several orders of magnitude. Multijet events with extremely rare but large stochastic fluctuations in the calorimetric measurements of jet energies can lead to values of α T slightly above 0.5. The edge at 0.5 sharpens with increasing H T for events containing at least three jets, primarily due to a corresponding increase in the average jet energy and consequently a (relative) improvement in the jet energy resolution.
For events containing at least two jets, thresholds on the minimum allowed α T values are applied independent of n jet and n b but dependent on H T , for events that satisfy 200 < H T < 800 GeV. The α T thresholds vary between 0.65 and 0.52 for, respectively, the regions 200 < H T < 250 GeV and 400 < H T < 800 GeV. No requirement on α T is made for the region H T > 800 GeV. The thresholds employed are summarised in Table 1. The α T thresholds are motivated both by the trigger conditions used to record the candidate signal events, described below, and by simulation-based studies and estimates of the multijet background derived from data.
An additional variable is based on the minimum azimuthal separation between a jet and the negative vector p T sum derived from all other jets in the event [28], This variable discriminates between final states with genuine p miss T , e.g. from the leptonic decay of the W boson, and energetic multijet events that have significant p miss T through jet energy mismeasurements or through the production of neutrinos, collinear with the axis of a jet, from semileptonic heavy-flavour decays. Multijet events populate the region Δφ * min < 0.5, with the multijet distribution peaking at a value of zero and falling approximately exponentially over five orders of magnitude to a single-event level at a value of Δφ * min ≈ 0.5, which is close to the distance parameter value of 0.4 used by the anti-k T jet clustering algorithm. Events with a genuine source of p miss T exhibit a long tail in Δφ * min with values as large as π .
A requirement of Δφ * min > 0.5 is sufficient to effectively suppress the multijet background to a low level while maintaining high efficiency for new-physics signatures. The combined rejection power of the α T and Δφ * min requirements for the region 200 < H T < 800 GeV is sufficient to suppress multijet events to the few percent level (and always <10%) with respect to all other SM backgrounds in all H T bins for all event categories of the signal region. For the region   Figure 1 shows the distributions of α T and Δφ * min observed in data for events that satisfy the full set of selection criteria used to define the signal region, summarised in Table 1, as well as the following modifications. The α T and Δφ * min distributions are constructed from events that satisfy n jet ≥ 2, p j 2 T > 100 GeV, and, respectively, H T > 300 or 800 GeV. In the case of Fig. 1 (top), the events with α T values greater than 0.55 must fulfill the full set of signal region criteria, including the Δφ * min > 0.5 requirement, while the events that satisfy α T < 0.55 are subject to the looser set of common preselection criteria defined in Table 1, excluding the H miss T > 130 GeV requirement. Hence, Fig. 1 (top) demonstrates the combined performance of several variables that are employed to suppress multijet events. For both distributions, the events are recorded with a set of inclusive trigger conditions that are independent of the α T and Δφ * min vari-ables. The distributions for the QCD multijet background are determined from simulation while all other SM backgrounds (vector boson production in association with jets, tt, and other residual contributions from rare SM processes) are estimated using a μ + jets data control sample, described in Sect. 6.2. The contribution from multijet events is observed to fall by more than five orders of magnitude for both variables.
The α T and Δφ * min requirements described above, in conjunction with the common preselection requirements H miss T > 130 GeV and H miss T /E miss T < 1.25, provide strong rejection power against the multijet background for events that satisfy n jet ≥ 2. For monojet events, a modification to the Δφ * min variable, which considers soft jets with p T > 25 GeV (Δφ * 25 min ), is utilised. No α T requirement is imposed, and Δφ * 25 min > 0.5 is sufficient to suppress contributions from multijet events to a negligible level. The aforementioned requirements complete the event selection criteria for the signal region.
Finally, Δφ * 25 min is also used as a control variable for events that satisfy n jet ≥ 2 to identify multijet contributions arising from instrumental effects, such as inefficient detector elements or detector noise. The axis of any jet that satisfies Δφ * 25 min < 0.5 is used to identify localised behaviour in the (η, φ) plane, which may be indicative of instrumental defects. No significant anomalies are observed in the sample of candidate signal events following the application of the dedicated vetoes described in Sect. 5.1.
Multiple trigger conditions are employed in combination to record candidate signal events. A set of trigger conditions utilise calculations of both H T and α T to record events with two or more jets. An event is recorded if it satisfies any of the following pairs of (H T [GeV], α T ) thresholds, (200, 0.57), (250, 0.55), (300, 0.53), (350, 0.52), or (400, 0.51), as well as a requirement on the mean value of the two highest p T jets, p T > 90 GeV. These requirements are collectively labelled as the "H T -α T " triggers. In addition, candidate signal events with one or more jets are also recorded if they satisfy the requirements H miss T > 90 GeV and E miss T > 90 GeV. Finally, for events that satisfy H T > 800 GeV, an additional trigger condition, defined by H T > 800 GeV, is employed in addition to the H T -α T trigger requirements to record events characterised by high activity in the calorimeters with high efficiency. The trigger-level jet energies are corrected to account for energy scale and pileup effects. The aforementioned triggers are employed in combination to provide efficiencies at or near 100% for all bins in the signal region.

Using H miss T templates
Following the event selection criteria described above, which provide a sample of candidate signal events with a negligible contribution from multijet events, further discriminating power is required to separate new-physics signatures from the remaining SM backgrounds, which are dominated by the production of tt or W(→ ν) + jets and Z(→ νν) + jets events. As discussed in Sect. 1, the production of coloured SUSY particles with decays to a weakly interacting LSP typically gives rise to a final state with multiple jets and large p miss T . The search therefore exploits the H miss T variable as an additional discriminant between new-physics and SM processes.
The search relies directly on simulation to determine a template for each (n jet , n b , H T ) bin that describes the expected distribution of events as a function of H miss T . These templates are used by the likelihood function as a model for the data, details of which can be found in Sect. 7. The templates are extensively validated against data in multiple control regions, and these studies are used to establish the uncertainty in the simulation modelling of the H miss T variable. The effects of theoretical and experimental uncertainties on the H miss T distributions are also studied. Further details can be found in Sect. 6.2.

Control regions
Four control regions in data are employed to estimate the background contributions from SM processes. The event selection criteria used to define the control regions comprise the common preselection requirements and additional sample-specific requirements, as summarised in Table 1. The first control region comprises a multijet-enriched sample of events, and is defined by the signal region selection criteria and the inverted requirement H miss T /E miss T > 1. 25. The events are recorded with the signal triggers described above, and the sample is used to estimate the multijet background in the signal region. Three additional control regions, defined by inverting one of the photon or lepton vetoes to select samples of γ + jets, μ + jets, or μμ + jets events, are used to estimate the background contributions from SM processes with final states containing genuine p miss T , which are primarily tt, W(→ ν) + jets, and Z(→ νν) + jets.
Additional kinematic requirements are employed to ensure the control samples are enriched in the same SM processes that contribute to background events in the signal region, and are depleted in contributions from multijet production or a wide variety of SUSY models (i.e. so-called signal contamination). The samples are defined, and their events are identically categorised and binned, such that the kinematic properties of events in the control regions and the candidate signal events resemble as closely as possible one another once the photon, muon, or dimuon system is ignored in the calculation of quantities such as H T and H miss T . The selections are summarised in Table 1 and described below.
The γ + jets event sample is defined by the common preselection requirements, but the photon veto is inverted and each event is required to contain a single isolated photon, as defined in Sect. 4, that satisfies p T > 200 GeV and |η| < 1.45 and is well separated from each jet j i in the event according to ΔR(γ , j i ) > 1.0. In addition, events must satisfy H T > 400 GeV, as well as the same H T -dependent α T requirements used to define the signal region. The events are recorded using a single-photon trigger condition and the selection criteria result in a trigger efficiency of 99%.
The μ + jets event sample is defined by the common preselection requirements, but the muon veto is inverted and each event is required to contain a single isolated muon, as defined in Sect. 4, that satisfies p T > 30 GeV and |η| < 2.1 and is well separated from each jet j i in the event according to ΔR(μ, j i ) > 0.5. The transverse mass formed by the muon p T and p miss T system must satisfy 30 < m T < 125 GeV to select a sample of events rich in W bosons, produced promptly or from the decay of top quarks. The μμ + jets sample uses a similar set of selection criteria as the μ + jets sample, but specifically requires two oppositely charged isolated muons that both satisfy p T > 30 GeV and |η| < 2.1 and are well separated from the jets in the event (ΔR(μ 1,2 , j i ) > 0.5). The muons are also required to have a dilepton invariant mass within a ±25 GeV window around the nominal mass of the Z boson [9]. For both the muon and dimuon samples, no requirement is made on α T in order to increase the statistical precision of the predictions from these samples. Both the μ + jets and μμ + jets samples are recorded using a trigger that requires an isolated muon. The selection criteria of the μ + jets and μμ + jets event samples are chosen so that the trigger is maximally efficient, with values of ∼90 and ∼99%, respectively.

Multijet background
The signal region is defined in a manner that suppresses the expected contribution from multijet production to a low level with respect to the total expected background from other SM processes for all signal region bins. This is achieved primarily through the application of very tight requirements on the variables α T and Δφ * min , as described in Sect. 5.3, as well as the requirement H miss T /E miss T < 1.25. In this section, we discuss these requirements further, and present the estimate of the multijet background.
The contamination from multijet events in the signal region is estimated using a multijet-enriched data sideband to the signal region, defined by the (inverted) requirement H miss T /E miss T > 1. 25. The observed counts in data, categorised according to n jet and H T , are corrected to account for contamination from nonmultijet SM processes, and the corrected counts N data (n jet , H T ) are assumed to arise solely from QCD multijet production. The nonmultijet processes, which comprise vector boson and tt production and residual contributions from other SM processes, are estimated using the μ + jets control region, as described in Sect. 6.2.
Independent ratios R QCD (n jet , H T ) of the number of multijet events that satisfy the requirement H miss T /E miss T < 1.25 to the number that fail this requirement are determined from simulation for events categorised according to n jet and H T , and inclusively with respect to n b and H miss T . The product of each ratio R QCD (n jet , H T ) and the corresponding corrected data count N data (n jet , H T ) provides the estimate of the multijet background P(n jet , H T ). The estimates as a function of n jet , H T , n b , and H miss T of the signal region are assumed to factorise as follows: where K n jet ,H T (n b , H miss T ) are multiplier terms that provide the estimated distribution of events as a function of n b and H miss T while preserving the normalisation P(n jet , H T ). The use of simulation to determine R QCD (n jet , H T ) is validated using a multijet-enriched data sideband defined by Δφ * min < 0.5. Each ratio R data (n jet , H T ) is constructed from data counts, corrected to account for contributions from nonmultijet processes, and compared with the corresponding ratio R QCD (n jet , H T ), determined from simulation, through the double ratio R data /R QCD , as shown in Fig. 2. The double ratios are statistically compatible with unity across the full phase space of the signal region, including the bins at high H T , which exhibit the highest statistical precision. In addition to statistical uncertainties as large as ∼100%, a systematic uncertainty of 100% in R QCD is assumed to adequately cover the observed level of agreement for the full signal region phase space.
The distribution of multijet events as a function of n b and H miss T , K n jet ,H T (n b , H miss T ), is assumed to be identical to the distribution expected for the nonmultijet backgrounds. This final assumption is based on studies in simulation and is a valid approximation given the magnitude of this background contribution, as well as the magnitude of the statistical and systematic uncertainties in the ratios R QCD (n jet , H T ), as described above.
6.2 Backgrounds with genuine E miss T Following the suppression of multijet events through the use of the α T and Δφ * min variables, the dominant nonmultijet backgrounds involve SM processes that produce highp T neutrinos in the final state. In events with few jets or few b quark jets, the associated production of W or Z bosons and jets, with the decays W ± → ν ( = e, μ, τ ) or Z → νν, Fig. 2 Validation of the ratio R QCD determined from simulation in bins of n jet and H T [GeV] by comparing with an equivalent ratio R data constructed from data in a multijet-enriched sideband to the signal region. A value of unity is expected for the double ratio R data /R QCD , and the grey shaded band represents the assumed systematic uncertainty of 100% in R QCD The veto of events containing at least one isolated track further suppresses these backgrounds, including those from single-prong τ -lepton decays. At higher jet or b-quark jet multiplicities, single top quark and tt production, followed by semileptonic top quark decay, also become an important source of background. The method to estimate the nonmultijet backgrounds in the signal region relies on the use of a transfer (T ) factor determined from simulation that is constructed per bin (in terms of n jet , n b , and H T ) per control region. Each T factor is defined as the ratio of the expected yields in the same (n jet , n b , H T ) bins of the signal region N SR MC and one of the control regions N CR MC . The T factors are used to extrapolate from the event yields observed in each bin of a data control sample N CR data to provide an estimate for the background, integrated over H miss T , from a particular SM process or processes in the corresponding bin of the signal region N SR data . The superscript SR or CR refers to, respectively, the process or processes being estimated and one of the μ + jets, μμ+jets, and γ +jets control regions, described in Sect. 5.5. The subscript refers to whether the counts are obtained from data, simulation ("MC"), or an estimate ("pred").
The method aims to minimise the effects of simulation mismodelling, as many systematic biases in the simulation are expected to largely cancel in the T factors, given that the events in any given (n jet , n b , H T ) bin of the control regions closely mirror those in the corresponding bin in the signal region in terms of the event energy scale, topology, and kinematics. In short, minimal extrapolations are made. Uncertainties in the T factors are determined from data, as described below.
Three independent estimates of the irreducible background of Z → νν + jets events are determined from the γ + jets, μμ + jets, and μ + jets data control samples. The γ + jets and Z → μμ + jets processes have similar kinematic properties when the photon or muons are ignored [87], albeit different acceptances. In addition, the γ + jets process has a larger production cross section than Z → νν + jets events. The μ + jets data sample is used to provide an estimate for both the Z → νν + jets background, as well as the other dominant SM processes, tt and W boson production (labelled collectively as W/tt). Residual contributions from all other SM relevant processes, such as single top quark, diboson, and Drell-Yan production, are also included as part of the W/tt estimate from the μ + jets sample. The definition of the various T factors used in the search are given below: The likelihood function, described in Sect. 7, encodes the estimate via the T factors of the W/tt background, as well as the three independent estimates of the Z → νν background, which are considered simultaneously.
Several sources of uncertainty in the T factors are evaluated. The most relevant effects are discussed below, and generally fall into one of two categories. The first category concerns uncertainties in the "scale factor" corrections applied to simulation, which are determined using inclusive data samples that are defined by loose selection criteria, to account for the mismodelling of theoretical and experimental parameters. The second category concerns "closure tests" in data that probe various aspects of the accuracy of the simulation to model correctly the T factors in the phase space of this search.
The uncertainties in the T factors are studied for variations in scale factors related to the jet energy scale (that result in uncertainties in the T factors as large as ∼15%), the efficiency and misidentification probability of b quark jets (up to 5%), and the efficiency to identify well-reconstructed, isolated leptons (up to ∼5%). A 5% uncertainty in the total inelastic cross section, σ in = 69.0 ± 3.5 mb [88], is assumed and propagated through to the reweighting procedure to account for differences between the simulated measured pileup, which results in changes of up to ∼10%. The modelling of the transverse momentum of top quarks ( p t T ) is evaluated by comparing the simulated and measured p T spectra of reconstructed top quarks in tt events [89]. Simulated events are reweighted according to scale factors that decrease from a value of ∼1.2 to ∼0.7, with uncertainties of ∼10-20%, within the range p t T < 400 GeV. The systematic uncertainties in T W/tt μ+jets arising from variations in the p t T scale factors are typically small ( 5%), due to the comparable phase space probed by the signal and control regions, while larger uncertainties ( 20%) in T Z→νν μ+jets are observed due to the potential for significant contamination from tt when using W(→ ν) + jets to predict Z(→ νν) + jets.
The aforementioned systematic uncertainties, resulting from variations in scale factors, are summarised in Table 3, along with representative magnitudes. These sources of uncertainty are each assumed to originate from a unique underlying source and so the effect of each source is varied assuming a fully correlated behaviour across the full phase space of the signal and control regions.
The second category of uncertainty is determined from sets of closure tests based on data control samples [31]. Each set uses the observed event counts in up to eight bins in H T for each of the nine n jet event categories in one of the three independent data control regions. These counts are used with the corresponding T factors, determined from simulation, to obtain a prediction N pred (n jet , H T ) of the observed yields N obs (n jet , H T ) in another control sample (or, in one case, n b event category).
Each set of tests targets a specific (potential) source of bias in the simulation modelling that may introduce an n jet -or H T -dependent source of systematic bias in the T factors [31]. Several sets of tests are performed. The Z/γ ratio determined from simulation is tested against the same ratio measured using Z(→ μμ) + jets events and the γ + jets sample. The W/Z ratio is also probed using the μ+jets and μμ+jets samples, which directly tests the simulation modelling of vector boson production, as well as the modelling of tt contamination in the μ+jets sample. A further set probes the modelling of the relative composition between W(→ ν) + jets and tt events using μ + jets events containing exactly zero or one more b-tagged jets, which represents a larger extrapolation in relative composition than used in the search. The effects of W polarisation are probed by using μ + jets events with a positively charged muon to predict those containing a negatively charged muon. Finally, the accuracy of the modelling of the efficiencies of the α T and Δφ * min requirements are estimated using the μ + jets sample.
For each set of tests, the level of closure, (N obs − N pred )/N obs , which considers only statistical uncertainties, is inspected to ensure no statistically significant biases are observed as a function of the nine n jet categories or the eight H T bins. In the absence of such a bias, the level of closure is recomputed by integrating over either all monojet and asymmetric n jet categories, or the symmetric n jet categories. The level of closure and its statistical uncertainty are combined in quadrature to determine additional contributions to the uncertainties in the T factors. These uncertainties are considered to be fully correlated between the monojet and asymmetric n jet categories or the symmetric n jet categories, and fully uncorrelated between these two regions in n jet and H T bins. If the closure tests use the μμ + jets sample, the level of closure is determined by additionally integrating over pairs of adjacent H T bins. These uncertainties, derived from the closure tests in data, are summarised in Table 3, along with representative magnitudes. These uncertainties are the dominant contribution to the total uncertainty in the T factors, due to the limited number of events in the data control regions.
As introduced in Sect. 5.4, templates are derived from simulation to predict the H miss T distributions of the background. The uncertainties in the T factors are used to constrain the normalisation of the H miss T templates. The uncertainties in the H miss T shape are discussed below. The accuracy to which the simulation describes the H miss T distributions is evaluated with respect to data in each (n jet , n b , H T ) bin in each of the μ + jets, μμ + jets, and γ + jets data control regions. The level of agreement between data and simulation, defined in terms of the ratio of observed and expected counts (from simulation) as a function of H miss T , is parameterised using an orthogonal first-order polynomial, f (x) = p 0 + p 1 (x − x), and described by two uncorrelated parameters, p 0 and p 1 . A binned likelihood fit is performed in each (n jet , n b , H T ) bin of each control region, and the best fit value p 1 and its uncertainty is used to determine the presence of biases dependent on H miss T . The pull of p 1 from a value of zero is defined as the best fit value over its standard deviation, considering only statistical uncertainties associated with the finite size of the data and simulated samples.
The lower bound of the final (open) bin in H miss T is not more than 800 GeV and is bounded from above by the upper bound of the H T bin in question. The lower bound of the final H miss T bin is merged with lower bins if fewer than ten events in the data control regions are observed. If a bin in (n jet , n b , H T ) contains fewer than ten events, the H miss T template is not used and the background estimates are determined inclu-sively with respect to H miss T . The merging of bins is typically only relevant for event categories that satisfy n b ≥ 2.
The presence of systematic biases is evaluated at a statistical level by considering the distribution of pulls obtained from each control region, which are consistent with statistical fluctuations, with no indication of trends across the full phase space of each control region. The p-values obtained from the fits are uniformly distributed.
The uncertainty in the H miss T modelling is extracted under the hypothesis of no bias. This is done using the maximum likelihood (ML) values of the fit parameters to determine the statistical precision to which this hypothesis can be confirmed. The quadrature sum of the ML value and its uncertainty for p 1 from each fit is used to define alternative templates that represent ±1σ variations to the nominal H miss T template. These alternative templates are encoded in the likelihood function, as described in Sect. 7. The observed variations are compatible with the expected values obtained from studies relying only on simulated event samples. The uncertainties in the final H miss T bin of the templates depend on the event category and H T bin, and are typically found to be in the range ∼10-100%.
The effect on the H miss T templates is determined under ±1σ variations in the jet energy scale, the efficiency and misidentification probability of b-quark jets, the efficiency to identify well-reconstructed, isolated leptons, the pileup reweighting, and the modelling of the top quark p T . These effects are easily covered by the uncertainties determined from data, as described above, across the full phase space of the control regions, which mirror closely that of the signal region.

Results
A model of the observations in all data samples, described by a likelihood function, is used to obtain a prediction of the SM backgrounds and to test for the presence of newphysics signals if the signal region is included in the ML fit. The observation in each bin defined by the n jet , n b , H T , and H miss T variables is modelled as a Poisson-distributed variable around the SM expectation and a potential signal contribution (assumed to be zero in the following discussion), where the SM expectation is the sum over the estimated contributions from all background processes according to the methods described in Sect. 6.
The nonmultijet backgrounds are related to the expected yields in the μ + jets, μμ + jets, and γ + jets control samples via the transfer factors derived from simulation, as described in Sect. 5.5. Estimates of the contribution from multijet events in the signal region are determined according to the method described in Sect. 6.1, and are included in the likelihood function.
The data are inspected to ascertain whether they are well described by the null (SM-only) hypothesis. This is done by considering the "pre-fit" SM background estimates, which are determined from observed data counts in the control regions only. The pre-fit result of this search is summarised in Figs. 3, 4 and 5 for, respectively, the monojet, asymmetric, and symmetric topologies. The figures also show the significance of deviations observed in data with respect to the pre-fit SM expectations expressed in terms of the total uncertainty in the SM expectations ("pull"). The data are well described by the background-only hypothesis. Figures 3, 4 and 5 also summarise the pulls from the post-fit result, which is based on a ML fit to observations in the signal region as well as the control regions.
A quantitative statement on the degree of compatibility between the observed yields and the SM expectations under the background-only hypothesis is obtained from a goodnessof-fit test based on a log likelihood ratio. The alternative hypothesis is defined by a "saturated" model [91], for which the background expectation is set equal to the observed number of events, and provides a reference for the largest value that the likelihood can take for any model for the given data set. Hence, this reference can be used as a reasonable normalisation for the maximum value observed for a more constraining model. A p-value of 0.20 is observed for the fit over the full signal region, and p-values in the range 0.04-1.00, Prefit Postfit The covariance and correlation matrices for the pre-fit SM expectations in all bins of the signal region, defined by n jet , n b , H T , and integrated over H miss T , are determined from 500 pseudo-experiments by sampling the pre-fit nuisance parameters under the background-only hypothesis. The SM expectations for different n jet and n b categories exhibit a nonnegligible level of covariance within the same H T bin, primarily as a result of the systematic uncertainties evaluated from closure tests, described in Sect. 6.2, that integrate yields over n jet and n b . Bins adjacent and next-to-adjacent in n jet and/or n b can have correlation coefficients in the range 0.2-0.4, and, infrequently, as large as ∼0.5. Otherwise, the correlation coefficients are <0.2. Anticorrelation coefficients are typically not larger than ∼0.2. Figure 6 shows the event yields observed in data and prefit SM expectations with their associated uncertainties as a function of H miss T for two categories of candidate signal events, which provide good sensitivity to models with highmass gluinos. For illustration only, the expected counts from benchmark signal models that assume the pair production and decay of gluinos, described further in Sect. 8, are also shown.

Specification for simplified models
The results of the search are used to constrain simplified SUSY models [92][93][94]. Each model assumes the pair production of gluinos or squarks and their subsequent prompt decays to SM particles and the LSP with a 100% branching fraction (unless indicated otherwise). The gluino decays contain intermediate on-shell SUSY particle states (such as the top squark or the chargino) for a subset of the models. All other SUSY particles are assumed to be too heavy (m g /m q = 10 TeV) to be produced directly. Three-body decays of gluinos are assumed to occur via off-shell squarks of light or heavy flavour. Off-shell decays are processed by pythia in a single three-or four-body step, without taking into account the width or polarisation of the parent: this is true for the top-squark four-body decay ( t → bff χ 0 1 ), as well as the three-body decay of the chargino ( χ ± 1 → ff χ 0 1 ), where f and f are fermions produced in the decay of an intermediate off-shell W boson.
Under the signal+background hypothesis, and in the presence of a nonzero signal contribution, a modified frequentist approach is used to determine upper limits at the 95% confidence level (CL) on the cross section, σ UL , to produce pairs of SUSY particles as a function of the parent SUSY particle and the LSP masses. The limits can be expressed in terms of the signal strength parameter, μ, which is determined relative to the theoretical cross section that is calculated at NLO + NLL accuracy. An Asimov data set [95] is used to determine the expected upper limit on the allowed cross section for a given model. The potential contributions from a new-physics signal to each of the signal and control regions are considered, even though the only significant contribution occurs in the signal region and not in the control regions (i.e. signal contamination). The approach is based on the one-sided (so called LHC-style) profile likelihood ratio as the test statistic [96] and the CL s criterion [97,98]. Asymptotic formulae [95] are utilised to approximate the distributions of the test statistics under the SM background-only and signal+background hypotheses.

Acceptances and uncertainties
The experimental acceptance times efficiency (A ε) and its uncertainty are evaluated independently for each model class as a function of (m SUSY , m LSP ). Table 5 summarises A ε for a number of benchmark models for which the search yields an expected exclusion (μ 1). For each topology, The lower panels indicate the significance of deviations (pulls) observed in data with respect to both the pre-fit SM expectations, expressed in terms of the total uncertainty in the SM expectations. The pulls cannot be considered independently due to inter-bin correlations typically two different pairs of parent SUSY particle and LSP masses (m SUSY , m LSP ) are chosen that are characterised by a large and a small (i.e. compressed) difference in parent SUSY particle and LSP masses. The four most sensitive event categories, defined in terms of n jet , are used to determine σ UL . The categories used per benchmark model are listed in Table 5, along with A ε determined for these four categories.
The effects of several sources of uncertainty on A ε, as well as the potential for event migration between bins of the signal region, are considered. The potential effect of each source of uncertainty is assessed by including in the likelihood function the alternative normalisations and shapes for the H miss T templates with respect to the nominal versions. The nominal and alternative templates are obtained from simulated event samples for the signal models, and the alternative templates effectively propagate the various input uncertainties to determine their effects on the n jet , n b , H T , and H miss T distributions for the signal model.
In addition to the uncertainty in the integrated luminosity of 2.7% [99], the following sources of uncertainty are dominant: the statistical uncertainty arising from the finite size of simulated signal samples, the modelling of ISR, the jet energy corrections (JEC) evaluated in simulation, and the modelling of scale factors applied to simulated event samples that correct for differences in the efficiency and misidentification probability of b-quark jets (SF b-tag ). The magnitude of each contribution depends on the model and the masses of the parent SUSY particle and LSP.
The A ε for models with small mass splittings (e.g. m q − m χ 0 1 m t ) is due largely to ISR, the modelling of which is evaluated by comparing the simulated and measured p T spectra of the system recoiling against the ISR jets in tt events, using the technique described in Ref. [100]. The uncertainty can be as large as ∼30%, and is the dominant systematic uncertainty for systems with a compressed mass spectrum. Uncertainties in the jet energy scale, as large as ∼40%, can also be dominant for models characterised by high jet multiplicities in the final state. The uncertainties in SF b-tag can be as large as ∼25%. Table 5 summarises these dominant contributions to the uncertainty in A ε for a range of benchmark models. Characteristic values for each model are expressed in terms of a range that is representative of the values across all bins of the signal region. The upper bound for each range may be subject to moderate statistical fluctuations.
Further uncertainties with subdominant contributions are considered on a similar footing. The uncertainties in the efficiency of identifying well-reconstructed, isolated leptons are considered, with a typical magnitude of ∼5% and treated as anticorrelated between the signal and control regions. The uncertainty of 5% in σ in is propagated through to the Gluino-mediated and direct production of light-flavour squarks T2qq_1fold pp → q q q → q χ 0 1 m q( q = uL) m uL Gluino-mediated production of off-shell third-generation squarks Natural gluino-mediated production of on-shell top squarks Direct production of on-shell third-generation squarks reweighting procedure to account for differences between the simulated and measured pileup. Finally, uncertainties in the simulation modelling of the efficiencies of the trigger strategy employed by the search are typically <10%. The choice of PDF set, or variations therein, predominantly affects A ε through changes in the p T spectrum of the system recoil, which is covered by the ISR uncertainty, hence no additional uncertainty is adopted. Uncertainties in A ε due to variations in the renormalisation and factorisation scales are determined to be ∼5%. In both cases, contributions to the uncertainty in the theoretical production cross section are considered.

Cross section and mass exclusions
Limits for each of the aforementioned benchmark models are summarised in Table 5, expressed in terms of μ. All benchmark models are expected to be excluded. The observed limits fluctuate around the expected μ values, with some models exhibiting a moderately weaker than expected limit. Figures 7 and 8 summarise the disfavoured regions of the mass parameter space for the fourteen classes of simplified models. These regions are derived by comparing the upper limits on the measured fiducial cross section, corrected for the experimental A ε, with the theoretical cross sections calculated at NLO+NLL accuracy in α s [33]. The former cross section value is determined as a function of m g or m q and m χ 0 1 , while the latter has a dependence solely on m g or m q . The exclusion of models is evaluated using observed data counts in the signal region (solid contours) and also expected counts based on an Asimov data set (dashed contours). Figure 7 (upper) shows the excluded mass parameter space for models that assume the gluino-mediated or direct production of light-flavour squarks. The excluded region extends to higher masses for the gluino-mediated production of lightflavour squarks (T1qqqq), with respect to direct pair production when assuming an eightfold degeneracy in mass (T2qq_8fold), due to a combination of a higher gluino pair production cross section and a final state characterised by higher jet multiplicities, which are exploited to provide better signal-to-background separation. The excluded mass region is significantly reduced when assuming only a single light squark (T2qq_1fold), with limits weakening due to the lower production cross section, compounded by the reduced signal-to-background ratios achieved in the core of distributions in the discriminating variables. Figure 7 (lower) shows the exclusion contours for models that assume the gluino-mediated pair production of off-shell third-generation squarks. For the topologies T1tttt and T1bbbb, each gluino is assumed to undergo a three-body decay via, respectively, an off-shell top or bottom squark to a quark-antiquark pair of the same flavour and the χ 0 1 . In the case of T1ttbb, each gluino is assumed to undergo a threebody decay to an on-shell chargino, χ ± 1 , a bottom quark, and The chargino decays promptly to the χ 0 1 and an off-shell W boson. The excluded mass regions differ significantly for these topologies, primarily due to the different number of (onshell) W bosons in their final states, resulting in the highest A ε for T1bbbb and lowest for T1tttt. Further, A ε has a strong dependence on jet multiplicity, which is highest for T1tttt, due to the Δφ * min variable. An additional feature for T1ttbb is the weakening of the mass limit at low values of m χ 0 1 , when m χ ± 1 = m χ 0 1 + 5 GeV m t . In this scenario, the χ ± 1 (and hence χ 0 1 ) is not highly Lorentz boosted relative to the top quark resulting from the three-body decay of the gluino. Hence, two χ 0 1 SUSY particles do not carry away significant p miss T , which is instead realised through W boson decays to neutrinos and "lost" leptons or τ leptons that decay to neutrinos and hadrons. The observed mass limits for these topologies are up to ∼2 standard deviations weaker than the expected limits. These differences are due to upward fluctuations in data for two contiguous bins that satisfy the requirements n jet ≥ 5, n b ≥ 2, and H T > 800 GeV. This region has the highest sensitivity to models involving gluino production and decays to third-generation quarks (via on-or offshell squarks). The observed counts are consistent with statistical fluctuations and the events do not exhibit anomalous nonphysical behaviour. The events are distributed in H miss T consistent with expectation, hence models characterised by high values of H miss T , such as T1bbbb with m g m χ 0 1 or m g ≈ m χ 0 1 , are less compatible with the data counts in this high-n jet , n b , and H T region. Figure 8 (upper) shows exclusion contours for models that assume gluino pair production, with each gluino decaying to a top quark and an on-shell top squark, the latter of which decays to SM particles and the LSP. As discussed earlier, these models can be considered as representations of a natural solution to the little hierarchy problem. Two different scenarios are considered for the decay of the top squarks.
Expected Observed Fig. 7 Observed and expected mass exclusions at 95% CL (indicated, respectively, by solid and dashed contours) for various classes of simplified models. (Top) Gluino-mediated or direct pair production of light-flavour squarks. The two scenarios involve, respectively, the decay g → qq χ 0 1 (T1qqqq) and q → q χ 0 1 , and the latter involves two assumptions on the mass degeneracy of the squarks (T2qq_8fold and T2qq_1fold). (Bottom) Three scenarios involving the gluino-mediated pair production of off-shell third-generation squarks: g → bb χ 0 1 (T1bbbb), g → t t * → tt χ 0 1 (T1tttt), and g → tb χ ± 1 → tbW * χ 0 1 (T1ttbb) The T5tttt_DM175 models assume a two-body decay to a top quark and the χ 0 1 , with the top squark mass defined relative to the   = 20 GeV (T5ttcc). Also shown, for comparison, is T1tttt. (Bottom) Six scenarios involving the direct pair production of third-generation squarks. The first scenario involves the pair production of bottom squarks, b → b χ 0 1 (T2bb). Two scenarios involve the decay of top squark pairs as follows: t → t χ 0 and branching fractions 50/50% (T2tb), or t → t χ 0 1 (T2tt). The final three scenarios consider top squark decays under the assumption 10 < m t − m χ 0 1 < 80 GeV: t → c χ 0 1 (T2cc), t → bW * χ 0 1 (T2tt_degen), and t → c χ 0 1 or t → bW * χ 0 1 with branching fractions 50/50% (T2tt_mixed). The grey shaded region denotes T2tt models that are not considered for interpretation χ 0 1 . The T2tt models assume top squarks are pair produced and each is assumed to undergo a two-or three-body decay to, respectively, a top quark and the   in this low-m t region due to the high levels of signal contamination found in the μ + jets control region for models that resemble the tt background in terms of their topological and kinematic properties. The T2tb models also assume the pair production of top squarks, with each undergoing a two-body decay to either a top quark and the χ 0 1 , or a bottom quark and the χ ± 1 , with equal branching fractions of 50%. As for the T1ttbb models, the chargino mass is defined relative to the neutralino mass via the expression m χ ± 1 −m χ 0 1 = 5 GeV, and the chargino decays promptly to the χ 0 1 and an off-shell W boson. The excluded mass regions differ significantly for the T2bb, T2tb, and T2tt topologies, in an analogous way to the T1bbbb, T1ttbb, and T1tttt models described above. The difference in the mass exclusions is due primarily to the different number of (on-shell) W bosons in the final states, which affects A ε through the presence of leptons from the decay of the W boson. An additional feature for T2tb is the weakening of the mass limit at low values of m χ 0 1 , when m χ ± 1 = m χ 0 1 + 5 GeV m t . Moderately weaker than expected mass limits are observed for all models involving two-body decays, which is traced to mild upward fluctuations in data for events satisfying n jet = 2, n b = 2, and 350 < H T < 500 GeV. Figure 8 (lower) also shows exclusion contours for models that assume the pair production of top squarks but a nearmass-degenerate system that satisfies 10 GeV < m t −m χ 0 1 < m W . Two decays of the top squark are considered. The T2cc and T2tt_degen models assume two-and fourbody decays of the top squark to, respectively, a charm quark and the χ 0 1 , or to bff χ 0 1 , where f andf are fermions produced in the decay of an intermediate off-shell W boson. A third class of models, T2tt_mixed, assumes both these decay modes with an equal branching fraction of 50%. For T2cc, the excluded mass region is relatively stable as a function of the mass splitting m = m t −m χ 0 1 , with t masses excluded up to 400 GeV. For T2tt_degen, the excluded mass region is strongly dependent on m, weakening considerably for increasing values of m due to the increased momentum phase space available to leptons produced in the four-body decay. The T2tt_mixed models exhibit an intermediate behaviour. Mass limits for all three model classes converge for the smallest mass splitting considered, m = 10 GeV, when the SM particles from the t decay are extremely soft and outside the experimental acceptance. An approximately contiguous mass exclusion limit is observed across the transition from the T2tt_degen four-body to the T2tt three-body decay of the t, as the top quark moves on-shell. The excluded mass region weakens further as m → m t . Table 6 summarises the strongest expected and observed mass limits for each class of simplified model.

Summary
An inclusive search for new-physics phenomena is reported, based on data from pp collisions at √ s = 13 TeV. The data are recorded with the CMS detector and correspond to an integrated luminosity of 2.3 ± 0.1 fb −1 . The final states analysed contain one or more jets with large transverse momenta and a significant imbalance of transverse momentum, as expected from the production of massive coloured SUSY particles, each decaying to SM particles and the lightest stable, weaklyinteracting, SUSY particle.
The sums of the standard model backgrounds are estimated from a simultaneous binned likelihood fit to the observed yields for samples of events categorised according to the number of reconstructed jets, the number of jets identified as originating from b quarks, and the scalar and the magnitude of the vector sums of the transverse momenta of jets. In addition to the signal region, μ + jets, μμ + jets, and γ + jets control regions are included in the likelihood fit. The observed yields are found to be in agreement with the expected contributions from standard model processes. The search result is interpreted in the mass parameter space of fourteen simplified SUSY models, which cover scenarios that involve the gluino-mediated or direct production of light-or heavy-flavour squarks, intermediate SUSY particle states, as well as natural and nearly mass-degenerate spectra.
The increase in the centre-of-mass energy of the LHC, from 8 to 13 TeV, provides a significant gain in sensitivity to heavy particle states such as gluinos. In the case of pair-produced gluinos, each decaying via an off-shell b squark to the b quark and the LSP, models with masses up to ∼1.6 and ∼1.0 TeV are excluded for, respectively, the gluino and LSP. These limits improve on those obtained at √ s = 8 TeV by, respectively, ∼250 and ∼300 GeV. In the case of direct pair production, models with masses up to ∼800 and ∼350 GeV are excluded for, respectively, the b squark and LSP. These mass limits are sensitive to the assumptions on the squark flavour and the presence of intermediate states such as charginos.
Finally, a comprehensive study of nearly mass-degenerate models involving top squark pair production is performed. The two decay modes of the top squark are the loop-induced two-body decay to the neutralino and one c quark, and the four-body decay to the neutralino, one b quark, and an offshell W boson. A third scenario is considered in which the two modes are simultaneously open, each with a branching fraction of 50%. Masses of the top squark and LSP up to, respectively, 400 and 360 GeV are excluded, depending on the decay modes considered.
In conclusion, the analysis provides sensitivity across a large region of the natural SUSY parameter space, as characterised by interpretations with several simplified models. In particular, these studies improve on existing limits for nearly mass-degenerate models involving the production of pairs of top squarks. 2012/07/E/ST2/01406; the Thalis and Aristeia programmes cofinanced by EU-ESF and the Greek NSRF; the National Priorities Research Program by Qatar National Research Fund; the Programa Clarín-COFUND del Principado de Asturias; the Rachadapisek Sompot Fund for Postdoctoral Fellowship, Chulalongkorn University and the Chulalongkorn Academic into Its 2nd Century Project Advancement Project (Thailand); and the Welch Foundation, contract C-1845.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .