Search for squarks and gluinos in events with hadronically decaying tau leptons, jets and missing transverse momentum in proton–proton collisions 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}=13$$\end{document}s=13 TeV recorded with the ATLAS detector

A search for supersymmetry in events with large missing transverse momentum, jets, and at least one hadronically decaying tau lepton has been performed using 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 of proton–proton collision data 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}=13{\mathrm {\ TeV}}$$\end{document}s=13TeV recorded by the ATLAS detector at the Large Hadron Collider in 2015. Two exclusive final states are considered, with either exactly one or at least two tau leptons. No excess over the Standard Model prediction is observed in the data. Results are interpreted in the context of gauge-mediated supersymmetry breaking and a simplified model of gluino pair production with tau-rich cascade decays, substantially improving on previous limits. In the GMSB model considered, supersymmetry-breaking scale (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}Λ) values below \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$92 {\mathrm {\ TeV}}$$\end{document}92TeV are excluded at the 95% confidence level, corresponding to gluino masses below \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2000 {\mathrm {\ GeV}}$$\end{document}2000GeV. For large values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tan \beta $$\end{document}tanβ, values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}Λ up to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$107 {\mathrm {\ TeV}}$$\end{document}107TeV and gluino masses up to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2300 {\mathrm {\ GeV}}$$\end{document}2300GeV are excluded. In the simplified model, gluino masses are excluded up to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1570 {\mathrm {\ GeV}}$$\end{document}1570GeV for neutralino masses around \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$100 {\mathrm {\ GeV}}$$\end{document}100GeV. Neutralino masses below \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$700 {\mathrm {\ GeV}}$$\end{document}700GeV are excluded for all gluino masses between 800 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1500 {\mathrm {\ GeV}}$$\end{document}1500GeV, while the strongest exclusion of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$750 {\mathrm {\ GeV}}$$\end{document}750GeV is achieved for gluino masses around \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1450 {\mathrm {\ GeV}}$$\end{document}1450GeV.


Introduction
Supersymmetry (SUSY) [1][2][3][4][5][6] introduces a symmetry between fermions and bosons, resulting in a SUSY partner (sparticle) for each Standard Model (SM) particle, with identical mass and quantum numbers, and a difference of half a unit of spin.Squarks ( q), gluinos (g), sleptons ( ˜ ) and sneutrinos (ν) are the superpartners of the quarks, gluons, charged leptons and neutrinos, respectively.The SUSY partners of the gauge and Higgs bosons are called gauginos and higgsinos, respectively.The charged electroweak gaugino and higgsino states mix to form charginos ( χ± i , i = 1,2), and the neutral states mix to form neutralinos ( χ0 j , j = 1,2,3,4).Finally, the gravitino ( G) is the SUSY partner of the graviton.As no supersymmetric particle has been observed, SUSY must be a broken symmetry.In general, the minimal supersymmetric Standard Model (MSSM) allows couplings which violate baryon-and lepton-number conservation, leading to, for example, a short proton lifetime.To ensure accordance with established measurements, R-parity [7][8][9][10][11] conservation is often assumed.In this scenario, sparticles are produced in pairs and decay through cascades involving SM particles and other sparticles until the lightest sparticle (LSP), which is stable, is produced.
Final states with tau leptons are of particular interest in SUSY searches, although they are experimentally challenging.Light sleptons could play a role in the co-annihilation of neutralinos in the early universe, and models with light scalar taus are consistent with dark-matter searches [12].Furthermore, should SUSY or any other physics beyond the Standard Model (BSM) be discovered, independent studies of all three lepton flavours are necessary to investigate the coupling structure of the new physics, especially with regard to lepton universality.If squarks and gluinos are within the reach of the Large Hadron Collider (LHC), their production may be among the dominant SUSY processes.This article reports on an inclusive search for squarks and gluinos produced via the strong interaction in events with jets, at least one hadronically decaying tau lepton, and large missing transverse momentum from undetected LSPs.Two distinct topologies are studied, with one tau lepton (1τ) or two or more tau leptons (2τ) in the final state.These mutually exclusive channels are optimised separately.The analysis is performed using 3.2 fb −1 of proton-proton (pp) collision data at √ s = 13 TeV recorded with the ATLAS detector at the LHC in 2015.Two SUSY scenarios are considered: a gauge-mediated SUSY-breaking (GMSB) model, and a simplified model of gluino pair production.Previous searches based on the same final state have been reported by the ATLAS [13,14] and CMS [15] collaborations.
In GMSB models [16][17][18], SUSY breaking is communicated from a hidden sector to the visible sector by a set of messenger fields that share the gauge interactions of the SM.SUSY is spontaneously broken in the messenger sector, leading to massive, non-degenerate messenger fields.Gauginos acquire their masses via one-loop diagrams involving messengers.Squarks and sleptons, which do not couple directly to the messenger sector, get their masses from two-loop diagrams involving messengers and SM gauge bosons, or messengers and gauginos.The free parameters of GMSB models are the SUSY-breaking mass scale in the messenger sector (Λ), the messenger mass scale (M mes ), the number of messenger multiplets (N 5 ) of the 5 + 5 representation of SU (5), the ratio of the two Higgs-doublet vacuum expectation values at the electroweak scale (tan β), the sign of the Higgsino mass term in the superpotential (sign(µ) = ±1), and a gravitino-mass scale factor (C grav ).Masses of gauginos and sfermions are proportional to Λ, and scale as N 5 and √ N 5 , respectively.The M mes scale is required to be larger than Λ, to avoid tachyonic messengers and charge-and colour-breaking vacua, and lower than the Planck mass to suppress flavour violation.The latter condition implies that the lightest supersymmetric particle is a very light gravitino.The C grav parameter, which results from the mechanism communicating SUSY breaking to the messengers, mainly affects the decay rate of the next-to-lightest supersymmetric particle (NLSP) into the LSP.As in previous ATLAS searches [13,14], the GMSB model is probed as a function of Λ and tan β, and the other parameters are set to M mes = 250 TeV, N 5 = 3, sign(µ) = 1 and C grav = 1.For this choice of parameters, the NLSP is the lightest scalar tau (τ 1 ) for large values of tan β, while for lower tan β values, the τ1 and the superpartners of the right-handed electron and muon (ẽ R , μR ) are almost degenerate in mass.The squark-antisquark production mechanism dominates at high values of Λ.A typical GMSB signal process is displayed in Fig. 1(a).The value of C grav corresponds to prompt decays of the NLSP.The region of small Λ and large tan β is unphysical since it leads to tachyonic states.
The other signal model studied in this analysis is a simplified model of gluino pair production [19] in an Rparity-conserving scenario.It is inspired by generic models such as the phenomenological MSSM [20,21] with dominant gluino pair production, light τ1 and a χ0 1 LSP.Gluinos are assumed to undergo a two-step cascade decay leading to tau-rich final states, as shown in Fig. 1(b).The two free parameters of the model are the masses of the gluino (m g) and the LSP (m χ0 1 ).Assumptions are made about the masses of other sparticles, namely the τ1 and ντ are mass-degenerate, and the χ0 2 and χ± 1 are also mass-degenerate, with Gluinos are assumed to decay to χ± 1 q q and χ0 2 q q with equal branching ratios, where q, q denote generic first-and second-generation quarks.Neutralinos χ0 2 are assumed to decay to ττ and ντ ν τ with equal probability, and charginos χ± 1 are assumed to decay to ντ τ and τν τ with equal probability.In the last step of the decay chain, τ and ντ are assumed to decay to τ χ0 1 and ν τ χ0 1 , respectively.All other SUSY particles are kinematically decoupled.The topology of signal events depends on the mass splitting between the gluino and the LSP.The sparticle decay widths are assumed to be small compared to sparticle masses, such that they play no role in the kinematics.

The ATLAS detector
The ATLAS experiment is described in detail in Ref. [22].It is a multi-purpose detector with a forwardbackward symmetric cylindrical geometry and a solid angle 1 coverage of nearly 4π.
The inner tracking detector (ID), covering the region |η| < 2.5, consists of a silicon pixel detector, a silicon microstrip detector and a transition radiation tracker.The innermost layer of the pixel detector, the insertable B-layer [23], was installed between Run 1 and Run 2 of the LHC.The inner detector is surrounded by a thin superconducting solenoid providing a 2 T magnetic field, and by a finely segmented lead/liquid-argon (LAr) electromagnetic calorimeter covering the region |η| < 3.2.An steel/scintillatortile hadronic calorimeter provides coverage in the central region |η| < 1.7.The end-cap and forward regions, covering the pseudorapidity range 1.5 < |η| < 4.9, are instrumented with electromagnetic and hadronic LAr calorimeters, with either steel, copper or tungsten as the absorber material.A muon spectrometer system incorporating large superconducting toroidal air-core magnets surrounds the calorimeters.Three layers of precision wire chambers provide muon tracking coverage in the range |η| < 2.7, while dedicated fast chambers are used for triggering in the region |η| < 2.4.
The trigger system, composed of two stages, was upgraded [24] before Run 2. The Level-1 trigger system, implemented with custom hardware, uses information from calorimeters and muon chambers to reduce the event rate from 40 MHz to a maximum of 100 kHz.The second stage, called the High-Level Trigger (HLT), reduces the data acquisition rate to about 1 kHz.The HLT is based on software and runs reconstruction algorithms similar to those used in the offline reconstruction.

Data and simulation samples
The data used in this analysis consist of pp collisions at a centre-of-mass energy of √ s = 13 TeV delivered by the LHC with a 25 ns bunch spacing and recorded by the ATLAS detector from August to November 2015.Data quality requirements are applied to ensure that all sub-detectors were operating normally, and that LHC beams were in stable-collision mode.The integrated luminosity of the resulting data set is 3.16 ± 0.07 fb −1 .Simulated Monte Carlo (MC) event samples are used to model both the SUSY signals and SM backgrounds, except multi-jet production, which is evaluated from data.All MC samples are generated at √ s = 13 TeV.In addition to the hard-scattering process, soft pp interactions (pile-up) are included in the simulation using the Pythia 8.186 [25] generator with the A2 [26] set of tuned parameters (tune) and MSTW2008LO [27] parton density function (PDF) set.Generated events are reweighted such that the average number of pp interactions per bunch crossing has the same distribution in data and simulation.For SM background samples, the interactions between generated particles and the detector material are simulated [28] using Geant4 [29] and a detailed description of the ATLAS detector.In the case of signal samples, a parameterised fast simulation [30] is used to describe the energy deposits in the calorimeters. 1ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point in the centre of the detector and the z-axis along the beam pipe.The x-axis points from the interaction point 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 beam pipe.The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2).
The W+jets and Z+jets processes are simulated with the Sherpa 2.1.1 [31] generator.Matrix elements (ME) are calculated for up to two partons at next-to-leading order (NLO) and up to four additional partons at leading order (LO) in perturbative QCD using the OpenLoops [32] and Comix [33] matrix element generators, respectively.The polarisation of tau leptons in W(τν)+jets and Z(ττ)+jets events is handled by the TauSpinner [34] program.The phase-space merging between the Sherpa parton shower (PS) [35] and matrix elements follows the ME+PS@NLO prescription [36].The CT10 [37] PDF set is used in conjunction with dedicated parton-shower tuning.Simulated samples are generated in bins of the transverse momentum (p T ) of the vector boson.The inclusive cross sections are normalised to a next-to-next-toleading-order (NNLO) calculation [38] in perturbative QCD based on the FEWZ program [39].
For the simulation of t t and single-top-quark production in the Wt-and s-channels, the Powheg-Box v2 [40] generator is used with the CT10 PDF set for the matrix elements calculation.Electroweak tchannel single-top-quark events are generated using the Powheg-Box v1 generator.This generator uses the four-flavour scheme for the NLO matrix element calculation together with the fixed four-flavour CT10f4 PDF set.For all top quark processes, top quark spin correlations are taken into account (for t-channel production, top quarks are decayed using MadSpin [41]).The parton shower, hadronisation, and underlying event are simulated using Pythia 6.428 [42] with the CTEQ6L1 [43] PDF set and the corresponding Perugia 2012 tune [44].Cross sections are calculated at NNLO in perturbative QCD with resummation of next-to-next-to-leading logarithmic (NNLL) soft gluon terms using the Top++ 2.0 program [45].
Diboson production is simulated using the Sherpa 2.1.1 generator with the CT10 PDF set.Processes with fully leptonic final states are calculated with up to one (4 , 2 + 2ν) or no partons (3 + 1ν) at NLO and up to three additional partons at LO. Diboson processes with one of the bosons decaying hadronically and the other leptonically are simulated with up to 1 (ZZ) or 0 (WW, WZ) partons at NLO and up to 3 additional partons at LO.The generator cross sections are used for these samples.
The simplified-model signal samples are generated using MG5_aMC@NLO v2.2.3 [46] interfaced to Pythia 8.186 with the A14 tune [47] for the modelling of the parton shower, hadronisation and underlying event.The ME calculation is performed at tree level and includes the emission of up to two additional partons.The PDF set used for the generation is NNPDF23LO [48].The ME-PS matching is done using the CKKW-L prescription [49], with a matching scale set to one quarter of the gluino mass.The GMSB signal samples are generated with the Herwig++ 2.7.1 [50] generator, with CTEQ6L1 PDFs and the UE-EE-5-CTEQ6L1 tune [51], using input files generated in the SLHA format with the SPheno v3.1.12[52] program.The parton shower evolution is performed using an algorithm described in Refs.[50,[53][54][55].Signal cross sections are calculated at NLO in the strong coupling constant, adding the resummation of soft gluon emission at next-to-leading-logarithm accuracy [56][57][58][59][60].The nominal cross section and the uncertainty are taken from an envelope of cross-section predictions using different PDF sets and factorisation and renormalisation scales, as described in Ref. [61].

Reconstruction of final-state objects
This analysis primarily requires the presence of jets, hadronically decaying tau leptons and missing transverse momentum in the final state.Jet b-tagging is used to separate the top quark background from vector bosons produced in association with jets (V+jets, where V = W, Z).Electrons and muons are vetoed in the 1τ channel, and muons are explicitly used in the 2τ channel for background modelling studies.
Primary vertices are reconstructed using inner-detector tracks with p T > 400 MeV that satisfy requirements on the number of hits in silicon tracking devices [62].Primary vertex candidates are required to have at least two associated tracks, and the candidate with the largest p 2 T is chosen as the primary vertex.
Jets are reconstructed with the anti-k t clustering algorithm [63] with a distance parameter R = 0.4.Clusters of topologically connected calorimeter cells [64] with energy above the noise threshold, calibrated at the electromagnetic energy scale, are used as input.The jet energy is calibrated using a set of global sequential calibrations [65,66].Energy from pile-up interactions is subtracted based on the jet area and the median energy density computed for each event [67].Jets are required to have p T > 20 GeV and |η| < 2.8.The discrimination between hard-interaction jets and pile-up jets is achieved by a jet-vertextagging algorithm [68].Events with jets originating from cosmic rays, beam background or detector noise are vetoed using the loose quality requirements defined in Ref. [69].Jets containing b-hadrons (b-jets) are identified using a multivariate algorithm exploiting the long lifetime, high decay multiplicity, hard fragmentation, and large mass of b-hadrons [70].The b-tagging algorithm identifies genuine b-jets with an efficiency of approximately 70% in simulated t t events.The rejection rates for c-jets, hadronically decaying tau leptons, and light-quark or gluon jets are approximately 8, 26 and 440, respectively [71].
Electron candidates are reconstructed from an isolated energy deposit in the electromagnetic calorimeter matched to an inner-detector track.They are required to have p T > 10 GeV, |η| < 2.47, and candidates reconstructed in the transition region between the barrel and end-cap calorimeters (1.37 < |η| < 1.52) are discarded.Electrons are required to satisfy a loose likelihood identification requirement [72,73] based on calorimeter shower shapes and track properties.The significance of the transverse impact parameter of the electron track is required to be less than five.
Muon candidates are reconstructed in the region |η| < 2.5 from muon spectrometer tracks matching ID tracks.Muons are required to have p T > 10 GeV and pass the medium identification requirements defined in Ref. [74], based on the number of hits in the ID and muon spectrometer, and the compatibility of the charge-to-momentum ratio measured in the two detector systems.
Hadronically decaying tau leptons are reconstructed [75] from anti-k t jets with E T ≥ 10 GeV and |η| < 2.5 calibrated with a local cluster weighting technique [76].Tau candidates are built from clusters of calorimeter cells within a cone of size ∆R ≡ √ (∆η) 2 + (∆φ) 2 = 0.2 centred on the jet axis.Tau energy scale corrections are applied to subtract the energy originating from pile-up interactions and correct for the calorimeter response.Tau leptons are required to have p T > 20 GeV, and candidates reconstructed within the transition region 1.37 < |η| < 1.52 are discarded.Tau leptons are required to have either one or three associated tracks, with a charge sum of ±1.A boosted-decision-tree discriminant is used to separate jets from tau leptons.It relies on track variables from the inner detector as well as longitudinal and transverse shower-shape variables from the calorimeters.The analysis makes use of loose and medium tau leptons, corresponding to identification working points with efficiencies of 60% and 55% for one-track tau leptons, respectively, and 50% and 40% for three-track tau leptons, respectively.Electrons mis-identified as 1track tau leptons are rejected by imposing a p T -and |η|-dependent requirement on the electron likelihood, which provides a constant efficiency of 95% for real tau leptons, with an inverse background efficiency (rejection factor) for electrons ranging from 30 to 150 depending on the |η| region.
The missing transverse momentum vector p miss T , whose magnitude is denoted by E miss T , is defined as the negative vector sum of the transverse momenta of all identified physics objects (electrons, muons, jets, and tau leptons) and an additional soft-term.Therefore, the E miss T calculation benefits from the dedicated calibration for each final-state object.The soft-term is constructed from all the tracks with p T > 400 MeV originating from the primary vertex which are not associated with any physics object.This track-based definition makes the soft-term largely insensitive to pile-up [77].
After object reconstruction, an overlap-removal procedure is applied to remove ambiguities in case the same object is reconstructed by several algorithms.Tau candidates are discarded if they are found to overlap with a light lepton (electron or muon).If an electron and a muon are reconstructed using the same inner-detector track, the electron is discarded.For overlapping jets and electrons, the electron is kept.Light leptons in the vicinity of a jet are considered to originate from secondary decays within the jet and are discarded.Finally, in case a jet is reconstructed as a tau lepton, the tau candidate is retained.
The successive steps of this procedure are summarised in Table 1.The final-state objects considered in the analysis are those surviving the overlap-removal algorithm.

Event selection
The trigger used in this analysis is the lowest-threshold missing transverse momentum trigger that was active without restrictions during the whole 2015 data-taking period.The efficiency of that trigger is measured using data collected by a set of single-jet triggers, with events containing at least two jets and one loose tau lepton.The trigger is found to have an efficiency greater than 99% when requiring E miss T > 180 GeV and a jet with p T > 120 GeV in the offline selection, which is referred to as trigger plateau conditions.
After trigger requirements, a pre-selection common to both channels is applied to ensure that only wellreconstructed events enter the analysis.Events containing no reconstructed primary vertex, mis-measured muon tracks, cosmic muon candidates, jets originating from calorimeter noise or reconstructed near inactive areas in the calorimeter, are vetoed.The presence of a second jet with p T > 20 GeV is required.To suppress the contribution from multi-jet events where large E miss T would arise from jet energy mismeasurement, a minimum angular separation in the transverse plane is imposed between either of the two leading jets and the missing transverse momentum, ∆φ(jet 1,2 , p miss T ) > 0.4.
As part of the pre-selection in the 1τ channel, events are required to contain exactly one medium tau lepton and no light lepton.The veto against electrons and muons was used in Run 1 to ensure the 1τ channel did not overlap with the eτ and µτ channels, which are not part of the present analysis.The lepton veto does not affect the expected sensitivity to the simplified model, and has been kept in the 1τ channel.Events with two or more loose tau leptons are rejected to make the 1τ and 2τ channels statistically independent.
In the 2τ channel, at least two loose tau leptons are required at pre-selection level.No veto is applied against light leptons, as this would cause a sizeable selection inefficiency for GMSB signals.
In each channel, signal regions (SRs) are defined for several signal scenarios.The following kinematic variables are found to provide discrimination between signal and background, or between backgrounds themselves: -the transverse mass of the system formed by p miss T and a lepton assumed to be massless, For events where the missing transverse momentum and the lepton originate from a W → ν decay, the m T distribution exhibits a Jacobian peak at the W boson mass.In the 1τ channel, the m τ T variable is used.In the 2τ channel, the m τ 1 T + m τ 2 T variable based on the two leading tau leptons is used, and m µ T is used for specific selections requiring a muon; -the scalar sum of the transverse momenta of all tau leptons and jets in the event, -the magnitude of the missing transverse momentum, E miss T ; -the effective mass, m eff = H T + E miss T ; -the m ττ T2 variable [78,79], also called stransverse mass, computed as where (a, b) refers to two invisible particles that are assumed to be produced with transverse momentum vectors p a,b T .In this calculation, (a, b) are assumed to be massless.The m ττ T2 distribution has a kinematic endpoint for processes where massive particles are pair-produced, each particle decaying to a tau lepton and an undetected particle.In cases where multiple tau leptons are produced in a decay chain, there is no a-priori way to select the pair leading to the desired characteristic.For events with more than two tau-lepton candidates, m ττ T2 is hence calculated using all possible tau-lepton pairs and the largest value is chosen; -the sum of the transverse masses of all jets and of the two leading tau leptons in the event, T , where m jet T is defined analogously to m T ; -the total number of jets, N jet ; -the number of b-tagged jets, N b-jet .These variables are also used for the selection of control regions (CRs) and validation regions (VRs) in the context of background modelling studies, as discussed in Section 6.
Figure 2 shows example kinematic distributions at pre-selection level.The dominant backgrounds are W(τν)+jets and t t production in the 1τ channel, with subdominant contributions from Z(νν)+jets and Z(ττ)+jets.In the 2τ channel, the pre-selection is dominated by t t, W(τν)+jets, and Z(ττ)+jets events.The multi-jet and diboson backgrounds do not contribute significantly.In the 1τ channel, three characteristic regions in the (m g, m χ0 1 ) parameter space of the simplified model are chosen as benchmark scenarios for the SR optimisation, with small (< 100 GeV), medium (500-900 GeV) and large (> 1200 GeV) mass splittings between the gluino and the LSP.In the following, the associated SRs are called Compressed, Medium-Mass and High-Mass SRs, respectively.The Compressed SR exploits topologies where a high-p T jet from initial-state radiation (ISR) recoils against the pair of gluinos.In this situation, the soft visible particles produced in the gluino decay receive a transverse Lorentz boost.While tau leptons and jets from gluino decays typically have low p T , such events have substantial E miss T since both LSPs tend to be emitted oppositely to the ISR jet in the transverse plane.In the case of large mass splitting, high-p T jets come mainly from gluino decays, and a high-p T requirement on the first two leading jets is effective in rejecting background without inducing a large inefficiency for the signal.A requirement on the transverse mass m τ T is applied in the Medium-Mass and High-Mass SRs to suppress W(τν)+jets events as well as semileptonic t t events with a tau lepton in the final state.The H T variable is also used in these two SRs, as H T increases for signal events with increasing mass splittings.No SR is defined for the GMSB model, as the expected sensitivity in the 1τ channel is significantly lower than that in the 2τ channel.
In the 2τ channel, two SRs are defined for the simplified model to cover small (< 900 GeV) and large (> 1200 GeV) mass-splitting scenarios.The Compressed SR imposes a requirement on m ττ T2 to exploit the fact that most of SM background contributions exhibit a kinematic endpoint around the W or Z boson mass, which is not the case for tau leptons produced in the cascade decay of gluinos.A requirement is also applied on m sum T to take advantage of the large E miss T and the large jet and tau lepton multiplicity that is expected for signal events.The High-Mass SR includes a requirement on H T , which is efficient for high-mass gluino signals.A requirement on m τ 1 T + m τ 2 T is also applied.In Z(ττ)+jets events, where E miss T originates from neutrinos from tau lepton decays, the trigger plateau requirement selects events with a high-p T Z boson recoiling against jets in the transverse plane.This topology leads to tau leptons with a small ∆φ separation, which results in low values of m τ 1 T + m τ 2 T given that tau neutrinos are themselves collimated with the visible decay products of tau leptons.For dileptonic t t events with two tau leptons and T > 120 GeV Tau leptons large genuine E miss T , and for W(τν)+jets events and semileptonic t t events with a high-p T jet mis-identified as a tau lepton, m τ 1 T + m τ 2 T can be larger, but even larger values are expected for high-mass gluino signals.A signal region is also defined for the GMSB model, and targets more specifically squark-antisquark production rather than gluino pair production in the region Λ 80 TeV, not excluded by Run 1 searches.Among the distinctive features that give large discrimination power to H T , decay chains are potentially longer than in the simplified model, and the almost-massless gravitino LSP leaves more phase space to other particles in the decay.Table 2 summarises the selection criteria for all the SRs of the 1τ and 2τ channels.

Background estimation
To predict the background contributions in the SRs, the normalisation of the dominant backgrounds is fitted to data in dedicated CRs.In each channel and for each SR, a simultaneous fit over the relevant CRs is performed using HistFitter [80] to extract these normalisation factors.Control regions are designed to have an enhanced contribution from a single background process, with contributions from other backgrounds as small as possible to reduce the uncertainties originating from correlations between CRs.Furthermore, CRs are defined in phase-space regions close to that of SRs, to avoid the extrapolation of the background normalisations over very different kinematic regimes.
A set of VRs is defined in intermediate phase-space regions between a SR and its associated CRs, where signal contributions are small.These VRs are not part of the fit; they are used to compare the fitted background predictions with the observed data in the vicinity of the SRs before unblinding those.

Vector-boson and top quark backgrounds
In both channels, the dominant backgrounds originate from SM processes involving the top quark or a massive vector boson and jets.These two backgrounds can be separated from each other by either requiring or vetoing the presence of a b-tagged jet in the event.In addition, a tau-lepton candidate can be either a genuine tau lepton (true tau lepton) from a W(τν) decay or a jet mis-identified as a tau lepton (fake tau lepton), which leads to two types of CRs.
In CRs targeting true tau-lepton contributions, the normalisation factor is used to absorb the theoretical uncertainties in cross-section computations, the experimental uncertainties in the integrated luminosity, and potential differences in the tau-lepton reconstruction and identification efficiencies between data and simulation.In the case of fake tau-lepton contributions, the normalisation factor combines several effects: the quark/gluon composition of jets mis-identified as tau leptons, which is process-dependent, the parton shower and hadronisation models of the generator, and the modelling in the simulation of jet shower shapes in the calorimeter, which mainly depends on the Geant4 hadronic interaction model and the modelling of the ATLAS detector.Other contributions affecting the background normalisation include the modelling of the kinematics and acceptance of background processes.These contributions are absorbed into the true and fake tau-lepton normalisation factors in the 1τ channel, while they are treated as a separate kinematic normalisation factor in the 2τ channel to avoid double-counting (true-or fake-tau normalisation factors are applied to each tau lepton, while the kinematic normalisation factor is applied once per event).
In the 1τ channel, four CRs are defined, with four associated normalisation factors.They target the top quark background (including t t and single-top-quark processes) with either a true or a fake tau lepton, and V+jets events with either a true or a fake tau lepton, respectively dominated by W(τν)+jets and Z(νν)+jets processes.The discrimination between true and fake tau-lepton contributions is achieved by a requirement on m τ T .A common set of four CRs is defined for the Medium-Mass and High-Mass SRs, due to the similarity of background compositions and event kinematics.These are separated from the SRs by an upper bound on m τ T .Another set of four CRs is defined for the Compressed SR, to study more specifically the background modelling for low-p T tau leptons.These CRs are separated from the SR by an upper bound on E miss T .The selection criteria defining these CRs are summarised in Tables 3 and 4. Figure 3 and fake tau-lepton selections.The separation of the VRs from the SRs is achieved by an inverted requirement on m τ T and E miss T , for true and fake tau-lepton VRs, respectively.The selection criteria defining all the VRs in the 1τ channel are listed in Tables 3 and 4.
In the 2τ channel, the dominant backgrounds are Z(ττ)+jets and dileptonic top quark contributions (including t t and single-top quark processes) with two true tau leptons, W(τν)+jets and semileptonic top quark contributions with one true and one fake tau lepton, and W( ν)+jets and top quark contributions with two fake tau leptons.Control regions are defined for W+jets and top quark backgrounds to extract normalisation factors related to the modelling of the process kinematics, and the modelling of real and fake tau leptons in the simulation.These CRs are separated from the SRs by replacing the requirement on the two tau leptons with requirements on different final-state objects, which stand in for true or fake tau leptons.To be independent from tau-lepton considerations, kinematic CRs are based on events with one muon, jets, E miss T , and without or with b-jets, to select W(µν)+jets and semileptonic top quark events with a final-state muon, respectively.The fake tau-lepton CRs use the same baseline selections as the kinematic CRs, but in addition, the presence of a loose tau-lepton candidate is required.Events with large m µ T values are discarded to suppress the dileptonic top quark background with a muon and a true tau lepton.The true tau-lepton CRs, which target W(τν)+jets and semileptonic top quark processes with a true tau lepton, are based on events with a loose tau lepton, jets and E miss T , without or with b-jets.Similarly, contributions from fake tau leptons are suppressed by a requirement on m τ T .A separate CR is designed to study Z(ττ)+jets events by inverting the m τ 1 T + m τ 2 T and H T requirements from the SRs.This selection requires two loose tau leptons of opposite electric charge.The selection criteria defining the various CRs are summarised in Table 5. Figure 4 illustrates the good background modelling in the CRs after the fit.
The VRs of the 2τ channel are presented in Table 6.For the Z(ττ)+jets background, the validation of the background extrapolation is performed from low-H T to high-H T , while keeping the upper bound on m τ 1 T + m τ 2 T which is effective in selecting Z(ττ)+jets events.The validity of the top quark and W(τν)+jets background predictions obtained with alternative object selections are checked for selections with two reconstructed tau leptons.High values of m τ 1 T + m τ 2 T are required to suppress Z(ττ)+jets events as in the SRs, while upper bounds on m T2 and H T ensure there is no overlap between SRs and VRs.The same set of CRs and VRs is used for the three SRs of the 2τ channel.
The resulting normalisation factors for both channels do not deviate from 1 by more than 25%, except for t t events with fake tau leptons in the 1τ channel, where the normalisation factor reaches 2 within large statistical uncertainty.The typical level of agreement between data and background distributions in the CRs after the fit can be seen in Figs. 3 and 4. Good modelling of kinematic distributions is observed in all CRs.The comparison between the number of observed events and the predicted background yields in the VRs is displayed in Fig. 5. Agreement within approximately one standard deviation is observed.T > 120 GeV Tau leptons  By construction, the total fitted background is equal to the number of observed events in each CR.The last bin includes overflow events.The shaded bands indicate the statistical uncertainties in the background predictions.

Multi-jet background
The multi-jet background contributes to the selection when two conditions are simultaneously fulfilled: jets have to be mis-identified as tau leptons, and large missing transverse momentum must arise from jet energy mis-measurement.This background is estimated from data, because final-state objects arising from mis-measurements are much more challenging to simulate than the reconstruction and identification of genuine objects.Moreover, the very large multi-jet production cross section at the LHC would imply simulating a prohibitively large number of multi-jet events.
The jet smearing method [81] employed in the 1τ channel proceeds in two steps.First, multi-jet events with well-measured jets are selected in a data sample collected by single-jet triggers.This is achieved by requiring E miss T /( E T ) 1/3 < 5 GeV 2/3 , where the objects entering the E T term are those entering T > 120 GeV Tau leptons Number of events  the E miss T calculation.Selected events are required to have at least two jets, no light lepton and exactly one tau candidate satisfying the medium identification criteria.The selection is dominated by multijet production, such that most tau candidates are jets mis-identified as tau leptons.In a second step, jet energies are smeared according to the p T -dependent jet energy resolution extracted from simulation.The smearing is performed multiple times for each event, leading to a large pseudo-data set where E miss T originates from resolution effects and which includes an adequate fraction of jets mis-identified as tau leptons.
This method cannot be used in the 2τ channel because of the limited number of events with well-measured jets that contain at least two loose tau candidates.Instead, a fake rate approach is adopted.The probability for jets to be mis-identified as tau leptons or muons, obtained from simulated dijet events, is applied to jets from an inclusive data sample collected by single-jet triggers and dominated by multi-jet events.

Systematic uncertainties
For all simulated processes, theoretical and experimental systematic uncertainties are considered.The former includes cross-section uncertainties, which are not relevant for the dominant backgrounds normalised to data, and generator modelling uncertainties.The latter refers to all the uncertainties related to the reconstruction, identification, calibration and corrections applied to jets, tau leptons, electrons, muons and missing transverse momentum.Specific uncertainties are evaluated for the multi-jet background, which is estimated from data.
Theoretical uncertainties are evaluated for all simulated samples.For backgrounds that are normalised in CRs, the uncertainty in the transfer factors, i.e. the ratio of the expected event yields in a SR or VR over the respective CR, is evaluated for all SRs and VRs.The difference between the nominal simulation and the systematically varied sample is used as an additional uncertainty.For backgrounds that are evaluated from simulation alone, i.e. the diboson background, a global normalisation uncertainty is added.Uncertainties for V+jets samples generated with Sherpa are estimated by up and down variations by factors of two and one-half in the renormalisation and factorisation scales, resummation scale (maximum scale of the additional emission to be resummed by the parton shower) and CKKW matching scale (matching between matrix elements and parton shower).The effect of scale variations is parameterised at generator level as a function of the vector boson p T and the number of particle jets.For the top quark background, the nominal predictions from Powheg-Box + Pythia6 are compared with predictions from alternative generators, and the differences are taken as systematic uncertainties.The MG5_aMC@NLO + Herwig++ generators are used to evaluate uncertainties in the modelling of the hard scattering, parton shower and hadronisation.The Powheg-Box + Herwig++ generators are used to compute a specific uncertainty in the parton shower and hadronisation models.An uncertainty in the ISR modelling is also assessed by varying the Powheg-Box parameter which controls the transverse momentum of the first additional parton emission beyond the Born configuration.In the case of the diboson background, a 6% uncertainty in the cross section due to scale and PDF uncertainties is considered.Uncertainties in signal cross sections are obtained by using different PDF sets and factorisation and renormalisation scales, as discussed in Section 3.
Systematic uncertainties affecting jets arise from the jet energy scale and jet energy resolution [66], as well as efficiency corrections for jet-vertex-tagging [68] and b-tagging [82].A set of p T -and η-dependent uncertainties in the jet energy scale and resolution is estimated by varying the conditions used in the simulation.Another set of uncertainties accounts for the modelling of the residual pile-up dependence.Additional uncertainties account for the jet flavour composition of samples that are used to derive in-situ energy scale corrections, where jets are calibrated against well-measured objects.A punch-through uncertainty for jets not entirely contained in the calorimeters, as well as a single-hadron response uncertainty, are also included for high-p T jets.An overall uncertainty in the jet energy resolution is applied to jets in the simulation as a Gaussian energy smearing.Systematic uncertainties affecting correctly identified tau leptons arise from the reconstruction, identification and tau-electron overlap-removal efficiencies, and the energy scale calibration [75].Most of the uncertainties are estimated by varying nominal parameters in the simulation: detector material, underlying event, hadronic shower model, pile-up and noise in the calorimeters.The uncertainty in the energy scale also includes non-closure of the calibration found in simulation, a single-pion response uncertainty, and an uncertainty in the in-situ energy calibration of data with respect to simulation derived in Z(ττ) events with a hadronically decaying tau lepton and a muon in the final state.In the case of signal samples, which undergo fast calorimeter simulation, a dedicated uncertainty takes into account the difference in performance between full and fast simulation.The effect of mis-identified tau leptons is largely constrained by the background estimation approaches.Uncertainties arise due to the extrapolation from the CRs to the VRs and SRs.These are considered as part of the theory uncertainties, which account for the impact of hadronisation on the mis-identification of jets as tau leptons.
Systematic uncertainties affecting electrons and muons are related to the energy or momentum calibration, as well as efficiency corrections for the reconstruction, identification and isolation requirements.These uncertainties have a negligible impact on the background predictions.
Systematic uncertainties in the missing transverse momentum originate from uncertainties in the energy or momentum calibration of jets, tau leptons, electrons, and muons, which are propagated to the E miss T calculation.Additional uncertainties are related to the calculation of the track-based soft-term.These uncertainties are derived by studying the p T balance in Z(µµ) events between the soft-term and the hardterm composed of all reconstructed objects.Soft-term uncertainties include scale uncertainties along the hard-term axis, and resolution uncertainties along and perpendicular to the hard-term axis [83].
A systematic uncertainty of the pile-up modelling is estimated by varying the distribution of the average number of interactions per bunch crossing in the simulation.The range of the variation is determined by studying the correlation in data and simulation between the average number of interactions and the number of reconstructed primary vertices.This uncertainty ranges from a few percent in the 1τ channel to about 15% in the poorly populated SRs of the 2τ channel.
The uncertainty in the integrated luminosity is ±2.1%.It is derived, following a methodology similar to that detailed in Ref. [84], from a calibration of the luminosity scale using x-y beam-separation scans performed in August 2015.This uncertainty only affects the diboson background prediction and the signal yields, as other backgrounds are normalised to the data.
The systematic uncertainty in the small multi-jet background contribution is estimated by considering alternative normalisation regions and different jet smearing parameters in the case of the 1τ channel.An uncertainty of the order of 70% is found for the 1τ channel, and a 100% uncertainty is assigned in the 2τ channel.
The influence of the main systematic uncertainties in the total background predictions in the SRs of the 1τ and 2τ channels are summarised in Table 7.The uncertainties reported in the table are derived assuming that no signal is present in the CRs.
The total uncertainties range between 13% and 72%.For all SRs, the largest uncertainties, between 7% and 60%, originate from the MC generator modelling of top quark events.Energy scale uncertainties affecting tau leptons and jets contribute significantly in all regions as well.Other uncertainties, e.g. in the b-tagging efficiency and the jet energy resolution, do not play a large role in most of the SRs.

Results
Kinematic distributions for extended SR selections of the 1τ and 2τ channels are shown in Figs. 6 and 7, respectively.These regions are defined by the same set of selection criteria as for the SRs, except that the criterion corresponding to the plotted variable is not applied.Data and background predictions, fitted to data in the control regions, are compared, and signal predictions are also shown for several benchmark models.Variables providing the most discrimination between signal and background are displayed: m τ T distributions for the 1τ channel, and m sum T , m τ 1 T + m τ 2 T and H T distributions for the 2τ channel.Reasonable agreement between data and background distributions is observed, given the low event yields remaining in data after these selections.The numbers of observed events and expected background events in all SRs are reported in Tables 8 and 9 for the 1τ and 2τ channels, respectively.No significant excess is observed in data over the SM predictions.Therefore, upper limits are set at the 95% confidence level (CL) on the number of hypothetical signal events, or equivalently, on the signal cross section.The one-sided profile-likelihood-ratio test statistic is used to assess the compatibility of the observed data with the background-only and signal-plus-background hypotheses.Systematic uncertainties are included in the likelihood function as nuisance parameters with Gaussian probability densities.Following the standards used for LHC analyses, p-values are computed according to the CL s prescription [85].
Model-independent upper limits are calculated for each SR, assuming no signal contribution in the CRs.The results are derived using profile-likelihood-ratio distributions obtained from pseudo-experiments.Upper limits on signal yields are converted into limits on the visible cross section (σ vis ) of BSM processes by dividing by the integrated luminosity of the data.The visible cross section is defined as the product of production cross section, acceptance and selection efficiency.Results are summarised at the bottom of Tables 8 and 9.The most stringent observed upper limits on the visible cross section are 1.17 fb for the High-Mass SR of the 1τ channel and 1.07 fb for the High-Mass and GMSB SRs of the 2τ channel.

ATLAS
Figure 8: Exclusion contours at the 95% confidence level for the simplified model of gluino pair production, based on results from the 1τ and 2τ channels.The red solid line and the blue dashed line correspond to the observed and median expected limits, respectively, for the combination of the two channels.The yellow band shows the onestandard-deviation spread of expected limits around the median.The effect of the signal cross-section uncertainty in the observed limits is shown as red dotted lines.Additionally, expected limits are shown for the 1τ and 2τ channels individually as dashed green and magenta lines, respectively.The previous ATLAS result [19] obtained with 20.3 fb −1 of 8 TeV data is shown as the grey filled area.
The results are interpreted in the context of the simplified model of gluino pair production and the GMSB model.In the case of model-dependent interpretations, the signal contribution in the control regions is included in the calculation of upper limits, and asymptotic properties of test-statistic distributions are used [86].Exclusion contours at the 95% CL are derived in the (m g, m χ0 1 ) plane for the simplified model and in the (Λ, tan β) plane for the GMSB model.Results are shown in Figs. 8 and 9.The solid lines and the dashed lines correspond to the observed and median expected limits, respectively.The band shows the one-standard-deviation spread of the expected limits around the median, which originates from statistical and systematic uncertainties in the background and signal.The theoretical uncertainty in the signal cross section is not included in the band.Its effect on the observed limits is shown separately as the dotted lines.For the simplified model, exclusion contours are shown for the 1τ and 2τ channels and their combination.The combination is performed by selecting, for each signal scenario, the SR with the lowest expected CL s value.The combination retains the Compressed SR of the 1τ channel in the region where the LSP mass is close to the gluino mass, and favours the High-Mass SR of the 2τ channel when the mass splitting is large.For the GMSB model, limits are shown for the GMSB SR of the 2τ channel.The stronger limits at high values of tan β are explained by the nature of the NSLP, which is the lightest scalar tau in this region.For both models, the exclusion limits obtained with 3.2 fb −1 of collision data at √ s = 13 TeV significantly improve upon the previous ATLAS results [14,19] established with 20.3 fb −1 of 8 TeV data.
[TeV] Λ  : Exclusion contours at the 95% confidence level for the gauge-mediated supersymmetry-breaking model, based on results from the 2τ channel.The red solid line and the blue dashed line correspond to the observed and median expected limits, respectively.The yellow band shows the one-standard-deviation spread of expected limits around the median.The effect of the signal cross-section uncertainty in the observed limits is shown as red dotted lines.The previous ATLAS result [14] obtained with 20.3 fb −1 of 8 TeV data is shown as the grey filled area.

Summary
A search for squarks and gluinos has been performed in events with hadronically decaying tau leptons, jets and missing transverse momentum, using 3.2 fb −1 of pp collision data at √ s = 13 TeV recorded by the ATLAS detector at the LHC in 2015.Two channels, with either one tau lepton or at least two tau leptons, are separately optimised.The numbers of observed events in the different signal regions are in agreement with the Standard Model predictions.Results are interpreted in the context of a gaugemediated supersymmetry breaking model and a simplified model of gluino pair production with taurich cascade decay.In the GMSB model, limits are set on the SUSY-breaking scale Λ as a function of tan β.Values of Λ below 92 TeV are excluded at the 95% CL, corresponding to gluino masses below 2000 GeV.A stronger exclusion is achieved for large values of tan β, where Λ and gluino mass values are excluded up to 107 TeV and 2300 GeV, respectively.In the simplified model, gluino masses are excluded up to 1570 GeV for neutralino masses around 100 GeV, neutralino masses up to 700 GeV are excluded for all gluino masses between 800 GeV and 1500 GeV, while the strongest neutralino-mass exclusion of 750 GeV is achieved for gluino masses around 1400 GeV.A dedicated signal region provides good sensitivity to scenarios with a small mass difference between the gluino and the neutralino LSP.

Figure 1 :
Figure1: Diagrams illustrating the two SUSY scenarios studied in this analysis.In the GMSB model, the scalar lepton ˜ is preferentially a scalar tau τ1 for high values of tan β.

Figure 2 :
Figure 2: Kinematic distributions at pre-selection level, for the 1τ and 2τ channels.The last bin includes overflow events.The shaded bands indicate the statistical uncertainties in the background predictions.The contribution labelled as "other" includes diboson and multi-jet events, and the V+jets processes not explicitly listed in the legend.Red arrows in the Data/SM ratio indicate bins where the corresponding entry falls outside the plotted range.
illustrates the good modelling of the background in the various CRs after the fit.Three types of VRs are used in the 1τ channel for the Medium-and High-Mass SRs, to validate the background extrapolation from low-H T to high-H T for selections based on true tau leptons, and the extrapolations along H T and m τ T for selections based on fake tau leptons.The separation of the VRs from the SRs is achieved by inverting the selections on m τ T or H T .For the Compressed SR, four VRs are used to validate the extrapolation of V+jets and top quark background predictions along E miss T , for both the true

Figure 3 :
Figure 3: Missing transverse momentum and leading-jet p T distributions in two different control regions of the 1τ channel after the fit, illustrating the overall background modelling in the CRs.By construction, the total fitted background is equal to the number of observed events in each CR.The last bin includes overflow events.The shaded bands indicate the statistical uncertainties in the background predictions.Red arrows in the Data/SM ratio indicate bins where the corresponding entry falls outside the plotted range.

Figure 4 :
Figure4: H T distribution in the top quark fake-tau CR and transverse momentum of the leading tau lepton in the Z(ττ)+jets CR of the 2τ channel after the fit, illustrating the overall background modelling in the control regions.By construction, the total fitted background is equal to the number of observed events in each CR.The last bin includes overflow events.The shaded bands indicate the statistical uncertainties in the background predictions.

Figure 5 :
Figure 5: Number of observed events, n obs , and predicted background yields after the fit, n pred , in the validation regions of the 1τ and 2τ channels.The background predictions are scaled using normalisation factors derived in the control regions.The shaded bands indicate the statistical uncertainties in the background predictions, and correspond to the σ tot uncertainties used in the lower part of the figure.

Figure 6 :
Figure 6: m τ T distributions for "extended SR selections" of the 1τ channel, for (a) the Compressed SR selection without the m τ T > 80 GeV requirement, (b) the Medium-Mass SR selection without the m τ T > 200 GeV requirement, and (c) the High-Mass SR selection without the m τ T > 200 GeV requirement.The last bin includes overflow events.The shaded bands indicate the statistical uncertainties in the background predictions.Red arrows in the Data/SM ratio indicate bins where the corresponding entry falls outside the plotted range.The signal region is indicated by the black arrow.Signal predictions are overlaid for several benchmark models, normalised to their predicted cross sections.For the simplified model, "LM" refers to a low mass splitting, or compressed scenario, with m g = 665 GeV and m χ0 1 = 585 GeV; "MM" stands for a medium mass splitting, with m g = 1145 GeV and m χ0 1 = 265 GeV; "HM" denotes a high mass splitting scenario, with m g = 1305 GeV and m χ0 1 = 105 GeV.

Figure 7 : 1 T + m τ 2 T
Figure 7: Kinematic distributions for "extended SR selections" of the 2τ channel, for (a) m sum T in the Compressed SR selection without the m sum T > 1400 GeV requirement, (b) m τ 1 T + m τ 2 T in the High-Mass SR selection without the m τ 1 T + m τ 2 T > 350 GeV requirement, and (c) H T in the GMSB SR selection without the H T > 1700 GeV requirement.The last bin includes overflow events.The shaded bands indicate the statistical uncertainties in the background predictions.The signal region is indicated by the black arrow.Signal predictions are overlaid for several benchmark models, normalised to their predicted cross sections.For the simplified model, "MM" refers to a medium mass splitting, with m g = 1145 GeV and m χ0 1 = 265 GeV; "HM" denotes a high mass splitting scenario, with m g = 1305 GeV and m χ0 1 = 105 GeV.The GMSB benchmark model corresponds to Λ = 90 TeV and tan β = 40.

Figure 9
Figure9: Exclusion contours at the 95% confidence level for the gauge-mediated supersymmetry-breaking model, based on results from the 2τ channel.The red solid line and the blue dashed line correspond to the observed and median expected limits, respectively.The yellow band shows the one-standard-deviation spread of expected limits around the median.The effect of the signal cross-section uncertainty in the observed limits is shown as red dotted lines.The previous ATLAS result[14] obtained with 20.3 fb −1 of 8 TeV data is shown as the grey filled area.

Table 1 :
Overview of the successive steps in the overlap removal algorithm.

Table 2 :
Selection criteria for the signal regions (SRs) of the 1τ channel (top) and the 2τ channel (bottom).

Table 3 :
Selection criteria defining the control and validation regions of the Compressed SR in the 1τ channel.The pre-selection criteria are also applied although not shown in the table.Unless mentioned, criteria apply to both the top quark and V+jets background regions.Symbols * indicate that criteria are only applied to regions targeting the top quark background, and symbols † denote criteria only applied to V+jets background regions.

Table 4 :
Selection criteria defining the control and validation regions of the Medium-Mass and High-Mass SRs in the 1τ channel.The pre-selection criteria are also applied although not shown in the table.Unless mentioned, criteria apply to both the top quark and V+jets background regions.Symbols * indicate that criteria are only applied to regions targeting the top quark background, and symbols † denote criteria only applied to V+jets background regions.

Table 5 :
Overview of the control regions used in the 2τ channel.

Table 6 :
Overview of the validation regions defined in the 2τ channel.

Table 7 :
Dominant systematic uncertainties in the total background predictions, in percent, for the signal regions of the 1τ (top) and 2τ (bottom) channels.The total systematic uncertainty accounts for other minor contributions not listed in this table, as well as correlations between the uncertainties.

Table 8 :
Number of observed events and predicted background yields in the three signal regions of the 1τ channel.The background prediction is scaled using normalisation factors derived in the control regions.All systematic and statistical uncertainties are included in the quoted uncertainties.The bottom part of the table shows the observed and expected model-independent upper limits at the 95% CL on the number of signal events, S 95 obs and S 95 exp , respectively, the corresponding observed upper limit on the visible cross section, σ vis 95 obs , and the CL b value, i.e. the confidence level observed for the background-only hypothesis.

Table 9 :
Number of observed events and predicted background yields in the three signal regions of the 2τ channel.The background prediction is scaled using normalisation factors derived in the control regions.All systematic and statistical uncertainties are included in the quoted uncertainties.The bottom part of the table shows the observed and expected model-independent upper limits at the 95% CL on the number of signal events, S

3
Department of Physics, University of Alberta, Edmonton AB, Canada ag Also at Section de Physique, Université de Genève, Geneva, Switzerland ah Also at Eotvos Lorand University, Budapest, Hungary ai Also at Departments of Physics & Astronomy and Chemistry, Stony Brook University, Stony Brook NY, United States of America a j Also at International School for Advanced Studies (SISSA), Trieste, Italy ak Also at Department of Physics and Astronomy, University of South Carolina, Columbia SC, United States of America al Also at School of Physics and Engineering, Sun Yat-sen University, Guangzhou, China am Also at Institute for Nuclear Research and Nuclear Energy (INRNE) of the Bulgarian Academy of Sciences, Sofia, Bulgaria an Also at Faculty of Physics, M.V.Lomonosov Moscow State University, Moscow, Russia ao Also at Institute of Physics, Academia Sinica, Taipei, Taiwan ap Also at National Research Nuclear University MEPhI, Moscow, Russia aq Also at Department of Physics, Stanford University, Stanford CA, United States of America ar Also at Institute for Particle and Nuclear Physics, Wigner Research Centre for Physics, Budapest, Hungary as Also at Flensburg University of Applied Sciences, Flensburg, Germany at Also at University of Malaya, Department of Physics, Kuala Lumpur, Malaysia au Also at CPPM, Aix-Marseille Université and CNRS/IN2P3, Marseille, France * Deceased