J H E P 0 6 ( 2 0 1 1 ) 0 7 7 Published for SISSA by Springer Received: April 16, 2011 Revised: May 23, 2011 Accepted: May 30, 2011 Published: June 17, 2011 Search for new physics with same-sign isolated dilepton events with jets and missing transverse energy at the LHC The CMS collaboration Abstract: The results of searches for new physics in events with two same-sign isolated leptons, hadronic jets, and missing transverse energy in the final state are presented. The searches use an integrated luminosity of 35 pb−1 of pp collision data at a centre-of-mass energy of 7 TeV collected by the CMS experiment at the LHC. The observed numbers of events agree with the standard model predictions, and no evidence for new physics is found. To facilitate the interpretation of our data in a broader range of new physics scenarios, information on our event selection, detector response, and efficiencies is provided. Keywords: Hadron-Hadron Scattering Open Access, Copyright CERN, for the benefit of the CMS collaboration doi:10.1007/JHEP06(2011)077 http://dx.doi.org/10.1007/JHEP06(2011)077 J H E P 0 6 ( 2 0 1 1 ) 0 7 7 Contents 1 Introduction 1 2 The CMS detector 2 3 Reconstruction of leptons, missing energy, and jets 3 4 Search regions 4 4.1 Searches using lepton triggers 5 4.2 Searches using hadronic triggers 6 5 Background estimation 7 5.1 Searches using lepton triggers 8 5.2 Search using hadronic triggers, electrons, and muons 10 5.3 Search using hadronic triggers and τh 15 5.4 Comparison of leptonic and HT -triggered analyses with electrons and muons in the final state 16 5.5 Electron charge mismeasurement 17 6 Signal acceptance and efficiency systematic uncertainties 18 7 Summary of results 19 8 Interpretation of results 20 9 Summary and conclusions 24 The CMS collaboration 28 1 Introduction Events with same-sign isolated lepton pairs from hadron collisions are very rare in the standard model (SM) but appear very naturally in many new physics scenarios. In par- ticular, they have been proposed as signatures of supersymmetry (SUSY) [1–3], universal extra dimensions [4], pair production of T5/3 (a fermionic partner of the top quark) [5], heavy Majorana neutrinos [6], and same-sign top-pair resonances as predicted in theories with warped extra dimensions [7]. In this paper we describe searches for new physics with same-sign isolated dileptons (ee, eµ, µµ, eτ , µτ , and ττ), missing transverse energy (Emiss T ), and hadronic jets. Our choice of signal regions is driven by two simple observa- tions. First, astrophysical evidence for dark matter [8] suggests that we concentrate on final states with Emiss T . Second, observable new physics signals with large cross sections are – 1 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 likely to be produced by strong interactions, and we thus expect significant hadronic activ- ity in conjunction with the two same-sign leptons. Beyond these simple guiding principles, our searches are as independent of detailed features of new physics models as possible. The results are based on a data sample corresponding to an integrated luminosity of 35 pb−1 collected in pp collisions at a centre-of-mass energy of 7 TeV by the Compact Muon Solenoid (CMS) experiment at the Large Hadron Collider (LHC) in 2010. This paper is organized as follows. The CMS detector is briefly described in section 2. The reconstruction of leptons, Emiss T , and jets at CMS is summarized in section 3. Section 4 describes our search regions. We perform separate searches based on leptonic and hadronic triggers in order to cover a wider region in the parameter space of new physics. Electron and muon triggers allow for searches that require less hadronic energy in the event, while hadronic triggers allow inclusion of lower transverse momentum (pT ) electrons and muons, as well as hadronic τ decays in the final state. The dominant backgrounds for all three searches are estimated from data, as discussed in section 5. Systematic uncertainties on the predicted number of signal events and results of these searches are discussed in sections 6 and 7. We conclude with a discussion on how to use our results to constrain a wide variety of new physics models in section 8. 2 The CMS detector A right-handed coordinate system is employed by the CMS experiment, with the origin at the nominal interaction point, the x-axis pointing to the centre of the LHC, and the y-axis pointing up (perpendicular to the LHC plane). The polar angle θ is measured from the positive z-axis and the azimuthal angle φ is measured in the xy plane. The pseudorapidity is defined as η = − ln [tan ( θ 2)]. The central feature of the CMS apparatus is a superconducting solenoid, of 6 m in- ternal diameter, 13 m in length, providing an axial field of 3.8 T. Within the field volume are several particle detection systems which each feature a cylindrical geometry, covering the full azimuthal range from 0 ≤ φ ≤ 2π. Silicon pixel and strip tracking detectors provide measurements of charged particle trajectories and extend to a pseudorapidity of |η| = 2.5. A homogeneous crystal electromagnetic calorimeter (ECAL) and a sampling brass/scintilator hadronic calorimeter (HCAL) surround the tracking volume and provide energy measurements of electrons, photons, and hadronic jets up to |η| = 3.0. An iron- quartz fiber hadronic calorimeter, which is also part of the HCAL system, is located in the forward region defined by 3.0 < |η| < 5.0. Muons are measured in the pseudorapidity range of |η| < 2.4, with detection planes made of three technologies: drift tubes, cathode strip chambers, and resistive plate chambers. These are instrumented outside of the mag- net coil within the steel return yoke. The CMS detector is nearly hermetic, allowing for energy balance measurements in the plane transverse to the beam direction. A two-tier trigger system is designed to select the most interesting pp collision events for use in physics analysis. A detailed description of the CMS detector can be found elsewhere [9]. – 2 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 3 Reconstruction of leptons, missing energy, and jets Muon candidates are required to be successfully reconstructed [10] using two algorithms, one in which tracks in the silicon detector are matched to consistent signals in the calorime- ters and muon system, and another in which a simultaneous global fit is performed to hits in the silicon tracker and muon system. The track associated with the muon candidate is required to have a minimum number of hits in the silicon tracker, have a high-quality global fit including a minimum number of hits in the muon detectors, and have calorimeter energy deposits consistent with originating from a minimum ionizing particle. Electron candidates are reconstructed [11] starting from a cluster of energy deposits in the ECAL, which is then matched to hits in the silicon tracker. A selection using electron identification variables based on shower shape and track-cluster matching is applied to the reconstructed candidates; the criteria are optimized in the context of the inclusive W → eν measurement [12] and are designed to maximally reject electron candidates from QCD multijet production while maintaining approximately 80% efficiency for electrons from the decay of W/Z bosons. Electron candidates within ∆R = √ ∆φ2 + ∆η2 < 0.1 of a muon are rejected to remove electron candidates due to muon bremsstrahlung and final- state radiation. Electron candidates originating from photon conversions are suppressed by looking for a partner track and requiring no missing hits for the track fit in the inner layers of the tracking detectors. Hadronic τ candidates (τh) are identified [13] starting with a hadronic jet clustered from the particles reconstructed using the particle-flow global-event reconstruction algorithm [14]. The highest-pT charged track within a cone of ∆R < 0.1 around the jet axis is required to have pT > 5 GeV. A variable size cone of ∆R < 5 GeV/pT is then defined around this track, and the boosted τ -decay products are expected to be confined within this narrow cone. Only τ candidates with one or three charged hadrons in this cone are selected. The discrimination between hadronic τ decays and generic QCD jets is based on an ensemble of five neural networks, each of which has been trained to identify one of the five main hadronic τ -decay modes using the kinematics of the reconstructed charged and neutral pions [15]. All lepton candidates are required to have |η| < 2.4, and be consistent with originating from the same interaction vertex. Charged leptons from the decay of W/Z bosons, as well as the new physics we are searching for, are expected to be isolated from other activity in the event. We calculate a relative measure of this isolation denoted as RelIso . This quantity is defined as the ratio of the scalar sum of transverse track momenta and trans- verse calorimeter energy deposits within a cone of ∆R < 0.3 around the lepton candidate direction at the origin, to the transverse momentum of the candidate. The contribution from the candidate itself is excluded. In order to suppress the background due to dileptons originating from the same jet, we require that selected dileptons have a minimum invariant mass of 5GeV. This helps to keep dileptons uncorrelated with respect to their RelIso observables, which is a feature we exploit in the analysis. We also remove events with a third lepton of opposite sign and same flavour as one of the two selected leptons if the invariant mass of the pair is between 76 – 3 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 � � Figure 1. An example of a process involving the production and decays of SUSY particles, which gives rise to two same-sign prompt leptons, jets, and missing transverse energy. and 106 GeV. This requirement further reduces an already small background contribution from WZ and ZZ production. Jets and Emiss T are reconstructed based on the particle-flow technique desribed in [14, 16]. For jet clustering, we use the anti-kT algorithm with the distance parameter R = 0.5 [17]. Jets are required to pass standard quality requirements [18] to remove those consistent with calorimeter noise. Jet energies are corrected for residual nonuniformity and nonlinearity of the detector response derived using collision data [19]. We require jets to have transverse energy above 30 GeV and to be within |η| < 2.5. We define the HT observable as the scalar sum of the pT of all such jets with ∆R > 0.4 to the nearest lepton passing all our requirements. 4 Search regions The searches discussed in this paper employ two different trigger strategies, electron and muon triggers in one case, and HT triggers in the other. The leptonic triggers allow for lower HT requirements, while the HT triggers allow for lower lepton pT , as well as final states with hadronic τ decays. The motivation for covering the widest possible phase space in this search can be illustrated by an example of a SUSY cascade, shown in figure 1, naturally giving rise to jets, Emiss T , and same-sign leptons: (gluinos/squarks) → (charged gaugino) → (lightest supersymmetric particle (LSP) neutralino). The mass difference between the gluino/squarks and the charged gaugino, typically arbitrary, defines the amount of hadronic activity one may expect in the event. The mass difference between the gaugino and a neutralino influences the lepton pT spectrum. Depending on the nature of the chargino and neutralino, their mass difference can be either arbitrary (e.g., wino and bino) or typically small (e.g., higgsinos). Moreover, there are a number of ways to generate a large production asymmetry between τ and e/µ leptons, which motivates us to look specifically for events with a τ . – 4 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 In the following we describe the search regions explored by each trigger strategy. As a new-physics reference point, we use LM0, a point in the constrained Minimal Super- symmetric Standard Model (CMSSM) [20] defined with the model parameters m0 = 200 GeV, m1/2 = 160 GeV, tanβ = 10, µ > 0, and A0 = −400 GeV. LM0 is one of the common CMSSM reference points used in CMS across many analyses. As the abbreviation suggests, LM0 provides squarks and gluinos with relatively low masses, and thus has a large pro- duction cross section. It is beyond the exclusion reach of the searches performed by LEP and Tevatron, but it has recently been excluded by searches [21–23] with ATLAS and CMS concurrently with this one. Nonetheless, we continue to use LM0 as it provides a common model for which to compare our sensitivity with that of other analyses. Aside from the LM0 point, several SM simulation samples are used to both validate and complement var- ious background estimation methods that are based on the data itself. These samples rely on either PYTHIA 6.4 [24] or MADGRAPH [25] for event generation and GEANT4 [26] for simulation of the CMS detector. Samples used include tt, single-t, γ+jets, W+jets, Z+jets, WW, WZ, ZZ, and QCD multijet production. Next-to-leading-order (NLO) cross sections are used for all samples except for QCD multijet production. 4.1 Searches using lepton triggers We start with a baseline selection inspired by our published tt → ℓ+ℓ− + X (ℓ = e or µ) cross section measurement [27]. Events are collected using single and dilepton triggers. The detailed implementation of these triggers evolved throughout the 2010 data-collecting period as the LHC instantaneous luminosity was increasing. Trigger efficiencies are measured from a pure lepton sample collected using Z → ℓ+ℓ− decays from data. The luminosity-averaged efficiency to trigger on events with two leptons with |η| < 2.4 and pT > 10 GeV, one of which also has pT > 20 GeV, is very high. For example, the trigger efficiency for an LM0 event passing the baseline selection described below is estimated to be (99 ± 1)%. One of the electrons and muons must have pT > 20 GeV and the second one must have pT > 10 GeV. Both leptons must be isolated. The isolation requirement is based on the RelIso variable introduced earlier. We require RelIso < 0.1 for leptons of pT > 20 GeV, and the isolation sum (i.e., the numerator of the RelIso expression) to be less than 2 GeV for pT < 20 GeV. We require the presence of at least two reconstructed jets, implying HT > 60 GeV. Finally, we require the missing transverse energy Emiss T > 30 GeV (ee and µµ) or Emiss T > 20 GeV (eµ). This defines our baseline selection. Following the guiding principles discussed in the introduction, we define two search re- gions. The first has high Emiss T (Emiss T > 80 GeV); the second has high HT (HT > 200 GeV). These Emiss T and HT values were chosen to obtain an SM background expectation in sim- ulation of 1/3 of an event in either of the two overlapping search regions. Figure 2 shows the HT versus Emiss T scatter plot for the baseline selection, indicating the Emiss T and HT requirements for the two search regions via horizontal and vertical lines, respectively. Figure 2 (left) shows three events (red dots) in the baseline region, one of which barely satisfies HT > 200 GeV, but fails the Emiss T > 80 GeV requirement. – 5 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 (GeV)TH 100 200 300 400 500 600 700 800 900 1000 ( G e V ) m is s T E 50 100 150 200 250 300 350 400 data, tight selection LM0 MC = 7 TeVsCMS, -1 = 35 pbintL (GeV) T H 100 200 300 400 500 600 700 800 900 1000 ( G e V ) m is s T E 50 100 150 200 250 300 350 400 data, loose selection = 7 TeVsCMS, -1 = 35 pbintL Figure 2. HT versus Emiss T scatter plots for baseline region. (Left) Overlay of the three observed events with the expected signal distribution for LM0. The three observed events all scatter in the lower left corner of the plot. (Right) Scatter plot of the background in data when only one of the two leptons is required to be isolated. In contrast, most of the signal from typical supersymmetry models tends to pass both of these requirements, as is visible in the LM0 expected signal distribution overlaid in figure 2 (left). Backgrounds to this analysis are dominated by events with jets mimicking leptons, as discussed in section 5. Requiring only one of the two leptons to be isolated thus allows us to increase the background statistics in order to display the expected distribution of SM background events in the (Emiss T , HT ) plane, as shown in figure 2 (right). Backgrounds clearly cluster at low Emiss T and low HT , with slightly more than half of the events failing both the Emiss T and HT selections. Moreover, comparing the left and right plots in figure 2 indicates that the lepton isolation requirement on both leptons versus only one lepton reduces the backgrounds by roughly a factor of ten. 4.2 Searches using hadronic triggers Hadronic triggers allow us to explore the phase space with low-pT electrons and muons, as well as final states with hadronic τ decays. We allow muons (electrons) with pT as low as 5 (10) GeV, and restrict ourselves to τh with visible transverse momentum > 15 GeV, where τh refers to hadronic τ candidates only. All leptons must be isolated with RelIso < 0.15. For the ee, eµ, and µµ final states, we require at least two jets, HT > 300 GeV, and Emiss T > 30 GeV. As backgrounds from QCD multijet production are significant for τh, we increase the Emiss T and HT requirements to Emiss T > 50 GeV and HT > 350 GeV in the eτh, µτh, and τhτh final states. Figure 3 shows the efficiency turn-on curves for the HT triggers used during three different data taking periods. The trigger thresholds were changing in order to cope with the increasing instantaneous luminosities over the 2010 running period. Roughly half of – 6 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 (GeV)TH 50 100 150 200 250 300 350 400 T ri g g e r e ff ic ie n c y 0 0.2 0.4 0.6 0.8 1 1 T H 2 T H 3 T H -1 = 35 pb int = 7 TeV, LsCMS Figure 3. HT Trigger efficiency as a function of the reconstructed HT for three data-collecting periods: 7 pb−1 with HT 1, 10 pb−1 with HT 2, and 18 pb−1 with HT 3. the integrated luminosity in 2010 was taken with the highest threshold trigger. This mea- surement indicates that at HT = 300 GeV the efficiency reaches (94 ± 5)% . These trigger turn-on curves are measured in data with events selected by muon triggers. 5 Background estimation Standard model sources of same-sign dilepton events with both leptons coming from a W or Z decay are very small in our data sample. Simulation-based predictions of the combined yields for qq → WZ and ZZ, double “W-strahlung” qq → q′q′W±W±, double parton scattering 2 × (qq → W±), ttW, and WWW comprise no more than a few percent of the total background in any of the final states considered. As these processes have never been measured in proton-proton collisions, and their background contributions are very small, we evaluate them using simulation, assigning a 50% systematic uncertainty. The background contribution from pp → Wγ, where the W decays leptonically and the photon converts in the detector material giving rise to an isolated electron, is also estimated from simulation and found to be negligible. All other backgrounds are evaluated from data, as discussed below. Backgrounds in all of our searches are dominated by one or two jets mimicking the lepton signature. Such lepton candidates can be genuine leptons from heavy-flavour decays, electrons from unidentified photon conversions, muons from meson decays in flight, hadrons reconstructed as leptons, or jet fluctuations leading to hadronic τ signatures. We will refer to all of these as ”fake leptons”. Leptons from W, Z, gauginos, etc., i.e., the signal we are searching for, will be referred to as ”prompt leptons”. The dominant background contribution is from events with one lepton, jets, and Emiss T — mostly tt̄ with one lepton from the W decay, and a second lepton from the decay of a heavy-flavour particle. These events contain one prompt and one fake lepton, and are estimated via two different techniques described in sections 5.1 and 5.2. While both techniques implement an extrapolation in lepton isolation, they differ in the assumptions made. Both techniques lead to consistent predictions as described in section 5.4, provid- – 7 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 ing additional confidence in the results. Backgrounds with two fake leptons are generally smaller, except in the final state with two hadronic τ leptons, where the dominant back- ground source is QCD multijet production. Contributions due to fake τh are estimated using an extrapolation from “loose” to “tight” τh identification, as described in section 5.3. For the ee and eµ final states, electron charge misreconstruction due to hard bremsstrahlung poses another potentially important background, as there are significant opposite-sign ee and eµ contributions, especially from tt̄, where both W’s from the top quarks decay leptonically. This is discussed in section 5.5. 5.1 Searches using lepton triggers Contributions from fake leptons are estimated using the so-called “tight-loose” (TL) method [27, 28]. In this method the probability ǫTL for a lepton passing loose selections to also pass the tight analysis selections is measured in QCD multijet events as a function of lepton pT and η. The key assumption of the method is that ǫTL is approximately universal, i.e., it is the same for all jets in all event samples. Tests of the validity of this assumption are described below. The main difference between the tight and loose lepton selections is that the re- quirement on the RelIso variable defined in section 3 is relaxed from RelIso < 0.1 to RelIso < 0.4. Other requirements that are relaxed are those on the distance of closest approach between the lepton track and the beamline (impact parameter) and, in the case of muons, the selection on the χ2 of the muon track fit. The quantity ǫTL is measured in a sample of lepton-trigger events with at least one jet satisfying pT > 40 GeV and well separated (∆R > 1) from the lepton candidate. We refer to this jet as the “away-jet”. We reduce the impact of electroweak background (W, Z, tt) by excluding events with Z → ℓℓ candidates, events with Emiss T > 20 GeV, and events where the transverse mass MT of the lepton and the Emiss T is greater than 25 GeV. Studies based on simulation indicate that this procedure results in an unbiased estimate of ǫTL up to lepton pT ≈ 40 GeV. At higher transverse lepton momenta the remaining electroweak contributions in the sample have a significant effect. Thus, ǫTL is measured only up to pT = 35 GeV. It is taken to be constant at higher transverse momenta, as suggested by simulation studies. The level of universality of ǫTL is tested with different jet samples. Two types of tests are relevant and both involve the parent jet from which the lepton originates. The first test explores the sensitivity to the jet’s pT , and the second test explores the sensitivity to the jet’s heavy-flavour content. Sensitivity to jet pT stems from the fact that the probability for a lepton of a given pT to pass the RelIso selection depends on the pT of the parton from which the lepton originates. To be explicit, a 10 GeV lepton originating from a 60 GeV b quark is less likely to pass our RelIso requirement than the same lepton originating from a 20 GeV b quark. The heavy flavour sensitivity can be traced to semileptonic decays, which are a source of leptons in bottom and charm jets, but not in light-quark and gluon jets. To test the jet pT dependence, we select loose leptons in events with the away-jet above a varying jet pT threshold. Since these events are mostly QCD dijets, the pT of the away-jet – 8 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 (GeV) T electron p 10 15 20 25 30 35 T L p ro b a b ili ty 0 0.1 0.2 0.3 0.4 0.5 > 20 GeV T jet p > 40 GeV T jet p > 60 GeV T jet p = 7 TeVsCMS, -1 = 35 pbintL (GeV) T muon p 10 15 20 25 30 35 T L p ro b a b ili ty 0 0.1 0.2 0.3 0.4 > 20 GeV T jet p > 40 GeV T jet p > 60 GeV T jet p = 7 TeVsCMS, -1 = 35 pbintL Figure 4. Electron (left) and muon (right) TL probability ǫTL computed from QCD multijet events with different requirements on the minimum pT of the away-jet. The probabilities shown are projections of the two-dimensional function ǫTL(η, pT ) onto the pT axis. is a good measure of the pT of the jet from which the lepton originates. We weight each event by ǫTL measured as described above, i.e., requiring the away-jet to have pT > 40 GeV. We then sum the weights and compare the sum to the number of observed leptons passing tight requirements. Varying the away-jet minimum pT requirement from 20 to 60 GeV, we find the observed yield to differ from the predicted yield by +54% (+49%) and −4% (−3%) for muons (electrons) in this test. The percentages here, as well as throughout this section refer to (observed − predicted)/predicted. This non-negligible jet pT dependence can also be seen in figure 4, where we show ǫTL calculated using different away-jet thresholds. To test the heavy-flavour dependence, we repeat the exercise requiring that the away-jet be above pT > 40 GeV and be b-tagged, i.e., a jet in which we find a secondary vertex well separated from the interaction point consistent with a b-hadron decay. By applying the b tag on the away-jet, the sample of jets from which the lepton originates is enriched in heavy flavours. Introducing this b tag we find the observed yield to differ from that predicted by −3% (−15%) for muons (electrons). We have thus shown that applying an ǫTL obtained without a b-tagging requirement to a sample with such a requirement leads to a modest difference between observed and predicted. This validates our assumption that the TL method is flavour universal. To predict the background from prompt lepton + jets events in a signal region, the TL probability is applied to a sample of dilepton events satisfying all the signal selection requirements, but where one of the leptons fails the tight selections and passes the loose ones. Each event is weighted by the factor ǫTL/(1 − ǫTL), where ǫTL is the tight-to-loose probability for the loose lepton in the event. The background contribution from this source is then estimated by summing the weights of all such events (S1). The sum S1 also includes the contribution from backgrounds with two fake leptons. However, these are double counted because in the case of two fake leptons passing loose – 9 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 requirements there are two combinations with one lepton passing the tight selections. The background contribution with two fake leptons is estimated separately by selecting events where both leptons pass the loose requirements but fail the tight requirements. Each event in this sample is weighted by the product of the two factors of ǫTL/(1−ǫTL) corresponding to the two leptons in the event, and the sum S2 of weights is used to estimate the background with two fake leptons. The total background from events with one or two fake leptons is then obtained as S1−S2. In kinematic regions of interest for this search, S2 is typically more than one order of magnitude smaller than S1, indicating that the main background contribution is from one prompt lepton and one fake lepton. The method has been tested on simulated tt̄ and W + jets events. In these tests, we use ǫTL measured from QCD simulation events to predict the number of same-sign dilepton events in these samples. In the tt̄ simulation sample, we find that the observed yield differs from the prediction for the baseline selection by -41% (fake muons) and -47% (fake electrons). Observed and predicted yields in this test are consistent with each other at the 5% confidence level (CL) because the simulation statistics are modest. The same level of agreement is also found for the two search regions. In the W + jets case, the ratio of predicted events to observed is 0.8 ± 0.4 for fake electrons; the statistics for fake muons are not sufficient to draw any definitive conclusions. Based on these studies, as well as the dependence on away-jet pT and heavy-flavour composition discussed above, we assign a ±50% systematic uncertainty on the ratio (observed − predicted)/predicted and, hence, on the estimation of backgrounds due to fake leptons. In addition to this systematic uncertainty, the method has significant statistical uncertainties based on the number of events in the samples to which the TL probability is applied. We find 6 (4) events in these samples for the Emiss T > 80 GeV (HT > 200 GeV) search regions. The resulting background estimates in the two regions are 1.1± 0.6 and 0.9± 0.6 events, respectively, including only statistical uncertainties. As an additional cross-check, we determine the background estimate and observed yields in the baseline region. We estimate 3.2 ± 0.9 ± 1.6 events from background due to fake leptons alone, and 3.4 ± 0.9 ± 1.6 after all backgrounds are taken into account. The uncertainties here are statistical and systematic, respectively. The composition of the total background is estimated to be 86% (7%) events with one (two) fake leptons, 3% due to charge misidentification, and 4% irreducible background for which both leptons are isolated leptons from leptonic W or Z decay. As mentioned in section 4.1, there are 3 events observed in the baseline region, in good agreement with the background estimate. Applying Emiss T > 80 GeV or HT > 200 GeV increases the fraction of events with one fake lepton, but statistical uncertainties on the individual components of the background estimate are too large to meaningfully quantify the change in the relative contributions of the components. 5.2 Search using hadronic triggers, electrons, and muons As described in section 4.2, hadronic triggers allow us to explore the phase space with low- pT leptons. However, lowering the lepton pT is expected to increase the relative contribution of events with two fake leptons. As shown below, this background now constitutes roughly 30% of the total background, as compared to only a few percent for the higher-pT thresholds – 10 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 in the search regions for the leptonic trigger analysis. This motivates the development of a method exclusively dedicated to predicting and understanding the QCD background with two fake leptons. At the same time, increasing the HT requirement to 300GeV, as driven by the hadronic trigger thresholds, reduces the expected W+jets background to only a few percent of the total background. The background with one fake lepton is now reduced to tt and single- t processes, where the fake lepton is due mostly to semileptonic b decays. Therefore, we tailor the method for estimating the background with one fake lepton to these expectations. The method is similar to the TL technique, but has a number of important differences that are discussed further below. For this analysis, the estimation of the background with fake leptons starts with an evaluation of background events with two fake leptons and then proceeds with an estimation of the contribution of events with only one such lepton. First, we define a preselection control sample of events from the HT -triggered data stream with at least two same-sign dileptons and with all event selection requirements applied, except for those related to Emiss T and isolation. We find 223 µµ, 6 ee, and 78 eµ events of this type. The large asymmetry between muons and electrons is mostly due to differences in the corresponding pT thresholds of 5 and 10 GeV, respectively. In addition, identification of electrons within jets is less efficient than that for muons. The preselection control sample is dominated by QCD multijet production. Studies based on simulation suggest that we should attribute about 10% of the preselection yields to tt contamination, while attributing a much smaller fraction to W+jets. The contribution from events with two fake leptons to the signal region is estimated by assuming that the three requirements, RelIso < 0.15 for each lepton and Emiss T > 30 GeV, are mutually independent and, hence, the total background-suppression efficiency can be written in the factorized form ǫtot = ǫℓ1 iso · ǫℓ2 iso · ǫMET. This assumption has been verified both in simulation and directly in data. With simulation it is straightforward to prove the principle in the nominal preselection region because we can safely measure the efficiencies in a dedicated QCD sample where we know all leptons can be considered background (i.e., no contamination from prompt leptons exists). In data the contribution from prompt leptons is non-negligible and therefore some extra selection requirements are necessary to isolate a QCD enriched control sample. We validate the factorized expression for ǫtot in two steps in data. First, we demon- strate that the selection requirement on RelIso is independent for each lepton. We begin by relaxing the HT selection to 200 GeV and add events collected with leptonic triggers to gain more statistics. We then require Emiss T < 20 GeV to suppress events with leptonic W decays. Figure 5 (left) shows that the single-muon efficiency can be squared to obtain the double-muon efficiency, thus validating the assumption that the RelIso observable is uncorrelated between the two fake leptons and the efficiencies can be factorized. In the second step, we demonstrate that the Emiss T and RelIso selection requirements are mutually independent. To accomplish this in data, we maintain HT above 300 GeV, but we include single-lepton events to increase statistics. To suppress the contributions from events with leptonic W decays, we modify the selection requirement on the lepton impact parame- – 11 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 Relative Isolation ­110 1 E ff ic ie n c y ­310 ­210 ­110 1 1 µ ε efficiency: µMeasured single­ 2 µ+ 1 µ ε efficiency: µMeasured double­ 2 1 µ ε efficiency: µPredicted double­ ­1 = 35 pb int = 7 TeV, LsCMS [GeV]miss TE 0 10 20 30 40 50 R e la ti v e I s o la ti o n E ff ic ie n c y 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 < 0.1 mm) 0 Muons (d < 0.1 mm) 0 Electrons (d > 0.1 mm) 0 Muons (d > 0.1 mm) 0 Electrons (d ­1 = 35 pb int = 7 TeV, LsCMS Figure 5. (Left) The lepton isolation efficiency for one (solid squares) and two (open squares) leptons as a function of the relative isolation parameter cut. Also shown is the predicted double- lepton efficiency if the two lepton efficiencies are assumed to be independent of each other. Only the dimuon sample is shown here. (Right) The lepton isolation efficiency as a function of the Emiss T cut for electrons and muons with different requirements on the lepton impact parameter. Details are given in the text. ter from the nominal d0 < 0.2 mm to d0 > 0.1 mm. Figure 5 (right) shows that the RelIso selection efficiency for muons and electrons remains constant as a function of the Emiss T selection requirement. The dashed lines represent the zeroth-order polynomial fits to the efficiency measurements made in the d0 control region for muons and electrons, respec- tively. For completeness, we also show the obvious bias arising when the impact parameter requirement is inverted to d0 < 0.1 mm to enrich the sample with leptonic W decays. It is important to note that no attempt is made to apply the RelIso selection efficiency measured in the control region defined by d0 > 0.1 mm to the above formula for ǫtot. This control region is only used to demonstrate the stability of the RelIso selection efficiency with respect to the Emiss T requirement. The actual values of ǫℓ1 iso and ǫℓ2 iso are measured in the nominal preselection region (d0 < 0.2 mm), where we assume this stability, and hence factorization, remains valid for events with two fake leptons. Having validated the selection factorization hypothesis, we proceed to measure the isolation and Emiss T selection efficiencies, one at a time, in the preselection control sample, where we obtain ǫµ iso = 0.036 ± 0.015, ǫe iso = 0.11 ± 0.08, ǫMET = 0.27 ± 0.03. Uncer- tainties quoted are statistical only. As before, we suppress leptonic W decays to reduce possible biases. We accomplish this by requiring either Emiss T < 20 GeV or RelIso > 0.2 when measuring ǫµ iso, ǫe iso, or ǫMET. The appropriate product of these efficiencies is then applied to the event counts observed in the preselection control sample, leading to the background estimate of 0.18 ± 0.12 ± 0.12 events. Uncertainties here are statistical and systematic, respectively. The systematic uncertainties quoted above have two dominant sources. One is due – 12 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 to limited statistics in simulation and data when validating that the three requirements are indeed independent. We take the statistical precision (25%) of this cross-check as our systematic uncertainty on the method. The other dominant source of systematic uncer- tainty can be attributed to the inability of the inverted selection requirements on Emiss T and RelIso to fully suppress contributions from leptonic W decays (e.g., tt, W+jets) while measuring ǫµ iso, ǫe iso, and ǫMET. Studies based on simulation suggest that the bias (over- estimate) can be as large as 60%, mostly via a bias in measuring ǫe iso. Conservatively, we do not correct for the possible bias, but take it as a systematic uncertainty. We thus arrive at a 65% systematic uncertainty on the estimate of backgrounds due to events with two fake leptons by adding these two effects in quadrature. In addition, we verified in simula- tion that the techniques used to suppress leptonic W decays (i.e., inversion of the Emiss T and RelIso requirements) do not alter or bias the selection efficiencies for fake leptons. It is worth mentioning that this method of evaluating background with two fake lep- tons does not require any reweighting of measured efficiencies. The average efficiencies are obtained from a QCD-dominated subset of the preselection sample, and then applied directly to the preselection sample as a whole to derive the prediction for the number of events with two fake leptons in the signal region. Next, we proceed with estimating the contribution of backgrounds with a single misidentified lepton. We start from a tight-loose control sample, to be further referred to as a sideband, and use the isolation selection efficiency for b jets, referred to as ǫ(b), to predict event counts in the signal region. The sideband control sample contains events passing all signal selection criteria, except one of the two leptons is now required to have RelIso > 0.15. To begin, we count the number of events in this sample: 11 (µµ), 2 (ee), 6 (eµ), and 5 (µe), the last lepton indicating which one in the pair is non-isolated. Then, we estimate the contribution of the background with two fake leptons to the sideband sample using the efficiencies quoted above from the factorization procedure. For example, for the dimuon channel, the contribution to the sideband from events with two fake leptons is Nµµ preselected · 2ǫµ iso(1 − ǫµ iso) · ǫMET. The resulting yield estimates for events with two fake leptons are 4.2 (µµ), 0.32 (ee), 2.3 (eµ), and 0.68 (µe). After subtracting this con- tribution, the remaining yields in the sideband are consistent with simulation predictions assuming that only tt (76%), single-t (7%), and W+jets (15%) contribute. This remaining sideband yield after subtraction is then scaled by an appropriate factor determined using the BTag-and-probe method [28], as described below. The BTag-and-probe method relies on the basic premise that events with one fake lepton can be attributed to tt production, with one prompt lepton from leptonic W decay and the second fake lepton from semi-leptonic b decay. The efficiencies ǫ (b) µ iso and ǫ (b) e iso are thus defined as the probabilities of a muon or electron from semi-leptonic b decay to pass the RelIso < 0.15 selection. These efficiencies can be measured in data using appropriately selected events from bb production. To determine ǫ (b) µ iso and ǫ (b) e iso, we select a bb enriched control sample by requiring one b-tagged away-jet and one lepton candidate. In addition, we require HT > 100 GeV to arrive at a b-quark pT spectrum similar to that expected for the tt background in the search. To reduce the bias from leptonic W or Z decays, we furthermore require Emiss T < 15 GeV and MT < 15 GeV, and veto events with two leptons – 13 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 Figure 6. Isolation variable distributions obtained with the BTag-and-probe method for muons (left) and electrons (right). Efficiencies for the RelIso < 0.15 (first bin in the distributions shown) are explicitly quoted. forming a mass within 7 GeV of the mass of the Z boson. Approximately 80% of the leptons in this sample are from semileptonic heavy-flavour decay. We find that the resulting bb̄ control sample differs sufficiently from the expected tt background in both lepton kinematics and jet multiplicity (Njets) to warrant corrections. We therefore measure the RelIso distribution in the bb control sample in data in bins of lepton pT and Njets, and reweight these distributions using event probabilities ω(pT , Njets) derived from a tt̄ simulation sample. The resulting reweighted RelIso distributions for these three samples are overlaid in figure 6 for muons (left) and electrons (right). The plots show distributions for tt̄ simulation (red crosses) after all selections except RelIso on one of the two leptons, reweighted bb̄ simulation (grey shade), and reweighted bb̄-enriched data (black dots). The agreement between the two simulation-based distributions validates the method. Agreement between data and simulation is observed but is not required for this method to be valid. The contents of the first bin of the two data plots are the above mentioned isolation selection efficiencies ǫ (b) µ iso and ǫ (b) e iso. We find ǫ (b) µ iso = 0.029+0.003 −0.002 and ǫ (b) e iso = 0.036+0.013 −0.008, with uncertainties due to statistics only. We probed four different potential sources of systematic uncertainties in the BTag-and- probe method, and added their contributions in quadrature to arrive at a total systematic uncertainty on the ǫ (b) iso efficiency of 54 (29)% for electrons (muons). The largest contribution to this uncertainty is due to the statistical precision with which the method in simulation is verified. Other potential sources of systematic errors were found to give subdominant contributions. We evaluated the sensitivity of the measured ǫ (b) iso to the choice of the HT requirement in the range from 100 to 150 GeV, when selecting the bb control sample in data. We allowed the W+jets contribution to the event count in the single fake control region to be twice as large as predicted from MC, with an ǫ (b) iso half/double of the one measured in the bb control sample. Finally, we used the MC to estimate the change in ǫ (b) iso for different b-purities of the bb control sample. We then multiply the sideband yield in each of the four channels µµ, ee, eµ, and µe by the appropriate factor ǫ (b) iso/(1 − ǫ (b) iso) to – 14 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 arrive at 0.52 ± 0.24 ± 0.26 as the estimate of the contribution of events with one fake lepton to the total background. The uncertainties quoted are statistical and systematic, respectively, taking correlations into account. While the BTag-and-probe technique described above is based on different assumptions than the TL method of section 5.1, we note that in their implementation the two techniques are quite similar. Both techniques use events in a RelIso sideband combined with a scale factor determined from an independent control sample to estimate the background in the signal region. The most notable differences are the requirement of the away-jet b-tag, which targets leptons from b decay in the BTag-and-probe method, the choice of variables used to parametrize ǫT/L in one case and ǫ (b) µ iso or ǫ (b) e iso in the other, and the size of the RelIso sideband used in the extrapolation. Combining the background estimates for events with one and two fake leptons and propagating all statistical and systematic uncertainties between channels, including their correlations, we arrive at a final estimate of background due to fake leptons of 0.70±0.23± 0.21 events, with the first (second) uncertainty being statistical (systematic). 5.3 Search using hadronic triggers and τh Simulation studies show clearly that the largest source of background for the τh channels is due to τh fake leptons. We estimate this background using the same ”Tight-Loose” (TL) method as was used for fake leptons in section 5.1, except that for the ”Loose” selection we loosen the τh identification instead of the isolation. To be specific, part of the discrimination between hadronic τ decays and generic QCD jets is based on five neural networks trained to identify different hadronic τ decay modes. The neural network requirements are used for the tight, but not the loose selection. As in section 5.1, in order to predict the number of events from fake τh, we measure the tight-to-loose ratio ǫTL in bins of η and pT . On average, ǫTL = 9.5 ± 0.5%, where the uncertainty is statistical only. This is measured using a single-τh control sample with HT > 300 GeV and Emiss T < 20 GeV. The HT requirement results in hadronic activity similar to our signal region, while the Emiss T requirement reduces contributions from W, Z plus jets, and tt, resulting in a control region that is dominated by QCD multijet production. The expected number of background events is estimated by selecting eτ , µτ , and ττ events where the τ candidates pass the loose selection but fail the tight selection. We find 1, 2, and 2 such events, respectively, in these three channels. For the eτ and µτ channels, these events are weighted by the corresponding factors of ǫTL/(1 − ǫTL), while for the ττ channel each event is weighted by the product of two such factors, one corresponding to each of the two τ leptons in the event. We perform two types of validation of this background estimate. First, we compare observation and prediction in the signal region for simulation. Second, we compare obser- vation and prediction in data after relaxing the HT selection from 350 GeV to 150 GeV, and removing the Emiss T requirement. In simulation the contribution of LM0 represents less than 4% of events with two same-sign isolated leptons and HT > 150 GeV. Table 1 presents both of these validations of the background estimation technique. We find good agreement between the observation and prediction in all channels in data and in the simulation. – 15 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 Simulation Data Only SM Relaxed selection Channel Observed Predicted Observed Predicted ττ 0.08±0.03 0.15±0.15 14 14.0±4.3±2.6 eτ 0.35±0.12 0.30±0.11 1 0.8±0.4±0.1 µτ 0.47±0.15 0.49±0.20 2 2.9±0.6 ±0.4 Table 1. Validation of the TL method. The number of observed events is compared to the number of predicted events in simulation (first two columns) and in a background-dominated control region with relaxed selection criteria (last two columns). The simulation is normalized to 35 pb−1. The first and second uncertainties in the number of predicted events in data are statistical and systematic, respectively. The largest source of systematic uncertainties in the prediction of background events is due to lack of statistics of simulated events to validate the method (30%). In addition, we find uncertainties of 18%, 8% and 7% in τhτh, eτh, and µτh, respectively, due to the correlation of ǫTL with HT . We measure ǫTL for HT > 150 GeV. We determine the systematic uncertainty as the difference in the number of predicted background events from the reference measurement at HT > 300 GeV with the measurement for HT > 150 GeV. From simulation studies it is found that an additional 10% systematic uncertainty must be added in quadrature in the eτh and µτh channels to account for neglecting background contributions from fake electrons and muons. Taking all of this into account, we arrive at an estimate of 0.28 ± 0.14 ± 0.09 events with fake leptons, where the uncertainties are statistical and systematic, respectively. 5.4 Comparison of leptonic and HT -triggered analyses with electrons and muons in the final state As discussed above, we utilize two different trigger strategies to define selections that cover the maximum phase space possible in our search for new physics. In addition, the fact that these two selections have an overlap allows us to perform direct comparisons and cross-checks that we present in this section. We start by defining an overlap preselection requiring one electron or muon of pT >20 GeV and a second with pT > 10 GeV, HT > 300 GeV, and no Emiss T or isolation requirements. We note that this corresponds to the preselection sample from section 5.2 with the pT of the leptons tightened to be consistent with section 5.1. Comparing yields from this selection for the two trigger strategies on an event-by-event basis, we find that all of the HT -triggered events are also present in the lepton-triggered sample. From this comparison, we calculate efficiencies for the HT and leptonic triggers to be (92± 4) % and (100+0 −2) %, respectively. While statistics are limited, this confirms the trigger efficiencies measured in independent data samples, as presented in section 4. Sections 5.1 and 5.2 introduced two alternative methods for estimating the background due to fake leptons. Here we compare the two independent predictions in the region of overlap for the two searches. To compare the two methods in a common signal region, we – 16 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 Figure 7. (Left) The probability to mismeasure the electron charge as a function of η in the pT range 10−100 GeV, as obtained from simulation. (Right) Same-sign ee invariant mass distribution in data compared with the Z → ee expectation from simulation. require Emiss T > 30 GeV and the lepton isolation described in section 5.1, in addition to the preselection defined above. The TL method predicts 0.68±0.39 based on a yield of 3 events that pass the loose selection. The second method introduced in section 5.2 results in an estimate of 0.27 ± 0.12 events based on 11 events in the RelIso sideband. Both of these uncertainties are statistical only. We thus conclude that the two trigger strategies lead to consistent results within the kinematic region where they overlap, and the two methods of estimating backgrounds due to fake leptons give consistent results in that region. 5.5 Electron charge mismeasurement A second potentially important source of background consists of opposite-sign dilepton events (e±e∓ or e±µ∓) where the sign of the charge of one of the electrons is mismeasured because of hard bremsstrahlung in the tracker volume. We measure the electron charge in three different ways. Two of the measurements are based on the reconstructed track from two separate tracking algorithms: the standard CMS track reconstruction algorithm [29, 30] and the Gaussian Sum Filter algorithm [31], optimized for the measurement of electron tracks that radiate in the tracker material. The third measurement is based on the relative position of the calorimeter cluster and the projection to the calorimeter of a line segment built out of hits in the pixel detector. To reduce the effect of charge mismeasurements, we require agreement among the three measurements. After this requirement, the probability of mismeasuring the charge of an electron in simulation is at the level of a few per mille, even in the |η| > 1 region where the amount of material is largest, as can be seen in figure 7 (left). To demonstrate our understanding of this probability, we show in figure 7 (right) the invariant-mass spectrum for same-sign ee events and our simulation-based prediction from Z → ee with one mismeasured charge. The Z sample shown uses the tight electron selection described in section 4.1, except with no jet or HT requirement. Instead, we require Emiss T < – 17 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 20 GeV and transverse mass < 25 GeV to reduce backgrounds from W+ jets. The highest- pT lepton has been used in the calculation of the transverse mass. Measurement of the electron momentum is dominated by the energy measurements in the calorimeter, while the measurement of its charge is dominated by measurements in the tracker. An electron with mismeasured charge will thus still have a correctly measured momentum, leading to the clear Z peak in the e±e± invariant-mass displayed in figure 7 (right). Normalized to 35 pb−1, the simulation predicts 7.4 ± 0.9 events in the Z mass region, with the quoted uncertainty due to statistics. In data, we observe 5 events in the same region. We predict the same-sign Z yield in simulation (6.37 ± 0.03) and data (4.9 ± 0.1) based on reweighting the opposite-sign Z → e±e∓ yield by a simulation-based parametrization of the probability for electron charge mismeasurement as a function of pT and η. For the leptonic trigger searches we estimate the number of background events due to charge misidentification by scaling the opposite-sign yields by the above probability function. We estimate the background due to electron charge misidentification as 0.012 ± 0.002 and 0.04 ± 0.01 for the Emiss T > 80 GeV and HT > 200 GeV regions, respectively. For the HT -triggered searches these backgrounds are further reduced since the opposite-sign yield is smaller given the tighter HT requirement. The resulting background prediction is 0.008 ± 0.005 (ee) and 0.004 ± 0.002 (eµ) events. For the search with τh in the final state, this background is negligible and ignored, as even in the opposite-sign eτh channel, background τh contributions dominate over those with a prompt τh. We assign a 50% systematic uncertainty on the estimated backgrounds due to electron charge mismeasurement. This is motivated by the statistics available in the doubly charged Z → e±e± signal region. 6 Signal acceptance and efficiency systematic uncertainties Electron and muon identification efficiencies above pT ≈ 20 GeV are known at the level of 3% per electron and 1.5% per muon, based on studies of large samples of Z → ee and µµ events in data and simulation. The uncertainties increase as the efficiencies themselves decrease towards lower pT , reaching 6% (8%) per muon (electron) at 5 (10) GeV. In addi- tion, there is a potential mismodelling of the lepton isolation efficiency between data and simulation that grows with the amount of hadronic activity per event. To assess this, we compare the isolation efficiency as a function of track multiplicity in data and simulation for Z → ee and µµ, and extrapolate to new physics signals with large hadronic activity using simulation, as discussed in more detail in section 8. Based on this, we assign an additional 5% systematic uncertainty per lepton. There is also a 1% (5%) uncertainty associated with the lepton (HT ) trigger efficiency. The efficiency of the hadronic τh selection is studied in data via the process Z → ττ , where one τ decays hadronically while the other decays into a muon [13]. The available statistics are an order of magnitude lower than the statistics available in Z → ee or µµ, at significantly lower purity. Accordingly, τh reconstruction versus pT can not be studied at the same level of detail in data as for the electron and muon reconstruction, and we depend – 18 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 to a greater extent on an accurate simulation than we do for electrons and muons. We assign an uncertainty of 30% [15] to the τh selection efficiency to account for the limited statistics available in data to validate the efficiency measured in simulation. An additional source of systematic uncertainty is associated with the current ≈ 5% uncertainty on the hadronic energy scale [19] at CMS. This scale uncertainty limits our understanding of the efficiency of the HT and Emiss T requirements. Clearly, final states where the typical HT and Emiss T are large compared to the selection values used in the analysis are less affected than those with smaller HT and Emiss T . We compute the systematic uncertainty due to this effect for the LM0 benchmark point with the four signal selections using the method of ref. [27]. We use the LM0 model as it is typical of the possible SUSY final states to which these analyses are sensitive. We find that the uncertainty varies between 1% at HT > 60 GeV and 7% at HT > 350 GeV, the values of HT used in the selections for the lepton-triggered baseline and τh search respectively. Uncertainties in the acceptance due to the modelling of initial- and final-state radiation and knowledge of the parton density functions (PDF) are estimated to be 2%. For the latter, we use the CTEQ6.6 [32] PDF and their uncertainties. Based on LM0 as a signal model, we arrive at total uncertainties on signal efficiencies of 12%, 15%, and 30% for the lepton triggered, HT triggered low pT , and HT triggered τh analyses, respectively. This includes a 4% luminosity systematic uncertainty [33]. In addition, to interpret these limits in terms of constraints on new physics models, one needs to take into account any model-dependent theoretical uncertainties. 7 Summary of results The results of our searches are summarized in table 2. The background (BG) predictions are given by the rows labelled ”predicted BG”. In addition to the background estimate from data, we also present an estimate of the background based on simulation in the rows labeled as ”MC”. While QCD multijet production samples are used for testing background estima- tion methods in our control regions, they are too statistically limited to provide meaningful estimates of yields in the signal regions listed in table 2, and are thus not included. All other SM simulation samples described in section 4 are included. Figure 8 summarizes the signal region yields and background composition in all four search regions presented in table 2. The lepton plus jets background where the second lepton candidate is a fake lepton from a jet clearly dominates all search regions. The low-pT -lepton analysis has a small, but non-negligible, background contribution from events with two fake leptons. Estimates for backgrounds due to events with one or two fake leptons were obtained directly from data in appropriately chosen control regions, as described in detail in sections 5.1, 5.2, and 5.3. In the ee and eµ final states, small additional background constributions are present due to the electron charge mismeasurement, as discussed in section 5.5. The remaining irreducible background from two prompt isolated same-sign leptons (WZ, ZZ, ttW, etc.) amounts to at most 10% of the total and is estimated based on theoretical cross section predictions and simulation. Uncertainties on the background prediction include statistical and systematic uncertainties added in quadrature. Contributions estimated with simulation are assigned a – 19 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 Search Region ee µµ eµ total 95% CL UL Yield Lepton Trigger Emiss T > 80 GeV MC 0.05 0.07 0.23 0.35 predicted BG 0.23+0.35 −0.23 0.23+0.26 −0.23 0.74 ± 0.55 1.2 ± 0.8 observed 0 0 0 0 3.1 HT > 200 GeV MC 0.04 0.10 0.17 0.32 predicted BG 0.71 ± 0.58 0.01+0.24 −0.01 0.25+0.27 −0.25 0.97 ± 0.74 observed 0 0 1 1 4.3 HT Trigger Low-pT MC 0.05 0.16 0.21 0.41 predicted BG 0.10 ± 0.07 0.30 ± 0.13 0.40 ± 0.18 0.80 ± 0.31 observed 1 0 0 1 4.4 eτh µτh τhτh total 95% CL UL Yield τh enriched MC 0.36 0.47 0.08 0.91 predicted BG 0.10 ± 0.10 0.17 ± 0.14 0.02 ± 0.01 0.29 ± 0.17 observed 0 0 0 0 3.4 Table 2. Observed and estimated background yields for all analyses. The rows labeled “predicted BG” refer to the sum of the data-driven estimates of the fake lepton contributions, and the residual contributions predicted by the simulation. The rows labeled “MC” refer to the background as predicted from the simulation alone. Rows labeled “observed” show the actual number of events seen in data. The last column (95% CL UL Yield) represents observed upper limits on event yields from new physics. 50% systematic uncertainty. Data-driven estimates are assigned a systematic uncertainty between 30% and 50% across the various signal regions and channels. The ee, eµ, and µµ channels have partially or fully correlated systematic uncertainties, as described in detail in section 5. We see no evidence of an event yield in excess of the background prediction and set 95% CL upper limits (UL) on the number of observed events using a Bayesian method [34] with a flat prior on the signal strength and log-normal priors for efficiency and background uncertainties. These include uncertainties on the signal efficiency of 12% , 15%, and 30% for the lepton triggered, HT triggered low-pT , and HT triggered τh analyses, respectively, as discussed in more detail in section 6. Based on the LM0 benchmark model, the simulation predicts 7.3, 9.6, 9.1, and 2.0 events for the four signal regions, respectively. These LM0 yields are based on the individual NLO cross sections for all production processes that contribute to the expected signal yield. 8 Interpretation of results One of the challenges of signature-based searches is to convey information in a form that can be used to test a variety of specific physics models. In this section we present additional information that can be used to confront models of new physics in an approximate way by – 20 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 0 1 2 3 H Trig taus Lept Trig E >80 E v e n ts observed bkgd with two fake leptons bkgd with one fake lepton opposite-sign leptons and charge mis-identification same-sign leptons Lept Trig H >200 H Trig leptons CMS s =7 TeV, L =35 pb -1 int T miss T T T Figure 8. A visual summary of the observed number of data events, the expected number of background events, and the composition of the background for the four search regions. generator-level simulation studies that compare the expected number of events in 35 pb−1 with our upper limits shown in table 2. The kinematic requirements described in section 4 are the first key ingredients of such studies. The HT variable can be approximated by defining it as the scalar sum of the pT of all final-state quarks (u, d, c, s, and b) and gluons with pT > 30 GeV produced in the hard-scattering process. The Emiss T can be defined as the magnitude of the vector sum of the transverse momentum over all non-interacting particles, e.g., neutrinos and LSP. The ratio of the mean detector responses for HT and Emiss T as defined above, to their true values are 0.94± 0.05, and 0.95± 0.05, respectively, where the uncertainties are dominated by the jet energy scale uncertainty. The resolution on these two quantities differs for the different selections. In addition, the Emiss T resolution depends on the total hadronic activity in the event. It ranges from about 7 to 25 GeV for events with HT in the range of 60 to 350 GeV. The HT resolution decreases from about 26% at 200GeV to 19% for 300 GeV and to 18% for 350 GeV. The HT resolution was measured in simulation using the LM0 reference model, while the Emiss T resolution was measured in data. Figure 9 shows the efficiency versus pT using the LM0 reference model for e, µ (left), and τh (right). Efficiencies here include reconstruction, isolation, and selection. We fit the curves in figure 9 to the functional form: efficiency(pT ) = ǫmax + A× (erf((pT − PTcut)/B) −1 ). We fix PTcut to 10, 5, 15 GeV, and find (A, B, ǫmax) of (0.40, 18, 0.66), (0.32, 18, 0.75), and (0.45, 31, 0.45) for e, µ, and τh, respectively. Lepton isolation efficiencies depend on the hadronic activity in the event, and in some extreme cases like significantly boosted top quarks, on the event topology. The number of charged particles at the generator level after fragmentation and hadronization is a good measure of the first of these two effects, and has been shown to agree reasonably well be- tween data and simulation [35]. We find that the isolation efficiency reduces roughly linearly – 21 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 / ndf 2χ 13.22 / 8 p0 0.0036± 0.7531 p1 0.0108± 0.3155 p2 0.99± 18.05 (GeV)gen T ,e pµ 20 40 60 80 100 120 140 160 180 200 E ff ic ie n c y 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 / ndf 2χ 13.22 / 8 p0 0.0036± 0.7531 p1 0.0108± 0.3155 p2 0.99± 18.05 / ndf 2χ 33.63 / 7 p0 0.0041± 0.6561 p1 0.0098± 0.4017 p2 0.82± 17.61 / ndf 2χ 33.63 / 7 p0 0.0041± 0.6561 p1 0.0098± 0.4017 p2 0.82± 17.61 CMSSM LM0 simulation µ efficiency fitµ e e efficiency fit -1 = 35 pb int = 7 TeV, LsCMS (GeV)gen T pτ 0 20 40 60 80 100 E ff ic ie n c y 0 0.1 0.2 0.3 0.4 0.5 -1 = 35 pb int = 7 TeV, LsCMS CMSSM LM0 simulation τ efficiency fitτ Figure 9. Electron, muon (left) and τh (right) selection efficiencies as a function of pT . The results of the fits described in the text are shown by the dotted lines. by 10% for every 15 charged particles with pT > 3 GeV within the detector acceptance. This was studied using Drell-Yan and LM0 simulation, and compared with Z → ℓ+ℓ− data. The linearity was thus shown to be valid within a range of charged multiplicities from about 10 to 40. The LM0 reference model shown in figure 9 has an average charged multiplicity of ∼ 25 for events that pass our selections. To arrive at the lepton efficiency for a new physics model with an average charged multiplicity of ∼ 40, one would take the efficiency parametrization depicted in figure 9 and multiply it by 0.9 to account for the 10% change in isolation efficiency due to the larger average charged multiplicity. The second effect, i.e., topologies as in boosted top quarks, is difficult to model at generator level, as the isolation efficiency may vary by an order of magnitude or more in extreme cases. Our results are thus not easily interpretable in new physics models with such characteristics without a detailed detector simulation. To show the level of precision obtainable based on such a simple efficiency model, figure 10 compares the exclusion line for CMSSM with tanβ = 3, A0 = 0 GeV, and µ > 0, as obtained from simulation (solid blue) with the corresponding curve (dashed black) for this simple efficiency model. For the exclusion plane, the SUSY particle spectrum is calculated using SoftSUSY [36] along with sparticle decay using SDECAY [37]. The signal events are generated with PYTHIA 6.4.22 [24] using CTEQ6m [38] PDF. The NLO cross sections are obtained with the cross section calculator Prospino [39] at each point in the (m0, m1/2) plane individually. When applying the simple efficiency model, we use the description above for the Emiss T and HT response and resolution, and the efficiency functions displayed in figure 9. No attempt was made here to correct for differences in hadronic event environment as those are small for the relevant (m0, m1/2) points. The width of the red shaded band around the blue line indicates uncertainties in the NLO cross section calculation. It includes variations – 22 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 Figure 10. Exclusion contour in the m0—m1/2 plane for CMSSM as described in the text. Comparing the width of the red shaded band (theoretical uncertainty) around the blue curve with the difference between the solid blue and dashed black curves shows that the imperfections in the simple efficiency model described in the text are small compared to the theoretical uncertainties. of the PDF and simultaneous variation by a factor of two of the renormalization and factor- ization scales. Both effects are added in quadrature. Figure 10 shows that the theoretical uncertainties are larger than the imperfections in the simple efficiency model. The specific limit shown here corresponds to the leptonic trigger result with Emiss T > 80 GeV for the purpose of illustration. The contour separates the bottom-left region where the expected event yield would be larger than the observed limit of 3.1 events (see table 2, first row, last column) from the top-right region where such expected yield would be lower. We choose this CMSSM model because it provides a common reference point to com- pare with previously published Tevatron results [40, 41]. The excluded regions from LEP are based on searches for sleptons and charginos [42–46]. Other final states, especially the all-hadronic [22], are better suited for this standard reference model, while the leptonic same-sign final state explored in this paper is more appropriate to constrain a wide variety of other new physics models [1–7]. As discussed in section 4, in a general supersymmetry context, one might expect gluino-gluino or gluino-squark production to lead to same-sign dilepton events via a decay chain involving a chargino. The salient, and very generic fea- ture here is one lepton per gluino with either sign being equally likely. As a result, 50% of the dilepton events will be same-sign. Moreover, these cascade decays are typically char- acterized by two mass difference scales that separately determine typical HT and lepton – 23 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 pT values. Different models of new physics may thus populate only one or the other of our different search regions. 9 Summary and conclusions Using two different trigger strategies, we have searched for new physics with same-sign dilepton events in the ee, µµ, eµ, eτh, µτh, and τhτh final states, and have seen no evidence for an excess over the background prediction. The τh leptons referred to here are recon- structed via their hadronic decays. The dominant background processes in all final states except τhτh involve events with one fake lepton. In the τhτh final state, events with two fake τh dominate. We have presented methods to derive background estimates from the data for all major background sources. We have set 95% CL upper limits on the number of signal events within |η| < 2.4 at 35 pb−1 in the range of 3.1 to 4.5 events, depending on signal region, and have presented details on signal efficiencies that can be used to confront a wide variety of models of new physics. Our analysis extends the region excluded by experiments at LEP and the Tevatron in the CMSSM model. Acknowledgments We wish to congratulate our colleagues in the CERN accelerator departments for the excel- lent performance of the LHC machine. We thank the technical and administrative staff at CERN and other CMS institutes, and acknowledge support from: FMSR (Austria); FNRS and FWO (Belgium); CNPq, CAPES, FAPERJ, and FAPESP (Brazil); MES (Bulgaria); CERN; CAS, MoST, and NSFC (China); COLCIENCIAS (Colombia); MSES (Croatia); RPF (Cyprus); Academy of Sciences and NICPB (Estonia); Academy of Finland, MEC, and HIP (Finland); CEA and CNRS/IN2P3 (France); BMBF, DFG, and HGF (Germany); GSRT (Greece); OTKA and NKTH (Hungary); DAE and DST (India); IPM (Iran); SFI (Ireland); INFN (Italy); NRF and WCU (Korea); LAS (Lithuania); CINVESTAV, CONA- CYT, SEP, and UASLP-FAI (Mexico); MSI (New Zealand); PAEC (Pakistan); SCSR (Poland); FCT (Portugal); JINR (Armenia, Belarus, Georgia, Ukraine, Uzbekistan); MST and MAE (Russia); MSTD (Serbia); MICINN and CPAN (Spain); Swiss Funding Agencies (Switzerland); NSC (Taipei); TUBITAK and TAEK (Turkey); STFC (United Kingdom); DOE and NSF (USA). Individuals have received support from the Marie-Curie programme and the European Research Council (European Union); the Leventis Foundation; the A. P. Sloan Foundation; the Alexander von Humboldt Foundation; the Associazione per lo Sviluppo Scientifico e Tecnologico del Piemonte (Italy); the Belgian Federal Science Policy Office; the Fonds pour la Formation à la Recherche dans l’Industrie et dans l’Agriculture (FRIA-Belgium); the Agentschap voor Innovatie door Wetenschap en Technologie (IWT-Belgium); and the Council of Science and Industrial Research, India. Open Access. This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited. – 24 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 References [1] R.M. Barnett, J.F. Gunion and H.E. Haber, Discovering supersymmetry with like sign dileptons, Phys. Lett. B 315 (1993) 349 [hep-ph/9306204] [SPIRES]. [2] M. Guchait and D.P. Roy, Like sign dilepton signature for gluino production at CERN LHC including top quark and Higgs boson effects, Phys. Rev. D 52 (1995) 133 [hep-ph/9412329] [SPIRES]. [3] H. Baer, C.-h. Chen, F. Paige and X. Tata, Signals for minimal supergravity at the CERN Large Hadron Collider II: multilepton channels, Phys. Rev. D 53 (1996) 6241 [hep-ph/9512383] [SPIRES]. [4] H.-C. Cheng, K.T. Matchev and M. Schmaltz, Bosonic supersymmetry? Getting fooled at the CERN LHC, Phys. Rev. D 66 (2002) 056006 [hep-ph/0205314] [SPIRES]. [5] R. Contino and G. Servant, Discovering the top partners at the LHC using same-sign dilepton final states, JHEP 06 (2008) 026 [arXiv:0801.1679] [SPIRES]. [6] F.M.L. Almeida, Jr., Y.D.A. Coutinho, J.A. Martins Simoes, P.P. Queiroz Filho and C.M. Porto, Same-sign dileptons as a signature for heavy Majorana neutrinos in hadron-hadron collisions, Phys. Lett. B 400 (1997) 331 [hep-ph/9703441] [SPIRES]. [7] Y. Bai and Z. Han, Top-antitop and top-top resonances in the dilepton channel at the CERN LHC, JHEP 04 (2009) 056 [arXiv:0809.4487] [SPIRES]. [8] G. Bertone, D. Hooper and J. Silk, Particle dark matter: evidence, candidates and constraints, Phys. Rept. 405 (2005) 279 [hep-ph/0404175] [SPIRES]. [9] CMS collaboration, R. Adolphi et al., The CMS experiment at the CERN LHC, 2008 JINST 3 S08004 [SPIRES]. [10] CMS collaboration, Performance of muon identification in pp collisions at √ s = 7TeV, CMS-PAS-MUO-10-002 (2010). [11] CMS collaboration, Electron reconstruction and identification at √ s = 7TeV, CMS-PAS-EGM-10-004. [12] CMS collaboration, V. Khachatryan et al., Measurements of inclusive W and Z cross sections in pp collisions at √ s = 7TeV, JHEP 01 (2011) 080 [arXiv:1012.2466] [SPIRES]. [13] CMS collaboration, Study of τ reconstruction algorithms using pp collisions data collected at√ s = 7TeV with CMS detector at LHC, CMS-PAS-PFT-10-004 (2010). [14] CMS collaboration, Particle-flow event reconstruction in CMS and performance for jets, taus, and Emiss T , CMS-PAS-PFT-09-001 (2009). [15] CMS collaboration, Performance of τ reconstruction algorithms in 2010 data collected with CMS, CMS-PAS-TAU-11-001 (2011). [16] CMS collaboration, Commissioning of the particle-flow reconstruction in minimum-bias and jet events from pp collisions at 7 TeV, CMS-PAS-PFT-10-002 (2010). [17] M. Cacciari, G.P. Salam and G. Soyez, The anti-kt jet clustering algorithm, JHEP 04 (2008) 063 [arXiv:0802.1189] [SPIRES]. [18] CMS collaboration, Jets in 0.9 and 2.36 TeV pp collisions, CMS-PAS-JME-10-001 (2010). [19] CMS collaboration, Jet performance in pp collisions at √ s = 7TeV, CMS-PAS-JME-10-003 (2010). – 25 – http://dx.doi.org/10.1016/0370-2693(93)91623-U http://arxiv.org/abs/hep-ph/9306204 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=HEP-PH/9306204 http://dx.doi.org/10.1103/PhysRevD.52.133 http://arxiv.org/abs/hep-ph/9412329 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=HEP-PH/9412329 http://dx.doi.org/10.1103/PhysRevD.53.6241 http://arxiv.org/abs/hep-ph/9512383 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=HEP-PH/9512383 http://dx.doi.org/10.1103/PhysRevD.66.056006 http://arxiv.org/abs/hep-ph/0205314 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=HEP-PH/0205314 http://dx.doi.org/10.1088/1126-6708/2008/06/026 http://arxiv.org/abs/0801.1679 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=0801.1679 http://dx.doi.org/10.1016/S0370-2693(97)00143-3 http://arxiv.org/abs/hep-ph/9703441 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=HEP-PH/9703441 http://dx.doi.org/10.1088/1126-6708/2009/04/056 http://arxiv.org/abs/0809.4487 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=0809.4487 http://dx.doi.org/10.1016/j.physrep.2004.08.031 http://arxiv.org/abs/hep-ph/0404175 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=HEP-PH/0404175 http://dx.doi.org/10.1088/1748-0221/3/08/S08004 http://dx.doi.org/10.1088/1748-0221/3/08/S08004 http://www-spires.slac.stanford.edu/spires/find/hep/www?j=JINST,3,S08004 http://cdsweb.cern.ch/record/1279140 http://cdsweb.cern.ch/record/1299116 http://dx.doi.org/10.1007/JHEP01(2011)080 http://arxiv.org/abs/1012.2466 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=1012.2466 http://cdsweb.cern.ch/record/1279358 http://cdsweb.cern.ch/record/1194487 http://cdsweb.cern.ch/record/1337004 http://cdsweb.cern.ch/record/1279341 http://dx.doi.org/10.1088/1126-6708/2008/04/063 http://dx.doi.org/10.1088/1126-6708/2008/04/063 http://arxiv.org/abs/0802.1189 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=0802.1189 http://cdsweb.cern.ch/record/1248210 http://cdsweb.cern.ch/record/1279362 J H E P 0 6 ( 2 0 1 1 ) 0 7 7 [20] G.L. Kane, C.F. Kolda, L. Roszkowski and J.D. Wells, Study of constrained minimal supersymmetry, Phys. Rev. D 49 (1994) 6173 [hep-ph/9312272] [SPIRES]. [21] ATLAS collaboration, G. Aad et al., Search for supersymmetry using final states with one lepton, jets and missing transverse momentum with the ATLAS detector in √ s = 7 TeV pp, Phys. Rev. Lett. 106 (2011) 131802 [arXiv:1102.2357] [SPIRES]. [22] CMS collaboration, V. Khachatryan et al., Search for supersymmetry in pp collisions at 7 TeV in events with jets and missing transverse energy, Phys. Lett. B 698 (2011) 196 [arXiv:1101.1628] [SPIRES]. [23] CMS collaboration, S. Chatrchyan et al., Search for physics beyond the standard model in opposite-sign dilepton events at √ s = 7 TeV, arXiv:1103.1348 [SPIRES]. [24] T. Sjöstrand, S. Mrenna and P.Z. Skands, PYTHIA 6.4 physics and manual, JHEP 05 (2006) 026 [hep-ph/0603175] [SPIRES]. [25] J. Alwall et al., MadGraph/MadEvent v4: the new web generation, JHEP 09 (2007) 028 [arXiv:0706.2334] [SPIRES]. [26] GEANT 4 collaboration, GEANT4 — A simulation toolkit, Nucl. Instrum. Meth. A 506 (2003) 250. [27] CMS collaboration, V. Khachatryan et al., First measurement of the cross section for top-quark pair production in proton-proton collisions at √ s = 7TeV, Phys. Lett. B 695 (2011) 424 [arXiv:1010.5994] [SPIRES]. [28] CMS collaboration, Performance of methods for data-driven background estimation in SUSY searches, CMS-PAS-SUS-10-001 (2010). [29] CMS collaboration, Tracking and vertexing results from first collisions, CMS-PAS-TRK-10-001 (2010). [30] W. Adam et al., Track reconstruction in the CMS tracker, CMS-NOTE-2006-041 (2006). [31] W. Adam, R. Fruhwirth, A. Strandlie and T. Todorov, Reconstruction of electrons with the Gaussian-sum filter in the CMS tracker at LHC, ECONF C 0303241 (2003) TULT009 [physics/0306087]. [32] H.-L. Lai et al., Uncertainty induced by QCD coupling in the CTEQ global analysis of parton distributions, Phys. Rev. D 82 (2010) 054021 [arXiv:1004.4624] [SPIRES]. [33] CMS collaboration, Measurement of CMS luminosity, CMS-PAS-EWK-10-004 (2010). [34] Particle Data Group collaboration, K. Nakamura et al., Review of particle physics, J. Phys. G 37 (2010) 075021 [SPIRES]. [35] CMS collaboration, V. Khachatryan et al., Charged particle multiplicities in pp interactions at √ s = 0.9, 2.36 and 7TeV, JHEP 01 (2011) 079 [arXiv:1011.5531] [SPIRES]. [36] B.C. Allanach, SOFTSUSY: a program for calculating supersymmetric spectra, Comput. Phys. Commun. 143 (2002) 305 [hep-ph/0104145] [SPIRES]. [37] M. Muhlleitner, A. Djouadi and Y. Mambrini, SDECAY: a Fortran code for the decays of the supersymmetric particles in the MSSM, Comput. Phys. Commun. 168 (2005) 46 [hep-ph/0311167] [SPIRES]. [38] P.M. Nadolsky et al., Implications of CTEQ global analysis for collider observables, Phys. Rev. D 78 (2008) 013004 [arXiv:0802.0007] [SPIRES]. – 26 – http://dx.doi.org/10.1103/PhysRevD.49.6173 http://arxiv.org/abs/hep-ph/9312272 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=HEP-PH/9312272 http://dx.doi.org/10.1103/PhysRevLett.106.131802 http://arxiv.org/abs/1102.2357 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=1102.2357 http://dx.doi.org/10.1016/j.physletb.2011.03.021 http://arxiv.org/abs/1101.1628 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=1101.1628 http://arxiv.org/abs/1103.1348 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=1103.1348 http://dx.doi.org/10.1088/1126-6708/2006/05/026 http://dx.doi.org/10.1088/1126-6708/2006/05/026 http://arxiv.org/abs/hep-ph/0603175 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=HEP-PH/0603175 http://dx.doi.org/10.1088/1126-6708/2007/09/028 http://arxiv.org/abs/0706.2334 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=0706.2334 http://dx.doi.org/10.1016/S0168-9002(03)01368-8 http://dx.doi.org/10.1016/S0168-9002(03)01368-8 http://dx.doi.org/10.1016/j.physletb.2010.11.058 http://dx.doi.org/10.1016/j.physletb.2010.11.058 http://arxiv.org/abs/1010.5994 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=1010.5994 http://cdsweb.cern.ch/record/1279147 http://cdsweb.cern.ch/record/1258204 http://cdsweb.cern.ch/record/934067 http://dx.doi.org/10.1088/0954-3899/31/9/N01 http://arxiv.org/abs/physics/0306087 http://dx.doi.org/10.1103/PhysRevD.82.054021 http://arxiv.org/abs/1004.4624 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=1004.4624 http://cdsweb.cern.ch/record/1279145 http://dx.doi.org/10.1088/0954-3899/37/7A/075021 http://dx.doi.org/10.1088/0954-3899/37/7A/075021 http://www-spires.slac.stanford.edu/spires/find/hep/www?j=JPHGB,G37,075021 http://dx.doi.org/10.1007/JHEP01(2011)079 http://arxiv.org/abs/1011.5531 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=1011.5531 http://dx.doi.org/10.1016/S0010-4655(01)00460-X http://dx.doi.org/10.1016/S0010-4655(01)00460-X http://arxiv.org/abs/hep-ph/0104145 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=HEP-PH/0104145 http://dx.doi.org/10.1016/j.cpc.2005.01.012 http://arxiv.org/abs/hep-ph/0311167 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=HEP-PH/0311167 http://dx.doi.org/10.1103/PhysRevD.78.013004 http://dx.doi.org/10.1103/PhysRevD.78.013004 http://arxiv.org/abs/0802.0007 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=0802.0007 J H E P 0 6 ( 2 0 1 1 ) 0 7 7 [39] W. Beenakker, R. Hopker, M. Spira and P.M. Zerwas, Squark and gluino production at hadron colliders, Nucl. Phys. B 492 (1997) 51 [hep-ph/9610490] [SPIRES]. [40] CDF collaboration, T. Aaltonen et al., Search for supersymmetry in pp̄ collisions at√ s = 1.96TeV using the trilepton signature of chargino-neutralino production, Phys. Rev. Lett. 101 (2008) 251801 [arXiv:0808.2446] [SPIRES]. [41] D0 collaboration, V.M. Abazov et al., Search for associated production of charginos and neutralinos in the trilepton final state using 2.3 fb−1 of data, Phys. Lett. B 680 (2009) 34 [arXiv:0901.0646] [SPIRES]. [42] ALEPH, DELPHI, L3 and OPAL collaboration, Joint SUSY working group, LEPSUSYWG-02-06-2. [43] ALEPH collaboration, A. Heister et al., Absolute mass lower limit for the lightest neutralino of the MSSM from e+e− data at s1/2 up to 209GeV, Phys. Lett. B 583 (2004) 247 [SPIRES]. [44] DELPHI collaboration, J. Abdallah et al., Searches for supersymmetric particles in e+e− collisions up to 208GeV and interpretation of the results within the MSSM, Eur. Phys. J. C 31 (2003) 421 [hep-ex/0311019] [SPIRES]. [45] L3 collaboration, P. Achard et al., Search for scalar leptons and scalar quarks at LEP, Phys. Lett. B 580 (2004) 37 [hep-ex/0310007] [SPIRES]. [46] OPAL collaboration, G. Abbiendi et al., Search for chargino and neutralino production at√ s = 192GeV to 209GeV at LEP, Eur. Phys. J. C 35 (2004) 1 [hep-ex/0401026] [SPIRES]. – 27 – http://dx.doi.org/10.1016/S0550-3213(97)00084-9 http://arxiv.org/abs/hep-ph/9610490 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=HEP-PH/9610490 http://dx.doi.org/10.1103/PhysRevLett.101.251801 http://dx.doi.org/10.1103/PhysRevLett.101.251801 http://arxiv.org/abs/0808.2446 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=0808.2446 http://dx.doi.org/10.1016/j.physletb.2009.08.011 http://arxiv.org/abs/0901.0646 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=0901.0646 http://dx.doi.org/10.1016/j.physletb.2003.12.066 http://www-spires.slac.stanford.edu/spires/find/hep/www?j=PHLTA,B583,247 http://dx.doi.org/10.1140/epjc/s2003-01355-5 http://dx.doi.org/10.1140/epjc/s2003-01355-5 http://arxiv.org/abs/hep-ex/0311019 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=HEP-EX/0311019 http://dx.doi.org/10.1016/j.physletb.2003.10.010 http://dx.doi.org/10.1016/j.physletb.2003.10.010 http://arxiv.org/abs/hep-ex/0310007 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=HEP-EX/0310007 http://dx.doi.org/10.1140/epjc/s2004-01780-x http://arxiv.org/abs/hep-ex/0401026 http://www-spires.slac.stanford.edu/spires/find/hep/www?eprint=HEP-EX/0401026 J H E P 0 6 ( 2 0 1 1 ) 0 7 7 The CMS collaboration Yerevan Physics Institute, Yerevan, Armenia S. Chatrchyan, V. Khachatryan, A.M. Sirunyan, A. Tumasyan Institut für Hochenergiephysik der OeAW, Wien, Austria W. Adam, T. Bergauer, M. Dragicevic, J. Erö, C. Fabjan, M. Friedl, R. Frühwirth, V.M. Ghete, J. Hammer1, S. Hänsel, M. Hoch, N. Hörmann, J. Hrubec, M. Jeitler, G. Kasieczka, W. Kiesenhofer, M. Krammer, D. Liko, I. Mikulec, M. Pernicka, H. Rohringer, R. Schöfbeck, J. Strauss, F. Teischinger, P. Wagner, W. Waltenberger, G. Walzel, E. Widl, C.-E. Wulz National Centre for Particle and High Energy Physics, Minsk, Belarus V. Mossolov, N. Shumeiko, J. Suarez Gonzalez Universiteit Antwerpen, Antwerpen, Belgium L. Benucci, E.A. De Wolf, X. Janssen, T. Maes, L. Mucibello, S. Ochesanu, B. Roland, R. Rougny, M. Selvaggi, H. Van Haevermaet, P. Van Mechelen, N. Van Remortel Vrije Universiteit Brussel, Brussel, Belgium F. Blekman, S. Blyweert, J. D’Hondt, O. Devroede, R. Gonzalez Suarez, A. Kalogeropoulos, J. Maes, M. Maes, W. Van Doninck, P. Van Mulders, G.P. Van Onsem, I. Villella Université Libre de Bruxelles, Bruxelles, Belgium O. Charaf, B. Clerbaux, G. De Lentdecker, V. Dero, A.P.R. Gay, G.H. Hammad, T. Hreus, P.E. Marage, L. Thomas, C. Vander Velde, P. Vanlaer Ghent University, Ghent, Belgium V. Adler, A. Cimmino, S. Costantini, M. Grunewald, B. Klein, J. Lellouch, A. Marinov, J. Mccartin, D. Ryckbosch, F. Thyssen, M. Tytgat, L. Vanelderen, P. Verwilligen, S. Walsh, N. Zaganidis Université Catholique de Louvain, Louvain-la-Neuve, Belgium S. Basegmez, G. Bruno, J. Caudron, L. Ceard, E. Cortina Gil, J. De Favereau De Jeneret, C. Delaere1, D. Favart, A. Giammanco, G. Grégoire, J. Hollar, V. Lemaitre, J. Liao, O. Militaru, S. Ovyn, D. Pagano, A. Pin, K. Piotrzkowski, N. Schul Université de Mons, Mons, Belgium N. Beliy, T. Caebergs, E. Daubie Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, Brazil G.A. Alves, D. De Jesus Damiao, M.E. Pol, M.H.G. Souza – 28 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 Universidade do Estado do Rio de Janeiro, Rio de Janeiro, Brazil W. Carvalho, E.M. Da Costa, C. De Oliveira Martins, S. Fonseca De Souza, L. Mundim, H. Nogima, V. Oguri, W.L. Prado Da Silva, A. Santoro, S.M. Silva Do Amaral, A. Sznajder, F. Torres Da Silva De Araujo Instituto de Fisica Teorica, Universidade Estadual Paulista, Sao Paulo, Brazil F.A. Dias, T.R. Fernandez Perez Tomei, E. M. Gregores2, C. Lagana, F. Marinho, P.G. Mercadante2, S.F. Novaes, Sandra S. Padula Institute for Nuclear Research and Nuclear Energy, Sofia, Bulgaria N. Darmenov1, L. Dimitrov, V. Genchev1, P. Iaydjiev1, S. Piperov, M. Rodozov, S. Stoykova, G. Sultanov, V. Tcholakov, R. Trayanov, I. Vankov University of Sofia, Sofia, Bulgaria A. Dimitrov, R. Hadjiiska, A. Karadzhinova, V. Kozhuharov, L. Litov, M. Mateev, B. Pavlov, P. Petkov Institute of High Energy Physics, Beijing, China J.G. Bian, G.M. Chen, H.S. Chen, C.H. Jiang, D. Liang, S. Liang, X. Meng, J. Tao, J. Wang, J. Wang, X. Wang, Z. Wang, H. Xiao, M. Xu, J. Zang, Z. Zhang State Key Lab. of Nucl. Phys. and Tech., Peking University, Beijing, China Y. Ban, S. Guo, Y. Guo, W. Li, Y. Mao, S.J. Qian, H. Teng, L. Zhang, B. Zhu, W. Zou Universidad de Los Andes, Bogota, Colombia A. Cabrera, B. Gomez Moreno, A.A. Ocampo Rios, A.F. Osorio Oliveros, J.C. Sanabria Technical University of Split, Split, Croatia N. Godinovic, D. Lelas, K. Lelas, R. Plestina3, D. Polic, I. Puljak University of Split, Split, Croatia Z. Antunovic, M. Dzelalija Institute Rudjer Boskovic, Zagreb, Croatia V. Brigljevic, S. Duric, K. Kadija, S. Morovic University of Cyprus, Nicosia, Cyprus A. Attikis, M. Galanti, J. Mousa, C. Nicolaou, F. Ptochos, P.A. Razis Charles University, Prague, Czech Republic M. Finger, M. Finger Jr. – 29 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 Academy of Scientific Research and Technology of the Arab Republic of Egypt, Egyptian Network of High Energy Physics, Cairo, Egypt Y. Assran4, S. Khalil5, M.A. Mahmoud6 National Institute of Chemical Physics and Biophysics, Tallinn, Estonia A. Hektor, M. Kadastik, M. Müntel, M. Raidal, L. Rebane Department of Physics, University of Helsinki, Helsinki, Finland V. Azzolini, P. Eerola, G. Fedi Helsinki Institute of Physics, Helsinki, Finland S. Czellar, J. Härkönen, A. Heikkinen, V. Karimäki, R. Kinnunen, M.J. Kortelainen, T. Lampén, K. Lassila-Perini, S. Lehti, T. Lindén, P. Luukka, T. Mäenpää, E. Tuominen, J. Tuominiemi, E. Tuovinen, D. Ungaro, L. Wendland Lappeenranta University of Technology, Lappeenranta, Finland K. Banzuzi, A. Korpela, T. Tuuva Laboratoire d’Annecy-le-Vieux de Physique des Particules, IN2P3-CNRS, Annecy-le-Vieux, France D. Sillou DSM/IRFU, CEA/Saclay, Gif-sur-Yvette, France M. Besancon, S. Choudhury, M. Dejardin, D. Denegri, B. Fabbro, J.L. Faure, F. Ferri, S. Ganjour, F.X. Gentit, A. Givernaud, P. Gras, G. Hamel de Monchenault, P. Jarry, E. Locci, J. Malcles, M. Marionneau, L. Millischer, J. Rander, A. Rosowsky, I. Shreyber, M. Titov, P. Verrecchia Laboratoire Leprince-Ringuet, Ecole Polytechnique, IN2P3-CNRS, Palaiseau, France S. Baffioni, F. Beaudette, L. Benhabib, L. Bianchini, M. Bluj7, C. Broutin, P. Busson, C. Charlot, T. Dahms, L. Dobrzynski, S. Elgammal, R. Granier de Cassagnac, M. Hague- nauer, P. Miné, C. Mironov, C. Ochando, P. Paganini, D. Sabes, R. Salerno, Y. Sirois, C. Thiebaux, B. Wyslouch8, A. Zabi Institut Pluridisciplinaire Hubert Curien, Université de Strasbourg, Univer- sité de Haute Alsace Mulhouse, CNRS/IN2P3, Strasbourg, France J.-L. Agram9, J. Andrea, D. Bloch, D. Bodin, J.-M. Brom, M. Cardaci, E.C. Chabert, C. Collard, E. Conte9, F. Drouhin9, C. Ferro, J.-C. Fontaine9, D. Gelé, U. Goerlach, S. Greder, P. Juillot, M. Karim9, A.-C. Le Bihan, Y. Mikami, P. Van Hove – 30 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 Centre de Calcul de l’Institut National de Physique Nucleaire et de Physique des Particules (IN2P3), Villeurbanne, France F. Fassi, D. Mercier Université de Lyon, Université Claude Bernard Lyon 1, CNRS-IN2P3, Institut de Physique Nucléaire de Lyon, Villeurbanne, France C. Baty, S. Beauceron, N. Beaupere, M. Bedjidian, O. Bondu, G. Boudoul, D. Boumediene, H. Brun, R. Chierici, D. Contardo, P. Depasse, H. El Mamouni, J. Fay, S. Gascon, B. Ille, T. Kurca, T. Le Grand, M. Lethuillier, L. Mirabito, S. Perries, V. Sordini, S. Tosi, Y. Tschudi, P. Verdier Institute of High Energy Physics and Informatization, Tbilisi State University, Tbilisi, Georgia D. Lomidze RWTH Aachen University, I. Physikalisches Institut, Aachen, Germany G. Anagnostou, M. Edelhoff, L. Feld, N. Heracleous, O. Hindrichs, R. Jussen, K. Klein, J. Merz, N. Mohr, A. Ostapchuk, A. Perieanu, F. Raupach, J. Sammet, S. Schael, D. Sprenger, H. Weber, M. Weber, B. Wittmer RWTH Aachen University, III. Physikalisches Institut A, Aachen, Germany M. Ata, W. Bender, E. Dietz-Laursonn, M. Erdmann, J. Frangenheim, T. Hebbeker, A. Hinzmann, K. Hoepfner, T. Klimkovich, D. Klingebiel, P. Kreuzer, D. Lanske†, C. Magass, M. Merschmeyer, A. Meyer, P. Papacz, H. Pieta, H. Reithler, S.A. Schmitz, L. Sonnenschein, J. Steggemann, D. Teyssier, M. Tonutti RWTH Aachen University, III. Physikalisches Institut B, Aachen, Germany M. Bontenackels, M. Davids, M. Duda, G. Flügge, H. Geenen, M. Giffels, W. Haj Ahmad, D. Heydhausen, T. Kress, Y. Kuessel, A. Linn, A. Nowack, L. Perchalla, O. Pooth, J. Rennefeld, P. Sauerland, A. Stahl, M. Thomas, D. Tornier, M.H. Zoeller Deutsches Elektronen-Synchrotron, Hamburg, Germany M. Aldaya Martin, W. Behrenhoff, U. Behrens, M. Bergholz10, K. Borras, A. Cakir, A. Campbell, E. Castro, D. Dammann, G. Eckerlin, D. Eckstein, A. Flossdorf, G. Flucke, A. Geiser, J. Hauk, H. Jung1, M. Kasemann, I. Katkov11, P. Katsas, C. Kleinwort, H. Kluge, A. Knutsson, M. Krämer, D. Krücker, E. Kuznetsova, W. Lange, W. Lohmann10, R. Mankel, M. Marienfeld, I.-A. Melzer-Pellmann, A.B. Meyer, J. Mnich, A. Mussgiller, J. Olzem, D. Pitzl, A. Raspereza, A. Raval, M. Rosin, R. Schmidt10, T. Schoerner-Sadenius, N. Sen, A. Spiridonov, M. Stein, J. Tomaszewska, R. Walsh, C. Wissing – 31 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 University of Hamburg, Hamburg, Germany C. Autermann, V. Blobel, S. Bobrovskyi, J. Draeger, H. Enderle, U. Gebbert, K. Kaschube, G. Kaussen, R. Klanner, J. Lange, B. Mura, S. Naumann-Emme, F. Nowak, N. Pietsch, C. Sander, H. Schettler, P. Schleper, M. Schröder, T. Schum, J. Schwandt, H. Stadie, G. Steinbrück, J. Thomsen Institut für Experimentelle Kernphysik, Karlsruhe, Germany C. Barth, J. Bauer, V. Buege, T. Chwalek, W. De Boer, A. Dierlamm, G. Dirkes, M. Feindt, J. Gruschke, C. Hackstein, F. Hartmann, M. Heinrich, H. Held, K.H. Hoffmann, S. Honc, J.R. Komaragiri, T. Kuhr, D. Martschei, S. Mueller, Th. Müller, M. Niegel, O. Oberst, A. Oehler, J. Ott, T. Peiffer, D. Piparo, G. Quast, K. Rabbertz, F. Ratnikov, N. Ratnikova, M. Renz, C. Saout, A. Scheurer, P. Schieferdecker, F.-P. Schilling, M. Schmanau, G. Schott, H.J. Simonis, F.M. Stober, D. Troendle, J. Wagner-Kuhr, T. Weiler, M. Zeise, V. Zhukov11, E.B. Ziebarth Institute of Nuclear Physics ”Demokritos”, Aghia Paraskevi, Greece G. Daskalakis, T. Geralis, K. Karafasoulis, S. Kesisoglou, A. Kyriakis, D. Loukas, I. Manolakos, A. Markou, C. Markou, C. Mavrommatis, E. Ntomari, E. Petrakou University of Athens, Athens, Greece L. Gouskos, T.J. Mertzimekis, A. Panagiotou, E. Stiliaris University of Ioánnina, Ioánnina, Greece I. Evangelou, C. Foudas, P. Kokkas, N. Manthos, I. Papadopoulos, V. Patras, F.A. Triantis KFKI Research Institute for Particle and Nuclear Physics, Budapest, Hungary A. Aranyi, G. Bencze, L. Boldizsar, C. Hajdu1, P. Hidas, D. Horvath12, A. Kapusi, K. Krajczar13, F. Sikler1, G.I. Veres13, G. Vesztergombi13 Institute of Nuclear Research ATOMKI, Debrecen, Hungary N. Beni, J. Molnar, J. Palinkas, Z. Szillasi, V. Veszpremi University of Debrecen, Debrecen, Hungary P. Raics, Z.L. Trocsanyi, B. Ujvari Panjab University, Chandigarh, India S. Bansal, S.B. Beri, V. Bhatnagar, N. Dhingra, R. Gupta, M. Jindal, M. Kaur, J.M. Kohli, M.Z. Mehta, N. Nishu, L.K. Saini, A. Sharma, A.P. Singh, J.B. Singh, S.P. Singh University of Delhi, Delhi, India S. Ahuja, S. Bhattacharya, B.C. Choudhary, P. Gupta, S. Jain, S. Jain, A. Kumar, K. Ranjan, R.K. Shivpuri – 32 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 Bhabha Atomic Research Centre, Mumbai, India R.K. Choudhury, D. Dutta, S. Kailas, V. Kumar, A.K. Mohanty1, L.M. Pant, P. Shukla Tata Institute of Fundamental Research - EHEP, Mumbai, India T. Aziz, M. Guchait14, A. Gurtu, M. Maity15, D. Majumder, G. Majumder, K. Mazumdar, G.B. Mohanty, A. Saha, K. Sudhakar, N. Wickramage Tata Institute of Fundamental Research - HECR, Mumbai, India S. Banerjee, S. Dugad, N.K. Mondal Institute for Research and Fundamental Sciences (IPM), Tehran, Iran H. Arfaei, H. Bakhshiansohi16, S.M. Etesami, A. Fahim16, M. Hashemi, A. Jafari16, M. Khakzad, A. Mohammadi17, M. Mohammadi Najafabadi, S. Paktinat Mehdiabadi, B. Safarzadeh, M. Zeinali18 INFN Sezione di Bari a, Università di Bari b, Politecnico di Bari c, Bari, Italy M. Abbresciaa,b, L. Barbonea,b, C. Calabriaa,b, A. Colaleoa, D. Creanzaa,c, N. De Filippisa,c,1, M. De Palmaa,b, L. Fiorea, G. Iasellia,c, L. Lusitoa,b, G. Maggia,c, M. Maggia, N. Mannaa,b, B. Marangellia,b, S. Mya,c, S. Nuzzoa,b, N. Pacificoa,b, G.A. Pierroa, A. Pompilia,b, G. Pugliesea,c, F. Romanoa,c, G. Rosellia,b, G. Selvaggia,b, L. Silvestrisa, R. Trentaduea, S. Tupputia,b, G. Zitoa INFN Sezione di Bologna a, Università di Bologna b, Bologna, Italy G. Abbiendia, A.C. Benvenutia, D. Bonacorsia, S. Braibant-Giacomellia,b, L. Brigliadoria, P. Capiluppia,b, A. Castroa,b, F.R. Cavalloa, M. Cuffiania,b, G.M. Dallavallea, F. Fabbria, A. Fanfania,b, D. Fasanellaa, P. Giacomellia, M. Giuntaa, S. Marcellinia, G. Masetti, M. Meneghellia,b, A. Montanaria, F.L. Navarriaa,b, F. Odoricia, A. Perrottaa, F. Primaveraa, A.M. Rossia,b, T. Rovellia,b, G. Sirolia,b, R. Travaglinia,b INFN Sezione di Catania a, Università di Catania b, Catania, Italy S. Albergoa,b, G. Cappelloa,b, M. Chiorbolia,b,1, S. Costaa,b, A. Tricomia,b, C. Tuvea INFN Sezione di Firenze a, Università di Firenze b, Firenze, Italy G. Barbaglia, V. Ciullia,b, C. Civininia, R. D’Alessandroa,b, E. Focardia,b, S. Frosalia,b, E. Galloa, S. Gonzia,b, P. Lenzia,b, M. Meschinia, S. Paolettia, G. Sguazzonia, A. Tropianoa,1 INFN Laboratori Nazionali di Frascati, Frascati, Italy L. Benussi, S. Bianco, S. Colafranceschi19, F. Fabbri, D. Piccolo INFN Sezione di Genova, Genova, Italy P. Fabbricatore, R. Musenich – 33 – J H E P 0 6 ( 2 0 1 1 ) 0 7 7 INFN Sezione di Milano-Biccoca a, Università di Milano-Bicocca b, Milano, Italy A. Benagliaa,b, F. De Guioa,b,1, L. Di Matteoa,b, A. Ghezzia,b, M. Malbertia,b, S. Malvezzia, A. Martellia,b, A. Massironia,b, D. Menascea, L. Moronia, M. Paganonia,b, D. Pedrinia, S. Ragazzia,b, N. Redaellia, S. Salaa, T. Tabarelli de Fatisa,b, V. Tancinia,b INFN Sezione di Napoli a, Università di Napoli ”Federico II” b, Napoli, Italy S. Buontempoa, C.A. Carrillo Montoyaa,1, N. Cavalloa,20, A. De Cosaa,b, F. Fabozzia,20, A.O.M. Iorioa,1, L. Listaa, M. Merolaa,b, P. Paoluccia INFN Sezione di Padova a, Università di Padova b, Università di Tren