Search for physics beyond the standard model in events with two leptons of same sign, missing transverse momentum, and jets in proton-proton collisions at sqrt(s) = 13 TeV

A data sample of events from proton-proton collisions with two isolated same-sign leptons, missing transverse momentum, and jets is studied in a search for signatures of new physics phenomena by the CMS Collaboration at the LHC. The data correspond to an integrated luminosity of 35.9 inverse femtobarns, and a center-of-mass energy of 13 TeV. The properties of the events are consistent with expectations from standard model processes, and no excess yield is observed. Exclusion limits at 95% confidence level are set on cross sections for the pair production of gluinos, squarks, and same-sign top quarks, as well as top-quark associated production of a heavy scalar or pseudoscalar boson decaying to top quarks, and on the standard model production of events with four top quarks. The observed lower mass limits are as high as 1500 GeV for gluinos, 830 GeV for bottom squarks. The excluded mass range for heavy (pseudo)scalar bosons is 350-360 (350-410) GeV. Additionally, model-independent limits in several topological regions are provided, allowing for further interpretations of the results.


Introduction
Final states with two leptons of same charge, denoted as same-sign (SS) dileptons, are produced rarely by standard model (SM) processes in proton-proton (pp) collisions. Because the SM rates of SS dileptons are low, studies of these final states provide excellent opportunities to search for manifestations of physics beyond the standard model (BSM). Over the last decades, a large number of new physics mechanisms have been proposed to extend the SM and address its shortcomings. Many of these can give rise to potentially large contributions to the SS dilepton signature, e.g., the production of supersymmetric (SUSY) particles [1,2], SS top quarks [3,4], scalar gluons (sgluons) [5,6], heavy scalar bosons of extended Higgs sectors [7,8], Majorana neutrinos [9], and vector-like quarks [10].
In the SUSY framework [11][12][13][14][15][16][17][18][19][20], the SS final state can appear in R-parity conserving models through gluino or squark pair production when the decay of each of the pair-produced particles yields one or more W bosons. For example, a pair of gluinos (which are Majorana particles) can give rise to SS charginos and up to four top quarks, yielding signatures with up to four W bosons, as well as jets, b quark jets, and large missing transverse momentum (E miss T ). Similar signatures can also result from the pair production of bottom squarks, subsequently decaying to charginos and top quarks.
While R-parity conserving SUSY models often lead to signatures with large E miss T , it is also interesting to study final states without significant E miss T beyond what is produced by the neutrinos from leptonic W boson decays. For example, some SM and BSM scenarios can lead to the production of SS or multiple top quark pairs, such as the associated production of a heavy (pseudo)scalar, which subsequently decays to a pair of top quarks. This scenario is realized in Type II two Higgs doublet models (2HDM) where associated production with a single top quark or a tt pair can in some cases provide a promising window to probe these heavy (pseudo)scalar bosons [21][22][23].
This paper extends the search for new physics presented in Ref. [24]. We consider final states with two leptons (electrons and muons) of same charge, two or more hadronic jets, and moderate E miss T . Compared to searches with zero or one lepton, this final state provides enhanced sensitivity to low-momentum leptons and SUSY models with compressed mass spectra. The results are based on an integrated luminosity corresponding to 35.9 fb −1 of √ s = 13 TeV protonproton collisions collected with the CMS detector at the CERN LHC. Previous LHC searches in the SS dilepton channel have been performed by the ATLAS [25-27] and CMS [24, 28-32] Collaborations. With respect to Ref. [24], the event categorization is extended to take advantage of the increased integrated luminosity, the estimate of rare SM backgrounds is improved, and the (pseudo)scalar boson interpretation is added.
The results of the search are interpreted in a number of specific BSM models discussed in Section 2. In addition, model-independent results are also provided in several kinematic regions to allow for further interpretations. These results are given as a function of hadronic activity and of E miss T , as well as in a set of inclusive regions with different topologies. The full analysis results are also summarized in a smaller set of exclusive regions to be used in combination with the background correlation matrix to facilitate their reinterpretation.

Background and signal simulation
Monte Carlo (MC) simulations are used to estimate SM background contributions and to estimate the acceptance of the event selection for BSM models. The MADGRAPH5 aMC@NLO 2 2 Background and signal simulation 2.2.2 [33][34][35] and POWHEG v2 [36,37] next-to-leading order (NLO) generators are used to simulate almost all SM background processes based on the NNPDF3.0 NLO [38] parton distribution functions (PDFs). New physics signal samples, as well as the same-sign W ± W ± process, are generated with MADGRAPH5 aMC@NLO at leading order (LO) precision, with up to two additional partons in the matrix element calculations, using the NNPDF3.0 LO [38] PDFs. Parton showering and hadronization, as well as the double-parton scattering production of W ± W ± , are described using the PYTHIA 8.205 generator [39] with the CUETP8M1 tune [40,41]. The GEANT4 package [42] is used to model the CMS detector response for background samples, while the CMS fast simulation package [43] is used for signal samples.
To improve on the MADGRAPH modeling of the multiplicity of additional jets from initialstate radiation (ISR), MADGRAPH tt MC events are reweighted based on the number of ISR jets (N ISR J ), so as to make the light-flavor jet multiplicity in dilepton tt events agree with the one observed in data. The same reweighting procedure is applied to SUSY MC events. The reweighting factors vary between 0.92 and 0.51 for N ISR J between 1 and 6. We take one half of the deviation from unity as the systematic uncertainty in these reweighting factors.
The new physics signal models probed by this search are shown in Figs. 1 and 2. In each of the simplified SUSY models [44,45] of Fig. 1, only two or three new particles have masses sufficiently low to be produced on-shell, and the branching fraction for the decays shown are assumed to be 100%. Gluino pair production models giving rise to signatures with up to four b quarks and up to four W bosons are shown in Figs. 1a-e. In these models, the gluino decays to the lightest squark ( g → qq), which in turn decays to same-flavor ( q → q χ 0 1 ) or differentflavor ( q → q χ ± 1 ) quarks. The chargino decays to a W boson and a neutralino ( χ ± 1 → W ± χ 0 1 ), where the χ 0 1 escapes detection and is taken to be the lightest SUSY particle (LSP). The first two scenarios considered in Figs. 1a and 1b include an off-shell third-generation squark ( t or b) leading to the three-body decay of the gluino, g → tt χ 0 1 (T1tttt) and g → tb χ + 1 (T5ttbbWW), resulting in events with four W bosons and four b quarks. In the T5ttbbWW model, the mass splitting between chargino and neutralino is set to m χ ± 1 − m χ 0 1 = 5 GeV, so that two of the W bosons are produced off-shell and can give rise to low transverse momentum (p T ) leptons. The next two models shown (Figs. 1c and d) include an on-shell top squark with different mass splitting between the t and the χ 0 1 , and consequently different decay modes: in the T5tttt model the mass splitting is equal to the top quark mass (m t − m χ 0 1 = m t ), favoring the t → t χ 0 1 decay, while in the T5ttcc model the mass splitting is only 20 GeV, favoring the flavor changing neutral current t → c χ 0 1 decay. In Fig. 1e, the decay proceeds through a virtual light-flavor squark, leading to a three-body decay to g → qq χ ± 1 , resulting in a signature with two W bosons and four light-flavor jets. The two W bosons can have the same charge, giving rise to SS dileptons. This model, T5qqqqWW, is studied as a function of the gluino and χ 0 1 mass, with two different assumptions for the chargino mass: ), producing mostly on-shell W bosons, and m χ ± 1 = m χ 0 1 + 20 GeV, producing off-shell W bosons. Finally, Fig. 1f shows a model of bottom squark production followed by the b → t χ ± 1 decay, resulting in two b quarks and four W bosons. This model, T6ttWW, is studied as a function of the b and χ ± 1 masses, keeping the χ 0 1 mass at 50 GeV, resulting in two of the W bosons being produced off-shell when the χ ± 1 and χ 0 1 masses are close. The production cross sections for SUSY models are calculated at NLO plus next-to-leading logarithmic (NLL) accuracy [46][47][48][49][50][51].
The processes shown in Fig. 2, ttH, tHq, and tWH, represent the top quark associated production of a scalar (H) or a pseudoscalar (A). The subsequent decay of the (pseudo)scalar to a pair of top quarks then gives rise to final states including a total of three or four top quarks. For the purpose of interpretation, we use LO cross sections for the production of a heavy Higgs boson

The CMS detector and event reconstruction
The central feature of the CMS detector is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter (HCAL), each composed of a barrel and two endcap sections. Forward calorimeters extend the pseudorapidity coverage provided by the barrel and endcap detectors. Muons are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid. A more detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [49].
Events of interest are selected using a two-tiered trigger system [50]. The first level (L1), composed of custom hardware processors, uses information from the calorimeters and muon detectors to select events at a rate of around 100 kHz within a time interval of less than 4 µs. The second level, known as the high-level trigger (HLT), consists of a farm of processors running a version of the full event reconstruction software optimized for fast processing, and reduces the event rate to less than 1 kHz before data storage.
Events are processed using the particle-flow (PF) algorithm [51,52], which reconstructs and identifies each individual particle with an optimized combination of information from the various elements of the CMS detector. The energy of photons is directly obtained from the ECAL measurement. The energy of electrons is determined from a combination of the elec- in the context of the Type II 2HDM of Ref. [23]. The mass of the new particle is varied in the range [350, 550] GeV, where the lower mass boundary is chosen in such a way as to allow the decay of the (pseudo)scalar into on-shell top quarks.

The CMS detector and event reconstruction
The central feature of the CMS detector is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter (HCAL), each composed of a barrel and two endcap sections. Forward calorimeters extend the pseudorapidity (η) coverage provided by the barrel and endcap detectors. Muons are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid. A more detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [52].
Events of interest are selected using a two-tiered trigger system [53]. The first level (L1), composed of custom hardware processors, uses information from the calorimeters and muon detectors to select events at a rate of around 100 kHz within a time interval of less than 4 µs. The second level, known as the high-level trigger (HLT), consists of a farm of processors running a version of the full event reconstruction software optimized for fast processing, and reduces the event rate to less than 1 kHz before data storage.
Events are processed using the particle-flow (PF) algorithm [54,55], which reconstructs and identifies each individual particle with an optimized combination of information from the various elements of the CMS detector. The energy of photons is directly obtained from the ECAL measurement. The energy of electrons is determined from a combination of the electron momentum at the primary interaction vertex as determined by the tracker, the energy of the corresponding ECAL cluster, and the energy sum of all bremsstrahlung photons spatially compatible with the electron track [56]. The energy of muons is obtained from the curvature of the corresponding track, combining information from the silicon tracker and the muon system [57]. The energy of charged hadrons is determined from a combination of their momentum measured in the tracker and the matching ECAL and HCAL energy deposits, corrected for the response function of the calorimeters to hadronic showers. Finally, the energy of neutral hadrons is obtained from the corresponding corrected ECAL and HCAL energy.
Hadronic jets are clustered from neutral PF candidates and charged PF candidates associated with the primary vertex, using the anti-k T algorithm [58,59] with a distance parameter R = √ (∆η) 2 + (∆φ) 2 of 0.4. Jet momentum is determined as the vectorial sum of all PF candidate momenta in the jet. An offset correction is applied to jet energies to take into account the contribution from additional proton-proton interactions (pileup) within the same or nearby bunch crossings. Jet energy corrections are derived from simulation, and are improved with in situ measurements of the energy balance in dijet and photon+jet events [60,61]. Additional selection criteria are applied to each event to remove spurious jet-like features originating from isolated noise patterns in certain HCAL regions. Jets originating from b quarks are identified (b tagged) using the medium working point of the combined secondary vertex algorithm CSVv2 [62]. The missing transverse momentum vector p miss T is defined as the projection on the plane perpendicular to the beams of the negative vector sum of the momenta of all reconstructed PF candidates in an event [63]. Its magnitude is referred to as E miss T . The sum of the transverse momenta of all jets in an event is referred to as H T .

Event selection and search strategy
The event selection and the definition of the signal regions (SRs) follow closely the analysis strategy established in Ref. [24]. With respect to the previous search, the general strategy has remained unchanged. We target, in a generic way, new physics signatures that result in SS dileptons, hadronic activity, and E miss T , by subdividing the event sample into several SRs sensitive to a variety of new physics models. The number of SRs was increased to take advantage of the larger integrated luminosity. Table 1 summarizes the basic kinematic requirements for jets and leptons (further details, including the lepton identification and isolation requirements, can be found in Ref. [24]). Table 1: Kinematic requirements for leptons and jets. Note that the p T thresholds to count jets and b-tagged jets are different.
Events are selected using triggers based on two sets of HLT algorithms, one simply requiring two leptons, and one additionally requiring H T > 300 GeV. The H T requirement allows for the lepton isolation requirement to be removed and for the lepton p T thresholds to be set to 8 GeV for both leptons, while in the pure dilepton trigger the leading and subleading leptons are required to have p T > 23 (17) GeV and p T > 12 (8) GeV, respectively, for electrons (muons). Based on these trigger requirements, leptons are classified as high (p T > 25 GeV) and low (10 < p T < 25 GeV) momentum, and three analysis regions are defined: high-high (HH), highlow (HL), and low-low (LL).
The baseline selection used in this analysis requires at least one SS lepton pair with an invariant mass above 8 GeV, at least two jets, and E miss T > 50 GeV. To reduce Drell-Yan backgrounds, events are rejected if an additional loose lepton forms an opposite-sign same-flavor pair with one of the two SS leptons, with an invariant mass less than 12 GeV or between 76 and 106 GeV. Events passing the baseline selection are then divided into SRs to separate the different background processes and to maximize the sensitivity to signatures with different jet multiplicity (N jets ), flavor (N b ), visible and invisible energy (H T and E miss T ), and lepton momentum spectra (the HH/HL/LL categories mentioned previously). The m min T variable is defined as the smallest of the transverse masses constructed between p miss T and each of the leptons. This variable features a cutoff near the W boson mass for processes with only one prompt lepton, so it is used to create SRs where the nonprompt lepton background is negligible. To further improve sensitivity, several regions are split according to the charge of the leptons (++ or −−), taking advantage of the charge asymmetry of SM backgrounds, such as ttW or WZ, with a single W boson produced in pp collisions. Only signal regions dominated by such backgrounds and with a sufficient predicted yield are split by charge. In the HH and HL categories, events in the tail regions H T > 1125 GeV or E miss T > 300 GeV are inclusive in N jets , N b , and m min T in order to ensure a reasonable yield of events in these SRs. The exclusive SRs resulting from this classification are defined in Tables 2-4. The lepton reconstruction and identification efficiency is in the range of 45-70% (70-90%) for electrons (muons) with p T > 25 GeV, increasing as a function of p T and converging to the maximum value for p T > 60 GeV. In the low-momentum regime, 15 < p T < 25 GeV for electrons and 10 < p T < 25 GeV for muons, the efficiencies are 40% for electrons and 55% for muons. The lepton trigger efficiency for electrons is in the range of 90-98%, converging to the maximum value for p T > 30 GeV, and around 92% for muons. The chosen b tagging working point results in approximately a 70% efficiency for tagging a b quark jet and a <1% mistagging rate for light-flavor jets in tt events [62]. 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. [60,61,64].

Backgrounds
Standard model background contributions arise from three sources: processes with prompt SS dileptons, mostly relevant in regions with high E miss T or H T ; events with a nonprompt lepton, dominating the overall final state; and opposite-sign dilepton events with a chargemisidentified lepton, the smallest contribution. In this paper we use the shorthand "nonprompt leptons" to refer to electrons or muons from the decays of heavy-or light-flavor hadrons, hadrons misidentified as leptons, or electrons from conversions of photons in jets.
Several categories of SM processes that result in the production of electroweak bosons can give rise to an SS dilepton final state. These include production of multiple bosons in the same event (prompt photons, W, Z, and Higgs bosons), as well as single-boson production in association with top quarks. Among these SM processes, the dominant ones are WZ, ttW, and ttZ production, followed by the W ± W ± process. The remaining SM processes are grouped into two categories, "Rare" (including ZZ, WWZ, WZZ, ZZZ, tWZ, tZq, as well as tttt and double parton scattering) and "X+γ" (including Wγ, Zγ, ttγ, and tγ). The expected yields from these SM backgrounds are estimated from simulation, accounting for both the theoretical and experimental uncertainties discussed in Section 6.
For the WZ and ttZ backgrounds, a three-lepton (3L) control region in data is used to scale the simulation, based on a template fit to the distribution of the number of b jets. The 3L control region requires at least two jets, E miss T > 30 GeV, and three leptons, two of which must form an opposite-sign same-flavor pair with an invariant mass within 15 GeV of the Z boson mass. In the fit to data, the normalization and shapes of all the components are allowed to vary according to experimental and theoretical uncertainties. The scale factors obtained from the fit in the phase space of the 3L control region are 1.26 ± 0.09 for the WZ process, and 1.14 ± 0.30 for the ttZ process.
The nonprompt lepton background, which is largest for regions with low m min T and low H T , is estimated by the "tight-to-loose" method, which was employed in several previous versions of the analysis [28-32], and significantly improved in the latest version [24] to account for the kinematics and flavor of the parent parton of the nonprompt lepton. The tight-to-loose method uses two control regions, the measurement region and the application region. The measurement region consists of a sample of single-lepton events enriched in nonprompt leptons by requirements on E miss T and transverse mass that suppress the W → ν contribution. This sample is used to extract the probability for a nonprompt lepton that satisfies the loose selection to also satisfy the tight selection. This probability ( TL ) is calculated as a function of lepton p corr T (defined below) and η, separately for electrons and muons, and separately for lepton triggers with and without an isolation requirement. The application region is a SS dilepton region where both of the leptons satisfy the loose selection but at least one of them fails the tight selection. This region is subsequently divided into a set of subregions with the exact same kinematic requirements as those in the SRs. Events in the subregions are weighted by a factor TL /(1 − TL ) for each lepton in the event failing the tight requirement. The nonprompt background in each SR is then estimated as the sum of the event weights in the corresponding subregion. The p corr T parametrization, where p corr T is defined as the lepton p T plus the energy in the isolation cone exceeding the isolation threshold value, is chosen because of its correlation with the parent parton p T , improving the stability of the TL values with respect to the sample kinematics. To improve the stability of the TL values with respect to the flavor of the parent parton, the loose electron selection is adopted. This selection increases the number of nonprompt electrons from the fragmentation and decay of light-flavor partons, resulting in TL values similar to those from heavy-flavor parent partons.
The prediction from the tight-to-loose method is cross-checked using an alternative method based on the same principle, similar to that described in Ref. [65]. In this cross-check, which aims to remove kinematic differences between measurement and application regions, the measurement region is obtained from SS dilepton events where one of the leptons fails the impact parameter requirement. With respect to the nominal method, the loose lepton definition is adapted to reduce the effect of the correlation between isolation and impact parameter. The predictions of the two methods are found to be consistent within systematic uncertainties.
Charge misidentification of electrons is a small background that can arise from severe bremsstrahlung in the tracker material. Simulation-based studies with tight leptons indicate that the muon charge misidentification probability is negligible, while for electrons it ranges between 10 −5 and 10 −3 . The charge misidentification background is estimated from data using an opposite-sign control region for each SS SR, scaling the control region yield by the charge misidentification probability measured in simulation. A low-E miss T control region, with e + e − pairs in the Z boson mass window, is used to cross-check the MC prediction for the misidentification probability, both inclusively and -where the number of events in data allows it -as a function of electron p T and η.

Systematic uncertainties
Several sources of systematic uncertainty affect the predicted yields for signal and background processes, as summarized in Table 5. Experimental uncertainties are based on measurements in data of the trigger efficiency, the lepton identification efficiency, the b tagging efficiency [62], the jet energy scale, and the integrated luminosity [66], as well as on the inelastic cross section value affecting the pileup rate. Theoretical uncertainties related to unknown higher-order effects are estimated by varying simultaneously the factorization and renormalization scales by a factor of two, while uncertainties in the PDFs are obtained using replicas of the NNPDF3.0 set [38].
Experimental and theoretical uncertainties affect both the overall yield (normalization) and the relative population (shape) across SRs, and they are taken into account for all signal samples as well as for the samples used to estimate the main prompt SS dilepton backgrounds: WZ, ttW, ttZ, W ± W ± . For the WZ and ttZ backgrounds, the control region fit results are used for the normalization, so these uncertainties are only taken into account for the shape of the backgrounds. For the smallest background samples, Rare and X+γ, a 50% uncertainty is assigned in place of the scale and PDF variations.
The normalization and the shapes of the nonprompt lepton and charge misidentification backgrounds are estimated from control regions in data. In addition to the statistical uncertainties from the control region yields, dedicated systematic uncertainties are associated with the methods used in this estimate. For the nonprompt lepton background, a 30% uncertainty (increased to 60% for electrons with p T > 50 GeV) accounts for the performance of the method in simulation and for the differences in the two alternative methods described in Section 5. In addition, the uncertainty in the prompt lepton yield in the measurement region, relevant when estimating TL for high-p T leptons, results in a 1-30% effect on the estimate. For the charge misidentification background, a 20% uncertainty is assigned to account for possible mismodeling of the charge misidentification rate in simulation.

Results and interpretation
A comparison between observed yields and the SM background prediction is shown in Fig. 3 for the kinematic variables used to define the analysis SRs: H T , E miss T , m min T , N jets , and N b . The distributions are shown after the baseline selection defined in Section 4. The full results of the search in each SR are shown in Fig. 4 and Table 6. The SM predictions are generally consistent with the data. The largest deviations are seen in HL SR 36 and 38, with a local significance, taking these regions individually or combining them with other regions adjacent in phase space, that does not exceed 2 standard deviations.
These results are used to probe the signal models discussed in Section 2: simplified SUSY models, (pseudo)scalar boson production, four top quark production, and SS top quark production. We also interpret the results as model-independent limits as a function of H T and E miss T . With the exception of the new (pseudo)scalar boson limits, the results can be compared to the previous version of the analysis [24], showing significant improvements due to the increase in the integrated luminosity and the optimization of SR definitions.
To obtain exclusion limits at the 95% confidence level (CL), the results from all SRs -including Table 5: Summary of the sources of uncertainty and their effect on the yields of different processes in the SRs. The first two groups list experimental and theoretical uncertainties assigned to processes estimated using simulation, while the last group lists uncertainties assigned to processes whose yield is estimated from data. The uncertainties in the first group also apply to signal samples. Reported values are representative for the most relevant signal regions.

Source
Typical signal and background uncertainties and their correlations -are combined using an asymptotic formulation of the modified frequentist CL s criterion [67][68][69][70]. When testing a model, all new particles not included in the specific model are considered too heavy to take part in the interaction. To convert cross section limits into mass limits, the signal cross sections specified in Section 2 are used.
The observed SUSY cross section limits as a function of the gluino and LSP masses, as well as the observed and expected mass limits for each simplified model, are shown in Fig. 5 for gluino pair production models with each gluino decaying through a chain containing off-or on-shell third-generation squarks. These models, which result in signatures with two or more b quarks and two or more W bosons in the final state, are introduced in Section 2 as T1tttt, T5ttbbWW, T5tttt, and T5ttcc. Figure 6 shows the limits for a model of gluino production followed by a decay through off-shell first-or second-generation squarks and a chargino. Two different assumptions are made on the chargino mass, taken to be between that of the gluino and the LSP. These T5qqqqWW models result in no b quarks and either on-shell or off-shell W bosons. Bottom squark pair production followed by a decay through a chargino, T6ttWW, resulting in two b quarks and four W bosons, is shown in Fig. 7. For all of the models probed, the observed limit agrees well with the expected one, extending the reach of the previous analysis by 200-300 GeV and reaching 1.5, 1.1, and 0.83 TeV for gluino, LSP, and bottom squark masses, respectively.
The observed and expected cross section limits on the production of a heavy scalar or a pseudoscalar boson in association with one or two top quarks, followed by its decay to top quarks, are shown in Fig. 8. The limits are compared with the total cross section of the processes described in Section 2. The observed limit, which agrees well with the expected one, excludes scalar (pseudoscalar) masses up to 360 (410) GeV.
The SM four top quark production, pp → tttt, is normally included among the rare SM backgrounds. When treating this process as signal, its observed (expected) cross section limit is determined to be 42 (27 +13 −8 ) fb at 95% CL, to be compared to the SM expectation of 9.2 +2.9 −2.4 fb [33]. This is a significant improvement with respect to the observed (expected) limits obtained    2.92 ± 0.50 3 SR47 1.78 ± 0.42 3 SR48 1.95 ± 0.39 5 SR49 1.23 ± 0.30 3 SR50 1.46 ± 0.31 0 SR51 0.74 ± 0.18 0  Figure 5: Exclusion regions at 95% CL in the m χ 0 1 versus m g plane for the T1tttt (a) and T5ttbbWW (b) models, with off-shell third-generation squarks, and the T5tttt (c) and T5ttcc (d) models, with on-shell third-generation squarks. For the T5ttbbWW model, m χ ± 1 = m χ 0 1 + 5 GeV, for the T5tttt model, m t − m χ 0 1 = m t , and for the T5ttcc model, m t − m χ 0 1 = 20 GeV and the decay proceeds through t → c χ 0 1 . The right-hand side color scale indicates the excluded cross section values for a given point in the SUSY particle mass plane. The solid, black curves represent the observed exclusion limits assuming the NLO+NLL cross sections [46][47][48][49][50][51] (thick line), or their variations of ±1 standard deviation (thin lines). The dashed, red curves show the expected limits with the corresponding ±1 and ±2 standard deviation experimental uncertainties. Excluded regions are to the left and below the limit curves.   −24 ) fb [71]. The results of the search are also used to set a limit on the production cross section for SS top quark pairs, σ(pp → tt) + σ(pp → tt). The observed (expected) limit, based on the kinematics of a SM tt sample and determined using the number of b jets distribution in the baseline region, is 1.2 (0.76 +0. 3 −0.2 ) pb at 95% CL, significantly improved with respect to the 1.7 (1.5 +0.7 −0.4 ) pb observed (expected) limit of the previous analysis [24].

Model-independent limits and additional results
The yields and background predictions can be used to test additional BSM physics scenarios. To facilitate such reinterpretations, we provide limits on the number of SS dilepton pairs as a function of the E miss T and H T thresholds in the kinematic tails, as well as results from a smaller number of inclusive and exclusive signal regions.
The E miss T and H T limits are based on combining HH tail SRs, specifically SR42-45 for high E miss T and SR46-51 for high H T , and employing the CL s criterion without the asymptotic formulation as a function of the minimum threshold of each kinematic variable. These limits are presented in Fig. 9 in terms of σA , the product of cross section, detector acceptance, and selection efficiency. Where no events are observed, the observed and expected limits reach 0.1 fb, to be compared with a limit of 1.3 fb obtained in the previous analysis [24].
Results are also provided in Table 7 for a small number of inclusive signal regions, designed based on different topologies and a small number of expected background events. The background expectation, the event count, and the expected BSM yield in any one of these regions can be used to constrain BSM hypotheses in a simple way.
In addition, we define a small number of exclusive signal regions based on integrating over the standard signal regions. Their definitions, as well as the expected and observed yields, are specified in Table 8, while the correlation matrix for the background predictions in these regions is given in Fig. 10. This information can be used to construct a simplified likelihood for

Summary
A sample of same-sign dilepton events produced in proton-proton collisions at 13 TeV, corresponding to an integrated luminosity of 35.9 fb −1 , has been studied to search for manifestations of physics beyond the standard model. The data are found to be consistent with the standard model expectations, and no excess event yield is observed. The results are interpreted as limits at 95% confidence level on cross sections for the production of new particles in simplified supersymmetric models. Using calculations for these cross sections as functions of particle masses, the limits are turned into lower mass limits that are as high as 1500 GeV for gluinos and  830 GeV for bottom squarks, depending on the details of the model. Limits are also provided on the production of heavy scalar (excluding the mass range 350-360 GeV) and pseudoscalar (350-410 GeV) bosons decaying to top quarks in the context of two Higgs doublet models, as well as on same-sign top quark pair production, and the standard model production of four top quarks. Finally, to facilitate further interpretations of the search, model-independent limits are provided as a function of H T and E miss T , together with the background prediction and data yields in a smaller set of signal regions.

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 centers 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 construction and operation of the LHC and the CMS detector provided by the following funding agencies: BMWFW and FWF (Aus [25] ATLAS Collaboration, "Search for gluinos in events with two same-sign leptons, jets and missing transverse momentum with the ATLAS detector in pp collisions at √ s = 7 TeV", Phys. Rev. Lett. 108 (2012) 241802, doi:10.1103/PhysRevLett.108.241802, arXiv:1203.5763.