Search for new physics in same-sign dilepton events 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\,\text {TeV} $$\end{document}s=13TeV

A search for new physics is performed using events with two isolated same-sign leptons, two or more jets, and missing transverse momentum. The results are based on a sample of proton–proton collisions at a center-of-mass energy of 13\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\text {TeV}$$\end{document}TeV recorded with the CMS detector at the LHC, corresponding to an integrated luminosity of 2.3 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{fb}^{1}$$\end{document}fb1. Multiple search regions are defined by classifying events in terms of missing transverse momentum, the scalar sum of jet transverse momenta, the transverse mass associated with a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {W}$$\end{document}W boson candidate, the number of jets, the number of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{b} $$\end{document}b quark jets, and the transverse momenta of the leptons in the event. The analysis is sensitive to a wide variety of possible signals beyond the standard model. No excess above the standard model background expectation is observed. Constraints are set on various supersymmetric models, with gluinos and bottom squarks excluded for masses up to 1300 and 680\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\text {GeV}$$\end{document}GeV, respectively, at the 95 % confidence level. Upper limits on the cross sections for the production of two top quark-antiquark pairs (119\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\text {fb}$$\end{document}fb) and two same-sign top quarks (1.7\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\text {pb}$$\end{document}pb) are also obtained. Selection efficiencies and model independent limits are provided to allow further interpretations of the results.


Introduction
Searches for new physics in final states with two leptons that have same-sign (SS) charges provide a powerful probe for searches of new physics, both because standard model (SM) processes with this signature are few and have low cross sections, and because this signature is produced in a large number of important new-physics scenarios. Examples of the latter include the production of supersymmetric (SUSY) particles [1,2], Majorana neutrinos [3], vector-like quarks [4], and SS top quark pairs [5,6]. In the SUSY framework [7][8][9][10][11][12][13][14][15], the SS signature can arise through gluino pair production. For example, the Majorana nature of the gluino allows gluino pairs to decay via SS charginos, yielding two SS W e-mail: cms-publication-committee-chair@cern.ch bosons. Gluino pair production can also yield four W bosons, e.g., from the decay of four top quarks, which may result in the SS dilepton final state. Alternatively, cascade decays of pair-produced squarks can lead to the SS dilepton signature. Searches for new physics in the SS channel have been previously performed at the CERN LHC by the ATLAS [16][17][18] and CMS [19][20][21][22][23] Collaborations. This paper describes a search for new physics in the final state with two or more leptons and including a SS pair (μ ± μ ± , μ ± e ± , or e ± e ± , where μ is a muon and e an electron). The analysis is based on proton-proton (pp) collision data at √ s = 13 TeV, corresponding to an integrated luminosity of 2.3 fb −1 collected with the CMS detector in 2015. The search strategy resembles that used in our analysis of 19.5 fb 1 of data collected at √ s = 8 TeV [23], which excluded gluino masses in the four top quark signature up to about 1050 GeV. We design an inclusive analysis sensitive to a wide range of new-physics processes produced via strong interactions and yielding undetected particles in the final state. The interpretations of the results consider R-parity conserving SUSY models [24], as well as cross section limits on the production of two top quark-antiquark (tt) pairs and of two SS top quarks. We also provide model independent limits to allow further interpretations of the results. With respect to Ref. [23], the kinematic regions are redefined and improvements in the event selection are implemented, both of which increase the sensitivity to new-physics scenarios at √ s = 13 TeV.

The CMS detector
The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the field volume are several particle detection systems. Charged-particle trajectories are measured with silicon pixel and strip trackers, covering 0 ≤ φ < 2π in azimuth and |η| < 2.5 in pseudorapidity, where η ≡ − ln[tan(θ/2)] and θ is the polar angle of the trajectory of the particle with respect to the counterclockwise beam direction. The transverse momentum, namely the component of the momentum p in the plane orthogonal to the beam, is defined as p T = p sin θ . Surrounding the silicon trackers, a lead tungstate crystal electromagnetic calorimeter and a brass and scintillator hadron calorimeter provide energy measurements of electrons, photons, and hadronic jets in the range |η| < 3.0. Muons are identified and measured within |η| < 2.4 by gas-ionization detectors embedded in the steel flux-return yoke of the solenoid. Forward calorimeters on each side of the interaction point encompass 3.0 < |η| < 5.0. The CMS trigger consists of a two-stage system. The first level of the CMS trigger system, composed of custom hardware processors, uses information from the calorimeters and muon detectors to select events in a fixed time interval of less than 4 µs. The high-level trigger (HLT) processor farm further decreases the event rate from around 100 kHz to less than 1 kHz, before data storage. A more detailed description of the CMS detector can be found in Ref. [25].

Event selection and Monte Carlo simulation
Events are selected with two sets of HLT algorithms. The first requires two very loosely isolated leptons, one satisfying p T > 17 GeV and the other satisfying p T > 8 GeV for a muon and 12 GeV for an electron. The isolation is evaluated with respect to nearby tracks for a muon and to both tracks and calorimetric objects for an electron. The second set of triggers selects events with lowered p T thresholds of 8 GeV and without a restriction on the isolation, but requiring a hadronic activity H HLT T > 300 GeV, where H HLT T is the scalar p T sum of all jets with p T > 40 GeV and |η| < 3.0 identified by the HLT. Typical trigger efficiencies for leptons satisfying the selection criteria described below are 94 % (98 %) per muon (electron), with 100 % efficiency for the H HLT T requirement. In the subsequent analysis, muon candidates are reconstructed by combining information from the silicon tracker and the muon spectrometer in a global fit [26]. A selection is performed using the quality of the geometrical matching between the tracker and muon system measurements. We select muons with well-determined charge by imposing an additional criterion: δp T (μ)/ p T (μ) < 0.2, where δp T (μ) is the uncertainty in the measurement of the muon p T from the global fit.
Electron candidates are reconstructed by combining clusters of energy in the electromagnetic calorimeter with tracks in the silicon tracker [27]. The identification is performed using a Boosted Decision Tree multivariate discriminant [28] based on shower shape and track quality variables. The nominal selection criteria are designed to provide a maximum rejection of electron candidates from multijet production while maintaining approximately 90 % efficiency for electrons from the decay of W or Z bosons. A relaxed selection on the multivariate discriminant is used to define "loose" criteria for electron identification. To improve the accuracy of the electron charge reconstruction, we require the position of the calorimeter deposit, relative to the linear projection of the deposits in the pixel detector to the inner calorimeter surface, to be consistent with the charge determination from the full track fit. Electrons originating from photon conversions are suppressed by rejecting candidates that are either without energy deposits in the innermost layers of the tracking system, or that are associated with a displaced vertex compatible with a photon conversion.
Lepton candidates are required to be consistent with originating from the collision vertex for which the summed p 2 T of the associated physics objects is the largest. The transverse (longitudinal) impact parameter of the leptons must not exceed 0.5 (1.0) mm with respect to this vertex, and they must fulfill the requirement |d 3D |/σ (d 3D ) < 4, where d 3D is the three-dimensional impact parameter with respect to the vertex, and σ (d 3D ) is its uncertainty from the track fit.
The charged leptons produced in decays of heavy particles, such as W and Z bosons or SUSY particles ("prompt" leptons), are typically spatially isolated from the hadronic activity in the event, while leptons produced in hadron decays or in photon conversions, as well as hadrons misidentified as leptons, are usually embedded in jets ("nonprompt" leptons). This distinction becomes less evident for systems with a high Lorentz boost, where decay products tend to overlap and jets may contribute to the energy deposition around prompt leptons. This problem is mitigated with an isolation definition constructed using the following three variables: • the mini-isolation variable (I mini ) [29], computed as the ratio of the scalar p T sum of charged hadrons, neutral hadrons, and photons within a cone of radius R ≡ ( η) 2 + ( φ) 2 around the lepton candidate direction at the vertex, to the transverse momentum of the lepton candidate ( p T ( )). The cone radius R depends on p T ( ) as: The varying isolation cone definition takes into account the increased collimation of the decay products of a hadron as its p T increases, and it reduces the inefficiency from accidental overlap between the lepton and jets in a busy event environment. The momentum estimate of each particle is performed by the particle-flow (PF) algorithm [30,31], which identifies individual particles through a combination of information from different detector components.
• the ratio of the p T of the lepton to that of the closest jet within a distance R = 0.4: where the definition of a jet is given below. In case of no jet within this distance, the value of p ratio T is set to 1. The p ratio T variable is a measure of the isolation in a larger cone and improves the performance of the isolation definition, especially for lowp T nonprompt leptons, which are more likely than highp T leptons to appear in a jet that is wider than the I mini cone.
• the p rel T variable [32], defined as the transverse momentum of the lepton relative to the residual momentum of the closest jet after lepton momentum subtraction: This variable allows the identification of leptons that accidentally overlap with jets.
A lepton is considered to be isolated if the following condition is satisfied: The values of I i , with i = 1, 2, 3, depend on the lepton flavor: because the probability to misidentify a lepton is higher for electrons, tighter isolation values are used in this case (see Table 1). In addition, a "loose" isolation criterion is defined as I mini < 0.4. Muons (electrons) are required to have p T > 10 (15) GeV and |η| < 2.4 (2.5); at least one SS lepton pair with an invariant mass above 8 GeV must be present in the event. In order to reduce backgrounds from inclusive production of the Z boson and from low-mass resonances decaying into lepton pairs, the SS pair is rejected if there is an additional lepton in the event that satisfies loose requirements and that forms an opposite-sign, same-flavor pair with an invariant mass less than 12 GeV or between 76 and 106 GeV with one of the two SS leptons.
Jets and missing transverse momentum (E miss T ) are reconstructed with the PF algorithm. We define E miss T as the magnitude of the vector sum of all PF candidate transverse To identify jets originating from b quarks, the combined secondary vertex algorithm CSVv2 [36] is used. Jets with p T > 25 GeV and |η| < 2.4 are considered as b-tagged if they satisfy the requirements of the medium working point of the algorithm. These requirements result in approximately a 70 % efficiency for tagging a b quark jet, and a less than 1 % mistagging rate for light-quark and gluon jets in tt events. The number of b-tagged jets in the event is denoted as N b .
Monte Carlo (MC) simulation, which includes the contribution of pileup, is used to estimate the background from SM processes with prompt SS leptons (see Sect. 5) and to calculate the efficiency for various new-physics scenarios. The SM background samples are produced with the Mad-Graph5_aMC@NLO 2.2.2 generator [37] at leading order (LO) or next-to-leading order (NLO) accuracy in perturbative quantum chromodynamics, with the exception of diboson samples, which are produced with the powheg v2 [38,39] generator. The NNPDF3.0LO [40] parton distribution functions (PDFs) are used for the simulated samples generated at LO, and the NNPDF3.0NLO [40] PDFs for the samples generated at NLO. Parton showering and hadronization are described using the pythia 8.205 generator [41] with the CUETP8M1 tune [42,43]. The CMS detector response for the background samples is modeled with the Geant4 package [44]. The signal samples are generated with Mad-Graph5_aMC@NLO at LO precision, including up to two additional partons in the matrix element calculations; parton showering and hadronization, as well as decays of SUSY particles, are simulated with pythia, while the detector simulation is performed with the CMS fast simulation package [45].

Search strategy
This analysis is designed as an inclusive search, sensitive to models matching two assumptions: a strong-interaction production mechanism, leading to relatively large hadronic T5qqqqWW Fig. 1 Diagrams illustrating the simplified SUSY models used in this analysis activity, and the presence of undetected particles in the final state, yielding sizable E miss T . In the process of defining the search strategy, R-parity conserving SUSY is taken as a guideline because of its rich variety of signatures. In this context, signal models that can lead to the experimental signature of SS lepton pairs differentiate themselves in the numbers of W bosons, b jets, and light-flavor jets produced in the decays of SUSY particles. In addition, the mass differences between the SUSY particles involved in the decay chains affect the energy spectra of the decay products, resulting in differences between the models in the distributions of kinematic quantities such as the p T of the leptons, H T , and E miss T . We consider SUSY scenarios in the context of simplified models of new-particle production [46,47]. Models with four W bosons and four b jets involve gluino pair production, followed by the decay of each gluino through a chain containing third-generation squarks. If the gluino is lighter than all squarks, and the top squark is the lightest squark, the gluino undergoes a three-body decay mediated by an off-shell top squark. If the dominant top squark decay is t 1 → t χ 0 1 , where χ 0 1 is the lightest neutralino, taken to be the stable, undetected, lightest SUSY particle (LSP), then the gluino threebody decay is g → tt χ 0 1 (T1tttt model in Fig. 1, upper left). If instead the dominant top squark decay is t 1 → b χ + 1 , the gluino three-body decay is g → tb χ + 1 (T5ttbbWW model in Fig. 1, upper middle); the latter signature can also arise if the bottom squark is the lightest squark and decays as b 1 → t χ − 1 . If the top squark is light enough to be on-shell and decays predominantly to a top quark and the LSP, gluinos decay through the chain g → t 1 t → tt χ 0 1 (T5tttt model in Fig. 1, upper right). If instead the top squark mainly decays to the charm quark and the LSP, gluinos decay as in the T5ttcc model ( Fig. 1, lower left); in this case only two W bosons and two b jets are produced.
Events with four W bosons and two b jets can arise from bottom squark pair production, where each bottom squark decays to a top quark and a chargino, and the chargino decays to an LSP and a (possibly off-shell) W boson (T6ttWW model in Fig. 1, lower middle).
Finally, SS lepton pairs can be produced in association with large values of H T , E miss T , and N jet , but without b jets. In particular, events with two W bosons and four light-flavor quark jets can arise from gluino pair production if each gluino decays to two light quarks and a chargino. The two charginos can have the same charge and each decay to a W boson and the LSP (T5qqqqWW model in Fig. 1, lower right). In the case that the difference in mass between the chargino and the LSP is small, the W bosons are off-shell and produce soft leptons.
To increase the sensitivity to new-physics scenarios, we categorize events based on their kinematic properties as follows. First, three exclusive lepton selections are defined: • high-high (HH) selection: two SS leptons, each with p T ≥ 25 GeV; • high-low (HL) selection: two SS leptons, one with p T ≥ 25 GeV and the other with 10 ≤ p T < 25 GeV; • low-low (LL) selection: two SS leptons, each with 10 ≤ p T < 25 GeV.
The high lepton p T threshold suppresses the contribution from nonprompt leptons; hence the SM background in the HH region arises primarily from events with prompt SS leptons. The nonprompt lepton background is largely contained in the HL region, where the highp T lepton is typically prompt and the lowp T lepton nonprompt. The LL region is characterized by a very small background since all processes where at least one lepton originates from an on-shell vec-tor boson are suppressed by the lowp T requirements, while events with two nonprompt leptons are suppressed by the kinematic requirements described below; the main residual contribution in this region is from nonprompt leptons. Second, search regions (SR) are introduced so that the analysis is sensitive to a variety of new-physics scenarios. SRs are defined separately for the HH, HL, and LL selections using the H T , E miss T , N jet , and N b variables: N jet and N b separate signal from background for scenarios with a large production of jets and/or b jets, while H T and E miss T increase sensitivity to models with different masses of SUSY particles. In addition, we make use of the M min T variable, defined as: where the transverse mass and φ ,E miss T is the azimuthal angle difference between the directions of the lepton and of the missing transverse momentum [48]. In the case of an SS lepton pair from tt or W+jets processes, where one lepton is prompt and the other nonprompt, this variable has a cutoff near the W boson mass; consequently, the nonprompt lepton background is suppressed for SRs requiring M min T > 120 GeV and is large for M min T < 120 GeV. In order to better characterize the background we use a fine SR binning in kinematic regions where SM processes are abundant (e.g., low M min T and low E miss T ), while, due to the low background, we use a coarser binning in regions with tight selections.
Finally, inclusive search regions in the HH and HL categories are defined in the tails of the E miss T and H T variables; the boundaries E miss T > 300 GeV and H T > 1125 GeV (for E miss T ≤ 300 GeV) are chosen so that each of these regions typically contains 1 background event.
A summary of the selection criteria is presented in Tables 2, 3 and 4. All SRs are non-overlapping. They are combined statistically to obtain the final results (Sect. 7).

Backgrounds
Backgrounds in the SS dilepton final state can be divided into three categories: Depending on the signal region, this background is dominated by tt and W+jets processes; it represents the largest background for regions with low M min T and low H T .
• SM processes with SS dileptons: Standard model processes that yield an SS lepton pair include multi-boson production (considering W, Z, H, and prompt γ ), single boson production in association with a tt pair, and doubleparton scattering. The dominant sources are WZ and ttW production, which contribute primarily to SRs with zero and ≥1 b jets, respectively. WZ events contribute to the background when the Z boson decays leptonically and is off-shell, when one of the Z-boson decay leptons is not identified, or when the Z boson decays to τ leptons that result in a semileptonic final state. SM processes with SS dileptons are the largest background in the signal regions defined by tight kinematic selections. • Charge misidentification: Charge misidentification arises from events with opposite-sign isolated leptons in which the charge of an electron is misidentified, mostly due to severe bremsstrahlung in the tracker material.
Overall, this is a small background.
The nonprompt lepton background is estimated from data using the "tight-to-loose" ratio method, which was employed in previous versions of the analysis [19-23] but has been improved for the current study. It is based on a control sample of events (application region) where one lepton fails the nominal (tight) selection but passes the loose requirements, defined by relaxing the isolation selection for muons, and both the isolation and identification requirements for electrons. Events in this control region are reweighted by the factor TL /(1 − TL ), where TL is the probability for a nonprompt lepton that satisfies the loose selection to also satisfy the tight selection [19]. Its value is measured in a multijetenriched data set (measurement region), using events from single-lepton triggers after applying a selection designed to suppress electroweak processes (Drell-Yan and W+jets) and after subtracting their residual contribution; this selection requires only one lepton in the event, E miss T < 20 GeV, and M T < 20 GeV. The measurement is made as a function of the lepton p T and η, separately for each lepton flavor (μ or e) and trigger (with or without isolation).
The method assumes that TL has the same value in the measurement and application regions. The main sources of discrepancy are identified as differences in the momentum spectrum and the flavor of the parton producing the nonprompt lepton. These two effects are mitigated in the following way. First, TL is parameterized as a function of p corr T , defined as the lepton p T plus the energy in the isolation cone exceeding the isolation threshold value -this quantity is highly correlated with the mother parton p T , and thus the parameterization is robust against mother parton p T variations. The second effect, i.e., flavor dependence, is relevant for electrons only: while nonprompt muons originate predominantly from heavy-flavor decays, nonprompt electrons receive sizable contributions from misidentified Table 2 SR definitions for the HH selection. The notation ( * ) indicates that, in order to avoid overlaps with SR31, an upper bound E miss T < 300 GeV is used for regions with H T > 300 GeV. All unlabeled region are included in the SR above them, for example the unlabeled regions between SR3 and SR9 are included in SR3, with the exception of the region to the right of SR31, which is included in SR31 hadrons and conversions. The effect of variations in the flavor composition is suppressed by adjusting the loose electron identification criteria so that the numerical value of TL for electrons from light flavors matches that for electrons from heavy flavors. The loose lepton selection is defined based on MC studies, but we verify that TL is not significantly different in data events with and without b jets.
As a cross-check of the prediction, an alternative TL measurement, similar to that described in Ref. [49], is performed in the dilepton control region where one of the leptons fails the impact parameter requirement. The predictions from the two methods are found to be consistent, both in MC samples and in data.
The background from SM processes with a prompt SS lepton pair is evaluated from simulation, accounting for both theoretical and experimental uncertainties. The WZ back-ground is normalized to data in a control region requiring at least two jets, no b jets, E miss T > 30 GeV, and three leptons, where two of the leptons form a same-flavor, opposite-sign pair with an invariant mass within 15 GeV of the Z boson mass; the measured normalization factor is found to be compatible with unity within about one standard deviation. The MC simulation of WZ production is used to relate the number of expected WZ events in the signal regions to the WZ event yield in the control region.
Finally, the charge misidentification background is estimated by reweighting events with opposite-sign lepton pairs by the charge misidentification probability. For electrons this probability is obtained from simulated tt events and from e ± e ± data in the Z mass window, and it lies in the range 10 −5 -10 −3 depending on the electron p T and η. Studies of simulated events indicate that the muon charge misidentification probability is negligible. indicates that, in order to avoid overlaps with SR25, an upper bound E miss T < 300 GeV is used for regions with H T > 300 GeV. All unlabeled region are included in the SR above them, for example the unlabeled regions between SR3 and SR7 are included in SR3, with the exception of the region to the right of SR25, which is included in SR25 Inclusive >120 SR8

Systematic uncertainties
Systematic uncertainties can affect both the overall normalization and the relative population of signal and background processes. A summary of their effects on the SR yields is given in Table 5. Experimental systematic uncertainties are mostly the consequence of differing event selection efficiencies in data and simulation. Lepton identification and trigger efficiencies are computed with the "tag-and-probe" technique [26,27] with an uncertainty of 2 and 4 %, respectively. For signal samples, additional uncertainties of 4-10 % are included to account for differences in the lepton efficiency between the fast and Geant4-based simulations. The jet energy scale uncertainty varies between 2 and 8 %, depending on the jet p T and η. Its impact is assessed by shifting the energy of each jet and propagating the variation to all dependent kinematic quantities (H T , E miss T , N jet , N b , and M min T ); correlation effects due to the migration of events from one SR to another are taken into account. These variations yield estimated uncertainties of 2-10 %. A similar approach is used to estimate the uncertainties associated with the b tagging efficiencies for light- Table 5 Summary of systematic uncertainties in the event yields in the SRs. The first six uncertainties are related to experimental factors for all processes whose yield is estimated from simulation; the next five are uncertainties in these yields related to the event simulation process itself. The last three uncertainties are related to background processes whose yield is estimated from data   [36], which are parameterized as a function of p T and η and are found to be of order 5 % for the highly populated SRs. The uncertainty in the modeling of pileup is 1-5 % depending on the SR. The uncertainty in the integrated luminosity is 2.7 % [50].
The background sources estimated from simulation are subject to theoretical uncertainties related to unknown higher-order effects and to uncertainties in the knowledge of the PDFs. The former are estimated by simultaneously varying the renormalization and factorization scales up and down by a factor of two. The effect on the overall cross section is found to be 13 % for ttW events and 11 % for ttZ events, while the effect on the acceptance for the various SRs amounts to 3-8 % depending on H T . The magnitude of the uncertainty related to the PDFs is obtained using variations of the NNPDF3.0 set [40]. The overall uncertainty is ∼4 % for the ttW and ttZ samples. Theoretical uncertainties are also considered for the remaining minor backgrounds estimated from simulation: a similar procedure is used for the W ± W ± process, leading to an overall uncertainty of 30 %, while a 50 % uncertainty is assigned to processes with a prompt γ and to the sum of the other rare processes. For all backgrounds estimated from simulation we account for the statistical uncertainty of the MC samples.
The remaining sources of uncertainty are those related to the methods that are used to estimate the nonprompt lepton, charge misidentification, and WZ backgrounds. An overall normalization uncertainty of 30 % is assigned to the nonprompt lepton background prediction. This uncertainty accounts for the performance of the method on simulated data and for the differences in the prediction from the two alternative procedures described in Sect. 5. An additional uncertainty is associated with the subtraction procedure to remove Drell-Yan and W+jets events from the measurement region;   the overall effect on the nonprompt lepton background yield is 1-20 %, depending on the SR considered, and is larger for highp T leptons. Finally, we account for the statistical uncertainty in the number of events observed in the application region. The background from charge misidentification is assigned a systematic uncertainty of 26 %, which corresponds to the difference between the e ± e ± event yield in the Z mass window in data and simulation.
The uncertainty in the WZ background is measured to be 30 % in the control region. It includes statistical uncertainties and systematic uncertainties due to non-WZ background subtraction. Using the same procedure as described above, uncertainties in the extrapolation from the control to the signal regions are assessed from the propagation of the uncertainty in the jet energy scale and in the b tagging efficiencies.

Results
Distributions of the five kinematic variables used to define the SRs, H T , E miss T , M min T , N jet , and N b , are shown in Fig. 2 after a baseline selection requiring a pair of SS leptons, two jets, and either E miss T > 30 GeV or H T > 500 GeV. The results are shown in comparison to the background prediction. The event yields in the SRs after the full selection are  Fig. 3 and in Table 6; no significant deviation from the SM background prediction is observed. The largest local significances are 2.2 and 1.8 standard deviations in HL SR8 and in HH SR10, respectively. The results of the search are used to constrain the benchmark SUSY models presented in Sect. 4. For each mass point in the SUSY particle mass spectrum, results from all SRs are combined to extract cross section exclusion limits at the 95 % confidence level (CL), using the asymptotic formulation of the modified frequentist CL s criterion [51-54]. Signal and background uncertainties are included as log-normal nuisance parameters and, when relevant, take into account correlation effects among different SRs and/or different processes. Exclusion contours make use of the cross section values calculated at NLO plus next-to-leading logarithmic (NLL) accuracy, assuming that all SUSY particles other than those included in the respective diagram are too heavy to participate in the interaction [55][56][57][58][59][60]. In general, the SR with the largest sensitivity is HH SR31, which requires E miss T > 300 GeV and is inclusive in the other variables. However, depending on the model and the region of parameter space, other SRs contribute significantly to the total sensitivity: for instance, a considerable contribution comes from HL SR25 in case of signal models with a soft lepton, from HH SR32 and HL SR26 in case of high H T , from HH SR3 and SR8 in case of no b jets, and from HH SR24 and SR26 in case of 2 or more b jets.
Results for models with gluinos decaying to virtual third generation squarks are shown in Fig. 4 as a function of the gluino and LSP masses. For the T1tttt model (Fig. 4-left), in the regions of the SUSY parameter space with a large mass difference between the gluino and the LSP, the results are rather stable with respect to LSP mass variations, and gluino masses up to 1300 GeV are excluded. Near the kinematic threshold m g − m χ 0 1 = 2(m W + m b ), the gluino mass limit becomes weaker and is reduced to 1050 GeV for an LSP mass of 800 GeV. Results for the T5ttbbWW model with nearly degenerate χ ± 1 and χ 0 1 masses are shown in Fig. 4-right; the limit on the gluino mass lies in the range 950-1100 GeV except for very small χ ± 1 and χ 0 1 masses, where the sensitivity increases because of the large Lorentz boost of the leptons from the χ ± 1 decay. Results for models with a gluino decaying to an onshell top squark are shown in Fig. 5 as a function of the gluino and LSP masses. For the T5tttt model (Fig. 5-top), for which we take m t 1 = m χ 0 1 + m t , similar exclusion curves are obtained as for the T1tttt model in Fig. 4-left because the production cross section and the final-state particles are the same. The limit becomes weaker when there is a small mass difference between the top squark and the LSP: for m t 1 − m χ 0 1 = 20 GeV, the limit on the gluino mass is 1140 GeV for small LSP masses and about 850 GeV for m χ 0 1 = 700 GeV (Fig. 5-bottom left). In the case of the T5ttcc model with the same SUSY particle mass values, the sensitivity is slightly reduced because of the smaller number of leptons and b jets in the final state (Fig. 5-bottom right). Figure 6 shows the results for b squark production in the T6ttWW model in the chargino ( χ ± 1 ) versus b squark mass Results for the T5qqqqWW model are shown in Fig. 7 as a function of the gluino and LSP masses, with two different assumptions for the chargino mass: it is either assumed to be the average of m g and m χ 0 1 , or it is set to m χ 0 1 + 20 GeV. In the first case (Fig. 7-left), the exclusion limit on gluino masses exceeds 1100 GeV for LSP masses up to 400 GeV; for larger LSP masses the limit is reduced to 830 GeV at m χ 0 1 = 700 GeV. In the second case (Fig. 7-right), due to the smaller mass difference, leptons in the final state are soft and thus the sensitivity is reduced.
The results of the search are also used to set 95 % CL upper limits on the double tt production cross section, whose SM value computed at NLO precision [37] is 9.1 fb. The upper limit on σ (pp → tttt) is found to be 119 fb, with an expected result of 102 +57 −35 fb. With the current integrated luminosity, the sensitivity to this signature is limited by the statistical precision.
Limits at the 95 % CL on the SS top quark pair production cross section are determined using events that satisfy the baseline selection categorized according to number of b jets (Fig. 2-bottom right); apart from the charge requirement, the detector acceptance and the selection efficiency for the signal are assumed to match those of SM tt events. The observed (expected) upper limit on σ (pp → tt)+σ (pp → tt) is 1.7 pb (1.5 +0.7 −0.4 pb). Finally, we report model independent limits on the product of cross section, detector acceptance, and selection efficiency, σA , for the production of an SS dilepton pair in the two inclusive HH regions, SR31 and SR32, using the CL s criterion without the asymptotic approximation. In SR31 the limit is computed as a function of the minimum threshold on E miss T for H T > 300 GeV, while in SR32 it is computed as a function of the H T threshold for 50 < E miss T < 300 GeV. The results are shown in Fig. 8, where, in regions with no observed events, the minimum limit value of 1.3 fb is obtained. These limits can be used to test additional BSM models, after accounting for the event selection efficiency. The lepton efficiency ranges between 70-85 % (45-70 %) for generated muons (electrons) with |η| < 2.4 and p T > 25 GeV, increasing as a function of p T and converging to the maximum value for p T > 60 GeV; the efficiencies of the H T and E miss T requirements are mostly determined by the jet energy and E miss T resolutions, which are discussed in Refs. [33,35].

Summary
The results of a search for new physics in same-sign dilepton events using the CMS detector at the LHC and based on a data sample of pp collisions at √ s = 13 TeV, corresponding to an integrated luminosity of 2. Acknowledgments We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centres and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construc-