Search for supersymmetry in final states with two same-sign or three leptons and jets using 36 fb$^{-1}$ of $\sqrt{s}=13$ TeV $pp$ collision data with the ATLAS detector

A search for strongly produced supersymmetric particles using signatures involving multiple energetic jets and either two isolated same-sign leptons ($e$ or $\mu$), or at least three isolated leptons, is presented. The analysis relies on the identification of $b$-jets and high missing transverse momentum to achieve good sensitivity. A data sample of proton--proton collisions at $\sqrt{s}= 13$ TeV recorded with the ATLAS detector at the Large Hadron Collider in 2015 and 2016, corresponding to a total integrated luminosity of 36.1 fb$^{-1}$, is used for the search. No significant excess over the Standard Model prediction is observed. The results are interpreted in several simplified supersymmetric models featuring $R$-parity conservation or $R$-parity violation, extending the exclusion limits from previous searches. In models considering gluino pair production, gluino masses are excluded up to 1.87 TeV at 95% confidence level. When bottom squarks are pair-produced and decay to a chargino and a top quark, models with bottom squark masses below 700 GeV and light neutralinos are excluded at 95% confidence level. In addition, model-independent limits are set on a possible contribution of new phenomena to the signal region yields.


Introduction
Supersymmetry (SUSY) [1][2][3][4][5][6] is one of the best-motivated extensions of the Standard Model (SM).A general review can be found in Ref. [7].In its minimal realization (the MSSM) [8,9] it predicts a new bosonic (fermionic) partner for each fundamental SM fermion (boson), as well as an additional Higgs doublet.If R-parity [10] is conserved (RPC) the lightest supersymmetric particle (LSP) is stable and can be the lightest neutralino1 χ0 1 .In many models, the LSP can be a dark-matter candidate [11,12] and produce signatures with large missing transverse momentum.If instead R-parity is violated (RPV), the LSP decay can generate events with high jet and lepton multiplicity.Both RPC and RPV scenarios can produce the final-state signatures considered in this article.
In order to address the SM hierarchy problem with SUSY models [13][14][15][16], TeV-scale masses are required [17,18] for the partners of the gluons (gluinos g) and of the top quarks (top squarks tL and tR ), due to the large top Yukawa coupling. 2The latter also favours significant tL -tR mixing, so that the mass eigenstate t1 is lighter than all the other squarks in many scenarios [19,20].Bottom squarks ( b1 ) may also be light, being bound to top squarks by S U(2) L invariance.This leads to potentially large production cross-sections for gluino pairs (gg), top-antitop squark pairs ( t1 t * 1 ) and bottom-antibottom squark pairs ( b1 b * 1 ) at the Large Hadron Collider (LHC) [21].Production of isolated leptons may arise in the cascade decays of those superpartners to SM quarks and neutralinos χ0 1 , via intermediate neutralinos χ0 2,3,4 or charginos χ± 1,2 that in turn lead to W, Z or Higgs bosons, or to lepton superpartners (sleptons, ˜ ).Light third-generation squarks would also enhance gluino decays to top or bottom quarks relative to the generic decays involving light-flavour squarks, favouring the production of heavy-flavour quarks and, in the case of top quarks, additional isolated leptons.This article presents a search for SUSY in final states with two leptons (electrons or muons) of the same electric charge, referred to as same-sign (SS) leptons, or three leptons (3L), jets and in some cases also missing transverse momentum, whose magnitude is referred to as E miss T .Only prompt decays of SUSY particles are considered.It is an extension of an earlier search performed by the ATLAS experiment [22] with √ s = 13 TeV data [23], and uses the data collected in proton-proton (pp) collisions during 2015 and 2016.Similar searches for SUSY in this topology were also performed by the CMS experiment at √ s = 13 TeV [24][25][26].While the same-sign or three-lepton signatures are present in many scenarios of physics beyond the SM (BSM), SM processes leading to such final states have very small cross-sections.Compared to other BSM searches, analyses based on these signatures therefore allow the use of looser kinematic requirements (for example, on E miss T or on the momentum of jets and leptons), preserving sensitivity to scenarios with small mass differences between the produced gluinos/squarks and the LSP, or in which R-parity is not conserved.This sensitivity to a wide range of BSM physics processes is illustrated by the interpretation of the results in the context of twelve different SUSY simplified models [27][28][29] that may lead to same-sign or three-lepton signatures.
In Figure 1(f), the W * labels indicate largely off-shell W bosons -the mass difference between χ± 1 and χ0 1 is around 1 GeV.subsequent electroweakino-mediated decays (Figures 1(e) and 1(f)).The former is characterized by final states with bottom squark pairs decaying to t tWW χ0 1 χ0 1 .The latter, addressed here by looking at a final state with three same-sign leptons, is a model that could explain the slight excess seen in same-sign lepton signatures during Run 1 [30].Finally, a full SUSY model with low fine-tuning, the non-universal Higgs model with two extra parameters (NUHM2) [31,32], is also considered.When the soft-SUSY-breaking electroweakino mass, m 1/2 , is in the range 300-800 GeV, the model mainly involves gluino pair production with gluinos decaying predominantly to t t χ0 1 and tb χ± 1 , giving rise to final states with two same-sign leptons and E miss T .In the case of non-zero RPV couplings in the baryonic sector (λ i jk ), as proposed in scenarios with minimal flavour violation [33][34][35], gluinos and squarks may decay directly to top quarks, leading to final states with same-sign leptons [36,37] and b-quarks (Figures 1(g) and 1(h)).Although these figures illustrate decay modes mediated by non-zero λ 313 (resp.λ 321 ) couplings, the exclusion limits set for these scenarios also hold for non-zero λ 323 (resp.λ 311 or λ 322 ), as these couplings lead to experimentally indistinguishable final states.Alternatively a gluino decaying to a neutralino LSP that further decays to SM particles via a non-zero RPV coupling in the leptonic sector, λ , or in the baryonic sector λ , is also possible (Figures 1(i) and 1(j)).Lower E miss T is expected in these scenarios, as there is no stable LSP, and the E miss T originates from neutrinos produced in the χ0 1 and top quark decays.Pair production of same-sign down squarkrights 3 (Figures 1(k) and 1(l)) is also considered.In all of these scenarios, antisquarks decay into the charge-conjugate final states of those indicated for the corresponding squarks, and gluinos decay with equal probabilities into the given final state or its charge conjugate.

ATLAS detector
The ATLAS experiment [22] is a multipurpose particle detector with a forward-backward symmetric cylindrical geometry and nearly 4π coverage in solid angle. 4The interaction point is surrounded by an inner detector (ID) for tracking, a calorimeter system, and a muon spectrometer (MS).The ID provides precision tracking of charged particles with pseudorapidities |η| < 2.5 and is surrounded by a superconducting solenoid providing a 2 T axial magnetic field.It consists of silicon pixel and silicon micro-strip detectors inside a transition radiation tracker.One significant upgrade for the √ s = 13 TeV running period is the presence of the insertable B-Layer [38], an additional pixel layer close to the interaction point, which provides high-resolution hits at small radius to improve the tracking and vertexing performance.In the pseudorapidity region |η| < 2.5, high-granularity lead/liquid-argon electromagnetic sampling calorimeters are used.A steel/scintillator tile calorimeter measures hadron energies for |η| < 1.7.The endcap and forward regions, spanning 1.5 < |η| < 4.9, are instrumented with liquid-argon calorimeters for both the electromagnetic and hadronic measurements.The MS consists of three large superconducting toroids with eight coils each and a system of trigger and precision-tracking chambers, which provide triggering and tracking capabilities in the ranges |η| < 2.4 and |η| < 2.7, respectively.A two-level trigger system is used to select events [39].The first-level trigger is implemented in hardware.This is followed by the software-based high-level trigger, which can run algorithms similar to those used in the offline reconstruction software, reducing the event rate to about 1 kHz.

Data set and simulated event samples
The data used in this analysis were collected during 2015 and 2016 with a peak instantaneous luminosity of L = 1.4 × 10 34 cm −2 s −1 .The mean number of pp interactions per bunch crossing (pile-up) in the data set is 24.After the application of beam, detector and data-quality requirements, the integrated luminosity considered corresponds to 36.1 fb −1 .The uncertainty in the combined 2015+2016 integrated luminosity is 3.2%.It is derived, following a methodology similar to that detailed in Ref. [40], from a preliminary calibration of the luminosity scale using x-y beam-separation scans performed in August 2015 and May 2016. 3 These RPV baryon-number-violating couplings only apply to S U(2) singlets. 4ATLAS 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 beam pipe.The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2).Rapidity is defined as y = 0.5 ln (E + p z )/(E − p z ) where E denotes the energy and p z is the component of the momentum along the beam direction.The transverse momentum p T , the transverse energy E T and the missing transverse momentum E miss T are defined in the x-y plane.
Monte Carlo (MC) simulated event samples are used to model the SUSY signals and to estimate the irreducible SM background with two same-sign and/or three "prompt" leptons.Prompt leptons are produced directly in the hard-scattering process, or in the subsequent decays of W, Z and H bosons or prompt τ leptons.The reducible background, mainly arising from t t production, is estimated from the data as described in Section 5.1.The MC samples were processed through a detailed ATLAS detector simulation [41] based on Geant4 [42] or a fast simulation using a parameterization of the calorimeter response and Geant4 for the ID and MS [43].To simulate the effects of additional pp collisions in the same and nearby bunch crossings, inelastic interactions were generated using the soft strong-interaction processes of Pythia 8.186 [44] with a set of tuned parameters referred to as the A2 tune [45] and the MSTW2008LO parton distribution function (PDF) set [46].These MC events were overlaid onto the simulated hard-scatter event and reweighted to match the pile-up conditions observed in the data.Table 1 presents, for all samples, the event generator, parton shower, cross-section normalization, PDF set and the set of tuned parameters for the modelling of the parton shower, hadronization and underlying event.In all MC samples, except those produced by the Sherpa event generator, the EvtGen v1.2.0 program [47] was used to model the properties of bottom and charm hadron decays.The SUSY signals from Figure 1 are defined by an effective Lagrangian describing the interactions of a small number of new particles [27][28][29].All SUSY particles not included in the decay of the pair-produced squarks and gluinos are effectively decoupled.These simplified models assume one production process and one decay channel with a 100% branching fraction.Apart from Figure 1(j), where events were generated with Herwig++ [51], all simplified models were generated from leading-order (LO) matrix elements with up to two extra partons in the matrix element (only up to one for the g → q q( /νν) χ0 1 model) using MG5_aMC@NLO 2.2.3 [48] interfaced to Pythia 8 with the A14 tune [50] for the modelling of the parton shower, hadronization and underlying event.Jet-parton matching was realized following the CKKW-L prescription [65], with a matching scale set to one quarter of the pair-produced superpartner mass.All signal models were generated with prompt decays of the SUSY particles.Signal cross-sections were calculated at next-to-leading order (NLO) in the strong coupling constant, adding the resummation of soft-gluon emission at next-to-leading-logarithmic accuracy (NLO+NLL) [52][53][54][55][56], except for the RPV models of Figures 1(k) and 1(l) and the NUHM2 model where NLO cross-sections were used [52,66].The nominal cross-sections and the uncertainties were taken from envelopes of cross-section predictions using different PDF sets and factorization and renormalization scales, as described in Refs.[21,57].Typical pair-production cross-sections are: 4.7 ± 1.2 fb for gluinos with a mass of 1.7 TeV, 28 ± 4 fb for bottom squarks with a mass of 800 GeV, and 15.0 ± 2.0 fb for down squark-rights with a mass of 800 GeV and a gluino mass of 2.0 TeV.
The two dominant irreducible background processes are t tV (with V being a W or Z/γ * boson) and diboson production with final states of four charged leptons , 5 three charged leptons and one neutrino, or two same-sign charged leptons and two neutrinos.The MC simulation samples for these are described in Refs.[67] and [62], respectively.For diboson production, the matrix elements contain the doubly resonant diboson processes and all other diagrams with four or six electroweak vertices, such as W ± W ± j j, with one (W ± W ± j j) or two (WZ, ZZ) extra partons.NLO cross-sections for t tW, t tZ/γ * (→ ) 6 and leptonic diboson processes are respectively 0.60 pb [60], 0.12 pb and 6.0 pb [62].The processes t tH and 4t, with NLO cross-sections of 507.1 fb [60] and 9.2 fb [48] respectively, are also considered.
Other background processes, with small cross-sections and no significant contribution to any of the signal regions, are grouped into a category labelled "rare".This category contains t tWW and t tWZ events generated with no extra parton in the matrix element, and tZ, tWZ, tt t, WH and ZH as well as triboson (WWW, WWZ, WZZ and ZZZ) production with fully leptonic decays, leading to up to six charged leptons.The processes WWW, WZZ and ZZZ were generated at NLO with additional LO matrix elements for up to two extra partons, while WWZ was generated at LO with up to two extra partons.

Event reconstruction and selection
Candidate events are required to have a reconstructed vertex [69] with at least two associated tracks with p T > 400 MeV.The vertex with the largest Σp 2 T of the associated tracks is chosen as the primary vertex of the event.
For the data-driven background estimations, two categories of electrons and muons are used: "candidate" and "signal" with the latter being a subset of the "candidate" leptons satisfying tighter selection criteria.Electron candidates are reconstructed from energy depositions in the electromagnetic calorimeter which were matched to an ID track and are required to have |η| < 2.47, p T > 10 GeV, and pass the "Loose" likelihood-based identification requirement [70].Candidates within the transition region between the barrel and endcap electromagnetic calorimeters, 1.37 < |η| < 1.52, are not considered.The track matched with the electron must have a significance of the transverse impact parameter d 0 with respect to the reconstructed primary vertex of |d 0 |/σ(d 0 ) < 5. Muon candidates are reconstructed in the region |η| < 2.5 from muon spectrometer tracks matching ID tracks.All muon candidates must have p T > 10 GeV and must pass the "Medium" identification requirements [71].
Jets are reconstructed with the anti-k t algorithm [72] with radius parameter R = 0.4, using three-dimensional topological energy clusters in the calorimeter [73] as input.All jets must have p T > 20 GeV and |η| < 2.8.For all jets the expected average energy contribution from pile-up is subtracted according to the jet area [74,75].Jets are then calibrated as described in Ref. [75].In order to reduce the effects of pileup, a significant fraction of the tracks in jets with p T < 60 GeV and |η| < 2.4 must originate from the primary vertex, as defined by the jet vertex tagger (JVT) [76].
Identification of jets containing b-hadrons (b-tagging) is performed with the MV2c10 algorithm, a multivariate discriminant making use of track impact parameters and reconstructed secondary vertices [77,78].A requirement is chosen corresponding to a 70% average efficiency for tagging b-jets in simulated t t events.The rejection factors for light-quark/gluon jets, c-quark jets and τ → ν + hadron decays in simulated t t events are approximately 380, 12 and 54, respectively [78,79].Jets with |η| < 2.5 which satisfy the b-tagging and JVT requirements are identified as b-jets.Correction factors and uncertainties determined from data for the b-tagging efficiencies and mis-tag rates are applied to the simulated samples [78].
After the object identification, overlaps between the different objects are resolved.Any jet within a distance ∆R y ≡ (∆y) 2 + (∆φ) 2 = 0.2 of a lepton candidate is discarded, unless the jet is b-tagged,7 in which case the lepton is discarded since it probably originated from a semileptonic b-hadron decay.Any remaining lepton within ∆R y = min{0.4,0.1 + 9.6 GeV/p T ( )} of a jet is discarded.In the case of muons, the muon is retained and the jet is discarded if the jet has fewer than three associated tracks.This reduces inefficiencies for high-energy muons undergoing significant energy loss in the calorimeter.
Signal electrons must satisfy the "Medium" likelihood-based identification requirement [70].In regions with large amounts of material in the tracker, an electron (positron) is more likely to emit a hard bremsstrahlung photon; if the photon subsequently converts to an asymmetric electron-positron pair, and the positron (electron) has high momentum and is reconstructed, the lepton charge can be misidentified (later referred to as "charge-flip").To reduce the impact of charge misidentification, signal electrons must satisfy |η| < 2.0.Furthermore, signal electrons that are likely to be reconstructed with an incorrect charge assignment are rejected using the electron cluster and track properties including the impact parameter, the curvature significance, the cluster width, and the quality of the matching between the cluster and its associated track, in terms of both energy and position.These variables, as well as the electron p T and η, are combined into a single classifier using a boosted decision tree (BDT) algorithm.A selection requirement on the BDT output is chosen to achieve a rejection factor of 7-8 for electrons with a wrong charge assignment while selecting correctly measured electrons with an efficiency of 97%.Correction factors to account for differences in the selection efficiency between data and MC simulation are applied to the selected electrons in MC simulation.These correction factors are determined using Z → ee events [80].
Signal muons must fulfil the requirement |d 0 |/σ(d 0 ) < 3. Tracks associated with the signal electrons or muons must have a longitudinal impact parameter z 0 with respect to the reconstructed primary vertex satisfying |z 0 sin θ| < 0.5 mm. Isolation requirements are applied to both the signal electrons and muons.The scalar sum of the p T of tracks within a variable-size cone around the lepton, excluding its own track, must be less than 6% of the lepton p T .
The track isolation cone size for electrons (muons) ∆R η ≡ (∆η) 2 + (∆φ) 2 is given by the smaller of ∆R η = 10 GeV/p T and ∆R η = 0.2 (0.3).In addition, in the case of electrons the calorimeter energy clusters in a cone of ∆R η = 0.2 around the electron (excluding the deposit from the electron itself) must be less than 6% of the electron p T .Simulated events are corrected to account for differences in the lepton trigger, reconstruction, identification and isolation efficiencies between data and MC simulation.
The missing transverse momentum is defined as the negative vector sum of the transverse momenta of all identified candidate objects (electrons, photons [81], muons and jets) and an additional soft term.The soft term is constructed from all tracks associated with the primary vertex but not with any physics object.In this way, the E miss T is adjusted for the best calibration of the jets and the other identified physics objects listed above, while maintaining approximate pile-up independence in the soft term [82,83].
Events are selected using a combination of dilepton and E miss T triggers, the latter being used only for events with E miss T > 250 GeV.The trigger-level requirements on E miss T and the leading and subleading lepton p T are looser than those applied offline to ensure that trigger efficiencies are constant in the relevant phase space.The event selection requires at least two signal leptons with p T > 20 GeV (apart from two signal regions where the lower bound on the subleading lepton p T is 10 GeV). 8If the event contains exactly two signal leptons, they must have the same electric charge.In order to reject detector noise and noncollision backgrounds (including those from cosmic rays, beam-gas and beam-halo interactions), events are discarded if they contain any jet not satisfying basic quality criteria [84,85].
To maximize the sensitivity to the signal models of Figure 1, 19 non-exclusive 9 signal regions (SRs) are defined in Table 2.The SRs are named in the form SNLMbX, where S indicates if the signal region is targeting an RPC or RPV model, N indicates the number of leptons required, M the number of b-jets required, and X indicates the severity of the E miss T or m eff requirements (Soft, Medium or Hard).All signal regions, except Rpv2L0b, allow any number of additional leptons in addition to a e ± e ± , e ± µ ± or µ ± µ ± pair.Signal regions with a three lepton selection can either require any lepton charge combination (Rpc3L0bH, Rpc3L0bS) or that all three leptons have the same charge (Rpc3LSS1b).The other requirements used to define the SRs are the number of signal leptons (N signal leptons ), number of b-jets with p T > 20 GeV (N b-jets ), number of jets with p T above 25, 40 or 50 GeV, regardless of their flavour (N jets ), E miss T , the effective mass (m eff ) and the charge of the signal leptons.The m eff variable is defined as the scalar sum of the p T of the signal leptons, jets and the E miss T .For SRs where the Z+jets background is important (Rpc3LSS1b, Rpv2L0b and Rpv2L2bH), events in which the invariant mass of two same-sign electrons is close to the Z boson mass are vetoed.For SRs targeting the production of down squark pairs (Rpv2L1bS, Rpv2L2bS, Rpv2L1bM), only events with at least two negatively charged leptons are considered, as the down squarks decay exclusively to top antiquarks.Finally, SRs targeting signal scenarios with lepton p T spectra softer than typical background processes impose an upper bound on the leptons' p T .The last column of Table 2 indicates the targeted signal model.The Rpc3L1b and Rpc3L1bH SRs are not motivated by a particular signal model and can be seen as a natural extension of the Rpc3L0b SRs with the same kinematic selections but requiring at least one b-jet.
The values of acceptance times efficiency of the SR selections for the RPC SUSY signal models, with masses near the exclusion limit, typically range between 0.5% and 7% for models with a light χ0 1 and between 0.5 and 2% for models with a heavy χ0 1 .For RPV SUSY signal models, these values are in the range 0.2-4%.To increase the signal efficiency for the SUSY models with low-energy leptons (Figure 1(b)), the p T threshold of leptons is relaxed from 20 GeV to 10 GeV in the SR definition.

Background estimation
Two main sources of SM background can be distinguished in this analysis.The first category is the reducible background, which includes events containing electrons with mismeasured charge, mainly from Signal region  the production of top quark pairs, and events containing at least one fake or non-prompt (FNP) lepton.The FNP lepton mainly originates from heavy-flavour hadron decays in events containing top quarks, or W or Z bosons.Hadrons misidentified as leptons, electrons from photon conversions and leptons from pion or kaon decays in flight are other possible sources.Data-driven methods used for the estimation of this reducible background in the signal and validation regions are described in Section 5.1.
The second background category is the irreducible background from events with two same-sign prompt leptons or at least three prompt leptons and is estimated using the MC simulation samples.Since diboson and t tV events are the main irreducible backgrounds in the signal regions, dedicated validation regions (VR) with an enhanced contribution from these processes, and small signal contamination, are defined to verify the background predictions from the simulation (Section 5.2).Section 5.3 discusses the systematic uncertainties considered when performing the background estimation in the signal and validation regions.

Reducible background estimation methods
Charge misidentification is only relevant for electrons.The contribution of charge-flip events to the SR/VR is estimated using the data.The electron charge-flip probability is extracted in a Z/γ * → ee data sample using a likelihood fit which takes as input the numbers of same-sign and opposite-sign electron pairs observed in a window of 10 GeV around the Z boson mass.The charge-flip probability is a free parameter of the fit and is extracted as a function of the electron p T and η.These probabilities are around 0.5% (1%) and 0.1% (0.2%) for the candidate and signal electrons for |η| < 1.37 (|η| > 1.52), respectively.The former is used only in the FNP lepton background estimation.The event yield of the charge-flip electron background in the signal or validation regions is obtained by multiplying the measured chargeflip probability with the number of events in data regions with the same kinematic requirements as the signal or validation regions but with opposite-sign lepton pairs.Two data-driven methods are used to estimate the FNP lepton background, referred to as the "matrix method" and the "MC template method".The estimates from these methods are combined to give the final estimate.These two methods are described below.
The first estimation of the FNP lepton background is performed with a matrix method similar to that described in Ref. [86].Two types of lepton identification criteria are defined: "tight", corresponding to the signal lepton criteria described in Section 4, and "loose", corresponding to candidate leptons after object overlap removal and the charge-flip BDT selection described also in Section 4. The matrix method relates the number of events containing prompt or FNP leptons to the number of observed events with tight or loose-not-tight leptons using the probability for loose prompt or FNP leptons to satisfy the tight criteria.The probability for loose prompt leptons to satisfy the tight selection criteria (ε) is obtained using a Z/γ * → data sample and is modelled as a function of the lepton p T and η.The efficiencies for electrons (muons) rise from 60% (80%) at low p T to almost 100% at p T above 50 GeV -apart from endcap electrons, for which they reach only 95%.The probability for loose FNP leptons to satisfy the tight selection criteria (FNP lepton rate, f ) is determined from data in SS control regions enriched in non-prompt leptons mostly originating from heavy-flavour hadron decays in single-lepton t t events.These regions contain events with at least one b-jet, one well-isolated muon (referred to as the "tag"), and an additional loose electron or muon which is used for the measurement.The rates f are measured as a function of p T after subtracting the small contribution from prompt-lepton processes predicted by simulation and the data-driven estimation of events with electron charge-flip. 10For electrons, and muons with |η| < 2.3, f is constant at around 10% for p T < 30 GeV (20% for muons with |η| > 2.3) and increases at higher p T .With these values of ε and f , the method has been demonstrated to correctly estimate the FNP lepton background.
The second method for FNP lepton estimation is the MC template method described in details in Refs.[86,87].It relies on the correct modelling of the kinematic distributions of the FNP leptons and chargeflipped electron processes in t t and V+jets samples.These samples were simulated with the Powheg-Box generator [88][89][90][91] and the parton shower and hadronization performed by either Pythia 6.428 [92] (t t) or Pythia 8.186 (V+jets).The FNP leptons are classified in five categories, namely electrons and muons originating from b-and light-quark jets as well as electrons from photon conversions.Normalization factors for each of the five sources are adjusted to match the observed data in dedicated control regions.Events are selected with at least two same-sign signal leptons, E miss T > 40 GeV, two or more jets, and are required not to belong to the SRs.They are further split into regions with or without b-jets and with different lepton flavours of the same-sign lepton pair, giving a total of six control regions.The global normalization factors applied to the MC samples for estimating the reducible background in each SR vary from 1.2 ± 1.1 to 2.9 ± 2.0, where the errors account for statistical uncertainties and uncertainties related to the choice of event generator (see Section 5.3).
Since the FNP lepton predictions from the MC template and matrix methods in the signal and validation regions are consistent with each other, a weighted average of the two results is used.With this approach, the combined estimate is always dominated by systematic uncertainties, which is not always the case when only the matrix method is used due to small number of events in the control regions.To check the validity and robustness of the FNP lepton estimate, the distributions of several discriminating variables in data are compared with the predicted background after various requirements on the number of jets and b-jets.Examples of such distributions are shown in Figure 2, and illustrate that the data are described by the prediction within uncertainties.The apparent disagreement for m eff above 1 TeV in Figure 2(d) is covered by the large theory uncertainty for the diboson background, which is not shown but amounts to about 30% for m eff above 1 TeV.The uncertainty bands include the statistical uncertainties for the background prediction as well as the systematic uncertainties for fake-or non-prompt-lepton backgrounds (using the matrix method) and charge-flip electrons.Not included are theoretical uncertainties in the irreducible background contributions.The rare category is defined in the text.

Validation of irreducible background estimates
Dedicated validation regions are defined to verify the estimate of the t tV, WZ and W ± W ± background in the signal regions.The corresponding selections are summarized in Table 3.The overlap with the signal regions is resolved by removing events that are selected in the signal regions.The purity of the targeted

All VRs
Veto events belonging to any SR ).The two leading-p T leptons are referred to as 1,2 with decreasing p T .Additional requirements are set on E miss T , m eff , the invariant mass of the two leading electrons m e ± e ± , the presence of SS leptons or a pair of same-flavour opposite-sign leptons (SFOS) and its invariant mass m SFOS .A minimum angular separation between the leptons and the jets (∆R η ( 1,2 , j)) and between the two leptons (∆R η ( 1 , 2 )) is imposed in the W ± W ± j j VR.For the two WZ VRs the selection also relies on the ratio of the E miss T in the event to the sum of p T of all signal leptons p T (E miss T / p T ).The ratio of the scalar sum of the p T of all b-jets to that of all jets in the event ( p b-jet background processes in these regions ranges from 35% to 65%.The expected signal contamination is generally below 5% for models near the limit of exclusion in t tZ, WZ and W ± W ± VRs and about 20% in the t tW VR.The observed yields, compared with the background predictions and uncertainties, are shown in Table 4.There is good agreement between data and the estimated background in all the validation regions.

Systematic uncertainties
Statistical uncertainties due to the number of data events in the loose and tight lepton control regions are considered in the FNP lepton background estimate.In the matrix method, the systematic uncertainties mainly come from potentially different compositions of b-jets, light-quark jets and photon conversions between the signal regions and the regions where the FNP lepton rates are measured.The uncertainty coming from the prompt-lepton contamination in the FNP lepton control regions is also considered.Overall, the uncertainty in the FNP lepton rate f amounts to 30% at low p T , and can reach 85% for muons with p T > 40 GeV, and 50% for electrons with p T > 20 GeV; these values are driven respectively by the dependency of the isolation of non-prompt muons on the kinematic properties of the jets which emit them, and the uncertainty in the proportion of non-prompt electrons from heavy-flavoured hadron decays with respect to other sources of FNP electrons (mainly converted photons).The uncertainties in the promptlepton efficiency ε are much smaller.The uncertainties in the FNP lepton background estimated with the matrix method in each VR and SR are then evaluated by propagating the f and ε uncertainties.In the MC template method, the systematic uncertainty is obtained by changing the generator from Powheg-Box to Sherpa and propagating uncertainties from the control region fit to the global normalization scale factors applied to the MC samples.The uncertainties in these scale factors are in the range 75-80%, depending Table 4: The numbers of observed data and expected background events in the validation regions.The rare category is defined in the text.Background categories with yields shown as "-" do not contribute to a given region (e.g.charge flips in three-lepton regions) or their estimates are below 0.01 events.The displayed yields include all statistical and systematic uncertainties described in Section 5.3.
on the SRs.When combining the results of the MC template method and the matrix method to obtain the final estimate, systematic uncertainties are propagated assuming conservatively a full correlation between the two methods.
The uncertainty in the electron charge-flip probability mainly originates from the number of events in the regions used in the charge-flip probability measurement and the uncertainty related to the background subtraction from the Z boson's mass peak.The relative error in the charge-flip rate is below 20% (30%) for signal (candidate) electrons with p T above 20 GeV.
The systematic uncertainties related to the estimated background from same-sign prompt leptons arise from the experimental uncertainties (jet energy scale calibration, jet energy resolution and b-tagging efficiency) as well as theoretical modelling and theoretical cross-section uncertainties.The statistical uncertainty of the simulated event samples is also taken into account.
The cross-sections used to normalize the MC samples are varied according to the uncertainty in the crosssection calculation, which is 13% for t tW, 12% for t tZ production [60], 6% for diboson production [62], 8% for t tH [60] and 30% for 4t [48].Additional uncertainties are assigned to some of these backgrounds to account for the theoretical modelling of the kinematic distributions in the MC simulation.For t tW and t tZ, the predictions from the MG5_aMC@NLO and Sherpa generators are compared, and the renormalization and factorization scales used to generate these samples are varied independently within a factor of two, leading to a 15-35% uncertainty in the expected SR yields for these processes.For diboson production, uncertainties are estimated by varying the QCD and matching scales, as well as the parton shower recoil scheme, leading to a 30-40% uncertainty for these processes after the SR selections.For t tH, 4t and rare production processes, a 50% uncertainty in their total contribution is assigned.

(b)
Figure 3: Comparison of (a) the observed and expected event yields in each signal region and (b) the relative uncertainties in the total background yield estimate.For the latter, "statistical uncertainty" corresponds to reducible and irreducible background statistical uncertainties.The background predictions correspond to those presented in Table 5 and the rare category is explained in the text.

Results and interpretation
Figure 3(a) shows the event yields for data and the expected background contributions in all signal regions.Detailed information about the yields can be found in Table 5.In all 19 SRs the number of observed data events is consistent with the expected background within the uncertainties.The contributions listed in the rare category are dominated by triboson, tWZ and t tWW production 11 : the triboson processes generally dominate in the SRs with no b-jets, while tWZ and t tWW dominate in the SRs with one and two b-jets, respectively.
Figure 3(b) summarizes the contributions from the different sources of systematic uncertainty to the total SM background predictions in the signal regions.The uncertainties amount to 25-50% of the total background depending on the signal region, dominated by systematic uncertainties coming from the reducible background or the theory.
In the absence of any significant deviation from the SM predictions, upper limits on possible BSM contributions to the signal regions are derived, as well as exclusion limits on the masses of SUSY particles in the benchmark scenarios of Figure 1.The HistFitter framework [93], which utilizes a profile-likelihood-ratio test [94], is used to establish 95% confidence intervals using the CL s prescription [95].The likelihood is built as the product of a Poisson probability density function describing the observed number of events in the signal region and, to constrain the nuisance parameters associated with the systematic uncertainties, Gaussian distributions whose widths correspond to the sizes of these uncertainties; Poisson distributions are used instead for MC simulation statistical uncertainties.Correlations of a given nuisance parameter between the backgrounds and the signal are taken into account when relevant.The hypothesis tests are performed for each of the signal regions independently.
Table 5 presents 95% confidence level (CL) observed (expected) model-independent upper limits on the number of BSM events, S 95  obs (S 95 exp ), that may contribute to the signal regions.Normalizing these by the integrated luminosity L of the data sample, they can be interpreted as upper limits on the visible BSM cross-section (σ vis ), defined as σ vis = σ prod × A × = S 95  obs /L, where σ prod is the production crosssection, A the acceptance and the reconstruction efficiency.The largest deviation of the data from the background prediction corresponds to an excess of 1.5 standard deviations in the Rpv2L1bM SR.
Exclusion limits at 95% CL are also set on the masses of the superpartners involved in the SUSY benchmark scenarios considered.Apart from the NUHM2 model, simplified models are used, corresponding to a single production mode and with 100% branching ratio to a specific decay chain, with the masses of the SUSY particles not involved in the process set to very high values.Figures 4, 5 and 6 show the exclusion limits in all the models considered in Figure 1 and the NUHM2 model.The assumptions about the decay chain considered for the different SUSY particles are stated above each figure.For each region of the signal parameter space, the SR with the best expected sensitivity is chosen.
For the RPC models, the limits set are compared with the existing limits set by other ATLAS SUSY searches [23,96].For the models shown in Figure 4, the mass limits on gluinos and bottom squarks are up to 400 GeV higher than the previous limits, reflecting the improvements in the signal region definitions as well as the increase in integrated luminosity.Gluinos with masses up to 1.75 TeV are excluded in scenarios with a light χ0 1 in Figure 4(a).This limit is extended to 1.87 TeV when χ0 2 and slepton masses are in-between the gluino and the χ0 1 masses (Figure 4(c)).More generally, gluino masses below 1.57TeV and bottom squarks with masses below 700 GeV are excluded in models with a massless LSP.The "compressed" regions, where SUSY particle masses are close to each other, are also better covered and LSP masses up to 1200 and 250 GeV are excluded in the gluino and bottom squark pairproduction models, respectively.Of particular interest is the observed exclusion of models producing gluino pairs with an off-shell top quark in the decay (Figure 1(b)), see Figure 4(a).In this case, models are excluded for mass differences between the gluino and neutralino of 205 GeV (only 35 GeV larger than the minimum mass difference for decays into two on-shell W bosons and two b-quarks) for a gluino mass below 0.9 TeV.The Rpc3LSS1b SR allows the exclusion of top squarks with masses below 700 GeV Table 5: Numbers of events observed in the signal regions compared with the expected backgrounds.The rare category is defined in the text.Background categories with yields shown as a "-" do not contribute to a given region (e.g.charge flips in three-lepton regions) or their estimates are below 0.01.The 95% confidence level (CL) upper limits are shown on the observed and expected numbers of BSM events, S 95 obs and S 95 exp (as well as the ±1σ excursions from the expected limit), respectively.The 95% CL upper limits on the visible cross-section (σ vis ) are also given.Finally, the p-values (p 0 ) give the probabilities to observe a deviation from the predicted background at least as large as that in the data.The number of equivalent Gaussian standard deviations (Z) is also shown when p 0 < 0.5.
[GeV] g m 800 1000 1200 1400 1600 1800 [GeV]      Figure 4: Observed and expected exclusion limits on the g, b1 , t1 and χ0 1 masses in the context of RPC SUSY scenarios with simplified mass spectra.The signal regions used to obtain the limits are specified in the subtitle of each scenario.All limits are computed at 95% CL.The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty.The contours of the band around the expected limit are the ±1σ results (±2σ is also considered in Figure (e), including all uncertainties except the theoretical uncertainties in the signal cross-section.In Figures (a)-(d), the diagonal line indicates the kinematic limit for the decays in each specified scenario and results are compared with the observed limits obtained by previous ATLAS searches [23,96].Figure 5: Observed and expected exclusion limits on the g, t1 , dR and χ0 1 masses in the context of RPV SUSY scenarios with simplified mass spectra featuring gg or dR dR pair production with exclusive decay modes.The signal regions used to obtain the limits are specified in the subtitle of each scenario.All limits are computed at 95% CL.The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty.The contours of the band around the expected limit are the ±1σ results, including all uncertainties except theoretical uncertainties in the signal cross-section (±2σ is also considered in Figures 5(e [31,32].The signal region Rpc2L2bH is used to obtain the limits.The contours of the green (yellow) band around the expected limit are the ±1σ (±2σ) results, including all uncertainties.The limits are computed at 95% CL.Finally, in the NUHM2 model with low fine-tuning, values of the parameter m 1/2 below 615 GeV are excluded, corresponding to gluino masses below 1500 GeV (Figure 6).

Conclusion
A search for supersymmetry in events with two same-sign leptons or at least three leptons, multiple jets, bjets and large E miss T and/or large m eff is presented.The analysis is performed with proton-proton collision data at √ s = 13 TeV collected in 2015 and 2016 with the ATLAS detector at the Large Hadron Collider corresponding to an integrated luminosity of 36.1 fb −1 .With no significant excess over the Standard Model prediction observed, results are interpreted in the framework of simplified models featuring gluino and squark production in R-parity-conserving and R-parity-violating scenarios.Lower limits on particle masses are derived at 95% confidence level.In the gg simplified RPC models considered, gluinos with masses up to 1.87 TeV are excluded in scenarios with a light χ0 1 .RPC models with bottom squark masses below 700 GeV are also excluded in a b1 b * 1 simplified model with b1 → tW − χ0 1 and a light χ0 1 .In RPV scenarios, masses of down squark-rights are probed up to m dR ≈ 500 GeV.All models with gluino masses below 1.3 TeV are excluded, greatly extending the previous exclusion limits obtained within this search.Model-independent limits on the cross-section of a possible signal contribution to the signal regions are set.

and a χ0 2 (
(a)) or off-shell (Figure 1(b)) top quarks, as well as on-shell light quarks.The latter are accompanied by a cascade decay involving a χ± 1 Figure 1(c)) or a χ0 2 and light sleptons (Figure 1(d)).The other two RPC scenarios target the direct production of third-generation squark pairs with
signal leptons ), the number of b-jets with p T > 20 GeV (N b-jets ), the number of jets (N jets ) above a certain p T threshold (p jet T ), E miss T , m eff and/or E miss T /m eff .The last column indicates the targeted signal model.The Rpc3L1b and Rpc3L1bH SRs are not motivated by a particular signal model and can be seen as a natural extension of the Rpc3L0b SRs with the same kinematic selections but requiring at least one b-jet.

Figure 2 :
Figure 2: Distributions of (a) the number of jets, (b) the number of b-tagged jets and (c), (d) the effective mass.The distributions are made after requiring at least two jets (p T > 40 GeV) and E miss T > 50 GeV, as well as at least two same-sign leptons (a, b, c) or three leptons (d).The uncertainty bands include the statistical uncertainties for the background prediction as well as the systematic uncertainties for fake-or non-prompt-lepton backgrounds (using the matrix method) and charge-flip electrons.Not included are theoretical uncertainties in the irreducible background contributions.The rare category is defined in the text.
Figure5: Observed and expected exclusion limits on the g, t1 , dR and χ0 1 masses in the context of RPV SUSY scenarios with simplified mass spectra featuring gg or dR dR pair production with exclusive decay modes.The signal regions used to obtain the limits are specified in the subtitle of each scenario.All limits are computed at 95% CL.The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty.The contours of the band around the expected limit are the ±1σ results, including all uncertainties except theoretical uncertainties in the signal cross-section (±2σ is also considered in Figures5(e) and 5(f)).In Figures5(a)-5(d), the diagonal line indicates the kinematic limit for the decays in each specified scenario.For Figures 5(e) and 5(f), theoretical production cross-sections are shown for two different gluino masses in red (1.4 TeV) and blue (2.0 TeV).

Figure 6 :
Figure6: Observed and expected exclusion limits as a function of m 1/2 in the NUHM2 model[31,32].The signal region Rpc2L2bH is used to obtain the limits.The contours of the green (yellow) band around the expected limit are the ±1σ (±2σ) results, including all uncertainties.The limits are computed at 95% CL.

→ χ± 1 W 1 (
when the top squark decays to a top quark and a cascade of electroweakinos χ0 2 ∓ → W * W ∓ χ0 see Figure4(e) for the conditions on the sparticle masses).For the RPV models with gluino pair production (Figures 5(a) -5(d)), a generic exclusion of gluinos with masses below 1.3 TeV is obtained.Weaker exclusion limits, typically around 500 GeV, are obtained in models with pair production of dR (Figures 5(e), 5(f)).

Table 2 :
Summary of the signal region definitions.Unless explicitly stated in the table, at least two signal leptons with p T >20 GeV and same charge (SS) are required in each signal region.Requirements are placed on the number of signal leptons (N

Table 3 :
Summary of the event selection in the validation regions (VRs).Requirements are placed on the number of signal leptons (N signal leptons ), the number of b-jets with p T > 20 GeV (N b-jets ) or the number of jets (N jets ) above a certain p T threshold (p jet T