EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) CERN-EP/2016-097 2017/01/04 CMS-SUS-14-003 Searches for R-parity-violating supersymmetry in pp collisions at √ s = 8 TeV in final states with 0–4 leptons The CMS Collaboration∗ Abstract Results are presented from searches for R-parity-violating supersymmetry in events produced in pp collisions at √ s = 8 TeV at the LHC. Final states with 0, 1, 2, or multi- ple leptons are considered independently. The analysis is performed on data collected by the CMS experiment corresponding to an integrated luminosity of 19.5 fb−1. No excesses of events above the standard model expectations are observed, and 95% con- fidence level limits are set on supersymmetric particle masses and production cross sections. The results are interpreted in models featuring R-parity-violating decays of the lightest supersymmetric particle, which in the studied scenarios can be either the gluino, a bottom squark, or a neutralino. In a gluino pair production model with baryon number violation, gluinos with a mass less than 0.98 and 1.03 TeV are ex- cluded, by analyses in a fully hadronic and one-lepton final state, respectively. An analysis in a dilepton final state is used to exclude bottom squarks with masses less than 307 GeV in a model considering bottom squark pair production. Multilepton fi- nal states are considered in the context of either strong or electroweak production of superpartners and are used to set limits on the masses of the lightest supersymmetric particles. These limits range from 300 to 900 GeV in models with leptonic and up to approximately 700 GeV in models with semileptonic R-parity-violating couplings. Published in Physical Review D as doi:10.1103/PhysRevD.94.112009. c© 2017 CERN for the benefit of the CMS Collaboration. CC-BY-3.0 license ∗See Appendix A for the list of collaboration members ar X iv :1 60 6. 08 07 6v 2 [h ep -ex ] 3 J an 20 17 11 Introduction Supersymmetry (SUSY) is an attractive extension of the standard model (SM) because it can solve the hierarchy problem and can ensure gauge coupling unification [1, 2]. The majority of the searches for SUSY focus on R-parity-conserving (RPC) models. The R-parity of a particle is defined by R = (−1)3B+L+2s, where B and L are its baryon and lepton numbers, respectively, and s is the particle spin [3]. In RPC SUSY, the lightest supersymmetric particle (LSP) is stable, which ensures proton stability and provides a dark matter candidate. All SM particles have R = +1; SUSY posits the existence of a superpartner with R = −1 corresponding to each SM particle. The most general gauge-invariant and renormalizable superpotential violates R-parity, and so R-parity violation is expected unless it is forbidden by some symmetry. Supersymmetric models with R-parity-violating (RPV) interactions can break baryon or lepton number conser- vation [4, 5]. The superpotential WRPV includes a bilinear term proportional to the coupling µ′i and three trilinear terms parameterized by the couplings λijk, λ′ijk, and λ ′′ ijk: WRPV = 1 2 λijkLiLjEk + λ′ijkLiQjDk + 1 2 λ′′ijkUiDjDk + µ ′ iHuLi, (1) where i, j, and k are generation indices; L, Q and Hu are the lepton, quark, and up-type Higgs SU(2)L doublet superfields, respectively; and E, D, and U are the charged lepton, down-type quark, and up-type quark SU(2)L singlet superfields, respectively. The third term violates baryon number conservation, while the other terms violate lepton number conservation. The final term, involving the lepton and up-type Higgs doublets, is also allowed in the superpoten- tial but the effects of this term are not considered in this analysis. Experimental bounds on leptonic, semileptonic, and hadronic RPV couplings are complemen- tary due to the strong constraint on the product of RPV couplings from nucleon stability mea- surements. For example, for squark masses of 1 TeV, stringent experimental limits on proton decay result in the constraint |λ′ijkλ∗i′ j′k′ | < O(10−9) for all generation indices [6]. Much stronger (by a factor of up to ≈ 1018 [4]) constraints are possible for couplings involving light genera- tions, and similar constraints exist for products of other RPV couplings. A subset of RPV scenarios focus on the RPV extension of the minimal supersymmetric model (MSSM) when the assumption of minimal flavor violation (MFV) is imposed [7]. Under this assumption, the only sources of R-parity violation are the SM Yukawa couplings, and the RPV couplings are therefore related to the components of the Cabibbo–Kobayashi–Maskawa matrix and the fermion masses. In some of these models, λ′′332 is the largest RPV coupling [8] and will be a focus of the searches involving hadronic R-parity violation. The missing transverse momentum vector ~pmissT is defined as the negative of the vector sum of momenta in the transverse direction. Its magnitude EmissT is often used in searches for RPC SUSY as the LSP is stable and leaves the detector undetected, leading to large values of EmissT . In the RPV models considered, the LSP decays promptly to SM particles and therefore no large EmissT is expected. Instead, we employ a variety of methods to search for the different types of RPV decays. We search for hadronic RPV SUSY, which arises when any of the λ′′ijk are nonzero, in events with zero or one lepton using the jet and b-tagged jet multiplicities of the event, and in dilep- ton events by means of a kinematic fit to reconstruct the bottom squark mass. We focus on couplings that involve top quarks, as motivated by MFV; the leptons in the final state are the result of leptonic W decay. The results of the fully hadronic and one-lepton analyses are inter- 2 2 CMS detector and reconstruction preted in a model in which the gluino decays via g˜ → t˜t, followed by decay of the top squark via a nonzero λ′′323 coupling: t˜ → bs (charge conjugate reactions are implied throughout this paper). Here the top squark is considered to be much heavier than the gluino, resulting in an effective three-body decay of the gluino. The analysis of the dilepton final state considers pair production of bottom squarks, which decay to a top quark and either a down or a strange quark. To search for leptonic and semileptonic RPV SUSY, which arise when λijk and λ′ijk respectively, are nonzero, we examine events with three or more leptons, binned in the multiplicity of recon- structed objects. Both strong and electroweak production of superpartners are considered. An analysis of a four-lepton final state targets production of squarks and gluinos in which the light- est neutralino decays to final states with electrons or muons. We also study final states with at least three leptons, which are sensitive to electroweak production of winos and Higgsinos. In all analyses considered, the LSP is assumed to decay promptly, meaning the decay vertex is indistinguishable from the primary interaction. This generally implies λ > 10−6 [3]. All RPV couplings are assumed to be zero, except for the specific coupling under study. In this paper, we present the results of these searches with interpretations in a variety of dif- ferent RPV models. The data set was recorded with the CMS detector at the CERN LHC in proton-proton collisions at a center-of-mass energy of 8 TeV and corresponds to an integrated luminosity of 19.3–19.5 fb−1. Searches for multijet resonances, a prominent signal when hadronic RPV SUSY is present, have been performed by CDF [9], ATLAS [10], and CMS [11–13]. The ATLAS Collaboration has also performed a search for RPV SUSY in high-multiplicity events [14]. Searches for RPV inter- actions in multilepton final states have been carried out at LEP [15–17], the Tevatron [18–20], HERA [21, 22], and the LHC [23–28]. The paper is organized as follows. Section 2 presents an overview of the CMS detector, and a description of simulated signal and background samples is given in Section 3. The analyses described in this paper use a common limit-setting procedure. This procedure, as well as the treatment of signal samples, is described in Section 4. The event selections that are common to all analyses in this paper are described in Section 5. The searches reported in this paper cover a wide range of signatures induced by the various RPV couplings. Sections 6 and 7 detail searches for hadronic R-parity violation in zero- and one-lepton final states, respectively. Section 8 describes a search for bottom squarks that decay to a top quark and either a down or a strange quark. Finally, searches for R-parity violation induced by leptonic and semileptonic RPV couplings in multilepton final states are described in Sections 9 and 10, respectively. The results of this paper are summarized in Section 11. 2 CMS detector and reconstruction A detailed description of the CMS detector, together with a definition of the coordinate system used, can be found in Ref. [29]. The central feature of the CMS apparatus is a superconduct- ing solenoid with an internal diameter of 6 m, which generates a 3.8 T uniform magnetic field along the axis of the LHC beams. A silicon pixel and strip tracker, a crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter (HCAL) are located within the magnet. Muons are identified and measured in gas-ionization detectors embedded in the outer steel magnetic flux-return yoke of the solenoid. The silicon tracker, the muon system, and the barrel and endcap calorimeters cover the pseudorapidity ranges |η| < 2.5, |η| < 2.4, 3and |η| < 3.0, respectively. The first level of the CMS trigger system, composed of custom hardware processors, uses information from the calorimeters and muon detectors to select the events most relevant for analysis in a fixed time interval of less than 4 µs. The high-level trigger processor farm further decreases the event rate from around 100 kHz to less than 1 kHz, before data storage. The particle-flow event algorithm [30–32] reconstructs and identifies individual particles with an optimized combination of information from the various elements of the CMS detector. The energy of photons is directly obtained from the ECAL measurement. The energy of electrons is determined from a combination of the electron momentum at the primary interaction vertex as determined by the tracker, the energy of the corresponding ECAL cluster, and the energy sum of all bremsstrahlung photons spatially compatible with originating from the electron track. The energy of muons is obtained from the curvature of the corresponding track. The energy of charged hadrons is determined from a combination of their momentum measured in the tracker and the matching ECAL and HCAL energy deposits. Finally, the energy of neutral hadrons is obtained from the corresponding corrected ECAL and HCAL energy. 3 Simulation The Monte Carlo (MC) simulation is used to estimate some of the SM backgrounds and to understand the efficiency of the signal models, including geometrical acceptance. The SM background samples are generated using MADGRAPH 5.1.3.30 [33], with parton showering and fragmentation modeled using PYTHIA (version 6.420) [34], and passed through a GEANT4- based [35] representation of the CMS detector. QCD multijet samples are generated with up to four partons in the matrix element and tt +jets events are generated with up to three extra partons in the matrix element. The parton shower is matched to the matrix element with the MLM prescription [36]. Signal samples [37] are generated with MADGRAPH and PYTHIA. Most of these samples are then passed through the CMS fast-simulation package [38]; the others are simulated with the same full simulation used for background processes. The CTEQ6L1 [39] set of parton distribution functions (PDFs) is used throughout. Background yields, when taken from simulation, are normalized to next-to-leading-order (NLO) or next-to-next-to-leading- order (NNLO) cross sections, when available. There is an additional uncertainty of 2.6% in these yields due to the imperfect knowledge of the integrated luminosity [40]. The modeling of multiple proton-proton primary interactions in a single bunch crossing, referred to as pileup, is corrected so that the pileup profile matches that of the data. In the data set used in this paper, the mean number of interactions per bunch crossing is 21. 4 Common procedures for signal samples and limits The analyses described in this paper use many shared procedures for the modeling of signal components and the setting of cross section limits, which we describe in this section. The analyses are interpreted in simplified models of SUSY [41–43]. In all of the interpretations, supersymmetric particles not explicitly considered are assumed to have very large masses so that their effect is negligible. However, the masses of intermediate states are assumed to be small enough that all supersymmetric particles decay promptly. The separate analyses probe different sectors of RPV parameter space, and the models used to interpret the results vary accordingly. The hadronic and one-lepton analyses are sensitive mainly to RPV SUSY in the hadronic sector and therefore assume that λ ′′ 332 6= 0 and all other 4 5 Object selection RPV couplings are zero; the experimental signature of a nonzero λ ′′ 331 coupling is identical as there is no discrimination between s and d quarks. Similarly, the two-lepton analysis assumes λ ′′ 332 6= 0 or λ ′′ 331 6= 0, with no other nonzero RPV couplings. As the multilepton searches analyze different lepton flavors separately, they are interpreted in terms of several models with nonzero lepton-flavor-violating couplings, λijk 6= 0 or λ′ijk 6= 0, for several values of i, j, and k. The uncertainty in the knowledge of the PDFs is obtained by applying the envelope prescrip- tion of the PDF4LHC working group [44, 45] with three different PDF sets (CTEQ6.6 [46], MSTW2008nlo68cl [47] and NNPDF2.0-100 [48]). Scale factors are applied to the fast-simulation samples, as a function of the transverse momentum pT and |η|, to reproduce the b jet identifica- tion and misidentification efficiencies obtained from a full simulation of the CMS detector. The uncertainty in the modeling of initial- and final-state radiation (ISR and FSR, respectively) is obtained from the discrepancies between data and simulation observed in Z+jets, diboson+jets, and tt+jets events as a function of the pT of the system recoiling against the ISR jets [49]. Limits are calculated using the CLs [50, 51] method. The LHC-style test statistic is used, within the formalism developed by the CMS and ATLAS collaborations in the context of the LHC Higgs Combination Group [52]. For each cross section σ being tested, the likelihood is pro- filed with respect to the nuisance parameters; that is, the nuisance parameters are treated as fit parameters subject to external constraints on their magnitude and distribution. We find the one-sided p-value of the observed data in the signal-plus-background hypothesis, denoted pσ. This is the fraction of pseudoexperiments with test statistic λp(σ) less than the value measured in data. We also generate pseudoexperiments with the signal cross section set to zero to con- struct the distribution of λp(σ) under the background-only hypothesis. From this distribution we obtain the p-value of data in the background-only hypothesis, denoted p0. Then CLs is de- fined as pσ/(1− p0). If CLs < 0.05, that value of σ is deemed to be excluded at a 95% confidence level (CL). The largest cross section not excluded corresponds to the CLs upper limit. Cross sections for SUSY signal processes, calculated at NLO with next-to-leading-log (NLL) resummation, are taken from the LHC SUSY Cross Sections Working Group [53–57]. To ac- count for theoretical uncertainties conservatively, mass exclusions are quoted using a signal production cross section that is reduced from the nominal value by the amount of the theoreti- cal uncertainty. 5 Object selection Electrons and muons are reconstructed using the tracker, calorimeter, and muon systems. De- tails of the reconstruction and identification for electrons and muons can be found in Ref. [58] and Ref. [59], respectively. In the leptonic analyses, we require that at least one electron or muon in each event has pT > 20 GeV. Additional electrons and muons must have pT > 10 GeV and all leptons must be within |η| < 2.4. The majority of hadronic decays of τ leptons (τh) yield either a single charged particle (one- prong) or three charged particles (three-prong), with or without additional electromagnetic en- ergy from neutral pion decays. We use one- and three-prong τh candidates with pT > 20 GeV, reconstructed with the “hadron plus strips” method [60], which has an efficiency of approxi- mately 70%. Leptons produced in τ decays are included with other electrons and muons. To ensure that electrons, muons, and τh candidates are isolated, we use the variable ET, cone, defined as the transverse energy in a cone of radius ∆R ≡ √ (∆η)2 + (∆φ)2 = 0.3 around the candidate, excluding the candidate itself. We remove energy from additional simultaneous 5proton-proton collisions by subtracting an average energy density computed on a per-event basis [32, 58]. For electrons and muons, we divide ET, cone by the lepton pT to obtain the relative isolation Irel = ET, cone/pT, which is required to be less than 0.15. We require ET, cone < 2 GeV for τh candidates. The difference in the reconstruction efficiencies of muons and electrons between data and simu- lations is estimated with a standard technique that uses dilepton decays of Z bosons. Scale fac- tors (SFs) are applied to simulation to match the data efficiencies and are pT- and η-dependent. The combined muon identification and isolation efficiency uncertainty is 11% at muon pT of 10 GeV and 0.2% at 100 GeV. The corresponding uncertainties for electrons are 14% and 0.6%. We reconstruct jets from particle flow (PF) candidates using the anti-kT algorithm [61] with a distance parameter of 0.5. Jets are required to have |η| < 2.5 and pT > 30 GeV and have ∆R > 0.3 with respect to any isolated electron, muon, or τh candidate. The jet energy scale (resolution) is corrected using pT- and η-dependent data-to-simulation scale (resolution) ratios [62]. Jet four-momenta are varied using the uncertainty on these correction factors to account for the uncertainty in the jet energy scale measurement. We account for any additional discrepancy in EmissT between simulation and data [32] arising from PF candidates that are not clustered into jets and find that this discrepancy results in a negligible systematic uncertainty. To determine if a jet originated from a bottom quark, we use the combined secondary-vertex (CSV) algorithm, which calculates a likelihood discriminant from the track impact parame- ter and secondary-vertex information [63]. A loose, medium, and tight discrimination selects b jets with average efficiencies of 85%, 70%, and 50%, c jets with average misidentification prob- abilities of 40%, 20%, and 5%, and light-parton jets (u, d, s, g) with average misidentification probabilities of 10%, 1.5%, and 0.1%, respectively. Scale factors, depending on pT and |η|, are measured in data control samples of tt and µ+jets events and are used to correct the tagging efficiencies obtained from simulation. A weight is applied to the response of the b-tagging al- gorithm for each jet that is matched to a bottom quark. A similar procedure is applied to model the mistag probability for jets originating from light partons and c quarks. The b-, c-, and light-parton-tagging efficiencies are varied separately within their statistical uncertainties, and data-to-simulation SFs are applied and varied within the measured uncertainties [63, 64]. The b and c quark SFs are treated as correlated, and the light-parton SFs are treated as uncorrelated with the heavy-flavor SF. 6 Fully hadronic final state Many signatures for physics beyond the standard model (BSM) result in long decay chains that produce high-multiplicity final states. Most searches for SUSY involve either leptonic final states or missing transverse momentum, but fully hadronic final states that do not result in missing transverse momentum have been explored less thoroughly. This section presents a search in a high-multiplicity, fully hadronic final state with no missing transverse momentum requirement. The multiplicity of b-tagged jets is used as a discriminating variable. Results are interpreted in a model in which pair-produced gluinos each decay via g˜ → tbs, which is allowed when λ′′332 6= 0, so that a top antisquark couples directly to b and s quarks. The top squark is assumed to be much heavier than the gluino, resulting in the three-body decay of the gluino shown in Fig. 1. All supersymmetric particles other than the top squark and the gluino are assumed to be decoupled. The top-squark mass and λ′′332 are assumed to take values such that the gluino decays promptly. Because the coupling λ′′332 involves heavy quarks, it is relatively unconstrained by measurements of nucleon stability or neutrino masses [4]. 6 6 Fully hadronic final state P1 P2 g˜ g˜ t b s s b t Figure 1: Diagram for pair production of gluinos that decay to tbs. 6.1 Event selection Events are selected by the trigger via a requirement on the sum of the pT of the jets in the event that varied between 650 and 750 GeV over the course of data taking. Substantial background suppression is achieved through the application of multiplicity re- quirements on the jets reconstructed in the event, together with pT threshold requirements. We require at least four jets with pT > 50 GeV, where at least one jet must additionally satisfy pT > 150 GeV. All jets with pT > 50 GeV are used to calculate the offline HT of the event, which is required to be greater than 1.0 TeV. With this selection, the trigger efficiency, measured with prescaled triggers that have lower HT requirements, is consistent with 100%; separate studies with leptonic triggers have shown that there is no source of inefficiency common to all triggers with HT requirements. The tight CSV requirement is used to identify jets arising from b quarks. At least two such b-tagged jets are required. Events are required to have no isolated muons or electrons with pT > 10 GeV. This requirement renders backgrounds due to feed-down from leptonic final states essentially negligible and ensures that this analysis is disjoint from the leptonic variant described in Section 7. The data are divided into bins of jet multiplicity Njet and the scalar sum of the transverse mo- menta of the jets, HT. The HT bins are 1.0 < HT < 1.75 TeV and HT > 1.75 TeV and the Njet bins are 4, 5, 6, 7, and ≥8. In each of these ten bins we fit the multiplicity of b-tagged jets, Nb. 6.2 Standard model background The dominant background in this analysis is from QCD multijet events (hereafter labeled as QCD), with contributions from tt becoming important at large values of Nb. Background sources other than multijet events are estimated directly from simulation. These backgrounds include tt production, hadronic decays of W and Z bosons, single top quark production, ZZ production, rare processes that include a tt pair (ttW, ttZ, ttH, and tttt), and leptonically decay- ing W bosons in which the lepton is not reconstructed correctly or the lepton is a hadronically decaying τ lepton. As the dominant background in this analysis arises from QCD multijet events, the modeling of this component is crucial. We proceed by deriving corrections from data to the simulated QCD background to predict the distribution of the number of b-tagged jets. There are three main concerns: the modeling of the Njet and HT distributions, the flavor composition of the QCD events, and the b quark production mechanisms. The theoretical uncertainty on the Nb distribution arising from mismodeling of the Njet and HT distributions is avoided by binning the sample in these two variables. The modeling of the flavor composition is corrected to match the data (Section 6.2.1). The modeling of the b quark production mechanism is also validated with the data (Section 6.2.2). With this procedure, we obtain an estimate of the QCD 6.2 Standard model background 7 CSV 0.9 0.92 0.94 0.96 0.98 1 En tri es / 0. 01 210 310 Data Total fit b jets c jets Light jets Non-QCD CMS = 8 TeVs -1L = 19.5 fb Figure 2: The CSV distribution in data for 4 ≤ Njet ≤ 5, HT > 1.1 TeV, and CSV > 0.9. The solid line is the result of a fit to the data with MC templates. Error bars reflect statistical uncertainties in the data (smaller than the marker) and MC samples. background from data in the variables Njet, HT, and Nb. Measurements of tt + jets events, which are the second-largest background to this analysis, have demonstrated that the jet multiplicity distribution is modeled within the uncertainties [65]. 6.2.1 Flavor composition correction Although this analysis will determine the overall normalization of the QCD multijet back- ground from data, an uncertainty arises from the poorly known flavor composition of this background. To ensure that the simulated QCD events have the appropriate flavor composi- tion, events are reweighted to match the flavor composition measured in data. The coefficients used in the reweighting procedure are derived from a fit to the distribution of the CSV dis- criminant. This fit is performed in a control region comprising events with only four or five reconstructed jets, to exclude a potential bias from signal contamination, and with the slightly tightened CSV discriminator requirement CSV > 0.9 (compared with the nominal 0.898). Ad- ditionally, to avoid bias due to the large weights arising from the low equivalent luminosity of the simulated QCD samples with HT < 1.0 TeV, the HT requirement is increased slightly to HT > 1.1 TeV; it has been verified that the flavor composition corrections are statistically compatible for requirements of 1.0 or 1.1 TeV. The fit of the distribution of the CSV discriminant is performed, including the statistical uncer- tainty in the MC prediction as nuisance parameters in the fit via the Barlow–Beeston method [66]. The overall QCD contribution is normalized to the data yield, less the expected non-QCD yield (obtained from simulation). Reconstructed jets are matched to the corresponding simulated jets, and templates for the CSV discriminant are formed for each flavor. The relative normal- ization of templates corresponding to jets matched at the generator level to bottom and charm partons are allowed to vary in the fit. The small contributions of non-QCD events (mainly tt) and light-parton jets are fixed in the fit, with the uncertainty in the light fraction considered as a systematic uncertainty. The fractions of bottom, charm, and light-parton jets prior to the fit are fb, fc, and flight, respec- 8 6 Fully hadronic final state tively. The fit provides new fractions f ′i defined as f ′i = ni nb + nc + nlight + nnon-QCD (2) where nb and nc are the fitted yields of bottom and charm jets, respectively; nlight and nnon-QCD are the fixed yields of light-parton and non-QCD jets, respectively. The index i corresponds to b, c, light-parton, and non-QCD events. The values of wb jeti = f ′ i / fi are listed in Table 1 (before reweighting) and the fit of the CSV distribution is shown in Fig. 2. The fit quality is good, with χ2/d.o.f. = 7.0/8 (where d.o.f. is the number of degrees of freedom in the fit), providing confidence in the modeling of the CSV distribution. For each simulated QCD multijet event, a weight is assigned based on the flavor fractions: wevent =∏ b jet wb jeti , (3) where wb jeti is a per-jet weight that is assigned to each b-tagged jet and depends on the flavor i of the parton matched to the reconstructed jet. This form of the per-event reweighting is mo- tivated by treating the corrections as independent corrections to the per-jet efficiency; an alter- nate reweighting procedure that reweights only bb pairs gives similar results. The reweighting procedure has a small effect on the Nb distributions, with no Nb bin changing by more than 1.2 standard deviations due to the reweighting. Although the fit models the data well, the good agreement between the model and the data could occur if mismodeled distributions accidentally have a linear combination that is con- sistent with the data. To eliminate this possibility, fits are performed with variations of the fit range. Even with an extreme variation in which the most sensitive region of the fit (CSV > 0.98) is removed, the fit results are still consistent with the nominal fit. There is no evidence for any systematic effect. As an additional cross-check, the fit is iterated: after the simulated events have been reweighted, the reweighted templates are fit to the data again. The resulting heavy-flavor weights are con- sistent with unity, as shown in Table 1. Table 1: The b jet weights wb jeti derived from the fit of the Njet = 4–5 control region before and after reweighting the QCD MC sample. The last column shows the result of the validation by iteration of the fit. Flavor Before reweighting After reweighting b 0.94± 0.03 1.02± 0.03 c 2.00± 0.43 0.84± 0.18 Light Fixed to 1.0 Fixed to 1.0 It is important to demonstrate that the weights derived in this fit of the Njet = 4–5 control region are applicable to the Njet ≥ 6 signal regions. This has been verified in two ways. First, the weights have been applied to the Njet = 6 region, which has a negligible signal contribution. The corrected predictions show good agreement with the data, as seen in Fig. 3. Second, as the expected signal yield is extremely small compared to the background in the region defined by HT > 1.1 TeV, Nb ≥ 2, and Njet ≥ 6, the CSV distribution can be fit directly in this region. The reweighting parameters resulting from this fit are all within one standard deviation of those from the low-Njet control region. The fit of the Njet ≥ 6 region is not used in the reweighting procedure because of the larger statistical uncertainty, as well as because of the potential bias arising from signal contamination. 6.2 Standard model background 9 bN 1.5 4.5 Ev en ts 1 10 210 310 410 CMS = 8 TeVs -1L = 19.5 fb QCD +jetstt Other Gluino (M=1 TeV) Data = 4jet N < 1.75 TeVT1.00 < H bN 2 3 4 re si du al s N or m al iz ed 3− 2− 1− 0 1 2 3 bN 1.5 5.5 Ev en ts 1 10 210 310 410 CMS = 8 TeVs -1L = 19.5 fb QCD +jetstt Other Gluino (M=1 TeV) Data = 5jet N < 1.75 TeVT1.00 < H bN 2 3 4 5 re si du al s N or m al iz ed 3− 2− 1− 0 1 2 3 bN 1.5 4.5 Ev en ts 1 10 210 310 410 CMS = 8 TeVs -1L = 19.5 fb QCD +jetstt Other Gluino (M=1 TeV) Data = 6jet N < 1.75 TeVT1.00 < H bN 2 3 4≥ re si du al s N or m al iz ed 3− 2− 1− 0 1 2 3 bN 1.5 3.5 Ev en ts 1 10 210 CMS = 8 TeVs -1L = 19.5 fb QCD +jetstt Other Gluino (M=1 TeV) Data = 4jet N > 1.75 TeVT H bN 2 3≥ re si du al s N or m al iz ed 3− 2− 1− 0 1 2 3 bN 1.5 4.5 Ev en ts 1 10 210 CMS = 8 TeVs -1L = 19.5 fb QCD +jetstt Other Gluino (M=1 TeV) Data = 5jet N > 1.75 TeVT H bN 2 3 4≥ re si du al s N or m al iz ed 3− 2− 1− 0 1 2 3 bN 1.5 3.5 Ev en ts 1 10 210 CMS = 8 TeVs -1L = 19.5 fb QCD +jetstt Other Gluino (M=1 TeV) Data = 6jet N > 1.75 TeVT H bN 2 3≥ re si du al s N or m al iz ed 3− 2− 1− 0 1 2 3 Figure 3: Distribution of Nb for data (dots with error bars) and corrected predictions. The upper (lower) row shows data in which events are required to have 1.00 < HT < 1.75 TeV (HT > 1.75 TeV). The jet multiplicity requirements are Njet = 4 (left), Njet = 5 (middle), and Njet = 6 (right). The hatched bands shows the statistical uncertainty of the simulated data. The bottom panels of each plot show the difference between the data and corrected prediction divided by the sum in quadrature of the statistical uncertainties associated with each. 6.2.2 Gluon splitting systematic uncertainty In QCD multijet events, jets containing b quarks are produced in three different ways: pair production (qq → bb), flavor excitation (bq → bq and charge conjugate), and gluon splitting (g → bb). The first two processes are important primarily in the initial hard scatter, with the second suppressed owing to the small intrinsic b-quark content of the proton. Pair production is known to be well modeled by the MADGRAPH generator, but the rate of gluon splitting is known to be off by up to a factor of two in the region of phase space dominated by the parton shower [67], necessitating an additional systematic uncertainty derived from data. The systematic uncertainty is obtained by constructing an alternative template to be used in the signal extraction fit described in Section 6.4. The alternative template for the QCD component differs from the nominal template by an additional g→ bb component, which may have a neg- ative normalization if the simulation overpredicts gluon splitting. The Nb distribution shape of this component is derived from g → bb simulated events. Its normalization is obtained by 10 6 Fully hadronic final state bbR∆ 0 1 2 3 4 5 6 En tri es / 0. 4 0 5 10 15 20 25 30 35 CMS = 8 TeVs -1L = 19.5 fb Data MC R < 1.6∆to data-MC in Gluon splitting scaled = 6jet N > 1.75 TeVT H Figure 4: Distribution of the angular distance between any two b-tagged jets for data (dots with error bars), uncorrected MC prediction (band), and generator-matched gluon splitting events scaled to the difference between data and simulation in ∆Rbb < 1.6 (histogram) for events with HT > 1.75 TeV and Njet = 6. The error bars and bands include the statistical uncertainty in the data and MC simulation, respectively. There are no entries with ∆Rbb < 0.4 as a consequence of the jet definition. comparing the ∆Rbb distributions for data and simulation, where ∆Rbb is the angular distance computed between any two b-tagged jets in the event. The simulated distribution, shown in Fig. 4, is normalized to data in the high-∆Rbb region, ∆Rbb > 2.4. The difference between data and simulation in ∆Rbb < 1.6 is assumed to arise entirely from g → bb events, and this differ- ence provides the normalization for the (negative) correction that we add to the QCD template to obtain an estimate of the systematic uncertainty of the g→ bb component. The assumption that the difference arises entirely from gluon splitting to bb provides a conservative estimate of the systematic uncertainty because, for example, events with gluon splitting to cc would gen- erate a smaller difference in the distribution of the multiplicity of b-tagged jets because of the smaller b-tagging efficiency for charm jets. The difference between data and simulation is de- termined separately for each (HT, Njet) bin. The constraint on the modeling of gluon splitting at low ∆Rbb is used as an uncertainty rather than a correction because in the large-Njet regions statistical fluctuations at low ∆Rbb are larger than the size of the correction. 6.3 Systematic uncertainties The systematic uncertainties due to the modeling of the QCD background constitute an impor- tant part of the total uncertainty. The uncertainty on the QCD flavor composition and gluon splitting are evaluated as described in Section 6.2.1 and Section 6.2.2, respectively. As tt is a subdominant background, the effect of its uncertainties are generally small. The tune of the underlying event, as well as variations of the renormalization, factorization, and match- ing scales are considered. Also, the inclusive tt production cross section is varied according to its NNLO+(next-to-NLL) uncertainty [68]. The top quark pT spectrum, which is not modeled well by simulation [69], is reweighted to agree with the data. The background contribution arising from ttbb production is doubled to match the data [65], with a 100% uncertainty. The cross sections of the remaining backgrounds are varied by 50%. The uncertainties in the jet pT scale, the jet resolution, and the b-tagging SFs for heavy-flavor and light-parton jets are 6.4 Control sample fit 11 evaluated as discussed in Section 5. The QCD MC simulation is affected by large statistical uncertainties, which are taken into ac- count by variations in which single bins of each Nb histogram are varied according to their statistical uncertainty. This is the largest systematic uncertainty in most Nb bins, and is about 40% in the most sensitive signal region (HT > 1.75 TeV and Njet ≥ 8). The systematic uncertain- ties are individually calculated for each (HT, Njet, Nb) bin. In general, the statistical uncertainty from data, the statistical uncertainty of the QCD MC simulation, and the sum of all other sys- tematic uncertainties are of similar magnitude in each bin, ranging from 1% to more than 50% across the bins. No uncertainty is assigned for trigger efficiency, which is consistent with 100% with per mille uncertainties. The signal samples are generated with a fast simulation. The efficiency for tagging b jets and the mistag rate for charm and light-flavor quarks is corrected to match the efficiency predicted by full simulation. Nuisance parameters parameterizing the uncertainty in these corrections for bottom jets, charm jets, and light-parton jets are considered separately and assumed to be mutually uncorrelated. Most signal systematic uncertainties are modeled as modification of the templates of the Nb distribution. The only exceptions are that of the luminosity uncertainty and the PDF uncer- tainty, which are modeled assuming a log-normal distribution for the corresponding nuisance parameter for each (Njet, HT) bin. 6.4 Control sample fit Signal-depleted control regions at low Njet (Njet = 4, 5, and 6) are studied before examining the signal region. For low jet multiplicities, tt backgrounds are less important, giving a largely pure sample of QCD events. A binned maximum likelihood fit of the Nb distributions is performed in which systematic uncertainties are profiled. Systematic uncertainties are included as shape uncertainties by in- terpolating between Nb histograms corresponding to ±1 standard deviation variations. As the HT and Njet dependence of the QCD contribution may not be modeled well, a separate normal- ization of the QCD contributions is allowed in each bin of HT and Njet. The likelihood used in the fit of the yields Nijk in the Nb distributions of the signal and control regions is L = ∏ i∈HT bins j∈Njet bins k∈Nb bins n∈syst P(Nijk|θn)Poisson ( Nijk|µsignalνijk,signal + µij,QCDνijk,QCD + νijk,other ) . (4) Here µsignal and µij,QCD are normalization constants. The parameters µsignal and µij,QCD do not have a dependence on k because the Nb input distribution is fixed for a given HT and Njet bin. The yields of signal, QCD background, and non-QCD background are relative to the nomi- nal values specified by νijk,signal, νijk,QCD, and νijk,other, respectively. In other words, there is a floating QCD normalization in each (HT, Njet) bin and fixed non-QCD background yields. The systematic uncertainties are included with nuisance parameters θn that can affect the inter- polation between the ±1 standard deviation templates; the parameters νijk,signal, νijk,QCD, and νijk,other are dependent on these nuisance parameters. These parameters are the same for all (HT, Njet) bins except for those associated with MC statistics, which have separate parameters for each (HT, Njet, Nb) bin. 12 7 Single-lepton final state In the control sample fit, the product over Njet bins is restricted to Njet = 4, 5, 6, and the signal yields are fixed to zero. The data and MC simulation inputs to this fit are shown in Fig. 3. All of the nuisance parameters are consistent within one standard deviation of their prefit un- certainties, except for a 1.4 standard deviation difference in the light jet fraction nuisance pa- rameter, which is however a subdominant uncertainty in the high-Njet signal region. 6.5 Results for fully hadronic final state The likelihood used in the fit of the signal region is that given by Eq. (4), with Njet = 6, 7, ≥8, and µsignal left free. Figure 5 shows a comparison of the data with the corrected simulation, where the QCD component has been scaled to the data yield minus the non-QCD background yields obtained from simulation. The yields corresponding to the Njet ≥ 8 and HT > 1.75 TeV region in Fig. 5 are shown in Table 2. Table 2: Summary of prefit expected background, expected signal for mg˜ = 1 TeV, and observed yields for Njet ≥ 8 and HT > 1.75 TeV. Uncertainties are statistical only. Background Nb = 2 Nb ≥ 3 QCD multijet 24.7 ± 3.8 1.4 ± 0.9 tt+jets 5.5 ± 0.6 0.7 ± 0.2 Other 0.6 ± 0.4 < 0.1 Total background 30.9 ± 3.9 2.2 ± 0.9 Data 31 2 Signal (mg˜ = 1.0 TeV) 22.8 ± 0.3 8.9 ± 0.2 At each gluino mass, the best fit returns zero signal events. The efficiency after applying all the selection criteria is shown in Fig. 6 and reaches a plateau of around 20% for mg˜ > 0.7 TeV. Figure 7 shows the expected and observed limits compared to the gluino pair production cross section. To summarize the fully hadronic search, the data in the signal regions are well described by the background predictions. The results are interpreted in terms of a specific model of RPV SUSY in which gluinos are pair produced and each gluino decays promptly via g˜→ tbs. Cross section limits are calculated and result in a 95%CL lower limit on the gluino mass of 0.98 TeV within this simplified model. 7 Single-lepton final state This section describes a search for g˜ → tbs events in which the top quark undergoes semilep- tonic decay. The strategy is similar to that of the previous section but with a requirement of at least one isolated charged lepton, which rejects most of the QCD background. We select final states with one lepton and multiple jets. Then we use the number of b-tagged jets as a discriminating variable to separate the signal from the background. 7.1 Event selection The analysis considers events selected by a trigger requiring an electron with pT > 27 GeV and |η| < 2.5 or a trigger requiring a muon with pT > 24 GeV and |η| < 2.1. Both triggers include loose isolation requirements. The offline selection raises the pT threshold to pT > 35 GeV for both muons and electrons while the |η| selection remains the same. 7.1 Event selection 13 bN 1.5 4.5 Ev en ts 1 10 210 310 410 CMS = 8 TeVs -1L = 19.5 fb QCD +jetstt Other Gluino (M=1 TeV) Data = 7jet N < 1.75 TeVT1.00 < H bN 2 3 4≥ re si du al s N or m al iz ed 3− 2− 1− 0 1 2 3 bN 1.5 4.5 Ev en ts 1 10 210 310 CMS = 8 TeVs -1L = 19.5 fb QCD +jetstt Other Gluino (M=1 TeV) Data 8≥ jet N < 1.75 TeVT1.00 < H bN 2 3 4≥ re si du al s N or m al iz ed 3− 2− 1− 0 1 2 3 bN 1.5 3.5 Ev en ts 1 10 210 CMS = 8 TeVs -1L = 19.5 fb QCD +jetstt Other Gluino (M=1 TeV) Data = 7jet N > 1.75 TeVT H bN 2 3≥ re si du al s N or m al iz ed 3− 2− 1− 0 1 2 3 bN 1.5 3.5 Ev en ts 1 10 210 CMS = 8 TeVs -1L = 19.5 fb QCD +jetstt Other Gluino (M=1 TeV) Data 8≥ jet N > 1.75 TeVT H bN 2 3≥ re si du al s N or m al iz ed 3− 2− 1− 0 1 2 3 Figure 5: Data (dots with error bars) and the corrected prediction of the Nb distribution in the high-Njet signal region. The hatched band shows the MC statistical uncertainty. The upper (lower) row shows events with 1.00 < HT < 1.75 TeV (HT > 1.75 TeV). The jet multiplicity requirements are Njet = 7 (left) and Njet ≥ 8 (right). The bottom panels of each plot show the difference between the data and corrected prediction divided by the sum in quadrature of the statistical uncertainties associated with each. 14 7 Single-lepton final state (GeV)g~m 400 600 800 1000 1200 Ef fic ie nc y 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 6≥ jetN = 6jetN = 7jetN 8≥ jetN CMS = 8 TeVs -1L = 19.5 fb Figure 6: Signal efficiencies as a function of mg˜ for Nb ≥ 2, HT > 1 TeV, and Njet ≥ 6, together with the breakdown among Njet bins. (GeV)g~m 600 800 1000 1200 (pb ) σ 3−10 2−10 1−10 1 10 σ 1±95% CL expected limit 95% CL observed limit cross sectiong~g~ g~ >> mq~m tbs→g~, g~g~→pp 0 leptons CMS = 8 TeVs -1L = 19.5 fb Figure 7: The 95% CL limit on the gluino pair production cross section as a function of mg˜ in the analysis of all-hadronic final states. The signal considered is pp → g˜g˜, followed by the decay g˜→ tbs. The red band shows the theoretical cross section and its uncertainty. The blue dashed lines show the uncertainty on the expected limit. We measure the trigger efficiency with respect to the offline selection in bins of η and pT of the lepton, and find it to vary from approximately 92% for |η| < 0.8 to 65% for 1.2 < |η| ≤ 2.4 in the case of muons. Electron efficiencies are within a few percent of those quoted for muons. For both lepton flavors, the variation of the efficiency over pT in a given η range is 1–2% for pT > 35 GeV. The baseline selection requires at least one lepton, and six jets with pT > 30 GeV. The medium working point of the CSV b tagging discriminator is used and at least one jet must pass this selection. Events with two identified leptons are allowed in the sample, but to avoid double counting we veto events in the electron sample if they also contain an identified muon. The electron and muon samples are distinguished to allow cross-checks; however most sys- tematic effects are correlated between the two samples and this is considered when fitting the 7.2 Standard model background 15 b-tagged jet multiplicity distribution to extract the signal. 7.2 Standard model background The selected events are divided into three signal regions, according to their jet multiplicity: 6, 7, and ≥8 jets. For each signal region, the signal and background yields are obtained by com- paring the observed multiplicity distribution of b-tagged jets with their respective background shapes. The main source of background to this search is the production of pairs of top quarks in as- sociation with jets. Additional contributions from single top quark, vector bosons, and QCD multijet production are relevant for low b jet multiplicities, becoming negligible for events with at least three b jets. Events with at least three b jets also have a small contamination, typically below 1%, from ttV events (where V = W or Z). The background from SM tttt production and ttH production is negligible, due to the small cross section for these processes. We study the b jet multiplicity of the background sources using MC simulation. Events are cor- rected for the different response of the b tagging algorithm in simulation and data as described in Section 5. We verified that the b tagging efficiency and the mistag rate vary negligibly in semileptonic tt events going from four to nine jets, always remaining within their uncertainty. Furthermore, we compared data and simulation in control regions with one or two leptons and four or five jets and found good agreement within the uncertainties in the b-tagging SFs. We account for residual small discrepancies by allowing for an MC mismodeling of SM events with four b quarks, as discussed in more detail in Section 7.3.1. The corrected b jet multiplicity provides the prediction for the SM background, which is com- pared with data in order to check for the presence of a signal. The final signal extraction fit obtains the background normalization from data using an extended likelihood. Therefore only the shape of the multiplicity distribution of b tagged jets is taken from simulation. This consid- erably reduces the systematic uncertainty in the background, since this shape is very weakly dependent on the jet energy scale and on the choices of matching and renormalization scales. 7.3 Systematic uncertainties 7.3.1 Background The background shape is affected by the jet energy scale uncertainty; the uncertainty in the b tagging efficiency SFs; the variation of renormalization, factorization, and matching scales; and the MC statistical uncertainty. Furthermore, we include a systematic effect parameterizing the mismodeling of the fraction of events with four bottom quarks. We evaluate the effect of jet energy, matching, renormalization, and factorization scales by repeating the selection procedure on MC simulated samples with the scales shifted up or down. The jet energy scale is varied as described in Section 5. The matching, renormalization, and factorization scales are fixed to factors of 2.0 and 0.5 with respect to the nominal scale, for positive and negative variations, respectively. The renormalization and factorization scales are varied simultaneously. The uncertainty from the b tagging SF is computed by comparing the b-tagged jet multiplicity distributions obtained by correcting the tagging efficiencies with SFs shifted by ±1 standard deviation. The uncertainties in the b jet and c jet SFs are taken to be correlated. The parameterization of the mismodeling of events with four b quarks is not as straightfor- ward as the computation of the other uncertainties. While the data-corrected MC distribution 16 7 Single-lepton final state is expected to account for events with multiple b-tagged jets originating from mistags, the con- tribution from gluon splitting to bb, and of SM four b quark events in general, is sensitive to the details of the MC modeling [65], as discussed in Section 6.2.2. We constrain the uncertainty in this contribution by studying the agreement between data and simulation in events with one identified electron, one identified muon, and associated jets. We consider separately events with four or five jets. Furthermore, we use single-lepton control regions by selecting events with one electron or one muon and four or five jets. These control regions provide a high- purity sample of tt+jets events, for which the signal contamination is expected to be negligible. Figure 8 shows that the largest difference between the prediction used in the analysis and the observed yield in the dileptonic control sample is an excess of less than one standard deviation in the three and four b jet bins for the four-jet sample. Total uncertainties are shown in this figure, including uncertainties that affect only the normalization of the background prediction and not the shape. The single-lepton control regions have similarly small discrepancies for Nb ≥ 3. These are the bins of the b tag multiplicity distribution that are most sensitive to the signal so we parameterize this effect and include it in the analysis as a systematic uncertainty. Ev en ts 0 100 200 300 400 500 600 Data +jetstt W+jets bWb Wtt Ztt Drell-Yan Single top QCD Sys uncertainty CMS = 8 TeVs -1L = 19.3 fb = 4jet channel, Nµe bN 1 2 3 4 D at a/ M C 0 0.5 1 1.5 2 2.5 Ev en ts 0 20 40 60 80 100 120 140 Data +jetstt W+jets bWb Wtt Ztt Drell-Yan Single top QCD Sys uncertainty CMS = 8 TeVs -1L = 19.3 fb = 5jet channel, Nµe bN 1 2 3 4 5 D at a/ M C 0 0.5 1 1.5 2 2.5 Figure 8: Distribution of the number of b-tagged jets for events with one electron, one muon, and Njets = 4 (left), or Njets = 5 (right) in data, compared to the background prediction from simulation corrected for the b tagging response. The hatched region represents the total uncer- tainty on the background yield. Using the data yields in the b tag multiplicity bins of the control samples defined above, we construct a system of three equations and three unknowns for each control region. In a sample with N events and J jets, we write the number of events with n b-tagged jets N(n, J) as a function of the b-tagging efficiency eb and the mistag rate emis. Assuming that the sample consists mainly of bb events, i.e. that b-tag multiplicities > 2 originate from mistagged jets, 7.3 Systematic uncertainties 17 one finds: N(n, J) =N 2 ∑ i=0 θ(n− i)θ(J − n− 2+ i) ( 2 i )( J − 2 n− i ) × (1− eb)2−ieib(1− emis)J−2−n+ien−imis , (5) where θ(m) is the Heaviside step function, where θ(m) = 1 for m ≥ 0 and θ(m) = 0 for m < 0 and the standard representation for the binomial coefficient is used. The index i runs over the number of true b quarks that are within the acceptance and tagged. Once a subset ∆N of the events is allowed to contain four real b quarks, Eq. 5 is modified as follows: N′(n, J) =(N − ∆N)N(n, J) N + ∆N 4 ∑ i=0 θ(n− i)θ(J − n− 4+ i) ( 4 i )( J − 4 n− i ) × (1− eb)4−ieib(1− emis)J−4−n+ien−imis . (6) Taking as input the yield observed for three values of n at a given J, one can solve a system of three equations and three unknowns and derive values for eb, emis, and ∆N. Rather than solving for ∆N, we introduce ∆ f4b = ∆N N ∣∣∣∣ data − ∆N N ∣∣∣∣ MC (7) and solve for each control region the set of equations e2mise 2 b(1− ∆ f4b) + e4b∆ f4b[ e2b(1− emis)2 + e2mis(1− eb)2 + 4emis(1− eb)eb(1− emis) ] (1− ∆ f4b) + 6e2b(1− eb)2∆ f4b = N4 N2[ 2e2bemis(1− emis) + 2eb(1− eb)e2mis ] (1− ∆ f4b) + 4e3b(1− eb)∆ f4b[ e2b(1− emis)2 + e2mis(1− eb)2 + 4emis(1− eb)eb(1− emis) ] (1− ∆ f4b) + 6e2b(1− eb)2∆ f4b = N3 N2[ e2b(1− emis)2 + e2mis(1− eb)2 + 4emis(1− eb)eb(1− emis) ] (1− ∆ f4b) + 6e2b(1− eb)2∆ f4b [2(1− eb)eb(1− emis)2 + 2(1− eb)2emis(1− emis)] (1− ∆ f4b) + 4eb(1− eb)3∆ f4b = N2 N1 . (8) Here Ni indicates the yield in the b tag multiplicity bin with Nb = i, and the unknowns are eb, emis, and ∆ f4b. We solve the system of equations numerically, neglecting terms quadratic in the mistag rate. Since ∆ f4b is common to all control regions, we use the resulting values of the average tagging efficiency and the average mistag rate determined in each control region to construct a global χ2: χ2 (∆ f4b) = ∑ i=1,...,Nb j∈CR ( Nijobs − NijMC − N4b(ejb, ejmis,∆ f4b) )2 σ2ij , (9) where the sum over j spans the different control regions, the index i gives the bin of the mul- tiplicity of b-tagged jets, and σij is the sum in quadrature of the statistical uncertainty in data and total uncertainty in simulation. Minimizing the χ2 results in an improved determination of ∆ f4b from the data in all control regions. 18 7 Single-lepton final state We associate a systematic uncertainty with the shape of the background, by determining ∆ f4b with the information from both the dilepton and single-lepton control regions, obtaining∆ f4b = −0.011± 0.049. The choice of combining the two control regions is justified by the fact that fit- ting them separately we obtain compatible results. We use ±1 standard deviation variations of ∆ f4b to construct two new background shapes in the 6 jets, 7 jets, and ≥8 jets signal regions. The difference between these shapes is used to evaluate a systematic uncertainty that is used in the signal extraction fit. The determination of the systematic uncertainty depends on the val- ues of the efficiencies used in the fit, on their uncertainty and on the choice of control regions. For this reason we compute the limit on the signal cross section for several different choices of control regions, tagging efficiencies, and mistag rates. The observed variations are below 10−3 of the cross section value obtained with the value of ∆ f4b from the combined fit. In the bins with fewer than three b-tagged jets the dominant sources of uncertainty come from the jet energy, renormalization, and factorization scales. For higher multiplicities of b-tagged jets, the uncertainty in the tagging SFs and the mismodeling of events with four b quarks be- come the main sources of uncertainty. 7.3.2 Signal The uncertainties in the signal efficiency and signal shape include the jet energy scale uncer- tainty, the b tagging SF uncertainty, the uncertainty in the PDFs, the uncertainty in the mea- sured integrated luminosity, the uncertainty in trigger and identification efficiencies for lep- tons, and the uncertainty in the MC modeling of ISR and FSR. The jet energy scale and b tagging SF uncertainties are computed in the same way as for the background. However, since the signal samples are processed through a fast rather than full detector simulation, we apply additional compensating SFs for the b-tagging efficiencies. The uncertainties in the reconstruction and trigger efficiencies of muons and electrons are estimated from Z→ `` events in bins of η and pT of the lepton; both are found to be always below 1%. The nominal efficiency from simulation is also corrected to reflect the lepton efficiency measured in data. The total efficiency for Njet ≥ 6 is 15.8% for mg˜ = 1.0 TeV. 7.4 Results for the single-lepton final state Figure 9 shows the b tag multiplicity for the 6 jets, 7 jets, and ≥8 jets signal regions, for events with at least one well-identified electron and for events with at least one well-identified muon. The hatched region represents the uncertainty propagated from the b-tagging correction fac- tors. We summarize the expected background, the expected signal for mg˜ = 1 TeV, and the observed yield in each signal region in Tables 3 and 4. The postfit uncertainties in the back- ground, shown in Tables 3 and 4, are considerably reduced with respect to the uncertainties in the prediction. The fit, described below, extracts the background normalization from data, reducing the total uncertainty to an uncertainty in the shape of the multiplicity distribution of b-tagged jets. Therefore uncertainties coming from jet energy, matching, renormalization, and factorization scales that affect mostly the total yield become almost negligible. The central values for the background prediction do not change significantly in the fit. The electron and muon samples, although presented separately to facilitate reinterpretations, are fitted simulta- neously; for the same reason, this fit does not assume a signal model and instead sets the signal yield to zero. Tables 3 and 4 also give the signal prediction for mg˜ = 1 TeV. Combining all the signal uncer- tainties gives a total uncertainty of ∼10–20% on the individual bins of b-tagged jet multiplicity. At high gluino masses (mg˜ ≥ 1 TeV) the uncertainty is dominated by the PDF uncertainty, 7.4 Results for the single-lepton final state 19 Ev en ts 0 500 1000 1500 2000 2500 Data +jetstt W+jets bWb Wtt Ztt Drell-Yan Single top QCD SF syst CMS = 8 TeVs -1L = 19.3 fb = 6jete channel, N bN 1 2 3 4 5≥ D at a/ M C 0 0.5 1 1.5 2 2.5 Ev en ts 0 100 200 300 400 500 Data +jetstt W+jets bWb Wtt Ztt Drell-Yan Single top QCD SF syst CMS = 8 TeVs -1L = 19.3 fb = 7jete channel, N bN 1 2 3 4 5≥ D at a/ M C 0 0.5 1 1.5 2 2.5 Ev en ts 0 10 20 30 40 50 60 70 80 90 100 Data +jetstt W+jets bWb Wtt Ztt Drell-Yan Single top QCD SF syst CMS = 8 TeVs -1L = 19.3 fb 8≥ jete channel, N bN 1 2 3 4 5≥ D at a/ M C 0 0.5 1 1.5 2 2.5 Ev en ts 0 500 1000 1500 2000 2500 3000 Data +jetstt W+jets bWb Wtt Ztt Drell-Yan Single top QCD SF syst CMS = 8 TeVs -1L = 19.3 fb = 6jet channel, Nµ bN 1 2 3 4 5≥ D at a/ M C 0 0.5 1 1.5 2 2.5 Ev en ts 0 100 200 300 400 500 600 Data +jetstt W+jets bWb Wtt Ztt Drell-Yan Single top QCD SF syst CMS = 8 TeVs -1L = 19.3 fb = 7jet channel, Nµ bN 1 2 3 4 5≥ D at a/ M C 0 0.5 1 1.5 2 2.5 Ev en ts 0 20 40 60 80 100 120 140 Data +jetstt W+jets bWb Wtt Ztt Drell-Yan Single top QCD SF syst CMS = 8 TeVs -1L = 19.3 fb 8≥ jet channel, Nµ bN 1 2 3 4 5≥ D at a/ M C 0 0.5 1 1.5 2 2.5 Figure 9: Distribution of the number of b-tagged jets for events with one electron (upper) or muon (lower) and 6 (left), 7 (middle) or ≥8 (right) jets, compared to the background prediction from simulation corrected for the b tagging response in data. The hatched region represents the uncertainty originating from the uncertainty in the b tagging correction factors. Most other uncertainties affect only the normalization and will cancel in the fit. while for lower masses the jet energy scale uncertainty constitutes the most important source of uncertainty, followed in magnitude by the uncertainty in the modeling of ISR and FSR. No sizable deviation from the expected SM yields is observed. We interpret the absence of an excess as an upper bound on the cross section for SUSY models predicting final states with one lepton and multiple b-tagged jets. The cross section limit is obtained with a maximum likelihood fit to the shape of the b-tagged jet multiplicity distribution and is converted to a bound on the mass of the gluino. In setting a limit on new physics, we consider the g˜ → tbs simplified model, which assumes λ′′332 6= 0. The main ingredient to the likelihood for a given multiplicity of b-tagged jets and lepton flavor in a given signal region (6 jets, 7 jets, ≥8 jets) is a Poisson function for n observed events, given an expected yield of eLσ+ B: P(n|eLσ+ B) = e −(eLσ+B) n! (eLσ+ B)n . (10) 20 7 Single-lepton final state Table 3: Summary of expected background, expected signal for mg˜ = 1 TeV, observed yields, and total background after the background-only fit for the electron samples considered in the analysis. The uncertainties given include all statistical and systematic uncertainties. e + 6 jets One b-tag Two b-tags Three b-tags Four b-tags Five b-tags Background prediction 2003±827 1701±762 281±130 27±17 8.0±6.8 Signal (mg˜ = 1 TeV) 1.9±0.3 2.9±0.5 1.9±0.3 0.41±0.10 0.03±0.01 Data 2128 1566 284 40 2 Background postfit 1967±54 1636±53 296.1±9.5 33.6±3.0 1.9±1.2 e+ 7 jets One b-tag Two b-tags Three b-tags Four b-tags Five b-tags Background prediction 373±200 352±199 67±39 8.7±6.3 1.1±1.1 Signal (mg˜ = 1 TeV) 2.0±0.3 3.4±0.5 2.7±0.4 0.86±0.15 0.07±0.02 Data 410 320 61 11 0 Background postfit 368±13 347±12 70.6±2.8 10.38±0.65 0.70±0.12 e+ 8 jets One b-tag Two b-tags Three b-tags Four b-tags Five b-tags Background prediction 73±51 70±49 18±15 2.7±2.1 0.47±0.38 Signal (mg˜ = 1 TeV) 2.4±0.4 4.9±0.8 4.7±0.7 2.0±0.3 0.23±0.04 Data 80 64 16 5 0 Background postfit 74.9±3.1 71.0±3.1 18.94±0.93 3.40±0.20 0.44±0.03 Table 4: Summary of the expected background, expected signal for mg˜ = 1 TeV, observed yields, and total background after the background-only fit for the muon samples. The un- certainties given include all statistical and systematic uncertainties. µ + 6 jets One b-tag Two b-tags Three b-tags Four b-tags Five b-tags Background prediction 2474±977 2002±801 322±152 30±29 7.7±6.5 Signal (mg˜ = 1 TeV) 3.0±0.5 4.6±0.7 2.8±0.4 0.6±0.1 0.04±0.03 Data 2585 1850 356 44 1 Background postfit 2425±60 1985±49 340±11 43.0±3.5 3.1±1.1 µ + 7 jets One b-tag Two b-tags Three b-tags Four b-tags Five b-tags Background prediction 493±203 448±180 88±39 10.7±7.2 0.9±2.8 Signal (mg˜ = 1 TeV) 3.0±0.5 5.0±0.7 3.9±0.6 1.1±0.2 0.09±0.04 Data 497 412 116 16 0 Background postfit 506±15 462±13 95.9±3.0 14.81±0.99 1.03±0.18 µ + 8 jets One b-tag Two b-tags Three b-tags Four b-tags Five b-tags Background prediction 112±47 104±46 26±12 4.3±2.1 0.39±0.75 Signal (mg˜ = 1 TeV) 3.7±0.6 7.0±1.0 6.4±0.9 2.5±0.4 0.33±0.06 Data 112 104 27 3 1 Background postfit 119.7±4.3 110.7±3.6 29.0±1.0 5.63±0.34 0.54±0.07 Here B is the expected background yield, e is the signal efficiency, L is the integrated luminosity of the data set, and σ is the cross section on which we want to set the limit. The extended likelihood is written as: L = e−(eLσ+B) n! (eLσ+ B)nPLN(e|e¯, δe)PLN(L|L¯, δL)PLN(B|B¯, δB). (11) We model the systematic uncertainty associated with the signal and the background prediction as log-normal functions PLN(x|x¯, δx) for the measured value x, given an expected value x¯ and an uncertainty δx. The full likelihood is obtained as the product of a set of likelihoods as the one in Eq. (11). The product runs over each multiplicity of b-tagged jets (one to five), lepton flavor (e, µ), and each of the three signal regions (6 jets, 7 jets,≥8 jets). The nuisance parameters 21 are taken to be fully correlated across the three signal regions with different jet multiplicities and the two lepton flavors: i.e. a common log-normal function for each nuisance multiplies the product of Poisson functions. For a given value of σ under test the likelihood is profiled with respect to the nuisance param- eters (L, e, and B). The result of this procedure is shown in Fig. 10 and results in a 95% CL lower limit on the gluino mass of 1.03 TeV when the gluino is assumed to decay exclusively to tbs. This is currently the strongest bound for this gluino decay mode. (GeV)g~m 400 600 800 1000 1200 (pb ) σ 3−10 2−10 1−10 1 10 95% CL expected limit 95% CL observed limit cross sectiong~g~ g~ >> mq~m tbs→g~, g~g~→pp 1 lepton CMS = 8 TeVs -1L = 19.3 fb Figure 10: The 95% CL limit on the gluino pair production cross section as a function of mg˜ in the one-lepton analysis. The signal considered is pp → g˜g˜, followed by the decay g˜ → tbs. The band shows the theoretical cross section and its uncertainty. The dashed lines show the uncertainty on the expected limit. 8 Dilepton final state We search for the RPV decays of the bottom squark (b˜) in an MSSM model featuring minimal flavor violation [8]. When the bottom squark is the LSP in this type of model, it can decay to a top quark and a down-type quark. We have chosen a model sensitive to the λ′′332 and λ′′331 hadronic RPV couplings, so the bottom squark decay of interest is to a top and either a strange or down quark. In contrast to Sections 6 and 7 which feature a gluino pair production model, this section focuses on a model with bottom squark pair production. We restrict this search to dilepton final states, where each top quark decays into a W boson, which in turn decays leptonically. A diagram of this process is shown in Fig. 11. To discriminate signal from background events, we analyze the reconstructed bottom squark mass and the transverse momenta of jets identified as coming from light quarks or gluons. The trigger requires two leptons, one of which has pT > 17 GeV and the other pT > 8 GeV. In the subsequent analysis at least two leptons passing identification and isolation criteria are required in each event. We require at least two selected jets to pass the loose b tagging selection, and in addition at least one of these must pass the medium b tagging selection. Additionally, we require that at least two jets fail the loose b tagging selection. This allows for unambiguous categorization of light-quark jets from the bottom squark decay and the b jets from the top decay. 22 8 Dilepton final state ¯˜ b b˜ t W+ t¯ W− p p s,d s¯, d¯ b b¯ `+ ν ν¯ `− Figure 11: Diagram for pair production of bottom squarks and their RPV decay. 8.1 Signal and background discrimination The dominant background to the signal originates from SM top quarks pair produced in asso- ciation with jets from ISR or FSR. Other SM processes account for a small (≈5%) contribution: single top quark production, diboson production, Drell–Yan production, and top quark pair production in association with vector bosons. Signal events contain a resonance that produces top quarks in association with a light-quark jet. The light-quark jet from this decay has a rela- tively high pT. We use both of these properties to discriminate between signal and background by the construction of a three-dimensional probability distribution over the reconstructed res- onance mass and the two pT values of the light-flavor jets. We associate the two highest pT non-b-tagged jets with light quarks from the bottom squark decays, the two highest pT b-tagged jets with bottom quarks from top quark decays, and the two highest pT leptons with leptons from W boson decays. In total, 6478 events pass the selec- tion requirements: 1723 in the ee channel, 1365 in the µµ channel, and 3390 in the eµ channel. Eleven events contain more than two leptons and are included in the analysis. 8.1.1 Light jet pT spectrum To model the light parton jet pT spectrum of SM processes, we assume that the light-parton jets are produced predominantly by ISR or FSR from tt events and therefore have a steeply falling pT spectrum. Signal events, on the other hand, are more likely to contain light-parton jets with relatively high pT. Letting p(1)T and p (2) T denote the transverse momenta of, respectively, the highest and second- highest pT light-parton jets in the event, we use simulated SM events to help choose the form of a two-dimensional probability density function ρSM2D (p (1) T , p (2) T ) with sufficient flexibility to fit the data well. This distribution is constructed as a sum of three 2D densities with signal region index i, ρ2Di (p (1) T , p (2) T ): ρSM2D (p (1) T , p (2) T ) = 2 3 ∑ i=1 fiρ2Di (p (1) T , p (2) T )θ(p (1) T − p(2)T ), (12) where ∑3i=1 fi = 1. The 2D distributions are multiplied by a Heaviside step function to enforce the ordering and a factor of 2 to normalize the function. For each component we find from simulation that the 2D density can be expressed as a product of 1D densities: ρ2Di (p (1) T , p (2) T ) = ρ jet i (p (1) T )ρ jet i (p (2) T ). (13) The 1D densities all have the same form: ρ jet i (pT) = λiα exp (−λipαT)pα−1T , (14) 8.1 Signal and background discrimination 23 where α is a parameter between 0 and 1 common to all components, while the λi are parameters differing in each component. This function has the steeply falling behavior that we expect from ISR or FSR, and potentially a longer tail than a pure exponential distribution, becoming identical to an exponential distribution for α = 1. The assumption that the density can be factorized is tested by comparing fit results from an SM MC simulation and found to work well. The parameters of these densities are determined by fitting to data, maximizing the likelihood function defined in Section 8.3. This function is sufficiently flexible to accommodate observed light-parton jet pT behavior while using a limited number of free parameters. Figure 12 shows the background-only hypothesis fits to the leading and subleading pT spectra for jets in data, illustrating good agreement. (GeV) T pLeading light jet 100 200 300 400 500 600 Ev en ts / ( 1 4.2 G eV ) 0 100 200 300 400 500 600 700 800 Data function Fitted /dof = 39.1/34 2χ 100 200 300 400 500 600 R at io 0 0.5 1 1.5 2 CMS L = 19.5 fb−1 √s = 8 TeV (GeV) T pSecond leading light jet 50 100 150 200 250 300 Ev en ts / ( 6 .75 G eV ) 0 200 400 600 800 1000 1200 1400 1600 1800 Data function Fitted /dof = 31.5/34 2χ 50 100 150 200 250 300 R at io 0 0.5 1 1.5 2 CMS L = 19.5 fb−1 √s = 8 TeV Figure 12: Background-only likelihood fits for the light-parton jet pT distributions, with signal cross section set to zero, for the (left) leading and the (right) second-leading light-quark jet. The line represents the fitted function and the points represent the data. The ratio of the data to the fitted function is also shown. To parameterize the signal distribution, we model jets as being sampled from the sum of two two-dimensional log-normal distributions and ordered by pT. The parameters of these two distributions and their relative fractions are determined by fitting to the signal simulation. We call this distribution ρsignal2D (p (1) T , p (2) T ). Figure 13 illustrates the difference between light-parton jet shapes. It shows the two-dimensional distribution of the background-only hypothesis fitted to the measured pT distribution and the signal distribution for an example mass point (mb˜ = 350 GeV). 8.1.2 Resonance reconstruction To extract the mass of the potential resonance, we reconstruct the bottom squark decay chain. First we reconstruct the top quark pair candidates using the lepton momenta, b-tagged jet mo- menta, and EmissT . Then these candidates are combined with the light-parton jets to reconstruct 24 8 Dilepton final state Pr ob ab ilit y de ns ity -510 -410 -310 (GeV) T pLeading light jet 50 100 150 200 250 300 350 400 (G eV ) Tp Se co nd le ad in g lig ht je t 50 100 150 200 250 300 350 400 CMS L = 19.5 fb−1 √s = 8 TeV Background-only fit to data Pr ob ab ilit y de ns ity -610 -510 -410 (GeV) T pLeading light jet 50 100 150 200 250 300 350 400 (G eV ) Tp Se co nd le ad in g lig ht je t 50 100 150 200 250 300 350 400 CMS MC Simulation √ s = 8 TeV 350 GeV b˜ signal point Figure 13: Two-dimensional light-parton jet pT distributions for (left) the background-only hy- pothesis fit to data and (right) the signal with mb˜ = 350 GeV. The scales are logarithmic and a line has been drawn along the diagonal to illustrate the ordering of the jets by pT. candidate bottom squark pairs. In tt decays resulting in two charged leptons, the momenta of the neutrinos can be determined by solving a quartic equation [70, 71], given the masses of the top quark and W boson. We assume that all EmissT in the event comes from neutrinos associated with leptonic W decays and fix all resonance masses to their on-shell values. We solve this quartic equation for both pairings of leptons and b-tagged jets yielding eight (possibly unphysical) solutions for the tt candidates in the event. If none of the eight solutions is physical, we vary the measured jet momenta within their resolutions using the procedure described in Ref. [72]. We sample 1000 times per event and choose the set of parameters with at least one physical solution for tt candidates and the smallest χ2 with respect to the measured jets. The resampling procedure improves the event selection efficiency by a factor of approximately 1.4 for both background and simulated signal events. The efficiency for finding real solutions in signal events that pass all other selection criteria is between 55% and 85%, depending on the particular signal model and mass point. The background efficiency is approximately 81%. The additional data events from the resampling procedure improve the signal sensitivity by 15%. The physical solutions for tt candidates are combined with the light-parton jets to form candi- date b˜ resonances, which can include up to 16 solutions. We select the pair with the smallest absolute value of the logarithm of the ratio of the two candidate resonance masses. For the remainder of the analysis, we substitute the average of these two masses, denoted by mtj, for the mass of each candidate within the chosen pair. It is bounded from below by the top quark mass. The resulting mtj distributions for background and signal are parameterized by different func- tions. The background distribution displays a turn-on behavior due to pT requirements on the leptons and jets, and a falling tail. It is modeled by the sum of a gamma distribution and a log-normal distribution, constrained to peak at the same value. There are four parameters in the background model that are left free in the fit: two widths, one peak position, and the rel- ative normalization of the two components. The signal is parameterized as the sum of two gamma distributions with parameters determined by a fit to the signal simulation. One of 8.2 Systematic uncertainties 25 the gamma functions models the correctly paired events and the other models the incorrectly paired events. Events that are correctly paired have a narrower mass width than those that are incorrectly paired, which also have a large tail extending to high mass. The success rate for cor- rectly matching the reconstructed and true objects increases from 30% to 50% with increasing bottom squark mass. In order to increase the sensitivity of the analysis, events are categorized into different regions according to p(2)T , as illustrated in Table 5. We label these as three signal-enhanced regions (denoted SR1–3) and one signal-depleted control region (denoted CR). The mass spectrum for each region is fit separately. We write the background distribution as ρSMmass(mtj|p(2)T ) and the signal distribution as ρsignalmass (mtj|p(2)T ). Table 5: Definition of signal and control regions in the dilepton analysis. Second-leading light-parton jet pT Region 30 < p(2)T < 50 GeV Control region (CR) 50 < p(2)T < 80 GeV Signal region 1 (SR1) 80 < p(2)T < 110 GeV Signal region 2 (SR2) p(2)T > 110 GeV Signal region 3 (SR3) Figure 14 shows the background-only hypothesis fits to the measured invariant mass spectra in the four light-parton jet regions (SR1–3 and CR). 8.2 Systematic uncertainties The fitted parameters of the background model are determined exclusively from data. Potential systematic uncertainties in the background model could arise from the choice of functional forms in Section 8.1.1. We perform studies to ensure that the model is sufficiently general to approximate the true shape of the distributions, with any differences negligible compared to the statistical uncertainties. The signal is modeled by a MC simulation and we consider several systematic uncertainties that can affect the light-parton jet pT spectrum and the mass reconstruction procedure. For each source of systematic uncertainty, we determine the change in the signal distribution parameters and encode this as a covariance matrix of all parameters. A multivariate normal distribution is constructed from the sum of all such covariance matrices and used to constrain the parameters of the signal distribution. Allowing the invariant mass parameters to vary within this con- straint has a negligible effect on the signal sensitivity. The final constraint function sets these parameters to their maximum likelihood values in signal simulation and does not allow them to vary. The light-parton jet parameters and selection efficiencies are allowed to vary within their uncertainties. The total selection efficiencies include the 10.9% branching fraction [73] of each of the top quarks to a leptonic final state and the production of tau leptons that decay to electrons or muons. We estimate that in bottom squark pair production with a dileptonic final state, about 60% of events produce leptons, b quarks, and light partons within the fiducial volume of the detector. About 85% of these events have the physics objects reconstructed and about 65% of those events have the jets from b quarks correctly tagged. The trigger efficiency is close to 95%. Approximately 90% of the remaining events pass all selection requirements. The total selec- tion efficiencies are approximately 2.0× 10−3, 4.0× 10−3, and 1.5× 10−3 for the ee, eµ, and µµ channels, respectively, and have a slight dependence on the b˜ mass. 26 8 Dilepton final state (GeV)t,jetm 200 400 600 800 1000 1200 Ev en ts / ( 2 5.7 G eV ) 0 100 200 300 400 500 Data function Fitted /dof = 39.9/36 2χ 200 400 600 800 1000 1200 R at io 0 0.5 1 1.5 2 CMS L = 19.5 fb−1 √s = 8 TeV (GeV)t,jetm 200 400 600 800 1000 1200 Ev en ts / ( 2 5.7 G eV ) 0 20 40 60 80 100 120 140 160 180 200 Data function Fitted /dof = 54.3/36 2χ 200 400 600 800 1000 1200 R at io 0 0.5 1 1.5 2 CMS L = 19.5 fb−1 √s = 8 TeV (GeV)t,jetm 200 400 600 800 1000 1200 Ev en ts / ( 2 5.7 G eV ) 0 10 20 30 40 50 60 Data function Fitted /dof = 51.0/36 2χ 200 400 600 800 1000 1200 R at io 0 0.5 1 1.5 2 CMS L = 19.5 fb−1 √s = 8 TeV (GeV)t,jetm 200 400 600 800 1000 1200 Ev en ts / ( 2 5.7 G eV ) 0 5 10 15 20 25 30 Data function Fitted /dof = 24.7/36 2χ 200 400 600 800 1000 1200 R at io 0 0.5 1 1.5 2 CMS L = 19.5 fb−1 √s = 8 TeV Figure 14: Reconstructed invariant mass distributions for data together with the result of the likelihood function maximization with signal cross section set to zero for the four light-parton jet regions defined in Table 5: CR (upper left), SR1 (upper right), SR2 (lower left) and SR3 (lower right). 8.3 The likelihood function 27 Potential discrepancies between lepton efficiencies in data and simulation are taken into ac- count varying lepton energies using scale factors, which are within 2% of unity when the pT is above 10 GeV [58, 59]. For electrons the energy is varied by 1.5% in the range 1.5 < |η| < 2.5 of the detector and 0.6% in the range |η| < 1.5; for muons this value is 0.2%. To match the efficiency for tagging b jets and the mistag rate for charm and light-flavor quarks between data and simulation, a scale factor is applied to the simulation. There is a correspond- ing uncertainty introduced by using these scale factors. We apply an additional systematic uncertainty to account for the change in efficiency and mistag rates of the b tagging algorithm when applying our specific selection criteria. Table 6 reports the uncertainty in the signal selection efficiency due to the different sources of systematic uncertainty for an example mass point. In combination these systematic uncertain- ties change the calculated upper limit on the cross section between 1% and 10%, depending on the mass point, compared to the upper limit calculated only with statistical uncertainties. The dominant systematic effect comes from variations of the jet energy scale. Table 6: Relative systematic uncertainty in the signal selection efficiency broken down by source of signal systematic uncertainty for a characteristic mass point. Signal simulation uncertainty mb˜ = 350 GeV Heavy-flavor SF for b tagging 4.9% Light-flavor SF for b tagging 4.7% Jet energy scale 4.6% Signal MC statistics 2.1% Jet energy resolution 1.8% Pileup 1.5% PDFs 1.0% MC b tagging efficiency for b jets 0.4% MC b tagging efficiency for c jets 0.3% MC b tagging efficiency for light-parton jets 0.5% Electron energy scale 0.2% Muon energy scale <0.1% Integrated luminosity 2.6% Total 9.2% 8.3 The likelihood function For both signal and background, we define a three-dimensional probability density function constructed using the two-dimensional light-parton jet distributions defined in Section 8.1.1 and the invariant mass distributions defined in Section 8.1.2. These three-dimensional distri- butions can be written as: ρSM3D (mtj, p (1) T , p (2) T ) = ρ SM mass(mtj|p(2)T )ρSM2D (p(1)T , p(2)T ) (15a) ρ signal 3D (mtj, p (1) T , p (2) T ) = ρ signal mass (mtj|p(2)T )ρsignal2D (p(1)T , p(2)T ), (15b) and the complete distribution is: ρtotal3D = (µ SMρSM3D + eLσsignalρsignal3D )/(µSM + eLσsignal). (16) Here µSM is the SM yield, e is the signal efficiency, L is the total integrated luminosity, and σsignal is the signal cross section. 28 8 Dilepton final state Constraints on the signal shape parameters are derived as described in Section 8.2. We write the constraint distribution as ρsyst. There are no constraints on the parameters describing the SM distribution. We construct an extended unbinned likelihood function from our data and these distributions as L(σsignal, θ) = ρsyst (µSM + eLσsignal)N exp (−µSM − eLσsignal) N! N ∏ i=0 ρtotal3D (mtj,i, p (1) T i, p (2) T i), (17) where N is the number of events in our sample; mtj,i, p (1) T i, and p (2) T i are, respectively, the mass, p(1)T , and p (2) T of the ith event; σ is the cross section for production of the bottom squark reso- nance pair; and θ is the set of all nuisance parameters included in ρsyst. 8.4 Results for the dilepton final state We observe consistency with the SM expectation and set limits on the cross section for each of the signal models. We construct unified intervals [74] on the signal cross section and observe only intervals with lower edges of zero. The upper edge of these intervals is considered as the upper limit. For each signal mass, pseudoexperiments generated using the background- only model are used to determine the distribution of upper limits in the absence of signal. The median of each distribution, along with intervals containing the central 68% and 95% of each distribution, is found. (GeV)b~m 200 250 300 350 400 450 500 550 600 650 jj) [pb ] t t → b~ b~ → (pp σ -310 -210 -110 1 10 210 310 410 SUSY cross section 95% CL observed limit 95% CL expected limit 68% band on exp. limit 95% band on exp. limit CMS L = 19.5 fb−1 √s = 8 TeV Figure 15: Observed and expected 95% CL upper limits on the production cross section of b˜ pairs, where the b˜ decays to a top and down (strange) quark via the RPV coupling λ′′331 (λ ′′ 332), as a function of b˜ mass derived using CLs intervals. The difference between observed and expected limits is correlated between neighboring signal points. The hatched region represents the theoretical uncertainty on the signal cross section. Figure 15 shows the observed and expected 95% CL upper limits on the production cross sec- tion of b˜ pairs using CLs intervals, assuming λ′′332 or λ′′331 is nonzero. The λ ′′ 332 and λ ′′ 331 couplings are from the RPV SUSY Lagrangian defined in Eq. (1) and control the branching fraction of the decay of a bottom squark to a top quark and a light down-type quark, s and d, respectively. We find that the median expected lower limit on the b˜ mass is 295 GeV with the central 68% of the limit distribution falling in the range 282–304 GeV. The measured 95% CL lower limit on the b˜ mass is 307 GeV. 29 Figure 16 shows the 350 GeV b˜ point where the invariant mass distributions peak in simi- lar locations to the background and the search sensitivity is mainly attributable to the two- dimensional light-parton jet distribution. The limits at large masses (mb˜ > 400 GeV) are correlated as the signal distributions are all concentrated in the SR3 region and the invariant mass shapes have a large overlap. (GeV)t,jetm 200 400 600 800 1000 1200 Ev en ts / ( 2 5.7 G eV ) 0 10 20 30 40 50 60 Data function Fitted component b~350 GeV /dof = 53.8/36 2χ 200 400 600 800 1000 1200 R at io 0 0.5 1 1.5 2 CMS L = 19.5 fb−1 √s = 8 TeV 80 GeV ≤ 2nd leading light jet pT < 110 GeV (GeV)t,jetm 200 400 600 800 1000 1200 Ev en ts / ( 2 5.7 G eV ) 0 5 10 15 20 25 30 Data function Fitted component b~350 GeV /dof = 27.5/36 2χ 200 400 600 800 1000 1200 R at io 0 0.5 1 1.5 2 CMS L = 19.5 fb−1 √s = 8 TeV 2nd leading light jet pT ≥ 110 GeV (GeV) T pLeading light jet 100 200 300 400 500 600 Ev en ts / ( 1 3.8 G eV ) 0 50 100 150 200 250 Data function Fitted component b~350 GeV /dof = 46.0/34 2χ 100 200 300 400 500 600 R at io 0 0.5 1 1.5 2 CMS L = 19.5 fb−1 √s = 8 TeV 2nd leading light jet pT ≥ 50 GeV Figure 16: Results of likelihood maximization with the signal cross section set to the calculated 95% CL upper limit of the 350 GeV b˜ point. The solid line shows the fitted function, the dashed line shows the signal component, and the points show the data. The ratio of the data to the fitted function is also shown. From left to right, these are the invariant mass distribution in SR2, the invariant mass distribution in SR3, and the leading light jet pT distribution for events with p(2)T > 50 GeV. 9 Four-lepton final state via strong production Multilepton final states, having clean identification criteria and well-understood background sources, present an important sector for searching for BSM physics. The SUSY models that contain leptonic or semileptonic RPV couplings produce multilepton events with very little EmissT . In this section we describe a search designed explicitly to look for evidence of nonzero leptonic RPV couplings λ122 and λ121 in models with a winolike neutralino and different types of strong SUSY production: gluino pairs, top-squark pairs and squark pairs composed of an equal mix- ture of up, down, strange, and charm. Examples of this production are shown in the diagrams in Fig. 17. In these events, the final state includes four light charged leptons (electrons and muons), and two neutrinos from the RPV decay of the neutralinos, as well as two or more jets from cascade decays of the strongly produced superpartners. The four leptons have a total charge of zero. The main SM contribution to the multilepton final state comes from processes with Z bosons decaying into lepton pairs. In contrast, the four leptons produced in RPV neutralino decays are not expected to accumulate at the Z boson mass in dilepton spectra. The leptons from both the SM and RPV decays are produced at the proton-proton interaction vertex; we refer to these leptons as prompt. We select events from a dilepton-triggered data set with exactly four isolated light leptons 30 9 Four-lepton final state via strong production P1 P2 q˜ ¯˜q χ˜01 χ˜01 q ν ` ` ν ` ` q¯ P1 P2 g˜ g˜ χ˜01 χ˜01 q q¯ ν ` ` ν ` ` q q¯ Figure 17: Diagrams of two simplified models with RPV [43]. Left: A simplified model fea- turing squark pair production, with q˜ → qχ˜01, and m(g˜)  m(q˜). Right: A simplified model featuring gluino pair production, with g˜ → qqχ˜01, and m(q˜)  m(g˜). In both models the neutralinos decay to two charged leptons and a neutrino via a RPV term. containing at least one opposite-sign, same-flavor (OSSF) lepton pair. For events with four isolated light leptons, the dilepton triggers are close to 100% efficient. We define M1 as the invariant mass of the OSSF lepton pair. If there are two OSSF lepton pairs in the event, M1 is the mass of the pair that is closest to the Z boson mass of 91 GeV. The invariant mass of the remaining lepton pair in the event is defined as M2. Then we classify each mass according to: • “below Z”: M < 75 GeV; • “at Z”: 75 < M < 105 GeV; • “above Z”: M > 105 GeV. This defines nine regions reflecting different kinds of resonant and nonresonant four-lepton production. Figure 18 presents the distribution of expected yields in the (M1, M2) plane for different SM processes and for an example RPV scenario. Based on the occupancy of different regions for typical ZZ production events, we define our signal region as “M1 above Z” or “M1 below Z and M2 above Z.” 9.1 Backgrounds The presence of four prompt leptons in the final state is a strong discriminant. The SM processes contributing to this signature are: • processes producing exactly four prompt leptons: ZZ(∗) and ttZ; • more rare processes producing four or more prompt leptons: ttWW, WWZ, WZZ, and ZZZ; • processes such as WZ+jets, ttW+jets, and ttZ that can produce three genuine prompt leptons and one candidate (nonprompt lepton) that arises primarily from the decay of a hadron of heavy flavor or from particle misidentification; • Drell–Yan production with two extra nonprompt leptons. We use MC simulated samples to estimate the backgrounds with four or more genuine prompt leptons. The primary background contribution is ZZ production. We normalize this sample to the ZZ production cross section measured by CMS [75], so our observation in the region with both M1 and M2 at Z is correlated with that measurement. Based on the consistency of our observations with the predictions in control regions, we assign a systematic uncertainty of 25% 9.1 Backgrounds 31 e xp ec te d ev en ts fo r 1 9. 5/ fb -410 -310 -210 -110 1 10 (GeV)1M 50 100 150 200 (G eV ) 2 M 50 100 150 200 = 8 TeVsCMS Simulation Z Z→pp e xp ec te d ev en ts fo r 1 9. 5/ fb -410 -310 -210 -110 (GeV)1M 50 100 150 200 (G eV ) 2 M 50 100 150 200 = 8 TeVsCMS Simulation Zt t →pp e xp ec te d ev en ts fo r 1 9. 5/ fb -410 -310 (GeV)1M 50 100 150 200 (G eV ) 2 M 50 100 150 200 = 8 TeVsCMS Simulation W Wt t →pp e xp ec te d ev en ts fo r 1 9. 5/ fb 200 400 600 800 (GeV)1M 50 100 150 200 (G eV ) 2 M 50 100 150 200 = 8 TeVsCMS Simulation 1 0χ∼q 1 0χ∼q → q~q~ →pp ν -l+ l→ 1 0χ∼ =425 GeVq~m =225 GeV 1 0χ∼m Figure 18: Expected shapes of the (M1, M2) distribution for the ZZ background (upper left), ttZ background (upper right), ttWW background (lower left), as well as for the model of Fig. 17 (left) with mq˜ = 425 GeV and mχ˜01 = 225 GeV (lower right). to the MC predictions of ZZ production contributions. The yields of the remaining background sources with four or more genuine leptons are estimated from MC simulation normalized to the theoretical cross sections, which are known to 5–10%. This uncertainty, combined with the 25% systematic uncertainty in the mismodeling observed in the ZZ sample, motivates the assignment of a conservative systematic uncertainty of 50% to background contributions of other rare processes. We estimate the contribution of nonprompt leptons with a sample of events containing an OSSF lepton pair with the requirement EmissT < 30 GeV, to increase the purity of the Drell-Yan control sample. The left plot in Fig. 19 presents the dilepton mass distribution of these events in which a single jet with pT > 30 GeV is present. The right plot presents the OSSF dilepton mass distri- bution for events in which a third isolated lepton is present. By comparing the yield of events in the Z peak in the left plot with that of events in the Z peak less the expected prompt lepton con- tribution in the right plot, we obtain a jet-to-lepton nonprompt rate of (0.103± 0.003 (stat))%. Figure 20 shows the contribution of events with three isolated leptons and one jet to the differ- ent regions of this search. These numbers are used as a reference for determining the contri- bution of nonprompt leptons to the four-lepton selection using the nonprompt rate evaluation of 0.103%. The jet momentum is used to calculate M2 and classify the events in Fig. 20. As the actual nonprompt lepton takes only a fraction of the heavy-flavor jet momentum, this proce- dure overestimates M2, which brings significant systematic uncertainty into this background estimation. We account for this by assigning a 100% systematic uncertainty to this contribution. Despite this large uncertainty, we tested that varying the estimation of this background down to zero keeps the exclusion result stable within 1%. 32 9 Four-lepton final state via strong production (GeV)llM 20 40 60 80 100 120 140 Ev en ts / 5 G eV 0 200 400 600 800 1000 1200 1400 1600 1800 310× -1 =19.5 fb = 8 TeV LsCMS 2l +jet MET<30 GeV Data Rare ZZ+jets WZ+jets +jetstt DY (GeV)llM 20 40 60 80 100 120 140 Ev en ts / 5 G eV 0 200 400 600 800 1000 1200 1400 1600 1800 -1 =19.5 fb = 8 TeV LsCMS 3l MET<30 GeV Data Rare ZZ+jets WZ+jets +jetstt DY Figure 19: Invariant mass of OSSF lepton pairs in events with EmissT < 30 GeV. Data are over- laid with contributions from different sources, predicted by simulation. Left: Events with two isolated leptons and exactly one jet with pT > 30 GeV. Right: Invariant mass of the two OSSF isolated leptons closest to the Z boson mass in events with three isolated leptons. all Z2M Z2M =Z1M Z1M =Z2M >Z1M >Z2M >Z1M Ev en ts / 1 G eV 0 500 1000 1500 2000 2500 3000 3500 -1 =19.5 fb = 8 TeV LsCMS 3l +jet Data Rare ZZ+jets WZ+jets +jetstt DY Figure 20: Number of events with three isolated leptons and one jet to different (M1, M2) re- gions. Data are overlaid with predicted contributions from different sources. 9.2 Observations Table 7 shows the observed number of events in different regions together with the expectations from SM processes. The observations are consistent with the SM background expectations. Based on the observations in the search region (“M1 above Z”, or “M1 below Z and M2 above Z”), the 95% CL upper limit on cross section times integrated luminosity times efficiency (σLε) for any physics process beyond the SM contributing to this search region is 3.4 events. The expected upper limit for this observation is 4.7 events. 9.3 Generalized efficiency 33 Table 7: Expected background contributions from different SM sources and experimentally observed events in all analysis regions. The ZZ prediction in the region with both M1 and M2 falling in the “at Z” region is based on simulation normalized to the CMS ZZ production cross section measurement and is therefore correlated with the observation in this analysis. The uncertainties take into account statistical and systematic contributions, combined quadratically. Mass ranges are given in GeV. M1 < 75 75 < M1 < 105 M1 > 105 ZZ 10±2 32±8 1.0 ±0.2 Rare 0.3±0.1 3±1 0.01±0.01 M2 < 75 Nonprompt 0.3±0.3 0.8±0.8 0.06±0.06 All backgrounds 10±2 35±8 1.0±0.2 Observed 14 30 1 ZZ 0.10±0.03 150 0.05±0.01 Rare 0.12±0.05 3±1 0.06±0.03 75 < M2 < 105 Nonprompt 0.3±0.3 1±1 0.05±0.05 All backgrounds 0.5±0.3 153 0.16±0.06 Observed 0 160 0 ZZ 0.8±0.2 15±4 0.30±0.07 Rare 0.3±0.1 3±1 0.12±0.05 M2 > 105 Nonprompt 0.4±0.4 0.7±0.7 0.05±0.05 All backgrounds 1.4±0.5 18±4 0.5±0.1 Observed 0 20 0 9.3 Generalized efficiency To understand how our results can be applied to a generic SUSY model with RPV, we study how the signal efficiency varies across about 7300 sets of RPC phenomenological MSSM (pMSSM) [76] model parameters, each one containing 10000 events, selected to fulfill different pre-CMS observations [77]. The RPV leptonic decay of the pair of neutralinos leads to four prompt leptons. The kinematics of these leptons are generally driven by the momentum distribution of the decaying neutrali- nos and by their masses. In most of the SUSY scenarios the lepton momentum is well above the required threshold, which results in high efficiency before the isolation requirement. The isolation efficiency depends on the hadronic activity in the event. We find that nearly all models have a four-lepton isolation efficiency in the range 50–100%. Therefore we consider an efficiency band spanning this range, using a 30% uncertainty when combining with other uncertainties. 9.4 Interpretations of the four-lepton results Once an upper limit on σLε is extracted from the observations, and the efficiency is evaluated, the corresponding limit on the cross section, σSUSYtotal can be calculated. The cross section exclusion limits for squark pair production with leptonic RPV SUSY (Fig. 17, left) are presented in Fig. 21. They are obtained using the exclusion limit on the number of events in the two signal regions. Here the bands show only the component of the uncertainty attributable to the lepton ID efficiencies. The edges of the bands correspond to the two effi- ciency profiles, that is, the width reflects the uncertainty introduced by changing the amount of hadronic activity in the signal model. The left plot of Fig. 21 presents results for neutralinos decaying exclusively to electrons or muons. The right plot takes the appropriate mixture of 34 9 Four-lepton final state via strong production electrons and muons in neutralino decays for λ121 6= 0 and λ122 6= 0 cases into account. The experimental observations described in Section 9.2, together with the pMSSM-based [76] efficiency estimation, drive the exclusion for the cross section of total leptonic RPV SUSY pro- duction, which is presented in Fig. 22. The bands in Fig. 22 correspond to varying the four- lepton isolation efficiency from 50% to 100%. Note that this band covers RPV models with a wide range of underlying RPC SUSY models. We demonstrate that the kinematic efficiency is controlled by the neutralino mass and only weakly dependent on the neutralino momentum. (GeV) 1 0χ∼ m 200 400 600 800 1000 1200 (fb ) SU SY to ta l σ -110 1 10 95% CL limits ν−e+e→0χObserved: ν−e+e→0χExpected: ν−µ+µ→0χObserved: ν−µ+µ→0χExpected: -1 =19.5 fb = 8 TeV LsCMS 1 0χ∼q 1 0χ∼q → q~q~ →pp (GeV) 1 0χ∼ m 200 400 600 800 1000 1200 (fb ) SU SY to ta l σ -110 1 10 95% CL limits 0≠121λObserved: 0≠121λExpected: 0≠122λObserved: 0≠122λExpected: -1 =19.5 fb = 8 TeV LsCMS ν −l+ l→ 1 0χ∼, 1 0χ∼q 1 0χ∼q → q~q~ →pp Figure 21: The 95% CL upper limit on the total cross section of the process sketched in Fig. 17 (left), as a function of the mass of the RPV decaying neutralino. The band corresponds to the component of uncertainty from the lepton efficiency. Results are shown for neutralinos decaying exclusively to electrons or muons (left) and decaying with the appropriate lepton flavor mixture corresponding to λ121 6= 0 and λ122 6= 0 RPV scenarios (right). (GeV) 1 0χ∼ m 200 400 600 800 1000 1200 (fb ) SU SY to ta l σ -110 1 10 95% CL limits ν−e+e→0χObserved: ν−e+e→0χExpected: ν−µ+µ→0χObserved: ν−µ+µ→0χExpected: -1 =19.5 fb = 8 TeV LsCMS [0.5, 1]∈ ISOε (GeV) 1 0χ∼ m 200 400 600 800 1000 1200 (fb ) SU SY to ta l σ -110 1 10 95% CL limits 0≠121λObserved: 0≠121λExpected: 0≠122λObserved: 0≠122λExpected: -1 =19.5 fb = 8 TeV LsCMS [0.5, 1]∈ ISOε Figure 22: The 95% CL upper limit on total cross sections for generic SUSY models. The band corresponds to the efficiency uncertainty from the various pMSSM models. Results are shown for neutralinos decaying exclusively to electrons or muons (left) and decaying with the appro- priate lepton flavor mixture corresponding to λ121 6= 0 and λ122 6= 0 RPV scenarios (right). To further convert the cross section limit into a mass exclusion, we consider several SUSY pro- duction mechanisms: gluino pair production, top-squark pair production, and first and second generation squark pair production. Using these total cross sections as a function of the mass of the corresponding SUSY particle, we convert the cross section limit bands in Figs. 21 and 22 into exclusion bands for the masses of the parent SUSY particles, as a function of the LSP mass. 35 These results are presented in Fig. 23. A 30% theoretical uncertainty in NLO+NLL calculations of SUSY production cross sections is included in the uncertainty band, along with the efficiency uncertainties. (GeV) 1 0χ∼ m 200 400 600 800 1000 1200 95 % C L m as s LL (T eV ) 0 0.5 1 1.5 2 0≠ 121λ exclusionq~m exclusiong~m exclusiont~m 0≠ 122λ exclusionq~m exclusiong~m exclusiont~m Expected exclusions -1 =19.5 fb = 8 TeV LsCMS ν −l+ l→ 1 0χ∼, 1 0χ∼q 1 0χ∼q → q~q~ →pp ν −l+ l→ 1 0χ∼, 1 0χ∼q q 1 0χ∼q q→ g~g~ →pp (GeV) 1 0χ∼ m 200 400 600 800 1000 1200 95 % C L m as s LL (T eV ) 0 0.5 1 1.5 2 0≠ 121λ exclusionq~m exclusiong~m exclusiont~m 0≠ 122λ exclusionq~m exclusiong~m exclusiont~m Expected exclusions -1 =19.5 fb = 8 TeV LsCMS ν −l+ l→LSP [0.5, 1]∈ ISOε Figure 23: Lower limits on the masses of intermediate-state particles for various SUSY pro- duction mechanisms as a function of the neutralino mass when the neutralino is the LSP. The limits shown in the left plot on the gluino mass, top-squark mass, and squark mass are for gluino pair production, top-squark pair production, and first and second generation squark pair production, respectively. The first and second generation squark mass exclusion is ob- tained from the cross section limit of Fig. 21 (right), combined with the corresponding cross section for this model, displayed in Fig. 17 (left). The other exclusions are obtained in a similar fashion. The limits shown in the right plot are obtained from the generic total RPV SUSY cross section limit from Fig. 22 (right). The bands include the uncertainty in SUSY production cross sections, along with the efficiency uncertainties. The models below the diagonal line were not investigated as the neutralino is no longer the LSP. The weaker limits at low neutralino mass are a consequence of the low efficiency in this region. For the cases λ121 6= 0 or λ122 6= 0, and λ sufficiently large to lead to prompt neutralino decays, this model is excluded at a 95% CL for gluino masses below about 1.4 TeV for a neutralino mass higher than 400 GeV. With only top-squark production, a top-squark mass below about 950 GeV is generally excluded. The squark mass exclusion depends significantly on the gluino mass in the model, as the squark pair production cross section increases as gluino mass decreases. For the 2.4 TeV gluino benchmark point, squarks with a mass below about 1.6 TeV are excluded. 10 Multilepton final state via electroweak production When leptonic or semileptonic RPV couplings are allowed, LSP pair production can result in a final state with three or more leptons. In this section we expand the results from a previous search [27] to this class of RPV models. Our focus here is on models with a Higgsino- or wino- like LSP, and we explore the leptonic RPV couplings λ122, λ123, and λ233 and the semileptonic couplings λ′131, λ ′ 233, λ ′ 331, and λ ′ 333. 10.1 Event selection We select events with three or more leptons (including up to one hadronically decaying tau lepton τh) that are accepted by a trigger requiring two light leptons, which may be electrons or muons. Any OSSF pair of electrons or muons must have an invariant mass m`` > 12 GeV, removing low-mass bound states and most γ∗ → `+`− events. 36 10 Multilepton final state via electroweak production We define 32 exclusive signal regions depending on the total number of leptons, the number of τh candidates in the event, the number of b-tagged jets in the event, the number of OSSF pairs, and in the cases where an OSSF pair is present, different bins in the dilepton invariant mass, which is calculated only for light leptons (electrons and muons). When there is more than one OSSF pair, if any of these dilepton masses are within 15 GeV of the Z mass, the event is called on-Z. The binning is shown in the first four columns of Table 8. Each of these 32 signal regions is divided into five different ST bins, where ST is defined as the sum of the missing transverse momentum and the scalar sum of jet and charged lepton pT: ST =[0–300], [300–600], [600–1000], [1000–1500], and [>1500] GeV. 10.2 Backgrounds Multilepton searches have two main sources of background, the first arising from processes that produce genuine multilepton events. The most abundant examples are WZ and ZZ pro- duction, but rare processes such as ttW and ttZ also contribute. Samples simulating WZ and ZZ production have been validated in control regions in data. For the rarer background processes, we rely solely on simulation. The second source of background originates from objects that are misclassified as prompt, iso- lated leptons, but are actually hadrons, leptons from a hadron decay, etc. Misidentified leptons are classified in three categories: misidentified light leptons, misidentified τh leptons, and light leptons originating from asymmetric internal conversions. The methods used to estimate these backgrounds are described in more detail in Ref. [24]. We estimate the contribution of misidentified light leptons by measuring the number of iso- lated tracks and applying a scale factor between isolated leptons and isolated tracks. These scale factors are measured in control regions that contain leptonically decaying Z bosons and a third, isolated track, as well as in control regions with opposite-sign, different-flavor leptons, an isolated track, and a b-tagged jet, which are dominated by tt production. The scale factor is then the probability for the third track to pass the lepton identification criteria. We find the scale factors to be (0.9± 0.2)% for electrons and (0.7± 0.2)% for muons. The scale factors are applied to the number of events in the sideband region with two light leptons and an isolated track. The scale factors depend on the heavy-flavor content in the different signal regions. We parameterize this dependence as a function of the impact parameter distribution of nonisolated tracks. The tt contribution is taken from simulation. The τh misidentification rate is measured in a jet-dominated data control sample by comparing the number of τh candidates in the isolated region defined by ET, cone < 2 GeV to the number of nonisolated τh candidates, which have 6 < ET, cone < 15 GeV. We measure the average misidentification rate as 15% with a systematic uncertainty of 30% based on the variation in different control samples. We apply this scale factor to the number of events in the sideband region with two light leptons and one nonisolated τh candidate. Another source of background leptons is internal conversions of a virtual photon to a dilepton pair. We measure the conversion factors of photons to light leptons in a control region with low EmissT and low hadronic activity. The ratio of the number of ` +`−`± candidates to the number of `+`−γ candidates in Z boson decays defines the conversion factor, which is 2.1% ± 1.0% (0.5%± 0.3%) for electrons (muons). 10.3 Systematic uncertainties We scale the WZ and ZZ simulation samples to match the data in control regions. The overall systematic uncertainty in WZ and ZZ contributions to the signal regions varies between 15% 10.4 Results and interpretations for the multilepton final state 37 and 30%, depending on the kinematics, and is the combination of the normalization uncertain- ties with resolution uncertainties [27]. An uncertainty of 50% for the tt background contribution is due to the low event counts in the isolation distributions in high-ST bins, which are used to validate the misidentification rate. We apply a 50% uncertainty to the normalization of all rare processes to cover cross section and PDF uncertainties. 10.4 Results and interpretations for the multilepton final state The observed and expected yields in the regions described above are shown in Table 8. We find good agreement between the SM predictions and observed data. Table 8: Expected and observed yields for three- and four-lepton events. The channels are split by the number of leptons (N`), the number of τh candidates (Nτ), whether the event contains b- tagged jets (Nb), the number of OSSF pairs (NOSSF), binning in the dilepton invariant mass (m`` of light leptons only), and the ST. Events are considered on-Z if 75 < m`` < 105 GeV. Expected yields are the sum of simulation and estimates of backgrounds from data in each channel. The channels are mutually exclusive. The uncertainties include statistical and systematic uncertain- ties. The ST and m`` values are given in GeV. Reproduced from Ref. [27]. N` Nτ Nb NOSSF m`` 0 < ST < 300 300 < ST < 600 600 < ST < 1000 1000 < ST < 1500 ST > 1500 obs exp obs exp obs exp obs exp obs exp 4 0 0 0 — 0 0.06 ± 0.06 0 0.09 ± 0.07 0 0.00 ± 0.03 0 0.00 ± 0.03 0 0.00 ± 0.03 4 0 1 0 — 0 0.00 ± 0.03 0 0.00 ± 0.03 0 0.06 ± 0.05 0 0.00 ± 0.03 0 0.00 ± 0.03 4 0 0 1 on-Z 2 3.1 ± 0.90 5 1.9 ± 0.48 0 0.44 ± 0.16 1 0.06 ± 0.06 0 0.00 ± 0.03 4 0 1 1 on-Z 2 0.07 ± 0.05 2 1.1 ± 0.53 0 0.57 ± 0.30 0 0.12 ± 0.09 0 0.02 ± 0.03 4 0 0 1 off-Z 2 0.48 ± 0.18 0 0.27 ± 0.11 0 0.07 ± 0.05 0 0.00 ± 0.02 0 0.00 ± 0.03 4 0 1 1 off-Z 0 0.04 ± 0.04 0 0.34 ± 0.17 0 0.06 ± 0.08 0 0.04 ± 0.04 0 0.00 ± 0.03 4 0 0 2 on-Z 135 120 ± 29 26 43 ± 10 4 6.0 ± 2.0 1 0.63 ± 0.26 0 0.06 ± 0.04 4 0 1 2 on-Z 1 1.0 ± 0.27 4 3.2 ± 1.1 1 1.1 ± 0.39 0 0.11 ± 0.06 0 0.04 ± 0.04 4 0 0 2 off-Z 7 8.3 ± 2.3 3 1.1 ± 0.30 0 0.11 ± 0.05 0 0.01 ± 0.02 0 0.00 ± 0.02 4 0 1 2 off-Z 0 0.18 ± 0.07 1 0.22 ± 0.11 0 0.15 ± 0.08 0 0.00 ± 0.03 0 0.00 ± 0.03 4 1 0 0 — 2 1.1 ± 0.46 1 0.54 ± 0.20 0 0.12 ± 0.12 0 0.00 ± 0.03 0 0.00 ± 0.03 4 1 1 0 — 0 0.26 ± 0.16 0 0.29 ± 0.13 0 0.13 ± 0.11 0 0.01 ± 0.02 0 0.00 ± 0.03 4 1 0 1 on-Z 43 42 ± 11 10 12 ± 3.1 0 1.8 ± 0.63 0 0.11 ± 0.07 0 0.02 ± 0.03 4 1 1 1 on-Z 2 1.0 ± 0.40 2 1.7 ± 0.5 0 0.78 ± 0.33 0 0.04 ± 0.04 0 0.01 ± 0.03 4 1 0 1 off-Z 18 8.4 ± 2.2 4 2.1 ± 0.52 2 0.48 ± 0.18 0 0.13 ± 0.08 0 0.01 ± 0.03 4 1 1 1 off-Z 1 0.64 ± 0.31 0 1.2 ± 0.44 0 0.30 ± 0.13 0 0.02 ± 0.03 0 0.00 ± 0.03 3 0 0 0 — 72 80 ± 23 32 27 ± 11 3 3.1 ± 1.00 0 0.22 ± 0.18 0 0.07 ± 0.06 3 0 1 0 — 37 33 ± 16 42 39 ± 19 2 5.0 ± 2.0 0 0.36 ± 0.14 0 0.06 ± 0.07 3 0 0 1 on-Z 4255 4400 ± 690 669 740 ± 170 106 110 ± 41 11 15 ± 6.9 3 1.3 ± 0.76 3 0 1 1 on-Z 140 150 ± 25 122 110 ± 25 16 25 ± 7.0 2 3.3 ± 1.2 1 0.32 ± 0.22 3 0 0 1 m`` < 75 617 640 ± 100 84 86 ± 21 14 11 ± 3.6 0 1.2 ± 0.39 1 0.12 ± 0.09 3 0 1 1 m`` < 75 62 74 ± 28 52 57 ± 23 4 8.3 ± 2.7 1 0.69 ± 0.28 0 0.08 ± 0.06 3 0 0 1 m`` > 105 180 200 ± 34 63 66 ± 12 13 10 ± 2.5 2 1.1 ± 0.40 0 0.16 ± 0.09 3 0 1 1 m`` > 105 17 17 ± 6.5 36 35 ± 14 7 7.4 ± 2.5 0 0.54 ± 0.23 0 0.08 ± 0.05 3 1 0 0 — 1194 1300 ± 330 289 290 ± 130 26 28 ± 12 2 2.6 ± 1.3 0 0.23 ± 0.20 3 1 1 0 — 316 330 ± 160 410 480 ± 240 46 58 ± 28 2 3.9 ± 2.0 0 0.46 ± 0.32 3 1 0 1 on-Z 49916 49000 ± 15000 2099 2700 ± 770 108 70 ± 17 9 6.0 ± 1.6 0 0.33 ± 0.18 3 1 1 1 on-Z 795 830 ± 230 325 280 ± 74 17 17 ± 4.8 1 1.8 ± 0.64 0 0.30 ± 0.14 3 1 0 1 m`` < 75 10173 9200 ± 2700 290 280 ± 72 21 11 ± 3.5 1 0.97 ± 0.44 0 0.04 ± 0.06 3 1 1 1 m`` < 75 297 290 ± 97 167 170 ± 87 14 12 ± 6.0 0 1.1 ± 0.74 0 0.06 ± 0.08 3 1 0 1 m`` > 105 1620 1700 ± 480 285 370 ± 96 21 23 ± 7.2 1 1.4 ± 0.61 0 0.22 ± 0.23 3 1 1 1 m`` > 105 97 79 ± 36 169 190 ± 94 23 28 ± 14 1 2.2 ± 1.3 0 0.20 ± 0.18 We explore SUSY models featuring production of electroweak boson superpartners, where these sparticles can either be wino- or Higgsino-like. We investigate cases where one leptonic or semileptonic coupling is nonzero. The selected models with leptonic RPV contain sleptons with masses of 1.5 TeV to mediate the decay. In the semileptonic models, these mediators are bottom and top squarks that also have masses of 1.5 TeV. In all cases, we assume that the decays 38 10 Multilepton final state via electroweak production of the neutralinos are prompt, and the difference between the neutralino and chargino masses is roughly 1 GeV. Figure 24 contains plots that show where we exclude wino- and Higgsino-like neutralinos when there are nonzero RPV couplings, which from top to bottom are λ122, λ123, and λ233. In the models with winolike neutralinos, the lower limit on the neutralino mass ranges from ap- proximately 700 to 875 GeV. In the case of Higgsino-like neutralinos, the range is much larger, from approximately 300 to 900 GeV. We show limits in models that have Higgsino-like neutralinos with semileptonic RPV couplings λ′131, λ ′ 233, λ ′ 331, and λ ′ 333 in Fig. 25. We also feature two different values of tan β (2 and 40), which changes the relative coupling between the Higgsino and the mediator particles, which in turn changes the branching fraction to multileptons. The sensitivity of this analysis is best for the couplings that produce a light lepton in the final state (the first and third rows of Fig. 25). (GeV)χ∼m 100 200 300 400 500 600 700 800 900 1000 B (p b) × σ -410 -310 -210 -110 1 10 210 CMS -1 = 8 TeV, L = 19.5 fbs 122λWino RPV Observed 95% CL limts Expected 95% CL limits σ1±Expected σ2±Expected Theory (NLO+NLL) Theoretical Uncertainty (GeV)χ∼m 100 200 300 400 500 600 700 800 900 1000 B (p b) × σ -410 -310 -210 -110 1 10 210 CMS -1 = 8 TeV, L = 19.5 fbs 122λHiggsino RPV Observed 95% CL limts Expected 95% CL limits σ1±Expected σ2±Expected Theory (NLO+NLL) Theoretical Uncertainty (GeV)χ∼m 100 200 300 400 500 600 700 800 900 1000 B (p b) × σ -410 -310 -210 -110 1 10 210 CMS -1 = 8 TeV, L = 19.5 fbs 123λWino RPV Observed 95% CL limts Expected 95% CL limits σ1±Expected σ2±Expected Theory (NLO+NLL) Theoretical Uncertainty (GeV)χ∼m 100 200 300 400 500 600 700 800 900 1000 B (p b) × σ -410 -310 -210 -110 1 10 210 CMS -1 = 8 TeV, L = 19.5 fbs 123λHiggsino RPV Observed 95% CL limts Expected 95% CL limits σ1±Expected σ2±Expected Theory (NLO+NLL) Theoretical Uncertainty (GeV)χ∼m 100 200 300 400 500 600 700 800 900 1000 B (p b) × σ -410 -310 -210 -110 1 10 210 CMS -1 = 8 TeV, L = 19.5 fbs 233λWino RPV Observed 95% CL limts Expected 95% CL limits σ1±Expected σ2±Expected Theory (NLO+NLL) Theoretical Uncertainty (GeV)χ∼m 100 200 300 400 500 600 700 800 900 1000 B (p b) × σ -410 -310 -210 -110 1 10 210 CMS -1 = 8 TeV, L = 19.5 fbs 233λHiggsino RPV Observed 95% CL limts Expected 95% CL limits σ1±Expected σ2±Expected Theory (NLO+NLL) Theoretical Uncertainty Figure 24: Upper 95% CL cross section times branching ratio limits as a function of the neu- tralino mass in models with wino production (left) and Higgsino production (right), assum- ing nonzero λijk couplings: λ122 (top), λ123 (middle), and λ233 (bottom). The decays proceed promptly through slepton mediators. 10.4 Results and interpretations for the multilepton final state 39 (GeV)χ∼m 100 200 300 400 500 600 700 800 900 1000 B (p b) × σ -410 -310 -210 -110 1 10 210 CMS -1 = 8 TeV, L = 19.5 fbs ) = 2β tan(131 'λHiggsino RPV Observed 95% CL limts Expected 95% CL limits σ1±Expected σ2±Expected Theory (NLO+NLL) Theoretical Uncertainty (GeV)χ∼m 100 200 300 400 500 600 700 800 900 1000 B (p b) × σ -410 -310 -210 -110 1 10 210 CMS -1 = 8 TeV, L = 19.5 fbs ) = 40β tan(131 'λHiggsino RPV Observed 95% CL limts Expected 95% CL limits σ1±Expected σ2±Expected Theory (NLO+NLL) Theoretical Uncertainty (GeV)χ∼m 100 200 300 400 500 600 700 800 900 1000 B (p b) × σ -410 -310 -210 -110 1 10 210 CMS -1 = 8 TeV, L = 19.5 fbs ) = 2β tan(331 'λHiggsino RPV Observed 95% CL limts Expected 95% CL limits σ1±Expected σ2±Expected Theory (NLO+NLL) Theoretical Uncertainty (GeV)χ∼m 100 200 300 400 500 600 700 800 900 1000 B (p b) × σ -410 -310 -210 -110 1 10 210 CMS -1 = 8 TeV, L = 19.5 fbs ) = 40β tan(331 'λHiggsino RPV Observed 95% CL limts Expected 95% CL limits σ1±Expected σ2±Expected Theory (NLO+NLL) Theoretical Uncertainty (GeV)χ∼m 100 200 300 400 500 600 700 800 900 1000 B (p b) × σ -410 -310 -210 -110 1 10 210 CMS -1 = 8 TeV, L = 19.5 fbs ) = 2β tan(233 'λHiggsino RPV Observed 95% CL limts Expected 95% CL limits σ1±Expected σ2±Expected Theory (NLO+NLL) Theoretical Uncertainty (GeV)χ∼m 100 200 300 400 500 600 700 800 900 1000 B (p b) × σ -410 -310 -210 -110 1 10 210 CMS -1 = 8 TeV, L = 19.5 fbs ) = 40β tan(233 'λHiggsino RPV Observed 95% CL limts Expected 95% CL limits σ1±Expected σ2±Expected Theory (NLO+NLL) Theoretical Uncertainty (GeV)χ∼m 100 200 300 400 500 600 700 800 900 1000 B (p b) × σ -410 -310 -210 -110 1 10 210 CMS -1 = 8 TeV, L = 19.5 fbs ) = 2β tan(333 'λHiggsino RPV Observed 95% CL limts Expected 95% CL limits σ1±Expected σ2±Expected Theory (NLO+NLL) Theoretical Uncertainty (GeV)χ∼m 100 200 300 400 500 600 700 800 900 1000 B (p b) × σ -410 -310 -210 -110 1 10 210 CMS -1 = 8 TeV, L = 19.5 fbs ) = 40β tan(333 'λHiggsino RPV Observed 95% CL limts Expected 95% CL limits σ1±Expected σ2±Expected Theory (NLO+NLL) Theoretical Uncertainty Figure 25: Upper 95% CL cross section times branching fraction limits as a function of the neutralino mass in models with Higgsino production and nonzero λ′ijk couplings with (left) tan β = 2 and (right) tan β = 40: λ′131 (top), λ ′ 331 (second row), λ ′ 233 (third row), and λ ′ 333 (bot- tom). The decays proceed promptly through bottom and top squark mediators. The couplings of the Higgsino to the mediator particles varies with tan β, affecting the branching fraction to multileptons and the acceptance. 40 11 Summary 11 Summary This paper explores a variety of final states where R-parity-violating (RPV) supersymmetry could appear. Using data samples corresponding to 19.5 fb−1 of proton-proton collision data collected with the CMS detector at √ s = 8 TeV, no discrepancies from the expectations of the standard model are found. Limits on the masses of supersymmetric particles are set at 95% confidence level in several models that exhibit different RPV couplings and contain different lightest supersymmetric particles (LSP). The consequences of minimal flavor violation are explored in analyses that consider pair pro- duction of either gluinos or bottom squarks. The b-tagged and total jet multiplicity distribu- tions are used to set limits on the mass of a gluino that decays to a top, a bottom, and a strange quark via the superpotential coupling λ′′332. Gluinos with masses less than 0.98 TeV are ex- cluded. Using a search region characterized by one lepton and high multiplicity of jets and b-tagged jets, gluinos with masses less than 1.03 TeV are excluded in the same model. Another model assumes pair production of a bottom squark LSP that decays via λ′′332 or λ′′331. Using the reconstructed resonance mass distribution in the dilepton final state, bottom squark production is excluded for masses less than 307 GeV. Multilepton final states are sensitive to models with a variety of different leptonic or semilep- tonic RPV couplings. Limits on the mediator masses are established in a search in the four- lepton channel for strong production of neutralinos and their decay via the leptonic couplings λ121 and λ122. In a model of squark pair production with a neutralino LSP, lower limits on the squark mass are set at about 1.6 TeV. In models that feature electroweak production via leptonic couplings with wino- or Higgsino-like LSPs, limits are set on the LSP mass that range from 300 to 900 GeV. In those with semileptonic couplings, there are a variety of scenarios. In the regions of highest sensitivity, lower limits on the Higgsino masses reach approximately 700 GeV. Acknowledgments We congratulate our colleagues in the CERN accelerator departments for the excellent perfor- mance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we grate- fully acknowledge the computing centers and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses. Fi- nally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: the Austrian Federal Min- istry of Science, Research and Economy and the Austrian Science Fund; the Belgian Fonds de la Recherche Scientifique, and Fonds voor Wetenschappelijk Onderzoek; the Brazilian Fund- ing Agencies (CNPq, CAPES, FAPERJ, and FAPESP); the Bulgarian Ministry of Education and Science; CERN; the Chinese Academy of Sciences, Ministry of Science and Technology, and Na- tional Natural Science Foundation of China; the Colombian Funding Agency (COLCIENCIAS); the Croatian Ministry of Science, Education and Sport, and the Croatian Science Foundation; the Research Promotion Foundation, Cyprus; the Secretariat for Higher Education, Science, Technology and Innovation, Ecuador; the Ministry of Education and Research, Estonian Re- search Council via IUT23-4 and IUT23-6 and European Regional Development Fund, Estonia; the Academy of Finland, Finnish Ministry of Education and Culture, and Helsinki Institute of Physics; the Institut National de Physique Nucle´aire et de Physique des Particules / CNRS, and Commissariat a` l’E´nergie Atomique et aux E´nergies Alternatives / CEA, France; the Bundes- ministerium fu¨r Bildung und Forschung, Deutsche Forschungsgemeinschaft, and Helmholtz- References 41 Gemeinschaft Deutscher Forschungszentren, Germany; the General Secretariat for Research and Technology, Greece; the National Scientific Research Foundation, and National Innova- tion Office, Hungary; the Department of Atomic Energy and the Department of Science and Technology, India; the Institute for Studies in Theoretical Physics and Mathematics, Iran; the Science Foundation, Ireland; the Istituto Nazionale di Fisica Nucleare, Italy; the Ministry of Science, ICT and Future Planning, and National Research Foundation (NRF), Republic of Ko- rea; the Lithuanian Academy of Sciences; the Ministry of Education, and University of Malaya (Malaysia); the Mexican Funding Agencies (BUAP, CINVESTAV, CONACYT, LNS, SEP, and UASLP-FAI); the Ministry of Business, Innovation and Employment, New Zealand; the Pak- istan Atomic Energy Commission; the Ministry of Science and Higher Education and the Na- tional Science Centre, Poland; the Fundac¸a˜o para a Cieˆncia e a Tecnologia, Portugal; JINR, Dubna; the Ministry of Education and Science of the Russian Federation, the Federal Agency of Atomic Energy of the Russian Federation, Russian Academy of Sciences, and the Russian Foundation for Basic Research; the Ministry of Education, Science and Technological Devel- opment of Serbia; the Secretarı´a de Estado de Investigacio´n, Desarrollo e Innovacio´n and Pro- grama Consolider-Ingenio 2010, Spain; the Swiss Funding Agencies (ETH Board, ETH Zurich, PSI, SNF, UniZH, Canton Zurich, and SER); the Ministry of Science and Technology, Taipei; the Thailand Center of Excellence in Physics, the Institute for the Promotion of Teaching Science and Technology of Thailand, Special Task Force for Activating Research and the National Sci- ence and Technology Development Agency of Thailand; the Scientific and Technical Research Council of Turkey, and Turkish Atomic Energy Authority; the National Academy of Sciences of Ukraine, and State Fund for Fundamental Researches, Ukraine; the Science and Technology Fa- cilities Council, United Kingdom; the U.S. Department of Energy, and the U.S. National Science Foundation. Individuals have received support from the Marie-Curie program and the European Research Council and EPLANET (European Union); the Leventis Foundation; the A. P. Sloan Foun- dation; the Alexander von Humboldt Foundation; the Belgian Federal Science Policy Office; the Fonds pour la Formation a` la Recherche dans l’Industrie et dans l’Agriculture (FRIA- Belgium); the Agentschap voor Innovatie door Wetenschap en Technologie (IWT-Belgium); the Ministry of Education, Youth and Sports (MEYS) of the Czech Republic; the Council of Science and Industrial Research, India; the HOMING PLUS program of the Foundation for Polish Science, cofinanced from European Union, Regional Development Fund, the Mobility Plus program of the Ministry of Science and Higher Education, the OPUS program Contract No. 2014/13/B/ST2/02543 and Contract No. Sonata-bis DEC-2012/07/E/ST2/01406 of the National Science Center (Poland); the Thalis and Aristeia programs cofinanced by EU-ESF and the Greek NSRF; the National Priorities Research Program by Qatar National Research Fund; the Programa Cları´n-COFUND del Principado de Asturias; the Rachadapisek Sompot Fund for Postdoctoral Fellowship, Chulalongkorn University and the Chulalongkorn Academic into Its 2nd Century Project Advancement Project (Thailand); and the Welch Foundation, Contract No. C-1845. References [1] H. P. Nilles, “Supersymmetry, supergravity and particle physics”, Phys. Rep. 110 (1984) 1, doi:10.1016/0370-1573(84)90008-5. [2] H. E. Haber and G. L. Kane, “The search for supersymmetry: Probing physics beyond the standard model”, Phys. Rep. 117 (1985) 75, doi:10.1016/0370-1573(85)90051-1. 42 References [3] G. R. Farrar and P. Fayet, “Phenomenology of the production, decay, and detection of new hadronic states associated with supersymmetry”, Phys. Lett. B 76 (1978) 575, doi:10.1016/0370-2693(78)90858-4. [4] R. Barbier et al., “R-parity violating supersymmetry”, Phys. Rep. 420 (2005) 1, doi:10.1016/j.physrep.2005.08.006, arXiv:hep-ph/0406039. [5] Particle Data Group, J. Beringer et al., “Review of particle physics (RPP)”, Phys. Rev. D 86 (2012) 010001, doi:10.1103/PhysRevD.86.010001. [6] A. Y. Smirnov and F. Vissani, “Upper bound on all products of R-parity violating couplings λ′ and λ′′ from proton decay”, Phys. Lett. B 380 (1996) 317, doi:10.1016/0370-2693(96)00495-9, arXiv:hep-ph/9601387. [7] E. Nikolidakis and C. Smith, “Minimal flavor violation, seesaw mechanism, and R parity”, Phys. Rev. D 77 (2008) 015021, doi:10.1103/PhysRevD.77.015021, arXiv:0710.3129. [8] C. Csa´ki, Y. Grossman, and B. Heidenreich, “Minimal flavor violation supersymmetry: A natural theory for R-parity violation”, Phys. Rev. D 85 (2012) 095009, doi:10.1103/PhysRevD.85.095009, arXiv:1111.1239. [9] CDF Collaboration, “First Search for Multijet Resonances in √ s = 1.96 TeV pp¯ Collisions”, Phys. Rev. Lett. 107 (2011) 042001, doi:10.1103/PhysRevLett.107.042001, arXiv:1105.2815. [10] ATLAS Collaboration, “Search for pair production of massive particles decaying into three quarks with the ATLAS detector in √ s = 7 TeV pp collisions at the LHC”, JHEP 12 (2012) 086, doi:10.1007/JHEP12(2012)086, arXiv:1210.4813. [11] CMS Collaboration, “Search for three-jet resonances in pp collisions at √ s = 7 TeV”, Phys. Lett. B 718 (2012) 329, doi:10.1016/j.physletb.2012.10.048, arXiv:1208.2931. [12] CMS Collaboration, “Search for Three-Jet Resonances in pp Collisions at √ s = 7 TeV”, Phys. Rev. Lett. 107 (2011) 101801, doi:10.1103/PhysRevLett.107.101801, arXiv:1107.3084. [13] CMS Collaboration, “Searches for light- and heavy-flavour three-jet resonances in pp collisions at √ s = 8 TeV”, Phys. Lett. B 730 (2014) 193, doi:10.1016/j.physletb.2014.01.049, arXiv:1311.1799. [14] ATLAS Collaboration, “Search for massive supersymmetric particles decaying to many jets using the ATLAS detector in pp collisions at √ s = 8 TeV”, Phys. Rev. D 91 (2015) 112016, doi:10.1103/PhysRevD.91.112016, arXiv:1502.05686. [15] ALEPH Collaboration, “Search for supersymmetric particles with R-parity violating decays in e+e− collisions at √ s up to 209 GeV”, Eur. Phys. J. C 31 (2003) 1, doi:10.1140/epjc/s2003-01311-5, arXiv:hep-ex/0210014. [16] DELPHI Collaboration, “Search for supersymmetric particles assuming R-parity nonconservation in e+e− collisions at √ s = 192 GeV to 208 GeV”, Eur. Phys. J. C 36 (2004) 1, doi:10.1140/epjc/s2004-01881-6, arXiv:hep-ex/0406009. References 43 [17] L3 Collaboration, “Search for R parity violating decays of supersymmetric particles in e+e− collisions at LEP”, Phys. Lett. B 524 (2002) 65, doi:10.1016/S0370-2693(01)01367-3, arXiv:hep-ex/0110057. [18] D0 Collaboration, “Search for Resonant Second Generation Slepton Production at the Fermilab Tevatron”, Phys. Rev. Lett. 97 (2006) 111801, doi:10.1103/PhysRevLett.97.111801, arXiv:hep-ex/0605010. [19] D0 Collaboration, “Search for R-parity violating supersymmetry via the LLE couplings λ121, λ122 or λ133 in pp¯ collisions at √ s = 1.96 TeV”, Phys. Lett. B 638 (2006) 441, doi:10.1016/j.physletb.2006.05.077, arXiv:hep-ex/0605005. [20] CDF Collaboration, “Search for Anomalous Production of Multilepton Events in pp¯ Collisions at √ s = 1.96 TeV”, Phys. Rev. Lett. 98 (2007) 131804, doi:10.1103/PhysRevLett.98.131804, arXiv:0706.4448. [21] H1 Collaboration, “Search for squarks in R-parity violating supersymmetry in ep collisions at HERA”, Eur. Phys. J. C 71 (2011) 1572, doi:10.1140/epjc/s10052-011-1572-y, arXiv:1011.6359. [22] ZEUS Collaboration, “Search for stop production in R-parity-violating supersymmetry at HERA”, Eur. Phys. J. C 50 (2007) 269, doi:10.1140/epjc/s10052-007-0240-8, arXiv:hep-ex/0611018. [23] CMS Collaboration, “Search for physics beyond the standard model using multilepton signatures in pp collisions at √ s = 7 TeV”, Phys. Lett. B 704 (2011) 411, doi:10.1016/j.physletb.2011.09.047, arXiv:1106.0933. [24] CMS Collaboration, “Search for anomalous production of multilepton events in pp collisions at √ s = 7 TeV”, JHEP 06 (2012) 169, doi:10.1007/JHEP06(2012)169, arXiv:1204.5341. [25] ATLAS Collaboration, “Search for R-parity-violating supersymmetry in events with four or more leptons in √ s = 7 TeV pp collisions with the ATLAS detector”, JHEP 12 (2012) 124, doi:10.1007/JHEP12(2012)124, arXiv:1210.4457. [26] ATLAS Collaboration, “Search for a heavy narrow resonance decaying to eµ, eτ, or µτ with the ATLAS detector in √ s = 7 TeV pp collisions at the LHC”, Phys. Lett. B 723 (2013) 15, doi:10.1016/j.physletb.2013.04.035, arXiv:1212.1272. [27] CMS Collaboration, “Search for Top Squarks in R-Parity-Violating Supersymmetry Using Three or More Leptons and b-Tagged Jets”, Phys. Rev. Lett. 111 (2013) 221801, doi:10.1103/PhysRevLett.111.221801, arXiv:1306.6643. [28] ATLAS Collaboration, “Search for supersymmetry in events with four or more leptons in√ s = 8 TeV pp collisions with the ATLAS detector”, Phys. Rev. D 90 (2014) 052001, doi:10.1103/PhysRevD.90.052001, arXiv:1405.5086. [29] CMS Collaboration, “The CMS experiment at the CERN LHC”, JINST 3 (2008) S08004, doi:10.1088/1748-0221/3/08/S08004. [30] CMS Collaboration, “Study of tau reconstruction algorithms using pp collisions data collected at √ s = 7 TeV”, CMS Physics Analysis Summary CMS-PAS-PFT-10-004, 2010. 44 References [31] CMS Collaboration, “CMS Strategies for tau reconstruction and identification using particle-flow techniques”, CMS Physics Analysis Summary CMS-PAS-PFT-08-001, 2009. [32] CMS Collaboration, “Commissioning of the Particle-Flow Reconstruction in Minimum-Bias and Jet Events from pp Collisions at 7 TeV”, CMS Physics Analysis Summary CMS-PAS-PFT-10-002, 2010. [33] F. Maltoni and T. Stelzer, “MadEvent: Automatic event generation with MadGraph”, JHEP 02 (2003) 027, doi:10.1088/1126-6708/2003/02/027, arXiv:hep-ph/0208156. [34] T. Sjo¨strand, S. Mrenna, and P. Z. Skands, “A brief introduction to PYTHIA 8.1”, Comput. Phys. Commun. 178 (2008) 852, doi:10.1016/j.cpc.2008.01.036, arXiv:0710.3820. [35] GEANT4 Collaboration, “GEANT4—A simulation toolkit”, Nucl. Instrum. Meth. A 506 (2003) 250, doi:10.1016/S0168-9002(03)01368-8. [36] M. L. Mangano, M. Moretti, F. Piccinini, and M. Treccani, “Matching matrix elements and shower evolution for top-quark production in hadronic collisions”, JHEP 01 (2007) 013, doi:10.1088/1126-6708/2007/01/013, arXiv:hep-ph/0611129. [37] J. A. Evans and Y. Kats, “LHC coverage of RPV MSSM with light stops”, JHEP 04 (2013) 028, doi:10.1007/JHEP04(2013)028, arXiv:1209.0764. [38] CMS Collaboration, “The fast simulation of the CMS detector at LHC”, J. Phys. Conf. Ser. 331 (2011) 032049, doi:10.1088/1742-6596/331/3/032049. [39] J. Pumplin et al., “New generation of parton distributions with uncertainties from global QCD analysis”, JHEP 07 (2002) 012, doi:10.1088/1126-6708/2002/07/012, arXiv:hep-ph/0201195. [40] CMS Collaboration, “CMS Luminosity Based on Pixel Cluster Counting - Summer 2013 Update”, CMS Physics Analysis Summary CMS-PAS-LUM-13-001, 2013. [41] N. Arkani-Hamed et al., “MARMOSET: The path from LHC data to the new standard model via on-shell effective theories”, (2007). arXiv:hep-ph/0703088. [42] J. Alwall, P. Schuster, and N. Toro, “Simplified models for a first characterization of new physics at the LHC”, Phys. Rev. D 79 (2009) 075020, doi:10.1103/PhysRevD.79.075020, arXiv:0810.3921. [43] LHC New Physics Working Group Collaboration, “Simplified Models for LHC New Physics Searches”, J. Phys. G 39 (2012) 105005, doi:10.1088/0954-3899/39/10/105005, arXiv:1105.2838. [44] S. Alekhin et al., “The PDF4LHC Working Group interim report”, (2011). arXiv:1101.0536. [45] M. Botje et al., “The PDF4LHC Working Group interim recommendations”, (2011). arXiv:1101.0538. [46] P. M. Nadolsky et al., “Implications of CTEQ global analysis for collider observables”, Phys. Rev. D 78 (2008) 013004, doi:10.1103/PhysRevD.78.013004, arXiv:0802.0007. References 45 [47] A. D. Martin, W. J. Stirling, R. S. Thorne, and G. Watt, “Parton distributions for the LHC”, Eur. Phys. J. C 63 (2009) 189, doi:10.1140/epjc/s10052-009-1072-5, arXiv:0901.0002. [48] R. D. Ball et al., “A first unbiased global NLO determination of parton distributions and their uncertainties”, Nucl. Phys. B 838 (2010) 136, doi:10.1016/j.nuclphysb.2010.05.008, arXiv:1002.4407. [49] CMS Collaboration, “Search for top-squark pair production in the single-lepton final state in pp collisions at √ s = 8 TeV”, Eur. Phys. J. C 73 (2013) 2677, doi:10.1140/epjc/s10052-013-2677-2, arXiv:1308.1586. [50] T. Junk, “Confidence level computation for combining searches with small statistics”, Nucl. Instrum. Meth. A 434 (1999) 435, doi:10.1016/S0168-9002(99)00498-2, arXiv:hep-ex/9902006. [51] A. L. Read, “Presentation of search results: the CLs technique”, J. Phys. G 28 (2002) 2693, doi:10.1088/0954-3899/28/10/313. [52] ATLAS and CMS Collaborations, “Procedure for the LHC Higgs boson search combination in Summer 2011”, Technical Report CMS-NOTE-2011-005, ATL-PHYS-PUB-2011-11, 2011. [53] W. Beenakker et al., “Production of Charginos, Neutralinos, and Sleptons at Hadron Colliders”, Phys. Rev. Lett. 83 (1999) 3780, doi:10.1103/PhysRevLett.83.3780, arXiv:hep-ph/9906298. [54] A. Kulesza and L. Motyka, “Threshold Resummation for Squark-Antisquark and Gluino-Pair Production at the LHC”, Phys. Rev. Lett. 102 (2009) 111802, doi:10.1103/PhysRevLett.102.111802, arXiv:0807.2405. [55] A. Kulesza and L. Motyka, “Soft gluon resummation for the production of gluino-gluino and squark-antisquark pairs at the LHC”, Phys. Rev. D 80 (2009) 095004, doi:10.1103/PhysRevD.80.095004, arXiv:0905.4749. [56] W. Beenakker et al., “Soft-gluon resummation for squark and gluino hadroproduction”, JHEP 12 (2009) 041, doi:10.1088/1126-6708/2009/12/041, arXiv:0909.4418. [57] W. Beenakker et al., “Squark and gluino hadroproduction”, Int. J. Mod. Phys. A 26 (2011) 2637, doi:10.1142/S0217751X11053560, arXiv:1105.1110. [58] CMS Collaboration, “Performance of electron reconstruction and selection with the CMS detector in proton-proton collisions at √ s = 8 TeV”, JINST 10 (2015) P06005, doi:10.1088/1748-0221/10/06/P06005, arXiv:1502.02701. [59] CMS Collaboration, “Performance of CMS muon reconstruction in pp collision events at√ s = 7 TeV”, JINST 7 (2012) P10002, doi:10.1088/1748-0221/7/10/P10002, arXiv:1206.4071. [60] CMS Collaboration, “Performance of τ-lepton reconstruction and identification in CMS”, JINST 7 (2012) P01001, doi:10.1088/1748-0221/7/01/P01001, arXiv:1109.6034. [61] M. Cacciari, G. P. Salam, and G. Soyez, “The anti-kt jet clustering algorithm”, JHEP 04 (2008) 063, doi:10.1088/1126-6708/2008/04/063, arXiv:0802.1189. 46 References [62] CMS Collaboration, “Determination of jet energy calibration and transverse momentum resolution in CMS”, J. Instrum. 6 (2011) P11002, doi:10.1088/1748-0221/6/11/P11002. [63] CMS Collaboration, “Identification of b-quark jets with the CMS experiment”, JINST 8 (2013) P04013, doi:10.1088/1748-0221/8/04/P04013, arXiv:1211.4462. [64] CMS Collaboration, “Performance of b-tagging at √ s = 8 TeV in multijet, tt and boosted topology events”, CMS Physics Analysis Summary CMS-PAS-BTV-13-001, 2013. [65] CMS Collaboration, “Measurement of tt production with additional jet activity, including b quark jets, in the dilepton decay channel using pp collisions at √ s = 8 TeV”, Eur. Phys. J. C 76 (2016) 379, doi:10.1140/epjc/s10052-016-4105-x, arXiv:1510.03072. [66] R. Barlow and C. Beeston, “Fitting using finite Monte Carlo samples”, Comput. Phys. Commun. 77 (1993) 219, doi:10.1016/0010-4655(93)90005-W. [67] CMS Collaboration, “Measurement of BB angular correlations based on secondary vertex reconstruction at √ s=7 TeV”, J. High Energy Phys. 03 (2011) 136, doi:10.1007/JHEP03(2011)136. [68] M. Czakon, P. Fiedler, and A. Mitov, “Total Top-Quark Pair-Production Cross Section at Hadron Colliders Through O(α4S)”, Phys. Rev. Lett. 110 (2013) 252004, doi:10.1103/PhysRevLett.110.252004, arXiv:1303.6254. [69] CMS Collaboration, “Measurement of the differential cross section for top quark pair production in pp collisions at √ s = 8 TeV”, Eur. Phys. J. C 75 (2015) 542, doi:10.1140/epjc/s10052-015-3709-x, arXiv:1505.04480. [70] L. Sonnenschein, “Analytical solution of tt¯ dilepton equations”, Phys. Rev. D 73 (2006) 054015, doi:10.1103/PhysRevD.73.054015, arXiv:hep-ph/0603011. Erratum: [doi:10.1103/PhysRevD.78.079902]. [71] H.-C. Cheng et al., “Mass determination in SUSY-like events with missing energy”, JHEP 12 (2007) 076, doi:10.1088/1126-6708/2007/12/076, arXiv:0707.0030. [72] CMS Collaboration, “Measurement of the top-quark mass in tt events with dilepton final states in pp collisions at √ s = 7 TeV”, Eur. Phys. J. C 72 (2012) 2202, doi:10.1140/epjc/s10052-012-2202-z, arXiv:1209.2393. [73] CMS Collaboration, “Measurement of the tt production cross section in the eµ channel in proton-proton collisions at √ s = 7 and 8 TeV”, JHEP 08 (2016) 029, doi:10.1007/JHEP08(2016)029, arXiv:1603.02303. [74] G. J. Feldman and R. D. Cousins, “Unified approach to the classical statistical analysis of small signals”, Phys. Rev. D 57 (1998) 3873, doi:10.1103/PhysRevD.57.3873, arXiv:physics/9711021. [75] CMS Collaboration, “Measurement of W+W− and ZZ production cross sections in pp collisions at √ s = 8 TeV”, Phys. Lett. B 721 (2013) 190, doi:10.1016/j.physletb.2013.03.027, arXiv:1301.4698. [76] MSSM Working Group Collaboration, “The minimal supersymmetric standard model: Group summary report”, (1998). arXiv:hep-ph/9901246. References 47 [77] CMS Collaboration, “Phenomenological MSSM interpretation of CMS searches in pp collisions at √ s = 7 and 8 TeV”, JHEP 10 (2016) 129, doi:10.1007/JHEP10(2016)129, arXiv:1606.03577. 48 References 49 A The CMS Collaboration Yerevan Physics Institute, Yerevan, Armenia V. Khachatryan, A.M. Sirunyan, A. Tumasyan Institut fu¨r Hochenergiephysik der OeAW, Wien, Austria W. Adam, E. Asilar, T. Bergauer, J. Brandstetter, E. Brondolin, M. Dragicevic, J. Ero¨, M. Flechl, M. Friedl, R. Fru¨hwirth1, V.M. Ghete, C. Hartl, N. Ho¨rmann, J. Hrubec, M. Jeitler1, A. Ko¨nig, M. Krammer1, I. Kra¨tschmer, D. Liko, T. Matsushita, I. Mikulec, D. Rabady, N. Rad, B. Rahbaran, H. Rohringer, J. Schieck1, J. Strauss, W. Treberer-Treberspurg, W. Waltenberger, C.-E. Wulz1 National Centre for Particle and High Energy Physics, Minsk, Belarus V. Mossolov, N. Shumeiko, J. Suarez Gonzalez Universiteit Antwerpen, Antwerpen, Belgium S. Alderweireldt, T. Cornelis, E.A. De Wolf, X. Janssen, A. Knutsson, J. Lauwers, S. Luyckx, M. Van De Klundert, H. Van Haevermaet, P. Van Mechelen, N. Van Remortel, A. Van Spilbeeck Vrije Universiteit Brussel, Brussel, Belgium S. Abu Zeid, F. Blekman, J. D’Hondt, N. Daci, I. De Bruyn, K. Deroover, N. Heracleous, J. Keaveney, S. Lowette, S. Moortgat, L. Moreels, A. Olbrechts, Q. Python, D. Strom, S. Tavernier, W. Van Doninck, P. Van Mulders, I. Van Parijs Universite´ Libre de Bruxelles, Bruxelles, Belgium H. Brun, C. Caillol, B. Clerbaux, G. De Lentdecker, G. Fasanella, L. Favart, R. Goldouzian, A. Grebenyuk, G. Karapostoli, T. Lenzi, A. Le´onard, T. Maerschalk, A. Marinov, A. Randle- conde, T. Seva, C. Vander Velde, P. Vanlaer, R. Yonamine, F. Zenoni, F. Zhang2 Ghent University, Ghent, Belgium L. Benucci, A. Cimmino, S. Crucy, D. Dobur, A. Fagot, G. Garcia, M. Gul, J. Mccartin, A.A. Ocampo Rios, D. Poyraz, D. Ryckbosch, S. Salva, R. Scho¨fbeck, M. Sigamani, M. Tytgat, W. Van Driessche, E. Yazgan, N. Zaganidis Universite´ Catholique de Louvain, Louvain-la-Neuve, Belgium C. Beluffi3, O. Bondu, S. Brochet, G. Bruno, A. Caudron, L. Ceard, S. De Visscher, C. Delaere, M. Delcourt, L. Forthomme, B. Francois, A. Giammanco, A. Jafari, P. Jez, M. Komm, V. Lemaitre, A. Magitteri, A. Mertens, M. Musich, C. Nuttens, K. Piotrzkowski, L. Quertenmont, M. Selvaggi, M. Vidal Marono, S. Wertz Universite´ de Mons, Mons, Belgium N. Beliy, G.H. Hammad Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, Brazil W.L. Alda´ Ju´nior, F.L. Alves, G.A. Alves, L. Brito, M. Correa Martins Junior, M. Hamer, C. Hensel, A. Moraes, M.E. Pol, P. Rebello Teles Universidade do Estado do Rio de Janeiro, Rio de Janeiro, Brazil E. Belchior Batista Das Chagas, W. Carvalho, J. Chinellato4, A. Custo´dio, E.M. Da Costa, D. De Jesus Damiao, C. De Oliveira Martins, S. Fonseca De Souza, L.M. Huertas Guativa, H. Malbouisson, D. Matos Figueiredo, C. Mora Herrera, L. Mundim, H. Nogima, W.L. Prado Da Silva, A. Santoro, A. Sznajder, E.J. Tonelli Manganote4, A. Vilela Pereira Universidade Estadual Paulista a, Universidade Federal do ABC b, Sa˜o Paulo, Brazil S. Ahujaa, C.A. Bernardesb, A. De Souza Santosb, S. Dograa, T.R. Fernandez Perez Tomeia, 50 A The CMS Collaboration E.M. Gregoresb, P.G. Mercadanteb, C.S. Moona,5, S.F. Novaesa, Sandra S. Padulaa, D. Romero Abadb, J.C. Ruiz Vargas Institute for Nuclear Research and Nuclear Energy, Sofia, Bulgaria A. Aleksandrov, R. Hadjiiska, P. Iaydjiev, M. Rodozov, S. Stoykova, G. Sultanov, M. Vutova University of Sofia, Sofia, Bulgaria A. Dimitrov, I. Glushkov, L. Litov, B. Pavlov, P. Petkov Beihang University, Beijing, China W. Fang6 Institute of High Energy Physics, Beijing, China M. Ahmad, J.G. Bian, G.M. Chen, H.S. Chen, M. Chen, T. Cheng, R. Du, C.H. Jiang, D. Leggat, R. Plestina7, F. Romeo, S.M. Shaheen, A. Spiezia, J. Tao, C. Wang, Z. Wang, H. Zhang State Key Laboratory of Nuclear Physics and Technology, Peking University, Beijing, China C. Asawatangtrakuldee, Y. Ban, Q. Li, S. Liu, Y. Mao, S.J. Qian, D. Wang, Z. Xu Universidad de Los Andes, Bogota, Colombia C. Avila, A. Cabrera, L.F. Chaparro Sierra, C. Florez, J.P. Gomez, B. Gomez Moreno, J.C. Sanabria University of Split, Faculty of Electrical Engineering, Mechanical Engineering and Naval Architecture, Split, Croatia N. Godinovic, D. Lelas, I. Puljak, P.M. Ribeiro Cipriano University of Split, Faculty of Science, Split, Croatia Z. Antunovic, M. Kovac Institute Rudjer Boskovic, Zagreb, Croatia V. Brigljevic, D. Ferencek, K. Kadija, J. Luetic, S. Micanovic, L. Sudic University of Cyprus, Nicosia, Cyprus A. Attikis, G. Mavromanolakis, J. Mousa, C. Nicolaou, F. Ptochos, P.A. Razis, H. Rykaczewski Charles University, Prague, Czech Republic M. Finger8, M. Finger Jr.8 Universidad San Francisco de Quito, Quito, Ecuador E. Carrera Jarrin Academy of Scientific Research and Technology of the Arab Republic of Egypt, Egyptian Network of High Energy Physics, Cairo, Egypt A.A. Abdelalim9,10, E. El-khateeb11, T. Elkafrawy11, M.A. Mahmoud12,13 National Institute of Chemical Physics and Biophysics, Tallinn, Estonia B. Calpas, M. Kadastik, M. Murumaa, L. Perrini, M. Raidal, A. Tiko, C. Veelken Department of Physics, University of Helsinki, Helsinki, Finland P. Eerola, J. Pekkanen, M. Voutilainen Helsinki Institute of Physics, Helsinki, Finland J. Ha¨rko¨nen, V. Karima¨ki, R. Kinnunen, T. Lampe´n, K. Lassila-Perini, S. Lehti, T. Linde´n, P. Luukka, T. Peltola, J. Tuominiemi, E. Tuovinen, L. Wendland Lappeenranta University of Technology, Lappeenranta, Finland J. Talvitie, T. Tuuva 51 DSM/IRFU, CEA/Saclay, Gif-sur-Yvette, France M. Besancon, F. Couderc, M. Dejardin, D. Denegri, B. Fabbro, J.L. Faure, C. Favaro, F. Ferri, S. Ganjour, A. Givernaud, P. Gras, G. Hamel de Monchenault, P. Jarry, E. Locci, M. Machet, J. Malcles, J. Rander, A. Rosowsky, M. Titov, A. Zghiche Laboratoire Leprince-Ringuet, Ecole Polytechnique, IN2P3-CNRS, Palaiseau, France A. Abdulsalam, I. Antropov, S. Baffioni, F. Beaudette, P. Busson, L. Cadamuro, E. Chapon, C. Charlot, O. Davignon, L. Dobrzynski, R. Granier de Cassagnac, M. Jo, S. Lisniak, P. Mine´, I.N. Naranjo, M. Nguyen, C. Ochando, G. Ortona, P. Paganini, P. Pigard, S. Regnard, R. Salerno, Y. Sirois, T. Strebler, Y. Yilmaz, A. Zabi Institut Pluridisciplinaire Hubert Curien, Universite´ de Strasbourg, Universite´ de Haute Alsace Mulhouse, CNRS/IN2P3, Strasbourg, France J.-L. Agram14, J. Andrea, A. Aubin, D. Bloch, J.-M. Brom, M. Buttignol, E.C. Chabert, N. Chanon, C. Collard, E. Conte14, X. Coubez, J.-C. Fontaine14, D. Gele´, U. Goerlach, C. Goetzmann, A.-C. Le Bihan, J.A. Merlin15, K. Skovpen, P. Van Hove Centre de Calcul de l’Institut National de Physique Nucleaire et de Physique des Particules, CNRS/IN2P3, Villeurbanne, France S. Gadrat Universite´ de Lyon, Universite´ Claude Bernard Lyon 1, CNRS-IN2P3, Institut de Physique Nucle´aire de Lyon, Villeurbanne, France S. Beauceron, C. Bernet, G. Boudoul, E. Bouvier, C.A. Carrillo Montoya, R. Chierici, D. Contardo, B. Courbon, P. Depasse, H. El Mamouni, J. Fan, J. Fay, S. Gascon, M. Gouzevitch, B. Ille, F. Lagarde, I.B. Laktineh, M. Lethuillier, L. Mirabito, A.L. Pequegnot, S. Perries, A. Popov16, J.D. Ruiz Alvarez, D. Sabes, V. Sordini, M. Vander Donckt, P. Verdier, S. Viret Georgian Technical University, Tbilisi, Georgia T. Toriashvili17 Tbilisi State University, Tbilisi, Georgia Z. Tsamalaidze8 RWTH Aachen University, I. Physikalisches Institut, Aachen, Germany C. Autermann, S. Beranek, L. Feld, A. Heister, M.K. Kiesel, K. Klein, M. Lipinski, A. Ostapchuk, M. Preuten, F. Raupach, S. Schael, C. Schomakers, J.F. Schulte, J. Schulz, T. Verlage, H. Weber, V. Zhukov16 RWTH Aachen University, III. Physikalisches Institut A, Aachen, Germany M. Ata, M. Brodski, E. Dietz-Laursonn, D. Duchardt, M. Endres, M. Erdmann, S. Erdweg, T. Esch, R. Fischer, A. Gu¨th, T. Hebbeker, C. Heidemann, K. Hoepfner, S. Knutzen, M. Merschmeyer, A. Meyer, P. Millet, S. Mukherjee, M. Olschewski, K. Padeken, P. Papacz, T. Pook, M. Radziej, H. Reithler, M. Rieger, F. Scheuch, L. Sonnenschein, D. Teyssier, S. Thu¨er RWTH Aachen University, III. Physikalisches Institut B, Aachen, Germany V. Cherepanov, Y. Erdogan, G. Flu¨gge, H. Geenen, M. Geisler, F. Hoehle, B. Kargoll, T. Kress, A. Ku¨nsken, J. Lingemann, A. Nehrkorn, A. Nowack, I.M. Nugent, C. Pistone, O. Pooth, A. Stahl15 Deutsches Elektronen-Synchrotron, Hamburg, Germany M. Aldaya Martin, I. Asin, K. Beernaert, O. Behnke, U. Behrens, K. Borras18, A. Campbell, P. Connor, C. Contreras-Campana, F. Costanza, C. Diez Pardos, G. Dolinska, S. Dooling, G. Eckerlin, D. Eckstein, T. Eichhorn, E. Gallo19, J. Garay Garcia, A. Geiser, A. Gizhko, 52 A The CMS Collaboration J.M. Grados Luyando, P. Gunnellini, A. Harb, J. Hauk, M. Hempel20, H. Jung, A. Kalogeropou- los, O. Karacheban20, M. Kasemann, J. Kieseler, C. Kleinwort, I. Korol, W. Lange, A. Lelek, J. Leonard, K. Lipka, A. Lobanov, W. Lohmann20, R. Mankel, I.-A. Melzer-Pellmann, A.B. Meyer, G. Mittag, J. Mnich, A. Mussgiller, E. Ntomari, D. Pitzl, R. Placakyte, A. Raspereza, B. Roland, M.O¨. Sahin, P. Saxena, T. Schoerner-Sadenius, C. Seitz, S. Spannagel, N. Stefaniuk, K.D. Trippkewitz, G.P. Van Onsem, R. Walsh, C. Wissing University of Hamburg, Hamburg, Germany V. Blobel, M. Centis Vignali, A.R. Draeger, T. Dreyer, J. Erfle, E. Garutti, K. Goebel, D. Gonzalez, M. Go¨rner, J. Haller, M. Hoffmann, R.S. Ho¨ing, A. Junkes, R. Klanner, R. Kogler, N. Kovalchuk, T. Lapsien, T. Lenz, I. Marchesini, D. Marconi, M. Meyer, M. Niedziela, D. Nowatschin, J. Ott, F. Pantaleo15, T. Peiffer, A. Perieanu, N. Pietsch, J. Poehlsen, C. Sander, C. Scharf, P. Schleper, E. Schlieckau, A. Schmidt, S. Schumann, J. Schwandt, H. Stadie, G. Steinbru¨ck, F.M. Stober, H. Tholen, D. Troendle, E. Usai, L. Vanelderen, A. Vanhoefer, B. Vormwald Institut fu¨r Experimentelle Kernphysik, Karlsruhe, Germany C. Barth, C. Baus, J. Berger, C. Bo¨ser, E. Butz, T. Chwalek, F. Colombo, W. De Boer, A. Descroix, A. Dierlamm, S. Fink, F. Frensch, R. Friese, M. Giffels, A. Gilbert, D. Haitz, F. Hartmann15, S.M. Heindl, U. Husemann, I. Katkov16, A. Kornmayer15, P. Lobelle Pardo, B. Maier, H. Mildner, M.U. Mozer, T. Mu¨ller, Th. Mu¨ller, M. Plagge, G. Quast, K. Rabbertz, S. Ro¨cker, F. Roscher, M. Schro¨der, G. Sieber, H.J. Simonis, R. Ulrich, J. Wagner-Kuhr, S. Wayand, M. Weber, T. Weiler, S. Williamson, C. Wo¨hrmann, R. Wolf Institute of Nuclear and Particle Physics (INPP), NCSR Demokritos, Aghia Paraskevi, Greece G. Anagnostou, G. Daskalakis, T. Geralis, V.A. Giakoumopoulou, A. Kyriakis, D. Loukas, A. Psallidas, I. Topsis-Giotis National and Kapodistrian University of Athens, Athens, Greece A. Agapitos, S. Kesisoglou, A. Panagiotou, N. Saoulidou, E. Tziaferi University of Ioa´nnina, Ioa´nnina, Greece I. Evangelou, G. Flouris, C. Foudas, P. Kokkas, N. Loukas, N. Manthos, I. Papadopoulos, E. Paradas, J. Strologas MTA-ELTE Lendu¨let CMS Particle and Nuclear Physics Group, Eo¨tvo¨s Lora´nd University N. Filipovic Wigner Research Centre for Physics, Budapest, Hungary G. Bencze, C. Hajdu, P. Hidas, D. Horvath21, F. Sikler, V. Veszpremi, G. Vesztergombi22, A.J. Zsigmond Institute of Nuclear Research ATOMKI, Debrecen, Hungary N. Beni, S. Czellar, J. Karancsi23, J. Molnar, Z. Szillasi University of Debrecen, Debrecen, Hungary M. Barto´k22, A. Makovec, P. Raics, Z.L. Trocsanyi, B. Ujvari National Institute of Science Education and Research, Bhubaneswar, India S. Choudhury24, P. Mal, K. Mandal, A. Nayak, D.K. Sahoo, N. Sahoo, S.K. Swain Panjab University, Chandigarh, India S. Bansal, S.B. Beri, V. Bhatnagar, R. Chawla, R. Gupta, U.Bhawandeep, A.K. Kalsi, A. Kaur, M. Kaur, R. Kumar, A. Mehta, M. Mittal, J.B. Singh, G. Walia 53 University of Delhi, Delhi, India Ashok Kumar, A. Bhardwaj, B.C. Choudhary, R.B. Garg, S. Keshri, A. Kumar, S. Malhotra, M. Naimuddin, N. Nishu, K. Ranjan, R. Sharma, V. Sharma Saha Institute of Nuclear Physics, Kolkata, India R. Bhattacharya, S. Bhattacharya, K. Chatterjee, S. Dey, S. Dutta, S. Ghosh, N. Majumdar, A. Modak, K. Mondal, S. Mukhopadhyay, S. Nandan, A. Purohit, A. Roy, D. Roy, S. Roy Chowdhury, S. Sarkar, M. Sharan Bhabha Atomic Research Centre, Mumbai, India R. Chudasama, D. Dutta, V. Jha, V. Kumar, A.K. Mohanty15, L.M. Pant, P. Shukla, A. Topkar Tata Institute of Fundamental Research, Mumbai, India T. Aziz, S. Banerjee, S. Bhowmik25, R.M. Chatterjee, R.K. Dewanjee, S. Dugad, S. Ganguly, S. Ghosh, M. Guchait, A. Gurtu26, Sa. Jain, G. Kole, S. Kumar, B. Mahakud, M. Maity25, G. Majumder, K. Mazumdar, S. Mitra, G.B. Mohanty, B. Parida, T. Sarkar25, N. Sur, B. Sutar, N. Wickramage27 Indian Institute of Science Education and Research (IISER), Pune, India S. Chauhan, S. Dube, A. Kapoor, K. Kothekar, A. Rane, S. Sharma Institute for Research in Fundamental Sciences (IPM), Tehran, Iran H. Bakhshiansohi, H. Behnamian, S.M. Etesami28, A. Fahim29, M. Khakzad, M. Mohammadi Najafabadi, M. Naseri, S. Paktinat Mehdiabadi, F. Rezaei Hosseinabadi, B. Safarzadeh30, M. Zeinali University College Dublin, Dublin, Ireland M. Felcini, M. Grunewald INFN Sezione di Bari a, Universita` di Bari b, Politecnico di Bari c, Bari, Italy M. Abbresciaa ,b, C. Calabriaa,b, C. Caputoa ,b, A. Colaleoa, D. Creanzaa ,c, L. Cristellaa,b, N. De Filippisa ,c, M. De Palmaa,b, L. Fiorea, G. Iasellia ,c, G. Maggia,c, M. Maggia, G. Minielloa ,b, S. Mya ,b, S. Nuzzoa,b, A. Pompilia ,b, G. Pugliesea,c, R. Radognaa ,b, A. Ranieria, G. Selvaggia ,b, L. Silvestrisa,15, R. Vendittia,b INFN Sezione di Bologna a, Universita` di Bologna b, Bologna, Italy G. Abbiendia, C. Battilana, D. Bonacorsia ,b, S. Braibant-Giacomellia,b, L. Brigliadoria ,b, R. Campaninia ,b, P. Capiluppia,b, A. Castroa ,b, F.R. Cavalloa, S.S. Chhibraa,b, G. Codispotia ,b, M. Cuffiania ,b, G.M. Dallavallea, F. Fabbria, A. Fanfania,b, D. Fasanellaa,b, P. Giacomellia, C. Grandia, L. Guiduccia ,b, S. Marcellinia, G. Masettia, A. Montanaria, F.L. Navarriaa ,b, A. Perrottaa, A.M. Rossia,b, T. Rovellia ,b, G.P. Sirolia ,b, N. Tosia ,b ,15 INFN Sezione di Catania a, Universita` di Catania b, Catania, Italy G. Cappellob, M. Chiorbolia,b, S. Costaa ,b, A. Di Mattiaa, F. Giordanoa ,b, R. Potenzaa ,b, A. Tricomia,b, C. Tuvea ,b INFN Sezione di Firenze a, Universita` di Firenze b, Firenze, Italy G. Barbaglia, V. Ciullia,b, C. Civininia, R. D’Alessandroa,b, E. Focardia,b, V. Goria ,b, P. Lenzia ,b, M. Meschinia, S. Paolettia, G. Sguazzonia, L. Viliania ,b ,15 INFN Laboratori Nazionali di Frascati, Frascati, Italy L. Benussi, S. Bianco, F. Fabbri, D. Piccolo, F. Primavera15 INFN Sezione di Genova a, Universita` di Genova b, Genova, Italy V. Calvellia ,b, F. Ferroa, M. Lo Veterea,b, M.R. Mongea ,b, E. Robuttia, S. Tosia ,b 54 A The CMS Collaboration INFN Sezione di Milano-Bicocca a, Universita` di Milano-Bicocca b, Milano, Italy L. Brianza, M.E. Dinardoa,b, S. Fiorendia ,b, S. Gennaia, A. Ghezzia,b, P. Govonia,b, S. Malvezzia, R.A. Manzonia,b,15, B. Marzocchia ,b, D. Menascea, L. Moronia, M. Paganonia ,b, D. Pedrinia, S. Pigazzini, S. Ragazzia,b, N. Redaellia, T. Tabarelli de Fatisa,b INFN Sezione di Napoli a, Universita` di Napoli ’Federico II’ b, Napoli, Italy, Universita` della Basilicata c, Potenza, Italy, Universita` G. Marconi d, Roma, Italy S. Buontempoa, N. Cavalloa ,c, S. Di Guidaa,d ,15, M. Espositoa ,b, F. Fabozzia ,c, A.O.M. Iorioa ,b, G. Lanzaa, L. Listaa, S. Meolaa,d ,15, M. Merolaa, P. Paoluccia ,15, C. Sciaccaa,b, F. Thyssen INFN Sezione di Padova a, Universita` di Padova b, Padova, Italy, Universita` di Trento c, Trento, Italy P. Azzia ,15, N. Bacchettaa, L. Benatoa,b, D. Biselloa,b, A. Bolettia ,b, R. Carlina,b, P. Checchiaa, M. Dall’Ossoa ,b, P. De Castro Manzanoa, T. Dorigoa, U. Dossellia, F. Gasparinia ,b, U. Gasparinia ,b, A. Gozzelinoa, S. Lacapraraa, M. Margonia ,b, A.T. Meneguzzoa ,b, F. Montecassianoa, M. Passaseoa, J. Pazzinia,b ,15, M. Pegoraroa, N. Pozzobona,b, P. Ronchesea ,b, F. Simonettoa ,b, E. Torassaa, M. Tosia,b, S. Vaninia,b, M. Zanetti, P. Zottoa ,b, A. Zucchettaa,b INFN Sezione di Pavia a, Universita` di Pavia b, Pavia, Italy A. Braghieria, A. Magnania ,b, P. Montagnaa,b, S.P. Rattia ,b, V. Rea, C. Riccardia ,b, P. Salvinia, I. Vaia,b, P. Vituloa ,b INFN Sezione di Perugia a, Universita` di Perugia b, Perugia, Italy L. Alunni Solestizia,b, G.M. Bileia, D. Ciangottinia,b, L. Fano`a,b, P. Laricciaa ,b, R. Leonardia ,b, G. Mantovania,b, M. Menichellia, A. Sahaa, A. Santocchiaa,b INFN Sezione di Pisa a, Universita` di Pisa b, Scuola Normale Superiore di Pisa c, Pisa, Italy K. Androsova ,31, P. Azzurria,15, G. Bagliesia, J. Bernardinia, T. Boccalia, R. Castaldia, M.A. Cioccia,31, R. Dell’Orsoa, S. Donatoa ,c, G. Fedi, A. Giassia, M.T. Grippoa,31, F. Ligabuea ,c, T. Lomtadzea, L. Martinia,b, A. Messineoa,b, F. Pallaa, A. Rizzia,b, A. Savoy-Navarroa ,32, P. Spagnoloa, R. Tenchinia, G. Tonellia,b, A. Venturia, P.G. Verdinia INFN Sezione di Roma a, Universita` di Roma b, Roma, Italy L. Baronea ,b, F. Cavallaria, G. D’imperioa,b,15, D. Del Rea,b,15, M. Diemoza, S. Gellia,b, C. Jordaa, E. Longoa,b, F. Margarolia,b, P. Meridiania, G. Organtinia ,b, R. Paramattia, F. Preiatoa ,b, S. Rahatloua,b, C. Rovellia, F. Santanastasioa,b INFN Sezione di Torino a, Universita` di Torino b, Torino, Italy, Universita` del Piemonte Orientale c, Novara, Italy N. Amapanea ,b, R. Arcidiaconoa ,c,15, S. Argiroa ,b, M. Arneodoa ,c, N. Bartosika, R. Bellana ,b, C. Biinoa, N. Cartigliaa, M. Costaa,b, R. Covarellia ,b, A. Deganoa,b, N. Demariaa, L. Fincoa ,b, B. Kiania ,b, C. Mariottia, S. Masellia, E. Migliorea,b, V. Monacoa,b, E. Monteila ,b, M.M. Obertinoa,b, L. Pachera ,b, N. Pastronea, M. Pelliccionia, G.L. Pinna Angionia,b, F. Raveraa ,b, A. Romeroa ,b, M. Ruspaa ,c, R. Sacchia ,b, V. Solaa, A. Solanoa,b, A. Staianoa, P. Traczyka ,b INFN Sezione di Trieste a, Universita` di Trieste b, Trieste, Italy S. Belfortea, V. Candelisea ,b, M. Casarsaa, F. Cossuttia, G. Della Riccaa ,b, C. La Licataa ,b, A. Schizzia ,b, A. Zanettia Kangwon National University, Chunchon, Korea S.K. Nam Kyungpook National University, Daegu, Korea 55 D.H. Kim, G.N. Kim, M.S. Kim, D.J. Kong, S. Lee, S.W. Lee, Y.D. Oh, A. Sakharov, D.C. Son, Y.C. Yang Chonbuk National University, Jeonju, Korea J.A. Brochero Cifuentes, H. Kim, T.J. Kim33 Chonnam National University, Institute for Universe and Elementary Particles, Kwangju, Korea S. Song Korea University, Seoul, Korea S. Cho, S. Choi, Y. Go, D. Gyun, B. Hong, Y. Jo, Y. Kim, B. Lee, K. Lee, K.S. Lee, S. Lee, J. Lim, S.K. Park, Y. Roh Seoul National University, Seoul, Korea H.D. Yoo University of Seoul, Seoul, Korea M. Choi, H. Kim, H. Kim, J.H. Kim, J.S.H. Lee, I.C. Park, G. Ryu, M.S. Ryu Sungkyunkwan University, Suwon, Korea Y. Choi, J. Goh, D. Kim, E. Kwon, J. Lee, I. Yu Vilnius University, Vilnius, Lithuania V. Dudenas, A. Juodagalvis, J. Vaitkus National Centre for Particle Physics, Universiti Malaya, Kuala Lumpur, Malaysia I. Ahmed, Z.A. Ibrahim, J.R. Komaragiri, M.A.B. Md Ali34, F. Mohamad Idris35, W.A.T. Wan Abdullah, M.N. Yusli, Z. Zolkapli Centro de Investigacion y de Estudios Avanzados del IPN, Mexico City, Mexico E. Casimiro Linares, H. Castilla-Valdez, E. De La Cruz-Burelo, I. Heredia-De La Cruz36, A. Hernandez-Almada, R. Lopez-Fernandez, J. Mejia Guisao, A. Sanchez-Hernandez Universidad Iberoamericana, Mexico City, Mexico S. Carrillo Moreno, F. Vazquez Valencia Benemerita Universidad Autonoma de Puebla, Puebla, Mexico I. Pedraza, H.A. Salazar Ibarguen, C. Uribe Estrada Universidad Auto´noma de San Luis Potosı´, San Luis Potosı´, Mexico A. Morelos Pineda University of Auckland, Auckland, New Zealand D. Krofcheck University of Canterbury, Christchurch, New Zealand P.H. Butler National Centre for Physics, Quaid-I-Azam University, Islamabad, Pakistan A. Ahmad, M. Ahmad, Q. Hassan, H.R. Hoorani, W.A. Khan, S. Qazi, M. Shoaib, M. Waqas National Centre for Nuclear Research, Swierk, Poland H. Bialkowska, M. Bluj, B. Boimska, T. Frueboes, M. Go´rski, M. Kazana, K. Nawrocki, K. Romanowska-Rybinska, M. Szleper, P. Zalewski 56 A The CMS Collaboration Institute of Experimental Physics, Faculty of Physics, University of Warsaw, Warsaw, Poland G. Brona, K. Bunkowski, A. Byszuk37, K. Doroba, A. Kalinowski, M. Konecki, J. Krolikowski, M. Misiura, M. Olszewski, M. Walczak Laborato´rio de Instrumentac¸a˜o e Fı´sica Experimental de Partı´culas, Lisboa, Portugal P. Bargassa, C. Beira˜o Da Cruz E Silva, A. Di Francesco, P. Faccioli, P.G. Ferreira Parracho, M. Gallinaro, J. Hollar, N. Leonardo, L. Lloret Iglesias, M.V. Nemallapudi, F. Nguyen, J. Rodrigues Antunes, J. Seixas, O. Toldaiev, D. Vadruccio, J. Varela, P. Vischia Joint Institute for Nuclear Research, Dubna, Russia P. Bunin, M. Gavrilenko, I. Golutvin, I. Gorbunov, A. Kamenev, V. Karjavin, A. Lanev, A. Malakhov, V. Matveev38,39, P. Moisenz, V. Palichik, V. Perelygin, M. Savina, S. Shmatov, S. Shulha, N. Skatchkov, V. Smirnov, N. Voytishin, A. Zarubin Petersburg Nuclear Physics Institute, Gatchina (St. Petersburg), Russia V. Golovtsov, Y. Ivanov, V. Kim40, E. Kuznetsova41, P. Levchenko, V. Murzin, V. Oreshkin, I. Smirnov, V. Sulimov, L. Uvarov, S. Vavilov, A. Vorobyev Institute for Nuclear Research, Moscow, Russia Yu. Andreev, A. Dermenev, S. Gninenko, N. Golubev, A. Karneyeu, M. Kirsanov, N. Krasnikov, A. Pashenkov, D. Tlisov, A. Toropin Institute for Theoretical and Experimental Physics, Moscow, Russia V. Epshteyn, V. Gavrilov, N. Lychkovskaya, V. Popov, I. Pozdnyakov, G. Safronov, A. Spiridonov, M. Toms, E. Vlasov, A. Zhokin National Research Nuclear University ’Moscow Engineering Physics Institute’ (MEPhI), Moscow, Russia M. Danilov, O. Markin, E. Popova, V. Rusinov, E. Tarkovskii P.N. Lebedev Physical Institute, Moscow, Russia V. Andreev, M. Azarkin39, I. Dremin39, M. Kirakosyan, A. Leonidov39, G. Mesyats, S.V. Rusakov Skobeltsyn Institute of Nuclear Physics, Lomonosov Moscow State University, Moscow, Russia A. Baskakov, A. Belyaev, E. Boos, M. Dubinin42, L. Dudko, A. Ershov, A. Gribushin, V. Klyukhin, O. Kodolova, I. Lokhtin, I. Miagkov, S. Obraztsov, S. Petrushanko, V. Savrin, A. Snigirev State Research Center of Russian Federation, Institute for High Energy Physics, Protvino, Russia I. Azhgirey, I. Bayshev, S. Bitioukov, V. Kachanov, A. Kalinin, D. Konstantinov, V. Krychkine, V. Petrov, R. Ryutin, A. Sobol, L. Tourtchanovitch, S. Troshin, N. Tyurin, A. Uzunian, A. Volkov University of Belgrade, Faculty of Physics and Vinca Institute of Nuclear Sciences, Belgrade, Serbia P. Adzic43, P. Cirkovic, D. Devetak, J. Milosevic, V. Rekovic Centro de Investigaciones Energe´ticas Medioambientales y Tecnolo´gicas (CIEMAT), Madrid, Spain J. Alcaraz Maestre, E. Calvo, M. Cerrada, M. Chamizo Llatas, N. Colino, B. De La Cruz, A. Delgado Peris, A. Escalante Del Valle, C. Fernandez Bedoya, J.P. Ferna´ndez Ramos, J. Flix, M.C. Fouz, P. Garcia-Abia, O. Gonzalez Lopez, S. Goy Lopez, J.M. Hernandez, M.I. Josa, E. Navarro De Martino, A. Pe´rez-Calero Yzquierdo, J. Puerta Pelayo, A. Quintario Olmeda, I. Redondo, L. Romero, M.S. Soares 57 Universidad Auto´noma de Madrid, Madrid, Spain J.F. de Troco´niz, M. Missiroli, D. Moran Universidad de Oviedo, Oviedo, Spain J. Cuevas, J. Fernandez Menendez, S. Folgueras, I. Gonzalez Caballero, E. Palencia Cortezon, J.M. Vizan Garcia Instituto de Fı´sica de Cantabria (IFCA), CSIC-Universidad de Cantabria, Santander, Spain I.J. Cabrillo, A. Calderon, J.R. Castin˜eiras De Saa, E. Curras, M. Fernandez, J. Garcia-Ferrero, G. Gomez, A. Lopez Virto, J. Marco, R. Marco, C. Martinez Rivero, F. Matorras, J. Piedra Gomez, T. Rodrigo, A.Y. Rodrı´guez-Marrero, A. Ruiz-Jimeno, L. Scodellaro, N. Trevisani, I. Vila, R. Vilar Cortabitarte CERN, European Organization for Nuclear Research, Geneva, Switzerland D. Abbaneo, E. Auffray, G. Auzinger, M. Bachtis, P. Baillon, A.H. Ball, D. Barney, A. Benaglia, L. Benhabib, G.M. Berruti, P. Bloch, A. Bocci, A. Bonato, C. Botta, H. Breuker, T. Camporesi, R. Castello, M. Cepeda, G. Cerminara, M. D’Alfonso, D. d’Enterria, A. Dabrowski, V. Daponte, A. David, M. De Gruttola, F. De Guio, A. De Roeck, E. Di Marco44, M. Dobson, M. Dordevic, B. Dorney, T. du Pree, D. Duggan, M. Du¨nser, N. Dupont, A. Elliott-Peisert, S. Fartoukh, G. Franzoni, J. Fulcher, W. Funk, D. Gigi, K. Gill, M. Girone, F. Glege, R. Guida, S. Gundacker, M. Guthoff, J. Hammer, P. Harris, J. Hegeman, V. Innocente, P. Janot, H. Kirschenmann, V. Knu¨nz, M.J. Kortelainen, K. Kousouris, P. Lecoq, C. Lourenc¸o, M.T. Lucchini, N. Magini, L. Malgeri, M. Mannelli, A. Martelli, L. Masetti, F. Meijers, S. Mersi, E. Meschi, F. Moortgat, S. Morovic, M. Mulders, H. Neugebauer, S. Orfanelli45, L. Orsini, L. Pape, E. Perez, M. Peruzzi, A. Petrilli, G. Petrucciani, A. Pfeiffer, M. Pierini, D. Piparo, A. Racz, T. Reis, G. Rolandi46, M. Rovere, M. Ruan, H. Sakulin, J.B. Sauvan, C. Scha¨fer, C. Schwick, M. Seidel, A. Sharma, P. Silva, M. Simon, P. Sphicas47, J. Steggemann, M. Stoye, Y. Takahashi, D. Treille, A. Triossi, A. Tsirou, V. Veckalns48, G.I. Veres22, N. Wardle, H.K. Wo¨hri, A. Zagozdzinska37, W.D. Zeuner Paul Scherrer Institut, Villigen, Switzerland W. Bertl, K. Deiters, W. Erdmann, R. Horisberger, Q. Ingram, H.C. Kaestli, D. Kotlinski, U. Langenegger, T. Rohe Institute for Particle Physics, ETH Zurich, Zurich, Switzerland F. Bachmair, L. Ba¨ni, L. Bianchini, B. Casal, G. Dissertori, M. Dittmar, M. Donega`, P. Eller, C. Grab, C. Heidegger, D. Hits, J. Hoss, G. Kasieczka, P. Lecomte†, W. Lustermann, B. Mangano, M. Marionneau, P. Martinez Ruiz del Arbol, M. Masciovecchio, M.T. Meinhard, D. Meister, F. Micheli, P. Musella, F. Nessi-Tedaldi, F. Pandolfi, J. Pata, F. Pauss, G. Perrin, L. Perrozzi, M. Quittnat, M. Rossini, M. Scho¨nenberger, A. Starodumov49, M. Takahashi, V.R. Tavolaro, K. Theofilatos, R. Wallny Universita¨t Zu¨rich, Zurich, Switzerland T.K. Aarrestad, C. Amsler50, L. Caminada, M.F. Canelli, V. Chiochia, A. De Cosa, C. Galloni, A. Hinzmann, T. Hreus, B. Kilminster, C. Lange, J. Ngadiuba, D. Pinna, G. Rauco, P. Robmann, D. Salerno, Y. Yang National Central University, Chung-Li, Taiwan K.H. Chen, T.H. Doan, Sh. Jain, R. Khurana, M. Konyushikhin, C.M. Kuo, W. Lin, Y.J. Lu, A. Pozdnyakov, S.S. Yu National Taiwan University (NTU), Taipei, Taiwan Arun Kumar, P. Chang, Y.H. Chang, Y.W. Chang, Y. Chao, K.F. Chen, P.H. Chen, C. Dietz, F. Fiori, W.-S. Hou, Y. Hsiung, Y.F. Liu, R.-S. Lu, M. Min˜ano Moya, J.f. Tsai, Y.M. Tzeng 58 A The CMS Collaboration Chulalongkorn University, Faculty of Science, Department of Physics, Bangkok, Thailand B. Asavapibhop, K. Kovitanggoon, G. Singh, N. Srimanobhas, N. Suwonjandee Cukurova University, Adana, Turkey A. Adiguzel, S. Cerci51, S. Damarseckin, Z.S. Demiroglu, C. Dozen, I. Dumanoglu, S. Girgis, G. Gokbulut, Y. Guler, E. Gurpinar, I. Hos, E.E. Kangal52, A. Kayis Topaksu, G. Onengut53, K. Ozdemir54, S. Ozturk55, B. Tali51, H. Topakli55, C. Zorbilmez Middle East Technical University, Physics Department, Ankara, Turkey B. Bilin, S. Bilmis, B. Isildak56, G. Karapinar57, M. Yalvac, M. Zeyrek Bogazici University, Istanbul, Turkey E. Gu¨lmez, M. Kaya58, O. Kaya59, E.A. Yetkin60, T. Yetkin61 Istanbul Technical University, Istanbul, Turkey A. Cakir, K. Cankocak, S. Sen62 Institute for Scintillation Materials of National Academy of Science of Ukraine, Kharkov, Ukraine B. Grynyov National Scientific Center, Kharkov Institute of Physics and Technology, Kharkov, Ukraine L. Levchuk, P. Sorokin University of Bristol, Bristol, United Kingdom R. Aggleton, F. Ball, L. Beck, J.J. Brooke, D. Burns, E. Clement, D. Cussans, H. Flacher, J. Goldstein, M. Grimes, G.P. Heath, H.F. Heath, J. Jacob, L. Kreczko, C. Lucas, Z. Meng, D.M. Newbold63, S. Paramesvaran, A. Poll, T. Sakuma, S. Seif El Nasr-storey, S. Senkin, D. Smith, V.J. Smith Rutherford Appleton Laboratory, Didcot, United Kingdom K.W. Bell, A. Belyaev64, C. Brew, R.M. Brown, L. Calligaris, D. Cieri, D.J.A. Cockerill, J.A. Coughlan, K. Harder, S. Harper, E. Olaiya, D. Petyt, C.H. Shepherd-Themistocleous, A. Thea, I.R. Tomalin, T. Williams, S.D. Worm Imperial College, London, United Kingdom M. Baber, R. Bainbridge, O. Buchmuller, A. Bundock, D. Burton, S. Casasso, M. Citron, D. Colling, L. Corpe, P. Dauncey, G. Davies, A. De Wit, M. Della Negra, P. Dunne, A. Elwood, D. Futyan, Y. Haddad, G. Hall, G. Iles, R. Lane, R. Lucas63, L. Lyons, A.-M. Magnan, S. Malik, L. Mastrolorenzo, J. Nash, A. Nikitenko49, J. Pela, B. Penning, M. Pesaresi, D.M. Raymond, A. Richards, A. Rose, C. Seez, A. Tapper, K. Uchida, M. Vazquez Acosta65, T. Virdee15, S.C. Zenz Brunel University, Uxbridge, United Kingdom J.E. Cole, P.R. Hobson, A. Khan, P. Kyberd, D. Leslie, I.D. Reid, P. Symonds, L. Teodorescu, M. Turner Baylor University, Waco, USA A. Borzou, K. Call, J. Dittmann, K. Hatakeyama, H. Liu, N. Pastika The University of Alabama, Tuscaloosa, USA O. Charaf, S.I. Cooper, C. Henderson, P. Rumerio Boston University, Boston, USA D. Arcaro, A. Avetisyan, T. Bose, D. Gastler, D. Rankin, C. Richardson, J. Rohlf, L. Sulak, D. Zou 59 Brown University, Providence, USA J. Alimena, G. Benelli, E. Berry, D. Cutts, A. Ferapontov, A. Garabedian, J. Hakala, U. Heintz, O. Jesus, E. Laird, G. Landsberg, Z. Mao, M. Narain, S. Piperov, S. Sagir, R. Syarif University of California, Davis, Davis, USA R. Breedon, G. Breto, M. Calderon De La Barca Sanchez, S. Chauhan, M. Chertok, J. Conway, R. Conway, P.T. Cox, R. Erbacher, C. Flores, G. Funk, M. Gardner, W. Ko, R. Lander, C. Mclean, M. Mulhearn, D. Pellett, J. Pilot, F. Ricci-Tam, S. Shalhout, J. Smith, M. Squires, D. Stolp, M. Tripathi, S. Wilbur, R. Yohay University of California, Los Angeles, USA R. Cousins, P. Everaerts, A. Florent, J. Hauser, M. Ignatenko, D. Saltzberg, E. Takasugi, V. Valuev, M. Weber University of California, Riverside, Riverside, USA K. Burt, R. Clare, J. Ellison, J.W. Gary, G. Hanson, J. Heilman, P. Jandir, E. Kennedy, F. Lacroix, O.R. Long, M. Malberti, M. Olmedo Negrete, M.I. Paneva, A. Shrinivas, H. Wei, S. Wimpenny, B. R. Yates University of California, San Diego, La Jolla, USA J.G. Branson, G.B. Cerati, S. Cittolin, R.T. D’Agnolo, M. Derdzinski, R. Gerosa, A. Holzner, R. Kelley, D. Klein, J. Letts, I. Macneill, D. Olivito, S. Padhi, M. Pieri, M. Sani, V. Sharma, S. Simon, M. Tadel, A. Vartak, S. Wasserbaech66, C. Welke, J. Wood, F. Wu¨rthwein, A. Yagil, G. Zevi Della Porta University of California, Santa Barbara, Santa Barbara, USA J. Bradmiller-Feld, C. Campagnari, A. Dishaw, V. Dutta, K. Flowers, M. Franco Sevilla, P. Geffert, C. George, F. Golf, L. Gouskos, J. Gran, J. Incandela, N. Mccoll, S.D. Mullin, J. Richman, D. Stuart, I. Suarez, C. West, J. Yoo California Institute of Technology, Pasadena, USA D. Anderson, A. Apresyan, J. Bendavid, A. Bornheim, J. Bunn, Y. Chen, J. Duarte, A. Mott, H.B. Newman, C. Pena, M. Spiropulu, J.R. Vlimant, S. Xie, R.Y. Zhu Carnegie Mellon University, Pittsburgh, USA M.B. Andrews, V. Azzolini, A. Calamba, B. Carlson, T. Ferguson, M. Paulini, J. Russ, M. Sun, H. Vogel, I. Vorobiev University of Colorado Boulder, Boulder, USA J.P. Cumalat, W.T. Ford, F. Jensen, A. Johnson, M. Krohn, T. Mulholland, K. Stenson, S.R. Wagner Cornell University, Ithaca, USA J. Alexander, A. Chatterjee, J. Chaves, J. Chu, S. Dittmer, N. Eggert, N. Mirman, G. Nicolas Kaufman, J.R. Patterson, A. Rinkevicius, A. Ryd, L. Skinnari, L. Soffi, W. Sun, S.M. Tan, W.D. Teo, J. Thom, J. Thompson, J. Tucker, Y. Weng, L. Winstrom, P. Wittich Fermi National Accelerator Laboratory, Batavia, USA S. Abdullin, M. Albrow, G. Apollinari, S. Banerjee, L.A.T. Bauerdick, A. Beretvas, J. Berryhill, P.C. Bhat, G. Bolla, K. Burkett, J.N. Butler, H.W.K. Cheung, F. Chlebana, S. Cihangir, M. Cremonesi, V.D. Elvira, I. Fisk, J. Freeman, E. Gottschalk, L. Gray, D. Green, S. Gru¨nendahl, O. Gutsche, D. Hare, R.M. Harris, S. Hasegawa, J. Hirschauer, Z. Hu, B. Jayatilaka, S. Jindariani, M. Johnson, U. Joshi, B. Klima, B. Kreis, S. Lammel, J. Lewis, J. Linacre, D. Lincoln, R. Lipton, T. Liu, R. Lopes De Sa´, J. Lykken, K. Maeshima, J.M. Marraffino, S. Maruyama, D. Mason, 60 A The CMS Collaboration P. McBride, P. Merkel, S. Mrenna, S. Nahn, C. Newman-Holmes†, V. O’Dell, K. Pedro, O. Prokofyev, G. Rakness, E. Sexton-Kennedy, A. Soha, W.J. Spalding, L. Spiegel, S. Stoynev, N. Strobbe, L. Taylor, S. Tkaczyk, N.V. Tran, L. Uplegger, E.W. Vaandering, C. Vernieri, M. Verzocchi, R. Vidal, M. Wang, H.A. Weber, A. Whitbeck University of Florida, Gainesville, USA D. Acosta, P. Avery, P. Bortignon, D. Bourilkov, A. Brinkerhoff, A. Carnes, M. Carver, D. Curry, S. Das, R.D. Field, I.K. Furic, J. Konigsberg, A. Korytov, K. Kotov, P. Ma, K. Matchev, H. Mei, P. Milenovic67, G. Mitselmakher, D. Rank, R. Rossin, L. Shchutska, D. Sperka, N. Terentyev, L. Thomas, J. Wang, S. Wang, J. Yelton Florida International University, Miami, USA S. Linn, P. Markowitz, G. Martinez, J.L. Rodriguez Florida State University, Tallahassee, USA A. Ackert, J.R. Adams, T. Adams, A. Askew, S. Bein, J. Bochenek, B. Diamond, J. Haas, S. Hagopian, V. Hagopian, K.F. Johnson, A. Khatiwada, H. Prosper, A. Santra, M. Weinberg Florida Institute of Technology, Melbourne, USA M.M. Baarmand, V. Bhopatkar, S. Colafranceschi68, M. Hohlmann, H. Kalakhety, D. Noonan, T. Roy, F. Yumiceva University of Illinois at Chicago (UIC), Chicago, USA M.R. Adams, L. Apanasevich, D. Berry, R.R. Betts, I. Bucinskaite, R. Cavanaugh, O. Evdokimov, L. Gauthier, C.E. Gerber, D.J. Hofman, P. Kurt, C. O’Brien, I.D. Sandoval Gonzalez, P. Turner, N. Varelas, Z. Wu, M. Zakaria, J. Zhang The University of Iowa, Iowa City, USA B. Bilki69, W. Clarida, K. Dilsiz, S. Durgut, R.P. Gandrajula, M. Haytmyradov, V. Khristenko, J.-P. Merlo, H. Mermerkaya70, A. Mestvirishvili, A. Moeller, J. Nachtman, H. Ogul, Y. Onel, F. Ozok71, A. Penzo, C. Snyder, E. Tiras, J. Wetzel, K. Yi Johns Hopkins University, Baltimore, USA I. Anderson, B. Blumenfeld, A. Cocoros, N. Eminizer, D. Fehling, L. Feng, A.V. Gritsan, P. Maksimovic, M. Osherson, J. Roskes, U. Sarica, M. Swartz, M. Xiao, Y. Xin, C. You The University of Kansas, Lawrence, USA P. Baringer, A. Bean, C. Bruner, J. Castle, R.P. Kenny III, A. Kropivnitskaya, D. Majumder, M. Malek, W. Mcbrayer, M. Murray, S. Sanders, R. Stringer, Q. Wang Kansas State University, Manhattan, USA A. Ivanov, K. Kaadze, S. Khalil, M. Makouski, Y. Maravin, A. Mohammadi, L.K. Saini, N. Skhirtladze, S. Toda Lawrence Livermore National Laboratory, Livermore, USA D. Lange, F. Rebassoo, D. Wright University of Maryland, College Park, USA C. Anelli, A. Baden, O. Baron, A. Belloni, B. Calvert, S.C. Eno, C. Ferraioli, J.A. Gomez, N.J. Hadley, S. Jabeen, R.G. Kellogg, T. Kolberg, J. Kunkle, Y. Lu, A.C. Mignerey, Y.H. Shin, A. Skuja, M.B. Tonjes, S.C. Tonwar Massachusetts Institute of Technology, Cambridge, USA A. Apyan, R. Barbieri, A. Baty, R. Bi, K. Bierwagen, S. Brandt, W. Busza, I.A. Cali, Z. Demiragli, L. Di Matteo, G. Gomez Ceballos, M. Goncharov, D. Gulhan, D. Hsu, Y. Iiyama, G.M. Innocenti, 61 M. Klute, D. Kovalskyi, K. Krajczar, Y.S. Lai, Y.-J. Lee, A. Levin, P.D. Luckey, A.C. Marini, C. Mcginn, C. Mironov, S. Narayanan, X. Niu, C. Paus, C. Roland, G. Roland, J. Salfeld-Nebgen, G.S.F. Stephans, K. Sumorok, K. Tatar, M. Varma, D. Velicanu, J. Veverka, J. Wang, T.W. Wang, B. Wyslouch, M. Yang, V. Zhukova University of Minnesota, Minneapolis, USA A.C. Benvenuti, B. Dahmes, A. Evans, A. Finkel, A. Gude, P. Hansen, S. Kalafut, S.C. Kao, K. Klapoetke, Y. Kubota, Z. Lesko, J. Mans, S. Nourbakhsh, N. Ruckstuhl, R. Rusack, N. Tambe, J. Turkewitz University of Mississippi, Oxford, USA J.G. Acosta, S. Oliveros University of Nebraska-Lincoln, Lincoln, USA E. Avdeeva, R. Bartek, K. Bloom, S. Bose, D.R. Claes, A. Dominguez, C. Fangmeier, R. Gonzalez Suarez, R. Kamalieddin, D. Knowlton, I. Kravchenko, F. Meier, J. Monroy, F. Ratnikov, J.E. Siado, G.R. Snow, B. Stieger State University of New York at Buffalo, Buffalo, USA M. Alyari, J. Dolen, J. George, A. Godshalk, C. Harrington, I. Iashvili, J. Kaisen, A. Kharchilava, A. Kumar, A. Parker, S. Rappoccio, B. Roozbahani Northeastern University, Boston, USA G. Alverson, E. Barberis, D. Baumgartel, M. Chasco, A. Hortiangtham, A. Massironi, D.M. Morse, D. Nash, T. Orimoto, R. Teixeira De Lima, D. Trocino, R.-J. Wang, D. Wood, J. Zhang Northwestern University, Evanston, USA S. Bhattacharya, K.A. Hahn, A. Kubik, J.F. Low, N. Mucia, N. Odell, B. Pollack, M.H. Schmitt, K. Sung, M. Trovato, M. Velasco University of Notre Dame, Notre Dame, USA N. Dev, M. Hildreth, C. Jessop, D.J. Karmgard, N. Kellams, K. Lannon, N. Marinelli, F. Meng, C. Mueller, Y. Musienko38, M. Planer, A. Reinsvold, R. Ruchti, N. Rupprecht, G. Smith, S. Taroni, N. Valls, M. Wayne, M. Wolf, A. Woodard The Ohio State University, Columbus, USA L. Antonelli, J. Brinson, B. Bylsma, L.S. Durkin, S. Flowers, A. Hart, C. Hill, R. Hughes, W. Ji, B. Liu, W. Luo, D. Puigh, M. Rodenburg, B.L. Winer, H.W. Wulsin Princeton University, Princeton, USA O. Driga, P. Elmer, J. Hardenbrook, P. Hebda, S.A. Koay, P. Lujan, D. Marlow, T. Medvedeva, M. Mooney, J. Olsen, C. Palmer, P. Piroue´, D. Stickland, C. Tully, A. Zuranski University of Puerto Rico, Mayaguez, USA S. Malik Purdue University, West Lafayette, USA A. Barker, V.E. Barnes, D. Benedetti, L. Gutay, M.K. Jha, M. Jones, A.W. Jung, K. Jung, D.H. Miller, N. Neumeister, B.C. Radburn-Smith, X. Shi, J. Sun, A. Svyatkovskiy, F. Wang, W. Xie, L. Xu Purdue University Calumet, Hammond, USA N. Parashar, J. Stupak 62 A The CMS Collaboration Rice University, Houston, USA A. Adair, B. Akgun, Z. Chen, K.M. Ecklund, F.J.M. Geurts, M. Guilbaud, W. Li, B. Michlin, M. Northup, B.P. Padley, R. Redjimi, J. Roberts, J. Rorie, Z. Tu, J. Zabel University of Rochester, Rochester, USA B. Betchart, A. Bodek, P. de Barbaro, R. Demina, Y.t. Duh, Y. Eshaq, T. Ferbel, M. Galanti, A. Garcia-Bellido, J. Han, O. Hindrichs, A. Khukhunaishvili, K.H. Lo, P. Tan, M. Verzetti Rutgers, The State University of New Jersey, Piscataway, USA J.P. Chou, E. Contreras-Campana, Y. Gershtein, T.A. Go´mez Espinosa, E. Halkiadakis, M. Heindl, D. Hidas, E. Hughes, S. Kaplan, R. Kunnawalkam Elayavalli, S. Kyriacou, A. Lath, K. Nash, S. Randall, H. Saka, S. Salur, S. Schnetzer, D. Sheffield, S. Somalwar, R. Stone, S. Thomas, P. Thomassen, M. Walker University of Tennessee, Knoxville, USA M. Foerster, J. Heideman, G. Riley, K. Rose, S. Spanier, K. Thapa Texas A&M University, College Station, USA O. Bouhali72, A. Castaneda Hernandez72, A. Celik, M. Dalchenko, M. De Mattia, A. Delgado, S. Dildick, R. Eusebi, J. Gilmore, T. Huang, T. Kamon73, V. Krutelyov, R. Mueller, I. Osipenkov, Y. Pakhotin, R. Patel, A. Perloff, L. Pernie`, D. Rathjens, A. Rose, A. Safonov, A. Tatarinov, K.A. Ulmer Texas Tech University, Lubbock, USA N. Akchurin, C. Cowden, J. Damgov, C. Dragoiu, P.R. Dudero, J. Faulkner, S. Kunori, K. Lamichhane, S.W. Lee, T. Libeiro, S. Undleeb, I. Volobouev, Z. Wang Vanderbilt University, Nashville, USA E. Appelt, A.G. Delannoy, S. Greene, A. Gurrola, R. Janjam, W. Johns, C. Maguire, Y. Mao, A. Melo, H. Ni, P. Sheldon, S. Tuo, J. Velkovska, Q. Xu University of Virginia, Charlottesville, USA M.W. Arenton, P. Barria, B. Cox, B. Francis, J. Goodell, R. Hirosky, A. Ledovskoy, H. Li, C. Neu, T. Sinthuprasith, X. Sun, Y. Wang, E. Wolfe, F. Xia Wayne State University, Detroit, USA C. Clarke, R. Harr, P.E. Karchin, C. Kottachchi Kankanamge Don, P. Lamichhane, J. Sturdy University of Wisconsin - Madison, Madison, WI, USA D.A. Belknap, D. Carlsmith, S. Dasu, L. Dodd, S. Duric, B. Gomber, M. Grothe, M. Herndon, A. Herve´, P. Klabbers, A. Lanaro, A. Levine, K. Long, R. Loveless, A. Mohapatra, I. Ojalvo, T. Perry, G.A. Pierro, G. Polese, T. Ruggles, T. Sarangi, A. Savin, A. Sharma, N. Smith, W.H. Smith, D. Taylor, P. Verwilligen, N. Woods †: Deceased 1: Also at Vienna University of Technology, Vienna, Austria 2: Also at State Key Laboratory of Nuclear Physics and Technology, Peking University, Beijing, China 3: Also at Institut Pluridisciplinaire Hubert Curien, Universite´ de Strasbourg, Universite´ de Haute Alsace Mulhouse, CNRS/IN2P3, Strasbourg, France 4: Also at Universidade Estadual de Campinas, Campinas, Brazil 5: Also at Centre National de la Recherche Scientifique (CNRS) - IN2P3, Paris, France 6: Also at Universite´ Libre de Bruxelles, Bruxelles, Belgium 7: Also at Laboratoire Leprince-Ringuet, Ecole Polytechnique, IN2P3-CNRS, Palaiseau, France 63 8: Also at Joint Institute for Nuclear Research, Dubna, Russia 9: Also at Helwan University, Cairo, Egypt 10: Now at Zewail City of Science and Technology, Zewail, Egypt 11: Also at Ain Shams University, Cairo, Egypt 12: Also at Fayoum University, El-Fayoum, Egypt 13: Now at British University in Egypt, Cairo, Egypt 14: Also at Universite´ de Haute Alsace, Mulhouse, France 15: Also at CERN, European Organization for Nuclear Research, Geneva, Switzerland 16: Also at Skobeltsyn Institute of Nuclear Physics, Lomonosov Moscow State University, Moscow, Russia 17: Also at Tbilisi State University, Tbilisi, Georgia 18: Also at RWTH Aachen University, III. Physikalisches Institut A, Aachen, Germany 19: Also at University of Hamburg, Hamburg, Germany 20: Also at Brandenburg University of Technology, Cottbus, Germany 21: Also at Institute of Nuclear Research ATOMKI, Debrecen, Hungary 22: Also at MTA-ELTE Lendu¨let CMS Particle and Nuclear Physics Group, Eo¨tvo¨s Lora´nd University, Budapest, Hungary 23: Also at University of Debrecen, Debrecen, Hungary 24: Also at Indian Institute of Science Education and Research, Bhopal, India 25: Also at University of Visva-Bharati, Santiniketan, India 26: Now at King Abdulaziz University, Jeddah, Saudi Arabia 27: Also at University of Ruhuna, Matara, Sri Lanka 28: Also at Isfahan University of Technology, Isfahan, Iran 29: Also at University of Tehran, Department of Engineering Science, Tehran, Iran 30: Also at Plasma Physics Research Center, Science and Research Branch, Islamic Azad University, Tehran, Iran 31: Also at Universita` degli Studi di Siena, Siena, Italy 32: Also at Purdue University, West Lafayette, USA 33: Now at Hanyang University, Seoul, Korea 34: Also at International Islamic University of Malaysia, Kuala Lumpur, Malaysia 35: Also at Malaysian Nuclear Agency, MOSTI, Kajang, Malaysia 36: Also at Consejo Nacional de Ciencia y Tecnologı´a, Mexico city, Mexico 37: Also at Warsaw University of Technology, Institute of Electronic Systems, Warsaw, Poland 38: Also at Institute for Nuclear Research, Moscow, Russia 39: Now at National Research Nuclear University ’Moscow Engineering Physics Institute’ (MEPhI), Moscow, Russia 40: Also at St. Petersburg State Polytechnical University, St. Petersburg, Russia 41: Also at University of Florida, Gainesville, USA 42: Also at California Institute of Technology, Pasadena, USA 43: Also at Faculty of Physics, University of Belgrade, Belgrade, Serbia 44: Also at INFN Sezione di Roma; Universita` di Roma, Roma, Italy 45: Also at National Technical University of Athens, Athens, Greece 46: Also at Scuola Normale e Sezione dell’INFN, Pisa, Italy 47: Also at National and Kapodistrian University of Athens, Athens, Greece 48: Also at Riga Technical University, Riga, Latvia 49: Also at Institute for Theoretical and Experimental Physics, Moscow, Russia 50: Also at Albert Einstein Center for Fundamental Physics, Bern, Switzerland 51: Also at Adiyaman University, Adiyaman, Turkey 52: Also at Mersin University, Mersin, Turkey 64 A The CMS Collaboration 53: Also at Cag University, Mersin, Turkey 54: Also at Piri Reis University, Istanbul, Turkey 55: Also at Gaziosmanpasa University, Tokat, Turkey 56: Also at Ozyegin University, Istanbul, Turkey 57: Also at Izmir Institute of Technology, Izmir, Turkey 58: Also at Marmara University, Istanbul, Turkey 59: Also at Kafkas University, Kars, Turkey 60: Also at Istanbul Bilgi University, Istanbul, Turkey 61: Also at Yildiz Technical University, Istanbul, Turkey 62: Also at Hacettepe University, Ankara, Turkey 63: Also at Rutherford Appleton Laboratory, Didcot, United Kingdom 64: Also at School of Physics and Astronomy, University of Southampton, Southampton, United Kingdom 65: Also at Instituto de Astrofı´sica de Canarias, La Laguna, Spain 66: Also at Utah Valley University, Orem, USA 67: Also at University of Belgrade, Faculty of Physics and Vinca Institute of Nuclear Sciences, Belgrade, Serbia 68: Also at Facolta` Ingegneria, Universita` di Roma, Roma, Italy 69: Also at Argonne National Laboratory, Argonne, USA 70: Also at Erzincan University, Erzincan, Turkey 71: Also at Mimar Sinan University, Istanbul, Istanbul, Turkey 72: Also at Texas A&M University at Qatar, Doha, Qatar 73: Also at Kyungpook National University, Daegu, Korea