1 Introduction

Supersymmetry (SUSY) [16] is a generalization of space-time symmetries that predicts new bosonic partners for the fermions and new fermionic partners for the bosons of the Standard Model (SM). If R-parity is conserved [7], SUSY particles (called sparticles) are produced in pairs and the lightest supersymmetric particle (LSP) is stable and represents a possible dark-matter candidate. The scalar partners of the left- and right-handed quarks, the squarks \(\tilde{q}_{\mathrm {L}} \) and \(\tilde{q}_{\mathrm {R}} \), mix to form two mass eigenstates \(\tilde{q}_1\) and \(\tilde{q}_2\) ordered by increasing mass. Superpartners of the charged and neutral electroweak and Higgs bosons also mix to produce charginos (\(\tilde{\chi }^\pm \)) and neutralinos (\(\tilde{\chi }^0 \)). Squarks and the fermionic partners of the gluons, the gluinos (\(\tilde{g} \)), could be produced in strong-interaction processes at the Large Hadron Collider (LHC) [8] and decay via cascades ending with the stable LSP, which escapes the detector unseen, producing substantial missing transverse momentum (\(\varvec{E}\mathrm {^{miss}_T}\)).

The production of gluinos and squarks is the primary target for early supersymmetry searches in proton–proton (pp) collisions at a centre-of-mass energy of 13 \({\mathrm{TeV}}\) at the LHC because of the large expected cross-sections predicted for the production of supersymmetric particles which participate to the strong interaction. This document presents a search for these particles in final states containing only hadronic jets and large missing transverse momentum. Interest in this final state is motivated by the large number of R-parity-conserving models [9, 10] in which squarks (including anti-squarks) and gluinos can be produced in pairs (\(\tilde{g} \tilde{g} \), \(\tilde{q} \tilde{q} \), \(\tilde{q} \tilde{g} \)) and can decay through \(\tilde{q} \rightarrow q \tilde{\chi }^0_1 \) and \(\tilde{g} \rightarrow q \bar{q} \tilde{\chi }^0_1 \) to the lightest neutralino, \(\tilde{\chi }^0_1 \), assumed to be the LSP. Additional decay modes can include the production of charginos via \(\tilde{q} \rightarrow q\tilde{\chi }^\pm \) (where \(\tilde{q} \) and q are of different flavour) and \(\tilde{g} \rightarrow q\bar{q}\tilde{\chi }^\pm \). Subsequent chargino decay to \(W^{\pm }\tilde{\chi }^0_1 \) can lead to still larger multiplicities of jets. The analysis presented here adopts the same analysis strategy as the previous ATLAS search designed for the analysis of the 7 \({\mathrm{TeV}}\) and 8 \({\mathrm{TeV}}\) data collected during Run 1 of the LHC, described in Refs. [1115]. The CMS Collaboration has set limits on similar models in Refs. [1621].

In this search, events with reconstructed electrons or muons are rejected to reduce the background from events with neutrinos (\(W \rightarrow e\nu ,\mu \nu \)) and to avoid any overlap with a complementary ATLAS search in final states with one lepton, jets and missing transverse momentum [22]. The selection criteria are optimized in the \((m_{\tilde{g}}, m_{\tilde{\chi }^0_1})\) and \((m_{\tilde{q}}, m_{\tilde{\chi }^0_1})\) planes, (where \(m_{\tilde{g}}\), \(m_{\tilde{q}}\) and \(m_{\tilde{\chi }^0_1}\) are the gluino, squark and the LSP masses, respectively) for simplified models [2325] in which all other supersymmetric particles are assigned masses beyond the reach of the LHC. Although interpreted in terms of SUSY models, the results of this analysis could also constrain any model of new physics that predicts the production of jets in association with missing transverse momentum.

2 The ATLAS detector and data samples

The ATLAS detector [26] is a multi-purpose detector with a forward-backward symmetric cylindrical geometry and nearly 4\(\pi \) coverage in solid angle.Footnote 1 The inner tracking detector (ID) consists of pixel and silicon microstrip detectors covering the pseudorapidity region \(|\eta |<2.5\), surrounded by a transition radiation tracker which improves electron identification over the region \(|\eta |<2.0\). The innermost pixel layer, the insertable B-layer [27], was added between Run 1 and Run 2 of the LHC, at a radius of 33 mm around a new, narrower and thinner, beam pipe. The ID is surrounded by a thin superconducting solenoid providing an axial 2 T magnetic field and by a fine-granularity lead/liquid-argon (LAr) electromagnetic calorimeter covering \(|\eta |<3.2\). A steel/scintillator-tile calorimeter provides hadronic coverage in the central pseudorapidity range (\(|\eta |<1.7\)). The endcap and forward regions (\(1.5<|\eta |<4.9\)) of the hadronic calorimeter are made of LAr active layers with either copper or tungsten as the absorber material. The muon spectrometer with an air-core toroid magnet system surrounds the calorimeters. Three layers of high-precision tracking chambers provide coverage in the range \(|\eta |<2.7\), while dedicated chambers allow triggering in the region \(|\eta |<2.4\).

The ATLAS trigger system [28] consists of two levels; the first level is a hardware-based system, while the second is a software-based system called the High-Level Trigger. The events used in this search were selected using a trigger logic that accepts events with a missing transverse momentum above 70 GeV, calculated using a sum over calorimeter cells. The trigger is 100 % efficient for the event selections considered in this analysis. Auxiliary data samples used to estimate the yields of background events were selected using triggers requiring at least one isolated electron (\(p_{\text {T}} >24~{\mathrm{GeV}}\)), muon (\(p_{\text {T}} >20~{\mathrm{GeV}}\)) or photon (\(p_{\text {T}} >120~{\mathrm{GeV}}\)). To increase the efficiency at high momenta, additional single-electron and single-muon triggers that do not require any isolation were included with thresholds of \(p_{\text {T}} = 60~{\mathrm{GeV}}\) and \(p_{\text {T}} = 50~{\mathrm{GeV}}\), respectively.

The dataset used in this analysis was collected in 2015 with the LHC colliding proton beams at a centre-of-mass energy of 13 \({\mathrm{TeV}}\), with 25 ns bunch spacing. The peak delivered instantaneous luminosity was \(L = 5.2 \times 10^{33} ~\mathrm{cm^{-2} s^{-1}}\) and the mean number of additional pp interactions per bunch crossing in the dataset was \(\left\langle \mu\right\rangle \) = 14. Application of beam, detector and data-quality criteria resulted in a total integrated luminosity of 3.2 \(\mathrm{fb}^{-1}\). The uncertainty in the integrated luminosity is ±5 %. It is derived, following a methodology similar to that detailed in Ref. [29], from a preliminary calibration of the luminosity scale using a pair of xy beam-separation scans performed in August 2015.

3 Monte Carlo simulated samples

Simulated Monte Carlo (MC) data samples are used to optimize the selections, estimate backgrounds and assess the sensitivity to specific SUSY signal models.

SUSY signals are described in this paper by simplified models. They are defined by an effective Lagrangian describing the interactions of a small number of new particles, typically assuming one production process and one decay channel with a 100 % branching fraction. Signal samples used to describe squark- and gluino-pair production, followed by the directFootnote 2 decays of squarks (\(\tilde{q} \rightarrow q\tilde{\chi }^0_1 \)) and direct (\(\tilde{g} \rightarrow q\bar{q} \tilde{\chi }^0_1 \)) or one-stepFootnote 3 (\(\tilde{g} \rightarrow q\bar{q}'W \tilde{\chi }^0_1 \)) decays of gluinos as shown in Fig. 1, are generated with up to two extra partons in the matrix element using MG5_aMC@NLO event generator [30] interfaced to Pythia 8.186 [31]. The CKKW-L merging scheme [32] is applied with a scale parameter that is set to a quarter of the mass of the gluino for \(\tilde{g} \tilde{g} \) production or of the squark for \(\tilde{q} \tilde{q} \) production. The A14 [33] set of tuned parameters (tune) is used for underlying event together with the NNPDF2.3LO [34] parton distribution function (PDF) set. The EvtGen v1.2.0 program [35] is used to describe the properties of the b- and c- hadron decays in the signal samples and the background samples except those produced with Sherpa [36]. The signal cross-sections are calculated at next-to-leading order (NLO) in the strong coupling constant, adding the resummation of soft gluon emission at next-to-leading-logarithmic accuracy (NLO+NLL) [3741]. The nominal cross-section is taken from an envelope of cross-section predictions using different PDF sets and factorization and renormalization scales, as described in Ref. [42], considering only light-flavour quarks (u, d, s, c). Cross-sections are evaluated assuming masses of 450 \({\mathrm{TeV}}\) for the light-flavour squarks in case of gluino- or gluinos in case of squark-pair production. The free parameters are \(m_{\tilde{\chi }^0_1}\) and \(m_{\tilde{q}}\) (\(m_{\tilde{g}}\)) for gluino-pair (squark-pair) production models.

Fig. 1
figure 1

The decay topologies of a squark-pair production and b, c gluino-pair production, in the simplified models with direct decays of squarks and direct or one-step decays of gluinos

A summary of the SM background processes together with the MC generators, cross-section calculation orders in \(\alpha _\mathrm{s}\), PDFs, parton shower and tunes used is given in Table 1.

Table 1 The Standard Model background Monte Carlo simulation samples used in this paper. The generators, the order in \(\alpha _\mathrm{s}\) of cross-section calculations used for yield normalization, PDF sets, parton showers and tunes used for the underlying event are shown

The production of \(\gamma \), W or Z bosons in association with jets [43] is simulated using the Sherpa 2.1.1 generator. For W or Z bosons, the matrix elements are calculated for up to two partons at NLO and up to two additional partons at leading order (LO) using the Comix [44] and OpenLoops [45] matrix-element generators, and merged with the Sherpa parton shower [46] using the ME+PS@NLO prescription [47]. Events containing a photon in association with jets are generated requiring a photon transverse momentum above 35 \({\mathrm{GeV}}\). For these events, matrix elements are calculated at LO with up to three or four partons depending on the \(p_{\text {T}} \) of the photon, and merged with the Sherpa parton shower using the ME+PS@LO prescription [48]. In both cases (W / Z+jets or \(\gamma \)+jets production), the CT10 PDF set [49] is used in conjunction with dedicated parton shower-tuning developed by the authors of Sherpa. The W / Z + jets events are normalized to their NNLO cross-sections [50]. For the \(\gamma \)+jets process the LO cross-section, taken directly from the Sherpa MC generator, is multiplied by a correction factor as described in Sect. 7.

For the generation of \(t\bar{t}\) and single-top processes in the Wt and s-channel [51] the Powheg-Box v2 [52] generator is used with the CT10 PDF set. The electroweak (EW) t-channel single-top events are generated using the Powheg-Box v1 generator. This generator uses the four-flavour scheme for the NLO matrix-element calculations together with the fixed four-flavour PDF set CT10f4 [49]. For this process, the decay of the top quark is simulated using MadSpin tool [53] preserving all spin correlations, while for all processes the parton shower, fragmentation, and the underlying event are generated using Pythia 6.428 [54] with the CTEQ6L1 [55] PDF set and the corresponding Perugia 2012 tune (P2012) [56]. The top quark mass is set to 172.5 \({\mathrm{GeV}}\). The \(h_\mathrm{damp}\) parameter, which controls the \(p_{\text {T}} \) of the first additional emission beyond the Born configuration, is set to the mass of the top quark. The main effect of this is to regulate the high-\(p_{\text {T}} \) emission against which the ttbar system recoils [51]. The \(t\bar{t}\) events are normalized to the NNLO+NNLL  [57, 58]. The s- and t-channel single-top events are normalized to the NLO cross-sections [59, 60], and the Wt-channel single-top events are normalized to the NNLO+NNLL [61, 62].

For the generation of \(t\bar{t}\) + EW processes (\(t\bar{t} + W/Z/WW\)) [63], the MG5_aMC@NLO [30] generator at LO interfaced to the Pythia 8.186 parton-shower model is used, with up to two (\(t\bar{t}+W\)), one (\(t\bar{t}+Z\)) or no (\(t\bar{t}+WW\)) extra partons included in the matrix element. The ATLAS underlying-event tune A14 is used together with the NNPDF2.3LO PDF set. The events are normalized to their respective NLO cross-sections [64, 65].

Diboson processes (WW, WZ, ZZ) [66] are simulated using the Sherpa 2.1.1 generator. For processes with four charged leptons (4\(\ell \)), three charged leptons and a neutrino (3\(\ell \)+1\(\nu \)) or two charged leptons and two neutrinos (2\(\ell \)+2\(\nu \)), the matrix elements contain all diagrams with four electroweak vertices, and are calculated for up to one (4\(\ell \), 2\(\ell \)+2\(\nu \)) or no partons (3\(\ell \)+1\(\nu \)) at NLO and up to three partons at LO using the Comix and OpenLoops matrix-element generators, and merged with the Sherpa parton shower using the ME+PS@NLO prescription. For processes in which one of the bosons decays hadronically and the other leptonically, matrix elements are calculated for up to one (ZZ) or no (WW, WZ) additional partons at NLO and for up to three additional partons at LO using the Comix and OpenLoops matrix-element generators, and merged with the Sherpa parton shower using the ME+PS@NLO prescription. In all cases, the CT10 PDF set is used in conjunction with a dedicated parton-shower tuning developed by the authors of Sherpa. The generator cross-sections are used in this case.

The multi-jet background is generated with Pythia 8.186 using the A14 underlying-event tune and the NNPDF2.3LO parton distribution functions.

For all Standard Model background samples the response of the detector to particles is modelled with a full ATLAS detector simulation [67] based on Geant4 [68]. Signal samples are prepared using a fast simulation based on a parameterization of the performance of the ATLAS electromagnetic and hadronic calorimeters [69] and on Geant4 elsewhere.

All simulated events are overlaid with multiple pp collisions simulated with the soft QCD processes of Pythia 8.186 using the A2 tune [33] and the MSTW2008LO parton distribution functions [70]. The simulations are not reweighted to match the distribution of the mean number of interactions observed in data. It was checked that the effect of such pile-up reweighting is completely negligible.

4 Object reconstruction and identification

The reconstructed primary vertex of the event is required to be consistent with the luminous region and to have at least two associated tracks with \(p_{\text {T}} > 400\) \({\mathrm{MeV}}\). When more than one such vertex is found, the vertex with the largest \(\sum p_{\text {T}} ^2\) of the associated tracks is chosen.

Jet candidates are reconstructed using the anti-\(k_{t}\) jet clustering algorithm [71, 72] with jet radius parameter of 0.4 and starting from clusters of calorimeter cells [73]. The jets are corrected for energy from pile-up using the method suggested in Ref. [74]: a contribution equal to the product of the jet area and the median energy density of the event is subtracted from the jet energy [75]. Further corrections, referred to as the jet energy scale corrections, are derived from MC simulation and data and used to calibrate on average the energies of jets to the scale of their constituent particles [76]. Only jet candidates with \(p_\mathrm {T}> 20~{\mathrm{GeV}}\) and \(|\eta |<2.8\) after all corrections are retained. An algorithm based on boosted decision trees, ‘MV2c20’ [77], is used to identify jets containing a b-hadron (b-jets), with an operating point corresponding to an efficiency of 77 % in simulated \(t\bar{t}\) events, along with a rejection factor of 140 for gluon and light-quark jets and of 4.5 for charm jets [77, 78]. Candidate b-tagged jets are required to have \(p_\mathrm {T}> 50~{\mathrm{GeV}}\) and \(|\eta |<2.5\). Events with jets originating from detector noise and non-collision background are rejected if the jets fail to satisfy the ‘LooseBad’ quality criteria, or if at least one of the two leading jets with \(p_{\text {T}} >100~{\mathrm{GeV}}\) fails to satisfy the ‘TightBad’ quality criteria, both described in Ref. [79]. These selections affect less than 1 % of the events used in the search.

Two different classes of reconstructed lepton candidates (electrons or muons) are used in this analysis. When selecting samples used for the search, events containing a ‘baseline’ electron or muon are rejected. The selections applied to identify baseline leptons are designed to maximize the efficiency with which W+jets and top quark background events are rejected. When selecting ‘control region’ samples for the purpose of estimating residual W+jets and top quark backgrounds, additional requirements are applied to leptons to ensure greater purity of the these backgrounds. These leptons are referred to as ‘high-purity’ leptons below and form a subset of the baseline leptons.

Baseline muon candidates are formed by combining information from the muon spectrometer and inner tracking detectors as described in Ref. [80] and are required to have \(p_\mathrm {T}> 10~ {\mathrm{GeV}}\) and \(|\eta | <2.7\). High-purity muon candidates must additionally have \(|\eta |<2.4\), the significance of the transverse impact parameter with respect to the primary vertex, \(|d_0^{\mathrm {PV}}|/\sigma (d_0^{\mathrm {PV}})<\) 3, the longitudinal impact parameter with respect to the primary vertex \(|z_0^{\mathrm {PV}} \mathrm {sin}(\theta )|<\) 0.5 mm, and to satisfy ‘GradientLoose’ isolation requirements described in Ref. [80] which rely on the use of tracking-based and calorimeter-based variables and implement a set of \(\eta \)- and \(p_{\text {T}} \)-dependent criteria. The leading muon is also required to have \(p_{\text {T}} > 25~{\mathrm{GeV}}\).

Baseline electron candidates are reconstructed from an isolated electromagnetic calorimeter energy deposit matched to an ID track and are required to have \(p_\mathrm {T}> 10~{\mathrm{GeV}}\), \(|\eta | < 2.47\), and to satisfy ‘Loose’ likelihood-based identification criteria described in Ref. [81]. High-purity electron candidates additionally must satisfy ‘Tight’ selection criteria described in Ref. [81], and the leading electron must have \(p_{\text {T}} >25~{\mathrm{GeV}}\). They are also required to have \(|d_0^{\mathrm {PV}}|/\sigma (d_0^{\mathrm {PV}})<\) 5, \(|z_0^{\mathrm {PV}} \mathrm {sin}(\theta )|<\) 0.5 mm, and to satisfy similar isolation requirements as those applied to high-purity muons.

After the selections described above, ambiguities between candidate jets with \(|\eta |<2.8\) and leptons are resolved as follows: first, any such jet candidate lying within a distance \(\Delta R\equiv \sqrt{(\Delta y)^2+(\Delta \phi )^2}=0.2\) of a baseline electron is discarded; then any baseline lepton candidate remaining within a distance \(\Delta R =0.4\) of any surviving jet candidate is discarded, except in the case where the lepton is a muon (which can radiate a photon and be misidentified as a jet) and the number of tracks associated with the jet is less than three.

Additional ambiguities between electrons and muons in a jet, originating from the decays of hadrons, are resolved to avoid double counting and/or remove non-isolated leptons: the electron is discarded if a baseline electron and a baseline muon share the same ID track. If two baseline electrons are within \(\Delta R\) = 0.05, the electron with the lowest \(p_{\text {T}} \) is discarded.

The measurement of the missing transverse momentum vector \(\varvec{E}\mathrm {^{miss}_T}\) (and its magnitude \(E_{\text {T}}^{\text {miss}} \)) is based on the calibrated transverse momenta of all electron, muon, photon and jet candidates and all tracks originating from the primary vertex and not associated with such objects [82].

Reconstructed photons, although not used in the main signal-event selection, are selected in the region used to constrain the Z+jets background, as explained in Sect. 7. Photon candidates are required to satisfy \(p_\mathrm {T}> 130~{\mathrm{GeV}}\) and \(|\eta | < 2.37\), to satisfy photon shower shape and electron rejection criteria [83], and to be isolated. Ambiguities between candidate jets and photons (when used in the event selection) are resolved by discarding any jet candidates lying within \(\Delta R\) = 0.4 of a photon candidate. Additional selections to remove ambiguities between electrons or muons and photons are applied such that the photon is discarded if it is within \(\Delta R\) = 0.4 of an electron or muon.

Corrections derived from data control samples are applied to account for differences between data and simulation for the lepton trigger and reconstruction efficiencies, the lepton momentum/energy scale and resolution, and for the efficiency and mis-tag rate of the b-tagging algorithm.

5 Analysis strategy and fit description

To search for a possible signal, selections are defined to enhance the signal relative to the SM background. These signal region (SR) selections are optimized to maximize the expected significance for each model considered using MC simulation for the signal and the SM backgrounds. To estimate the SM backgrounds in a consistent and robust fashion, corresponding control regions (CRs) are defined for each of the signal regions. They are chosen to be non-overlapping with the SR selections in order to provide independent data samples enriched in particular background sources, and are used to normalize the background MC simulation. The CR selections are optimized to have negligible SUSY signal contamination for the models near the previously excluded boundary [14], while minimizing the systematic uncertainties arising from the extrapolation of the CR event yields to estimate backgrounds in the SR. Cross-checks of the background estimates are performed with data in several validation regions (VRs) selected with requirements such that these regions do not overlap with the CR and SR selections, again with a low expected signal contamination.

To extract the final results, three different classes of likelihood fit are employed: background-only, model-independent and model-dependent fits [84]. A background-only fit is used to estimate the background yields in each SR. The fit is performed using as constraints only the observed event yields from the CRs associated with the SR, but not the SR itself. It is assumed that signal events from physics beyond the Standard Model (BSM) do not contribute to these yields. The scale factors (\(\mu _{W\mathrm +jets}\), \(\mu _{Z\mathrm +jets}\), \(\mu _\mathrm{Top}\), \(\mu _{\text {Multi-jet}}\)) are fitted in each CR attached to a SR. The expected background in the SR is based on the yields predicted by simulation, corrected by the scale factors derived from the fit. The systematic uncertainties and the MC statistical uncertainties in the expected values are included in the fit as nuisance parameters which are constrained by Gaussian distributions with widths corresponding to the sizes of the uncertainties considered and by Poisson distributions, respectively. The background-only fit is also used to estimate the background event yields in the VRs.

If no excess is observed, a model-independent fit is used to set upper limits on the number of BSM signal events in each SR. This fit proceeds in the same way as the background-only fit, except that the number of events observed in the SR is added as an input to the fit, and the BSM signal strength, constrained to be non-negative, is added as a free parameter. The observed and expected upper limits at 95 % confidence level (CL) on the number of events from BSM phenomena for each signal region (\(S_\mathrm{obs}^{95}\) and \(S_\mathrm{exp}^{95}\)) are derived using the \(CL_\mathrm{s}\) prescription [85], neglecting any possible signal contamination in the control regions. These limits, when normalized by the integrated luminosity of the data sample, may be interpreted as upper limits on the visible cross-section of BSM physics (\(\langle \epsilon \sigma \rangle _\mathrm{obs}^{95}\)), where the visible cross-section is defined as the product of production cross-section, acceptance and efficiency. The model-independent fit is also used to compute the one-sided p-value (\(p_0\)) of the background-only hypothesis, which quantifies the statistical significance of an excess.

Finally, model-dependent fits are used to set exclusion limits on the signal cross-sections for specific SUSY models. Such a fit proceeds in the same way as the model-independent fit, except that both the yield in the signal region and the signal contamination in the CRs are taken into account. Correlations between signal and background systematic uncertainties are taken into account where appropriate. Signal-yield systematic uncertainties due to detector effects and the theoretical uncertainties in the signal acceptance are included in the fit.

6 Event selection and signal regions definitions

Due to the high mass scale expected for the SUSY models considered in this study, the ‘effective mass’, \(m_{\mathrm {eff}}\), is a powerful discriminant between the signal and most SM backgrounds. When selecting events with at least \(N_\mathrm{j}\) jets, \(m_{\mathrm {eff}}(N_\mathrm{j})\) is defined to be the scalar sum of the transverse momenta of the leading \(N_\mathrm{j}\) jets and \(E_{\text {T}}^{\text {miss}}\) . Requirements placed on \(m_{\mathrm {eff}}\) and \(E_{\text {T}}^{\text {miss}}\) form the basis of this search by strongly suppressing the multi-jet background where jet energy mismeasurement generates missing transverse momentum. The final signal selection uses requirements on both \(m_{\mathrm {eff}}(\mathrm{incl.})\), which sums over all jets with \(p_\mathrm {T}>50~{\mathrm{GeV}}\) and \(E_{\text {T}}^{\text {miss}}\) , which is required to be larger than 200 \({\mathrm{GeV}}\).

Following the object reconstruction described in Sect. 4, events are discarded if a baseline electron or muon with \(p_{\text {T}} >10~{\mathrm{GeV}}\) remains, or if they contain a jet failing to satisfy quality selection criteria designed to suppress detector noise and non-collision backgrounds (described in Sect. 4). Events are also rejected if no jets with \(p_{\text {T}} >50~{\mathrm{GeV}}\) are found. Reconstructed photons and hadronically decaying \(\tau \)-leptons are not used in SR selections.

In order to maximize the sensitivity in the \((m_{\tilde{g}},m_{\tilde{q}})\) plane, a variety of signal regions are defined. Squarks typically generate at least one jet in their decays, for instance through \(\tilde{q} \rightarrow q \tilde{\chi }^0_1 \), while gluinos typically generate at least two jets, for instance through \(\tilde{g} \rightarrow q \bar{q} \tilde{\chi }^0_1 \). Processes contributing to \(\tilde{q} \tilde{q} \) and \(\tilde{g} \tilde{g} \) final states therefore lead to events containing at least two or four jets, respectively. Decays of heavy SUSY and SM particles produced in longer \(\tilde{q} \) and \(\tilde{g} \) decay cascades (e.g. \(\tilde{\chi }^\pm _1 \rightarrow qq'\tilde{\chi }^0_1 \)) tend to further increase the jet multiplicity in the final state.

Seven inclusive SRs characterized by increasing minimum jet multiplicity from two to six, are defined in Table 2. Some of them require the same jet-multiplicity, but are distinguished by increasing background rejection, ranging from ‘loose’ (labelled ‘l’) to ‘tight’ (labelled ‘t’).

In each region, different thresholds are applied on jet momenta and on \(\Delta \phi (\text {jet},\varvec{E}\mathrm {^{miss}_T})_\mathrm {min}\), which is defined to be the smallest azimuthal separation between \(\varvec{E}\mathrm {^{miss}_T}\) and the momenta of any of the reconstructed jets with \(p_\mathrm {T}>50~{\mathrm{GeV}}\). Requirements on \(\Delta \phi (\text {jet},\varvec{E}\mathrm {^{miss}_T})_\mathrm {min}\) and \(E_{\text {T}}^{\text {miss}}/m_{\mathrm {eff}}(N_\mathrm{j})\) are designed to reduce the background from multi-jet processes. For the SRs which are optimized for squark-pair (gluino-pair) production followed by the direct decay of squarks (gluinos), the selection requires \(\Delta \phi (\text {jet},\varvec{E}\mathrm {^{miss}_T})_\mathrm {min}>0.8\) (\(\Delta \phi (\text {jet},\varvec{E}\mathrm {^{miss}_T})_\mathrm {min}>0.4\)) using up to three leading jets present in the event. For the SRs requiring at least four jets in the final state, an additional requirement \(\Delta \phi (\text {jet},\varvec{E}\mathrm {^{miss}_T})_\mathrm {min}>0.2\) is placed on all jets. Signal region 2jm makes use of the presence of jets due to initial-state radiation by requiring a higher \(p_{\text {T}} \) threshold for the most energetic jet in the event, and is optimized to target models with small mass differences between the SUSY particles (compressed scenarios).

In the 2-jet SRs the requirement on \(E_{\text {T}}^{\text {miss}}/m_{\mathrm {eff}}(N_\mathrm{j})\) is replaced by a requirement on \(E_{\text {T}}^{\text {miss}}/\sqrt{H_\mathrm{T}}\) (where \(H_\mathrm{T}\) is defined as the scalar sum of the transverse momenta of all jets), which was found to lead to enhanced sensitivity to models characterized by \(\tilde{q} \tilde{q} \) production. In the other regions, additional suppression of background processes is based on the aplanarity variable, which is defined as \(A = 3/2 \lambda _3\), where \(\lambda _3\) is the smallest eigenvalue of the normalized momentum tensor of the jets [86].

Table 2 Selection criteria and targeted signal model used to define each of the signal regions in the analysis. Each SR is labelled with the inclusive jet multiplicity considered (‘2j’, ‘4j’ etc.) together with the degree of background rejection. The latter is denoted by labels ‘l’ (‘loose’), ‘m’ (‘medium’) and ‘t’ (‘tight’). The \(E_{\text {T}}^{\text {miss}}/m_{\mathrm {eff}}(N_\mathrm{j})\) cut in any \(N_\mathrm{j}\)-jet channel uses a value of \(m_{\mathrm {eff}}\) constructed from only the leading \(N_\mathrm{j}\) jets (\(m_{\mathrm {eff}}(N_\mathrm{j})\)). However, the final \(m_{\mathrm {eff}}(\mathrm{incl.})\) selection, which is used to define the signal regions, includes all jets with \(p_{\text {T}} >50~{\mathrm{GeV}}\)

7 Background estimation and validation

Standard Model background processes contribute to the event counts in the signal regions. The dominant sources are: \(Z+\)jets, \(W+\)jets, top quark pairs, single top quarks, dibosons and multi-jet production. Diboson production is estimated with MC simulated data normalized to NLO cross-section predictions, as described in Sect. 3. Most of the W+jets background is composed of \(W\rightarrow \tau \nu \) events in which the \(\tau \)-lepton decays to hadrons, with additional contributions from \(W\rightarrow e\nu , \mu \nu \) events in which no baseline electron or muon is reconstructed. The largest part of the Z+jets background comes from the irreducible component in which \(Z\rightarrow \nu \bar{\nu }\) decays generate large \(E_{\text {T}}^{\text {miss}} \). Top quark pair production followed by semileptonic decays, in particular \(t \bar{t} \rightarrow b \bar{b} \tau \nu q q' \) (with the \(\tau \)-lepton decaying to hadrons), as well as single-top-quark events, can also generate large \(E_{\text {T}}^{\text {miss}} \) and satisfy the jet and lepton-veto requirements. The multi-jet background in the signal regions is due to missing transverse momentum from misreconstruction of jet energies in the calorimeters, as well as neutrino production in semileptonic decays of heavy-flavour hadrons. After applying the requirements based on \(\Delta \phi (\text {jet},\varvec{E}\mathrm {^{miss}_T})_\mathrm {min}\) and \(E_{\text {T}}^{\text {miss}}/m_{\mathrm {eff}}(N_\mathrm{j})\) listed in Table 2 the remaining multi-jet background is negligible.

Table 3 Control regions used in the analysis. Also listed are the main targeted background in the SR in each case, the process used to model the background, and the main CR requirement(s) used to select this process. The transverse momenta of high-purity leptons (photons) used to select CR events must exceed 25 (130) \({\mathrm{GeV}}\). The jet \(p_{\text {T}} \) thresholds and \(m_{\mathrm {eff}}(\mathrm{incl.})\) selections match those used in the corresponding SRs

In order to estimate the backgrounds in a consistent and robust fashion, four control regions are defined for each of the seven signal regions, giving 28 CRs in total. The CR selections are optimized to maintain adequate statistical precision while minimizing the systematic uncertainties arising from the extrapolation of the CR event yield to estimate the background in the SR. This latter requirement is addressed through the use of CR jet \(p_{\text {T}} \) thresholds and \(m_{\mathrm {eff}}\)(incl.) selections which match those used in the SR. The CR definitions are listed in Table 3.

The CR\(\gamma \) region is used to estimate the contribution of \(Z(\rightarrow \nu \bar{\nu })\)+jets background events to each SR by selecting a sample of \(\gamma \)+jets events with \(p_{\text {T}} (\gamma )>130~{\mathrm{GeV}}\) and then treating the reconstructed photon as contributing to \(E_{\text {T}}^{\text {miss}} \). For \(p_{\text {T}} (\gamma )\) significantly larger than \(m_Z\) the kinematic properties of such events strongly resemble those of Z+jets events [13]. In order to reduce the theoretical uncertainties associated with the \(Z/\gamma ^*\)+jets background expectations in SRs arising from the use of LO \(\gamma \)+jets cross-sections, a correction factor is applied to the CR\(\gamma \) events. This correction factor, \(\kappa =1.5 \pm 0.1\), is determined by comparing CR\(\gamma \) observations with those in a highly populated auxiliary control region defined by selecting events with two electrons or muons for which the invariant mass lies within 25 \({\mathrm{GeV}}\) of the mass of the Z boson, satisfying \(200~{\mathrm{GeV}}< | \varvec{E}\mathrm {^{miss}_T} + \varvec{p}\mathrm {_T}(\ell \bar{\ell }) | < 300~{\mathrm{GeV}}\), together with at least two jets.

The CRW and CRT regions aim to select samples rich in \(W(\rightarrow \ell \nu )\)+jets and semileptonic \(t\bar{t}\) background events respectively. Consequently, they differ in their number of b-jets (zero or greater or equal to one respectively) but apply the same selection requirements on the transverse mass \(m_\mathrm{T}\) formed by the \(E_\mathrm{T}^{\mathrm {miss}}\) and a high-purity lepton with \(p_\mathrm{T}\) > 25 \({\mathrm{GeV}}\). These samples are used to estimate respectively the W+jets and combined \(t\bar{t}\) and single-top background populations, treating the lepton as a jet with the same momentum to model background events in which a hadronically decaying \(\tau \)-lepton is produced or events in which no baseline electron or muon is reconstructed because it is outside the detector acceptance or below the required \(p_{\text {T}} \) threshold. The CRW and CRT selections omit the SR selection requirements on \(\Delta \phi (\text {jet},\varvec{E}\mathrm {^{miss}_T})_\mathrm {min}\) or \(E_{\text {T}}^{\text {miss}}/m_{\mathrm {eff}}(N_\mathrm{j})\) (\(E_{\text {T}}^{\text {miss}}/\sqrt{H_\mathrm{T}}\) where appropriate) in order to increase the number of CR data events without significantly increasing the theoretical uncertainties associated with the background estimation procedure.

The CRQ region uses reversed selection requirements on \(\Delta \phi (\text {jet},\varvec{E}\mathrm {^{miss}_T})_\mathrm {min}\) and on \(E_{\text {T}}^{\text {miss}}/m_{\mathrm {eff}}(N_\mathrm{j})\) (or \(E_{\text {T}}^{\text {miss}}/\sqrt{H_\mathrm{T}}\) where appropriate) to produce samples enriched in multi-jet background events.

As an example, the \( m_{\mathrm {eff}}(\mathrm{incl.})\) distributions in control regions associated with SR 4jt are shown in Fig. 2. In all CRs, the data are consistent with the pre-fit MC background prediction within uncertainties, although the overall normalization is lower by approximately one standard deviation.

Fig. 2
figure 2

Observed \(m_{\mathrm {eff}}(\mathrm{incl.})\) distributions in control regions a CR\(\gamma \), b CRW and c CRT after selecting events with \(E_{\text {T}}^{\text {miss}} > 200 ~{\mathrm{GeV}}\) and at least four energetic jets with the corresponding transverse momenta as indicated in Table 2 for SR 4jt. No selection requirements on \(\Delta \phi (\text {jet},\varvec{E}\mathrm {^{miss}_T})_\mathrm {min}\) or \(E_{\text {T}}^{\text {miss}}/m_{\mathrm {eff}}(N_\mathrm{j})\) are applied in these distributions. The arrows indicate the values at which the requirements on \(m_{\mathrm {eff}}(\mathrm{incl.})\) are applied. The histograms denote the pre-fit MC background expectations, normalized to cross-section times integrated luminosity. The last bin includes the overflow. In the lower panels the hatched (red) error bands denote the combined experimental, MC statistical and theoretical modelling uncertainties

The background estimation procedure is validated by comparing the numbers of events observed in the VRs to the corresponding SM background expectations obtained from the background-only fits. Several VR samples are selected with requirements distinct from those used in the CRs, which maintain a low probability of signal contamination.

The CR\(\gamma \) estimates of the \(Z(\rightarrow \nu \bar{\nu })\)+jets background are validated using the samples of \(Z(\rightarrow \ell \bar{\ell })\)+jets events selected by requiring high-purity lepton pairs of opposite sign and identical flavour for which the dilepton invariant mass lies within 25 \({\mathrm{GeV}}\) of the mass of the Z boson (VRZ). In VRZ, the leptons are treated as contributing to \(E_{\text {T}}^{\text {miss}} \).

The CRW and CRT estimates of the W+jets and top quark background are validated with the same CRW and CRT selections, but reinstating the requirement on \(\Delta \phi (\text {jet},\varvec{E}\mathrm {^{miss}_T})_\mathrm {min}\) and \(E_{\text {T}}^{\text {miss}}/m_{\mathrm {eff}}(N_\mathrm{j})\) (or \(E_{\text {T}}^{\text {miss}}/\sqrt{H_\mathrm{T}}\) as appropriate), and treating the lepton either as a jet (VRW, VRT) or as contributing to \(E_{\text {T}}^{\text {miss}} \) (VRW\(\nu \), VRT\(\nu \)).

The CRQ estimates of the multi-jet background are validated with VRs for which the CRQ selection is applied, but with the SR \(E_{\text {T}}^{\text {miss}}/m_{\mathrm {eff}}(N_\mathrm{j})\) (\(E_{\text {T}}^{\text {miss}}/\sqrt{H_\mathrm{T}}\)) requirement reinstated (VRQa), or with a requirement of an intermediate value of \(\Delta \phi (\text {jet},\varvec{E}\mathrm {^{miss}_T})_\mathrm {min}\) applied (VRQb).

The results of the validation procedure are shown in Fig. 3. The entries in the matrix are the differences between the numbers of observed and expected events expressed as fractions of the one-standard deviation \((1\sigma )\) uncertainties on the latter. Most VR observations lie within \(1\sigma \) of the background expectations, with the largest discrepancy out of 49 VRs being \(-1.5\sigma \) the CRQb associated with the SR 4jt.

Fig. 3
figure 3

Differences between the numbers of observed events in data and the SM background predictions for each VR, expressed as a fraction of the total uncertainty which combines the uncertainty on the background expectations, and the expected statistical uncertainty of the test obtained from the number of expected events

8 Systematic uncertainties

Systematic uncertainties in background estimates arise from the use of extrapolation factors which relate observations in the control regions to background expectations in the signal regions, and from the MC modelling of minor backgrounds.

The overall background uncertainties, detailed in Table 4, range from 8 % in SR 2jl to 29 % in SR 6jt. In SR 2jl the loose selection minimizes theoretical uncertainties and the impact of statistical fluctuations in the CRs, while the opposite is true in SR 6jt.

For the backgrounds estimated with MC simulation-derived extrapolation factors, the primary common sources of systematic uncertainty are the jet energy scale (JES) calibration, jet energy resolution (JER), theoretical uncertainties, and limited event yields in the MC samples and data CRs. Correlations between uncertainties (for instance between JES or JER uncertainties in CRs and SRs) are taken into account where appropriate.

The JES uncertainty was measured using the techniques described in Refs. [76, 87, 88]. The JER uncertainty is estimated using the methods discussed in Refs. [76, 89]. An additional uncertainty in the modelling of energy not associated with reconstructed objects, used in the calculation of \(E_{\text {T}}^{\text {miss}} \) and measured with unassociated charged tracks, is also included. The combined JES, JER and \(E_{\text {T}}^{\text {miss}} \) (Jet/\(E_{\text {T}}^{\text {miss}} \)) uncertainty ranges from <1 % of the expected background in 2-jet SRs to 5 % in SR 6jt.

Uncertainties arising from theoretical modelling of background processes are evaluated by comparing samples produced with different MC generators. The W / Z+jets events generated with Sherpa are compared to events generated with MG5_aMC@NLO at leading order and interfaced to the Pythia 8.186 parton shower model. Uncertainties in the modelling of top quark pair production are estimated by comparing Powheg-Box to aMc@Nlo [90], and by accounting for different generator and radiation tunes. Uncertainties associated with PDF modelling of top quark pair production are found to be negligible. Uncertainties in diboson production due to PDF, renormalization, factorization and resummation scale uncertainties (estimated by increasing and decreasing the scales used in the MC generators by a factor of two) are accounted for by applying a uniform 50 % uncertainty in all SRs, and are the dominant source of uncertainty in SRs 2jl and 2jm. Uncertainties associated with the modelling of Z+jets production are largest in the SRs with tight selection cuts (up to 14 %). The statistical uncertainty arising from the use of MC samples is largest (8 %) in SR 6jt. The uncertainties arising from the data-driven correction procedure applied to events selected in the CR\(\gamma \) region, described in Sect. 7, are included in Table 4 under ‘CR\(\gamma \) corr. factor’ and reach a value of 4 % in most of the SRs. The impact of lepton reconstruction uncertainties, and of the uncertainties related to the b-tag/b-veto efficiency, on the overall background uncertainty are found to be negligible for all SRs. The total background uncertainties for all SRs, broken down into the main contributing sources, are summarized in Table 4.

Table 4 Breakdown of the dominant systematic uncertainties in the background estimates. The individual uncertainties can be correlated, and do not necessarily add in quadrature to the total background uncertainty. \(\Delta \mu \) uncertainties are the result of the control region statistical uncertainties and the systematic uncertainties entering a specific control region. In brackets, uncertainties are given relative to the expected total background yield, also presented in the Table. Empty cells (indicated by a ‘–’) correspond to uncertainties lower than 1 per mil

9 Results, interpretation and limits

The number of events observed in the data and the number of SM events expected to enter each of the signal regions, determined using the background-only fit, are shown in Table 5 and Fig. 4. The pre-fit background expectations are also shown in Table 5 for comparison. The normalisation factors extracted simultaneously through the fit range for the different signal regions between 0.7 and 1.2 for W+jets, 0.4 and 0.8 for \(t\bar{t} \)(+EW) + single top, and 1.0 and 1.6 for \(Z/\gamma ^*\)+jets backgrounds.

Table 5 Numbers of events observed in the signal regions used in the analysis compared with background expectations obtained from the fits described in the text. No signal contribution is considered in the CRs for the fit. Empty cells (indicated by a ‘–’) correspond to estimates lower than 0.01. The p-values (\(p_{0}\)) give the probabilities of the observations being consistent with the estimated backgrounds. For an observed number of events lower than expected, the p value is truncated at 0.5. Between parentheses, p-values are also given as the number of equivalent Gaussian standard deviations (Z). Also shown are 95 % CL upper limits on the visible cross-section (\(\langle \epsilon \sigma \rangle _\mathrm{obs}^{95}\)), the visible number of signal events (\(S_\mathrm{obs}^{95}\) ) and the number of signal events (\(S_\mathrm{exp}^{95}\)) given the expected number of background events (and \(\pm 1\sigma \) excursions of the expectation)
Fig. 4
figure 4

Comparison of the observed and expected event yields as a function of signal region. The background expectations are those obtained from the background-only fits, presented in Table 5

Distributions of \(m_{\mathrm {eff}}(\mathrm{incl.})\) obtained before the final selections on this quantity (but after applying all other selections), for data and the different MC samples normalized with the theoretical cross-sections, i.e. before applying the normalization from the CR fit, are shown in Figs. 5 and 6. Examples of typical expected SUSY signals are shown for illustration. These signals correspond to the processes to which each SR is primarily sensitive – \(\tilde{q} \tilde{q} \) production for the lower jet-multiplicity SRs and \(\tilde{g} \tilde{g} \) production for the higher jet-multiplicity SRs. In these figures, data and background distributions largely agree within uncertainties. The differences seen in the lower regions of \(m_{\mathrm {eff}}(\mathrm{incl.})\) distribution (1.2 – 2.0 \({\mathrm{TeV}}\)) in Fig. 6 do not affect the background expectations in the signal regions since the backgrounds are normalized using control regions (Table 3) with the same \(m_{\mathrm {eff}}(\mathrm{incl.})\) selections. The fit to the CRs for each SR compensates for the differences related to the overall normalization of the background seen in Figs. 5 and 6, leading to the good agreements between data and post-fit expectations in the SRs observed in Table 5 and Fig. 4.

Fig. 5
figure 5

Observed \(m_{\mathrm {eff}}(\mathrm{incl.})\) distributions for the a 2jl, b 2jm, c 2jt signal regions. The histograms denote the MC background expectations prior to the fits described in the text, normalized to cross-section times integrated luminosity. The last bin includes the overflow. In the lower panels the hatched (red) error bands denote the combined experimental, MC statistical and theoretical modelling uncertainties. The arrows indicate the values at which the requirements on \(m_{\mathrm {eff}}(\mathrm{incl.})\) are applied. Expected distributions for benchmark model points, normalized to NLO+NLL cross-section (Sect. 3) times integrated luminosity, are also shown for comparison (masses in \({\mathrm{GeV}}\))

Fig. 6
figure 6

Observed \(m_{\mathrm {eff}}(\mathrm{incl.})\) distributions for the a 4jt, b 5j, c 6jm and d 6jt signal regions. The histograms denote the MC background expectations prior to the fits described in the text, normalized to cross-section times integrated luminosity. The last bin includes the overflow. In the lower panels the hatched (red) error bands denote the combined experimental, MC statistical and theoretical modelling uncertainties. The arrows indicate the values at which the requirements on \(m_{\mathrm {eff}}(\mathrm{incl.})\) are applied. Expected distributions for benchmark model points, normalized to NLO+NLL cross-section (Sect. 3) times integrated luminosity, are also shown for comparison (masses in \({\mathrm{GeV}}\))

In the absence of a statistically significant excess, limits are set on contributions to the SRs from BSM physics. Upper limits at 95 % CL on the number of BSM signal events in each SR and the corresponding visible BSM cross-section are derived from the model-independent fits described in Sect. 5 using the \(CL_\mathrm{s}\) prescription. Limits are evaluated using MC pseudo-experiments. The results are presented in Table 5.

The model-dependent fits in all the SRs are then used to set limits on specific classes of SUSY models, using the result from the SR with the best expected sensitivity at each point in each model parameter space. ‘Observed limits’ are calculated from the observed SR event yields for the nominal signal cross-section. ‘Expected limits’ are calculated by setting the nominal event yield in each SR to the corresponding mean expected background.

In Fig. 7, limits are shown for two classes of simplified models in which only direct production of light-flavour squark or gluino pairs are considered. In these simplified model scenarios, the upper limit of the excluded light-flavour squark mass region is 1.03 \({\mathrm{TeV}}\) assuming massless \(\tilde{\chi }^0_1 \), as obtained from the signal region 2jt. The corresponding limit on the gluino mass is 1.51 \({\mathrm{TeV}}\)  if the \(\tilde{\chi }^0_1 \) is massless, as obtained from the signal region 4jt. The best sensitivity in the region of parameter space where the mass difference between the squark (gluino) and the lightest neutralino is small is obtained from the signal region 2jm.

Fig. 7
figure 7

Exclusion limits for direct production of a light-flavour squark pairs with decoupled gluinos and b gluino pairs with decoupled squarks. Gluinos (light-flavour squarks) are required to decay to two quarks (one quark) and a neutralino LSP. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The blue dashed lines show the expected limits at 95 % CL, with the light (yellow) bands indicating the \(1\sigma \) excursions due to experimental and background-only theoretical uncertainties. Observed limits are indicated by medium dark (maroon) curves where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalization and factorization scale and PDF uncertainties. Results are compared with the observed limits obtained by the previous ATLAS search [15]. The black stars indicate the benchmark models used in Figs. 5 and 6

In Fig. 8, limits are shown for pair-produced gluinos each decaying via an intermediate \(\tilde{\chi }^\pm _1 \) to two quarks, a W boson and a \(\tilde{\chi }^0_1 \). Results are presented for simplified models in which the mass of the chargino \(\tilde{\chi }^\pm _1 \) is fixed to \(m(\tilde{\chi }^\pm _1)=(m(\tilde{g})+m(\tilde{\chi }^0_1))/2\). For a \(\tilde{\chi }^0_1 \) mass of \(\sim \) 200 \({\mathrm{GeV}}\), the lower limit on the gluino mass, obtained from the signal region 4jt, extends up to 1.5 \({\mathrm{TeV}}\) in this model. In the region of parameter space where the mass difference between the gluino and the lightest neutralino is small, the best sensitivity is obtained from the signal region 2jm. Results are compared with the observed limits obtained from the statistical combination of the search with no lepton and the search with one isolated lepton, high-\(p_{\text {T}} \) jets and missing transverse momentum performed at ATLAS [15] using the 8 \({\mathrm{TeV}}\) data. Statistical combinations of these two searches, designed to be statistically independent in their signal and control region definitions, are performed in order to increase the exclusion reach in models in which at least two analyses obtain comparable sensitivities, and still provide the strongest exclusion limits in the region of parameter space in which the mass of gluino is between 700 and 1100 \({\mathrm{GeV}}\)  and the \(\tilde{\chi }^0_1 \) mass is above \(\sim \) 500 \({\mathrm{GeV}}\).

Fig. 8
figure 8

Exclusion limits for pair-produced gluinos each decaying via an intermediate \(\tilde{\chi }^\pm _1 \) to two quarks, a W boson and a \(\tilde{\chi }^0_1 \) for models with a fixed \(m(\tilde{\chi }^\pm _1)=(m(\tilde{g})+m(\tilde{\chi }^0_1))/2\) and varying values of \(m(\tilde{g})\) and \(m({\tilde{\chi }^0_1})\). Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The blue dashed lines show the expected limits at 95 % CL, with the light (yellow) bands indicating the \(1\sigma \) excursions due to experimental and background-only theoretical uncertainties. Observed limits are indicated by medium dark (maroon) curves where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalization and factorization scale and PDF uncertainties. Results are compared with the observed limits obtained from the statistical combination of the search with no lepton and the search with one isolated lepton, high-\(p_{\text {T}} \) jets and missing transverse momentum performed at ATLAS [15]. The black stars indicate the benchmark models used in Fig. 6

10 Conclusion

This paper reports a search for squarks and gluinos in final states containing high-\(p_{\text {T}}\) jets, large missing transverse momentum but no electrons or muons, based on a 3.2 \(\mathrm{fb}^{-1}\) dataset of \(\sqrt{s}=13 {\mathrm{TeV}}\) proton–proton collisions recorded by the ATLAS experiment at the LHC in 2015. Good agreement is seen between the numbers of events observed in the data and the numbers of events expected from SM processes.

Results are interpreted in terms of simplified models with only light-flavour squarks, or gluinos, together with a neutralino LSP, with the masses of all the other SUSY particles set beyond the reach of the LHC. For a massless lightest neutralino, gluino masses below 1.51 \({\mathrm{TeV}}\) are excluded at the 95 % confidence level in a simplified model with only gluinos and the lightest neutralino. For a simplified model involving the strong production of squarks of the first and second generations, with decays to a massless lightest neutralino, squark masses below 1.03 \({\mathrm{TeV}}\) are excluded, assuming mass-degenerate squarks. In simplified models with pair-produced gluinos, each decaying via an intermediate \(\tilde{\chi }^\pm _1 \) to two quarks, a W boson and a \(\tilde{\chi }^0_1 \), gluino masses below 1.5 \({\mathrm{TeV}}\) are excluded for \(\tilde{\chi }^0_1 \) masses of \(\sim \) 200 \({\mathrm{GeV}}\). These results substantially extend the region of supersymmetric parameter space excluded by previous LHC searches.