Search for gluinos in events with an isolated lepton, jets and missing transverse momentum at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}$$\end{document}s = 13 Te V with the ATLAS detector

The results of a search for gluinos in final states with an isolated electron or muon, multiple jets and large missing transverse momentum using proton–proton collision data at a centre-of-mass energy of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s} = 13 \mathrm{{\ Te V}}$$\end{document}s=13TeV are presented. The dataset used was recorded in 2015 by the ATLAS experiment at the Large Hadron Collider and corresponds to an integrated luminosity of 3.2 fb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{-1}$$\end{document}-1. Six signal selections are defined that best exploit the signal characteristics. The data agree with the Standard Model background expectation in all six signal selections, and the largest deviation is a 2.1 standard deviation excess. The results are interpreted in a simplified model where pair-produced gluinos decay via the lightest chargino to the lightest neutralino. In this model, gluinos are excluded up to masses of approximately 1.6 Te V depending on the mass spectrum of the simplified model, thus surpassing the limits of previous searches.


Introduction
Supersymmetry (SUSY) [1][2][3][4][5][6] is a theoretical framework of physics beyond the Standard Model (SM) that predicts for each SM particle the existence of a supersymmetric partner differing by half a unit of spin. The partner particles of the SM fermions (quarks and leptons) are the scalar squarks (q) and sleptons (˜ ). In the boson sector, the supersymmetric partner of the gluon is the fermionic gluino (g), whereas the supersymmetric partners of the Higgs (higgsinos) and the electroweak gauge bosons (winos and bino) mix to form charged and neutral mass eigenstates called charginos (χ ± 1,2 ) and neutralinos (χ 0 1,2,3,4 ). In the Minimal Supersymmetric extension of the Standard Model (MSSM) [7,8] two scalar Higgs doublets along with their Higgsino partners are predicted. SUSY addresses the SM hierarchy problem [9][10][11][12] provided that the masses of at least some of the supersyme-mail: atlas.publications@cern.ch metric particles (most notably the higgsinos, the top squarks and the gluinos) are near the TeV scale.
In R-parity-conserving SUSY [13], gluinos might be pairproduced at the Large Hadron Collider (LHC) via the strong interaction and decay either directly or via intermediate states to the lightest supersymmetric particle (LSP). The LSP is stable and is assumed to be only weakly interacting, making it a candidate for dark matter [14,15].
This paper considers a SUSY-inspired model where pairproduced gluinos decay via the lightest chargino (χ ± 1 ) to the LSP, which is assumed to be the lightest neutralino (χ 0 1 ). The three-body decay of the gluino to the chargino proceeds viag → qq χ ± 1 . The chargino decays to the LSP by emitting an on-or off-shell W boson, depending on the available phase space. In the MSSM this decay chain is realised when the gluino decays via a virtual squark that is the partner of the left-handed SM quark, to the chargino with a dominant wino component. In the MSSM the mass of the chargino is independent of the mass of the gluino.
The experimental signature characterising this search consists of a lepton (electron or muon), several jets, and missing transverse momentum ( p miss T with magnitude E miss T ) from the undetectable neutralinos and neutrino(s). The analysis is based on two complementary sets of search channels. The first set requires a low transverse momentum ( p T ) lepton (7/6 < p T (e/μ) < 35 GeV), and is referred to as the softlepton channel, while the second set requires a highp T lepton ( p T (e/μ) > 35 GeV) and is referred to as the hard-lepton channel. The two sets target SUSY models with small and large mass differences between the predicted supersymmetric particles, respectively. The search uses the ATLAS data collected in proton-proton LHC collisions in 2015 corresponding to an integrated luminosity of 3.2 fb −1 at a centreof-mass energy of 13 TeV.
The analysis extends previous ATLAS searches with similar event selections [16] which were performed with data collected during the first data-taking campaign between 2010 and 2012 (LHC Run 1) at a centre-of-mass energy of up to 8 TeV. The results of all Run 1 ATLAS searches targeting squark and gluino pair production are summarised in Ref. [17]. The CMS Collaboration has performed similar searches for gluinos with decays via intermediate supersymmetric particles in Run 1 [18,19] and Run 2 [20].
This paper is structured as follows. After a brief description of the ATLAS detector in Sect. 2, the simulated data samples for the background and signal processes used in the analysis as well as the dataset and the trigger strategy are detailed in Sect. 3. The reconstructed objects and quantities used in the analysis are described in Sect. 4 and the event selection is presented in Sect. 5. The background estimation and the systematic uncertainties associated with the expected event yields are discussed in Sects. 6 and 7, respectively, while details of the statistical interpretation are given in Sect. 8. Finally, the results of the analysis are presented in Sect. 9 and are followed by a conclusion.

ATLAS detector
ATLAS [21] is a general-purpose detector with a forwardbackward symmetric design that provides almost full solid angle coverage around the interaction point. 1 The main components of ATLAS are the inner detector (ID), which is surrounded by a superconducting solenoid providing a 2 T axial magnetic field, the calorimeter system, and the muon spectrometer (MS), which is immersed in a magnetic field generated by three large superconducting toroidal magnets. The ID provides track reconstruction within |η| < 2.5, employing pixel detectors close to the beam pipe, silicon microstrip detectors at intermediate radii, and a straw-tube tracker with particle identification capabilities based on transition radiation at radii up to 1080 mm. The innermost pixel detector layer, the insertable B-layer [22], was added during the shutdown between LHC Run 1 and Run 2, at a radius of 33 mm around a new, narrower, beam pipe. The calorimeters cover |η| < 4.9, the forward region (3.2 < |η| < 4.9) being instrumented with a liquid-argon (LAr) calorimeter for both the electromagnetic and the hadronic measurements. In the central region, a lead/LAr electromagnetic calorimeter covers |η| < 3.2, while the hadronic calorimeter uses two different detector technologies, with scintillator tiles (|η| < 1.7) or 1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upward. Cylindrical coordinates (r, φ) are used in the transverse plane, φ being the azimuthal angle around the z-axis. The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2) and the rapidity is defined as y = 0.5 ln[(E + p z )/(E − p z )], where E is the energy and p z the longitudinal momentum of the object of interest. LAr (1.5 < |η| < 3.2) as the active medium. The MS consists of three layers of precision tracking chambers providing coverage over |η| < 2.7, while dedicated fast chambers allow triggering over |η| < 2.4. The ATLAS trigger system (developed from Ref. [23]) consists of a hardware-based first-level trigger and a software-based high-level trigger.

Simulated event samples and data samples
The signal model considered in this search is a simplified model [24-26] that has been used in previous similar ATLAS searches [16]. In this model, exclusive pair-production of gluinos is assumed. The gluinos decay via an intermediate chargino, here the lightest charginoχ ± 1 , into the lightest supersymmetric particle, the lightest neutralinoχ 0 1 . The branching ratio of each supersymmetric particle decay considered is assumed to be 100 %. Other supersymmetric particles not entering the decay chain described are not considered in this simplified model and their masses are set to high values. The gluino decay is assumed to proceed only via virtual first-and second-generation quarks, hence no bottom or top quarks are produced in the simplified model. The free parameters of the model are the masses of the gluino (mg), the chargino (mχ± 1 ), and the neutralino (mχ0 1 ). Two types of scenarios are considered: in the first type, the mass of the neutralino is fixed to 60 GeV, and the sensitivity is assessed as a function of the gluino mass and a mass-ratio parameter defined as x = (mχ± 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-leadinglogarithmic accuracy (NLL) [33][34][35][36][37]. The nominal crosssection and its uncertainty are taken from an envelope of cross-section predictions using different PDF sets and factorisation and renormalisation scales [38,39].
The simulated event samples for the SM backgrounds are summarised in Table 1 The response of the detector to particles is modelled with a full ATLAS detector simulation [58] using Geant4 [59], or using a fast simulation [60] based on a parameterisation of the performance of the electromagnetic and hadronic calorimeters and on Geant4 elsewhere. All background (signal) samples are prepared using the full (fast) detector simulation.
All simulated samples are generated with a varying number of minimum-bias interactions (simulated using Pythia 8 with the MSTW2008LO PDF set [61] and the A2 tune [62]) overlaid on the hard-scattering event to model the multiple proton-proton interactions in the same and the nearby bunch crossings. Corrections are applied to the simulated samples to account for differences between data and simulation for trigger, identification and reconstruction efficiencies.
The proton-proton data analysed in this paper were collected by ATLAS in 2015 at a centre-of-mass energy of 13 TeV. During this period the instantaneous luminosity of the LHC reached 5.0 × 10 33 cm −2 s −1 with a mean number of additional pp interactions per bunch crossing of approximately 14. After application of data-quality requirements related to the beam and detector conditions, the total integrated luminosity amounts to 3.2 fb −1 , with an associated uncertainty of ±5 %. These values are derived following the same methodology as the one detailed in Ref. [63].
The data are collected using an E miss T trigger with a threshold of 70 GeV. This trigger is close to fully efficient after applying the requirement on the offline E miss T to be larger than 200 GeV.

Object reconstruction and identification
The reconstructed primary vertex of an event is required to be consistent with the interaction region and to have at least two associated tracks with p T > 400 MeV. When more than one such vertex is found, the vertex with the largest p 2 T of the associated tracks is chosen.
In the analysis, a distinction is made between preselected reconstructed objects, which fulfil a set of basic criteria and are used in the E miss T computation, and signal objects that enter the various control, validation and signal regions and are subject to more stringent requirements.
Jets are reconstructed from topological clusters in the calorimeters using the anti-k t algorithm with a radius parameter R = 0.4 [64,65]. Prior to jet reconstruction, clusters are calibrated to the electromagnetic scale response. Additional correction factors derived from simulation and data are applied to the measured jet energy to calibrate it to the particle level [65]. To mitigate the contributions from pileup, the median energy density of all the jets in the event, multiplied by the jet area, is subtracted from the reconstructed jet energy [66,67]. Preselected jets are required to have p T > 20 GeV and |η| < 4.5. The contamination from cosmic rays, other sources of non-collision background and detector noise is suppressed using dedicated jet-quality criteria [68]: when such criteria are not fulfilled, the event is rejected.
Electron candidates are reconstructed using ID tracks matched to energy clusters in the electromagnetic calorimeter. They are identified according to the likelihood-based loose criteria [69]. Preselected electrons in the soft-lepton (hard-lepton) channel must satisfy p T > 7(10) GeV and |η| < 2.47. When the angular separation R = ( y) 2 + ( φ) 2 between an electron candidate and a preselected jet amounts to 0.2 < R(e, jet) < 0.4, the jet is retained and the electron is rejected to remove electrons originating from b-hadron decays. Since all electrons are also reconstructed as jets, if R(e, jet) < 0.2 the electron is kept and the jet is discarded. Finally, electron candidates with a R(e, μ) < 0.01 with respect to a preselected muon (defined below) are rejected and the muon is kept to suppress the contribution of electron candidates from muon bremsstrahlung and subsequent photon conversion.
Muon candidates are reconstructed by combining tracks formed in the ID and the MS sub-systems. The medium identification criteria are applied, which offer good efficiency and purity for the selected muons [70]. Preselected muons in the soft-lepton (hard-lepton) channel are required to have p T > 6(10) GeV and |η| < 2.4. Muons with an angular separation of R(μ, jet) < 0.4 with respect to the closest preselected jet are rejected, after the electron-jet overlap ambiguities are resolved. However, if the number of tracks with p T > 500 MeV associated with the jet is less than three the jet is discarded and the muon kept.
The E miss T is calculated as the magnitude of the negative vector sum of the transverse momenta of identified and calibrated muons, electrons, jets and photons, in addition to the soft-track term. The soft-track term is defined as the vectorial sum of the p T of all reconstructed tracks associated with the primary vertex that are not associated with the identified objects entering explicitly the E miss T computation [71,72]. Signal jets are required to have p T > 25 GeV and |η| < 2.8. A likelihood discriminant, the jet-vertex tagger (JVT), is used to remove the residual contamination of pile-up jets. The JVT is constructed from track-based variables that are sensitive to the vertex of origin of the jet [73]. Jets with p T < 50 GeV, |η| < 2.4 and JVT score less than 0.64 are rejected.
Signal jets containing b-hadrons are identified using the MV2c20 algorithm [74] and are hereafter referred to as btagged jets. The MV2c20 algorithm uses as input the impact parameters of all associated tracks and any reconstructed secondary vertex. The requirement chosen in the analysis provides an inclusive b-tagging efficiency of 77 % in simulated tt events, along with a rejection factor of 140 for gluon and light-quark jets and of 4.5 for charm jets [74,75].
Signal muons and electrons in the soft-lepton and hardlepton channels are subject to an additional p T < 35 GeV or p T ≥ 35 GeV requirement, respectively. Electrons must satisfy likelihood-based tight criteria which are defined in Ref.
[69]. In both channels, signal leptons must satisfy the GradientLoose [70] isolation requirements, which rely on the use of tracking-based and calorimeter-based variables and implement a set of ηand p T -dependent criteria. The efficiency for prompt leptons with transverse momentum < 40 GeV to satisfy the GradientLoose requirements is measured to be about 95 % in Z → events, progressively rising up to 99 % at 100 GeV [70].
To enforce compatibility with the primary vertex, the distance |z 0 ·sin(θ )| is required to be less than 0.5 mm for signal lepton tracks, where z 0 is the longitudinal impact parameter with respect to the primary-vertex position. Moreover, in the transverse plane the distance of closest approach of the lepton track to the proton beam line, divided by the corresponding uncertainty, must be less than three for muons and less than five for electrons.
Reconstruction, identification and isolation efficiencies in simulation, when applicable, are calibrated to data for all reconstructed objects.

Event selection
Events selected by the trigger are further required to have a reconstructed primary vertex. An event is rejected if it contains a preselected jet which fails to satisfy the quality criteria designed to suppress non-collision backgrounds and detector noise [68]. Exactly one signal lepton is required in both the soft-and the hard-lepton channels. Any event with additional preselected leptons is vetoed to suppress the dilepton tt, single-top (W t-channel) and diboson backgrounds.
A dedicated optimisation study was performed to design signal region (SR) selection criteria and to maximise the signal sensitivity. Four hard-lepton signal regions and two softlepton signal regions are defined, targeting different mass hierarchy scenarios in the simplified model. The selection criteria used to define the signal regions are summarised in Table 2 for the soft-lepton channel and in Table 3 for the hard-lepton channel.
The observables defined below are used in the event selection.
The transverse mass (m T ) of the lepton and the p miss T is defined as where φ( p T , p miss T ) is the azimuthal angle between the lepton and the missing transverse momentum. This is used in the soft-lepton 2-jet signal region and all hard-lepton signal regions to reject W +jets and semileptonic tt events.
The inclusive effective mass (m inc eff ) is the scalar sum of the p T of the signal lepton and jets and the E miss T : where the index j runs over all the signal jets in the event with p T > 30 GeV. The inclusive effective mass provides good discrimination against SM backgrounds, without being too sensitive to the details of the SUSY cascade decay chain.

The transverse momentum scalar sum (H T ) is defined as
where the index j runs over all the signal jets in the event.
The H T variable is used to define the soft-lepton 5-jet signal region, as the many energetic jets in the signal model render this variable useful to separate signal from background.
The ratio E miss T /m inc eff is used in both the soft-and the hard-lepton channels; it provides good discrimination power between signal and background with fake E miss T due to instrumental effects.
Additional suppression of background processes is based on the aplanarity variable, which is defined as A = 3 2 λ 3 , where λ 3 is the smallest eigenvalue of the normalised momentum tensor of the jets [76]. Typical measured values lie in the range 0 A < 0.3, with values near zero indicating relatively planar background-like events.
The hard-lepton 5-jet region targets scenarios with high gluino masses and lowχ 0 1 masses in models with the chargino mass mχ± 1 chosen such that the mass-ratio parameter x = 1/2. Tight requirements on m T and m inc eff are applied. For the same set of models, the hard-lepton 6-jet region is designed to provide sensitivity to scenarios where the mass separation between the gluino and the neutralino is smaller. For this reason, the requirements on m T and m inc eff are relaxed with respect to the hard-lepton 5-jet region. Two distinct hard-lepton 4-jet Jet aplanarity ->0.02 Table 3 Overview of the selection criteria for the hard-lepton signal regions. The symbol p T refers to signal leptons. The mass-ratio parameter x used in the signal region labels is defined in Sect. 3 regions are used, both designed to target models where the neutralino mass is fixed to 60 GeV, while the gluino mass and the mass-ratio x vary. The 4-jet high-x region is designed for regions of the parameter space where the W boson produced in the chargino decay is significantly boosted, leading to highp T leptons. The main characteristics of signal events in the phase-space of this model are large m T values and relatively soft jets emitted from the gluino decays. In the 4-jet low-x region, the W boson tends to be virtual while the jets from the gluino decays tend to have high p T due to the large gluino-chargino mass difference. Therefore, the m T requirement is relaxed and more stringent jet p T requirements are imposed.
The soft-lepton channels focus on models with compressed mass spectra. The soft-lepton 2-jet region provides sensitivity to scenarios characterised by a relatively heavy neutralino and a small mass separation between the gluino, the chargino and the neutralino. Due to the small mass separation, most of the decay products tend to be low p T , or soft. Thus, a highp T initial-state radiation (ISR) jet recoiling against the rest of the event is required, in order to enhance the kinematic properties of the signal and to provide separation with respect to the backgrounds. The soft-lepton 5-jet region is designed to be sensitive to the configurations in parameter space with a large mass gap between the gluino and chargino and a small separation between mχ± 1 and mχ0 1 .
As a consequence, several energetic jets from the decay of the two gluinos to the charginos are expected, while the virtual W bosons produced in the decay of the charginos result in lowp T jets and leptons.

Background estimation
The two dominant background processes in final states with one isolated lepton, multiple jets and large missing transverse momentum are tt and W +jets. The differential distributions arising from these two background processes as predicted from simulation are simultaneously normalised to the number of data events observed in dedicated control regions (CR), through the fitting procedure explained in Sect. 8. The simulation is then used to extrapolate the measured background rates to the corresponding signal region.
The control regions are designed to have high purity in the process of interest, a small contamination from the signal model and enough events to result in a small statistical uncertainty in the background prediction. Moreover, they are designed to have kinematic properties resembling as closely as possible those of the signal regions, in order to provide good estimates of the kinematics of background processes there. This procedure limits the impact of potentially large systematic uncertainties in the expected yields.
Additional sources of background events are single-top events (s-channel, t-channel and associated production with a W boson), Z +jets and diboson processes (W W , W Z, Z Z, W γ , Z γ ), and tt production in association with a W or a Z boson. Their contributions are estimated entirely using simulated event samples normalised to the most accurate theoretical cross-sections available.
The contribution from multi-jet processes with a misidentified lepton is found to be negligible once lepton isolation criteria and a stringent E miss T requirement are imposed. A data-driven matrix method, following the implementation described in Ref. [16], confirms this background is consistent with zero. This is mainly a result of the improved lepton reconstruction and identification and the higher threshold on E miss T with respect to the previous searches performed in this final state [16]. As this background is found to be negligible it is ignored in all aspects of the analysis. Figure 2 visualises the criteria that define the control regions in the soft-lepton and hard-lepton channels. Based on these, separate control regions are defined to extract the normalisation factors for tt and W +jets by requiring at least one, or no, b-tagged signal jets, respectively. The crosscontamination between these two types of control regions is accounted for in the fit. Figure 3 shows the E miss T distribution in selected softlepton and hard-lepton control regions. The normalisation of the W +jets and tt simulations are adjusted to match the observed number of data events in the control region, so that the plots illustrate the modelling of the shape of each variable's distribution. In general, good agreement between data and background simulations is found within the uncertainties in all the control regions used in the analysis.
To gain confidence in the extrapolation from control to signal regions using simulated event samples, the results of the simultaneous fit are cross-checked in validation regions which are disjoint with both the control and the signal regions. The validation regions are designed to be kinematically close to the signal regions, as shown in Fig. 2, while expecting only a small contamination from the signal in the models considered in this search. The validation regions are not used to constrain parameters in the fit, but they provide a statistically independent cross-check of the extrapolation.
This analysis uses two validation regions per signal region. In the hard-lepton channel, one of the validation regions is used to test the extrapolation to larger m T values, while the other validation region tests the extrapolation to larger aplanarity values or, in the case of the 4-jet high-x selection, to larger values in E miss T /m incl eff . In the soft-lepton channel, the validation regions are used to test the extrapolation to larger E miss T , m T or H T values.

Systematic uncertainties
Two categories of systematic uncertainties have an impact on the results presented here: uncertainties arising from experimental effects and uncertainties associated with theoretical predictions and modelling. Their effects are evaluated for all signal samples and background processes. Since the normalisation of the dominant background processes is extracted in dedicated control regions, the systematic uncertainties only affect the extrapolation to the signal regions in these cases. Among the dominant experimental systematic uncertainties are the jet energy scale (JES) and resolution (JER) and the muon momentum resolution. The jet uncertainties are derived as a function of p T and η of the jet, as well as of the pile-up conditions and the jet flavour composition of the selected jet sample. They are determined using a combination of simulated samples and studies of data, such as measurements of the jet balance in dijet, Z +jet and γ +jet events [77]. The J/ψ → + − , W ± → ± ν and Z → + − decays in data and simulation are exploited to estimate the uncertainties in lepton reconstruction, identification, momentum/energy scale and resolution and isolation criteria [69,70]. In particular, muon momentum resolution and scale calibrations are derived for simulation from a template fit that compares the invariant mass of Z → μμ and J/ψ → μμ candidates in data and simulation. The corresponding uncertainties are computed from variations of several fit parameters, following the procedure described in Ref. [78].
The simulation is reweighted to match the distribution of the average number of proton-proton interactions per bunch crossing observed in data. In the signal regions characterised by a higher jet multiplicity, the uncertainty arising from this reweighting also becomes relevant.
The systematic uncertainties related to the modelling of E miss T in the simulation are estimated by propagating the uncertainties on the energy and momentum scale of each of the objects entering the calculation, as well as the uncertainties on the soft term resolution and scale.
Different uncertainties in the theoretical modelling of the SM production processes are considered, as described in the following.
For tt, single-top and W/Z +jets samples, the uncertainties related to the choice of QCD renormalisation and factorisation scales are assessed by varying the corresponding generator parameters up and down by a factor of two around their nominal values. Uncertainties in the resummation scale and the matching scale between matrix elements and parton shower are evaluated for the W +jets samples by varying up For tt and single-top production, specific samples with an increased and decreased amount of initial-and final-state radiation are compared to the nominal sample. The relative difference in the extrapolation factors (tt) or expected rates (single top) is assigned as an uncertainty. Moreover, the uncertainty associated with the parton-shower modelling An uncertainty arising from the choice of parton level generator is estimated for tt, diboson and W/Z +jets processes. In the former case, the predictions from powheg-box are compared to aMc@NLO 3 [80]; for dibosons, Sherpa is compared to powheg-box; for W/Z +jets, Sherpa is compared to Madgraph [81].
An uncertainty of 5 % in the inclusive Z +jets cross-section is assumed [82]. Uncertainties in the inclusive single-top cross-sections are assigned as 3.7 % (s-channel, top), 4.7 % (s-channel, anti-top), 4 % (t-channel, top), 5 % (t-channel, anti-top) and 5.3 % (W t-channel) [83]. Samples using dia- 2 The Herwig++ version 2.7.1 used. 3 The aMc@NLO version 2.1.1 is used. gram subtraction and diagram removal schemes are compared for assessing the sensitivity to the treatment of interference effects between single-top and tt production at NLO.
An overall systematic uncertainty of 6 % in the inclusive cross-section is assigned to the small contribution from W W , W Z, Z Z, W γ and Z γ processes, which are estimated entirely from simulation. The uncertainty accounts for missing higher-order corrections, for the uncertainty in the value of the strong coupling constant and for the uncertainties on the PDF sets. The uncertainties associated with the resummation, factorisation and renormalisation scales are computed by varying the corresponding Sherpa parameters.
For the very small contributions of tt + W/Z /W W , an uncertainty of 30 % is assigned.
Among the main systematic uncertainties on the total background predictions in the various signal regions are the ones associated with the finite size of the MC samples, which range from 11 % in the hard-lepton 6-jet SR to 33 % in the hard-lepton 4-jet high-x SR. Moreover, the uncertainties associated with the normalisation of the tt background, ranging from 7 % in the soft-lepton 5-jet SR to 21 % in the hard-lepton 4-jet low-x SR. Further important uncertainties are the theoretical uncertainties associated with the singletop background in the hard-lepton regions, which amount to 3 % in the 4-jet high-x signal region and increase to as much as 34 % in the 4-jet low-x signal region, and the theoretical uncertainties on the W +jets background in the soft-lepton regions (up to 11 %).
The theoretical systematic uncertainty affecting the modelling of ISR can become sizeable in the simplified signal models used in this analysis, especially when the SUSY particles' mass splitting becomes small. Variations of a factor of Number of events  two in the following Madgraph and Pythia parameters are used to estimate these uncertainties: the renormalisation and factorisation scales, the initial-and the final-state radiation scales, as well as the Madgraph jet matching scale. The overall uncertainties range from about 5 % for signal models with large mass differences between the gluino, the chargino and the neutralino, to 25 % for models with very compressed mass spectra.

Statistical analysis
The final results are based on a profile likelihood method [84] using the HistFitter framework [85]. To obtain a set of background predictions that is independent of the observation in the signal regions, the fit can be configured to use only the control regions to constrain the fit parameters; this is referred to as the background-only fit. For each signal region a background-only fit is performed, based on the following inputs: • the observed number of events in each of the control regions associated with the signal region, together with the number of events expected from simulation; • the extrapolation factors, including uncertainties, from control regions to the signal region, as obtained from simulation, for the W +jets and the tt backgrounds; • the yields of the smaller backgrounds such as the single top, tt+V , Z +jets and diboson backgrounds as obtained from simulation, including uncertainties.
Using this information a likelihood is constructed for every background-only fit. It consists of a product of Poisson probability density functions for every region and of constraint terms for systematic uncertainties as described below.
Multiple parameters are included in each likelihood: two normalisation parameters describing the normalisation of the W +jets and tt backgrounds and nuisance parameters associated with the systematic uncertainties (as described in Sect. 7) or the statistical uncertainties in simulated event yields. The nuisance parameters associated with the systematic uncertainties are constrained by Gaussian functions with their widths corresponding to the size of the uncertainty, while the statistical uncertainties are constrained by Poisson functions. The parameters are correlated between the control regions and the signal region.
In the fit, the likelihood is maximised by adjusting normalisation and nuisance parameters. The normalisation scale factor of the tt background is fitted to values between 0.34 +0.28 −0.25 (4-jet high-x control regions) and 0.92 +0.14 −0.12 (5-jet soft-lepton control regions), the normalisation of the W +jets background to values between 0.72 +0.31 −0.33 (6-jet control regions) and 1.00 ± 0.04 (2-jet soft-lepton control regions). Previous analyses [16] also found normalisation factors considerably smaller than unity for these background processes in similarly extreme regions of phase space. The fit introduces correlations between the normalisation parameters associated with the tt and the W +jets backgrounds and the nuisance parameters associated with systematic uncertainties. The uncertainty in the total background estimate may thus be smaller or larger than the sum in quadrature of the individual uncertainties.

Results
The results of the background-only fit described in Sect. 8 in the validation and signal regions are shown in Fig. 4 and are further detailed for the signal regions in Table 4. Good agreement between predicted and observed event yields is seen in all validation regions. Figure 5 shows the m T , E miss T and E miss T /m incl eff distributions before applying the requirement on the plotted variable in the signal regions. The predicted background yields and the observed number of events agree in all signal regions. The largest deviation, 2.1 standard deviations, is observed in the 6-jet hard-lepton signal region. This excess arises only from the muon channel, in which 8 events are observed, while 2.5 ± 0.7 events are predicted (local significance of 2.6 standard deviations). Table 5 The columns show from left to right: the name of the respective signal region; the 95 % confidence level (CL) upper limits on the visible cross-section ( σ 95 obs ) and on the number of signal events (S 95 obs ); the 95 % CL upper limit on the number of signal events (S 95 exp ), given the expected number (and ±1σ variations on the expectation) of background events; the two-sided C L b value, i.e. the confidence level observed for the background-only hypothesis and the one-sided discovery p-value ( p(s = 0)). The discovery p values are capped to 0.5 in the case of observing less events than the fitted background estimates The electron channel shows good agreement, with 2 events observed and 1.9 ± 0.6 predicted. Model-independent upper limits and discovery p values [85] in the signal regions are calculated in a modified fit configuration with respect to the background-only fit. The only region considered in these fits is the respective signal region. Control regions are not explicitly included and thus any signal contamination in the control regions is not taken into account, thus giving conservative limits. These fits use the background estimates as derived in the background-only fits as input and allow for a non-negative signal contribution in the signal region. An additional normalisation parameter for the signal contribution is included.
Observed and expected upper limits at 95 % confidence level (CL) on the number of events signifying new phenomena beyond the SM (S 95 obs and S 95 exp , respectively) are derived based on the C L s prescription [86] and are shown in Table 5 together with the upper limits on the visible beyond the SM cross-section (σ vis , defined as the product of acceptance, selection efficiency and production cross-section). The latter is calculated by dividing the observed upper limit on the beyond-SM events by the integrated luminosity of 3.2 fb −1 . The table also gives the background-only confidence level C L b . Table 5 also shows the discovery p values, giving the probability for the background-only assumption to produce event yields greater or equal to the observed data. The C L b and p values use different definitions of test statistics in their calculation, the former with the signal-strength parameter set to one and the latter to zero.
Model-dependent limits are calculated in a modified fit configuration with respect to the background-only fit. A sig- = 60 GeV models. The red solid line corresponds to the observed limit with the red dotted lines indicating the ±1σ variation of this limit due to the effect of theoretical scale and PDF uncertainties in the signal cross-section. The dark grey dashed line indicates the expected limit with the yellow band representing the ±1σ variation of the median expected limit due to the experimental and theoretical uncertainties. The exclusion limits at 95 % CL by previous ATLAS analyses [17] are shown as the grey area nal contribution is allowed and considered in all control and signal regions, with a non-negative signal-strength normalisation parameter included. For the signal processes, uncertainties due to detector effects and theoretical modelling are considered. The signal regions are explicitly used in the fit to constrain the likelihood parameters. Figure 6 shows the combined 95 % CL exclusion limits in the simplified models with gluino production using for each model point the signal region with the best expected sensitivity. Gluino masses up to 1.6 TeV are excluded for scenarios with large mass differences between the gluino and the neutralino and x = (mχ± In the same scenario and for models with a small mass difference between the gluino and the neutralino, gluino masses up to 870 GeV are excluded. The signal regions address very different sets of models and are complementary to each other. In the case of the hard-lepton 6-jet signal region (covering the central part in the (mg, mχ0 1 ) mass plane in Fig. 6), the observed exclusion limit is considerably weaker than the expected one due to the excess seen in this region.

Conclusion
A search for gluinos in events with one isolated lepton, jets and missing transverse momentum is presented. The analysis uses 3.2 fb −1 of proton-proton collision data collected by the ATLAS experiment in 2015 at √ s = 13 TeV at the LHC. Six signal regions requiring at least two to six jets are used to cover a broad spectrum of the targeted SUSY model parameter space. While four signal regions are based on highp T lepton selections and target models with large mass differences between the supersymmetric particles, two dedicated lowp T lepton regions are designed to enhance the sensitivity to models with compressed mass spectra.
The observed data agree with the Standard Model background prediction in the signal regions. The largest deviation is a 2.1 standard deviation excess in a channel requiring a highp T lepton and six jets. For all signal regions, limits on the visible cross-section are derived in models of new physics within the kinematic requirements of this search. In addition, exclusion limits are placed on models with gluino production and subsequent decays via an intermediate chargino to the lightest neutralino. The exclusion limits of previous searches conducted in LHC Run 1 are significantly extended. Gluino masses up to 1.6 TeV are excluded for low neutralino masses ( 300 GeV) and chargino masses of ∼850 GeV.